Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.

Messages - Imane

Pages: [1]

Could you please help me to extract data of the projected density of states for spin polarized calculations (spin up and spin down). I used this script already posted here:

dos = nlread('', DensityOfStates)[0]
energies = dos.energies()
DOS_s = open('file_s.dat','w')
projection_list_1 = ProjectionList(atoms=[0,1,2,3])
dos_s = dos.evaluate(Spin.All,projection_list_1)
dos_s = dos_s[0].inUnitsOf(eV**-1)
for i in range(0,len(dos_s)):
     DOS_s.write("%16.6f %16.6f \n" %(energies, 2*dos_s))

and I change Spin.All by Spin.Up/Spin.down, but it doesn't work.

Thank you in advance

Pages: [1]