Spin orbit using initial CHGCAR file from previoulsy converged spin-polarized calculation

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
chelman
Newbie
Newbie
Posts: 13
Joined: Fri Mar 11, 2011 4:58 pm
License Nr.: 5-793
Location: Buenos Aires, Argentina

Spin orbit using initial CHGCAR file from previoulsy converged spin-polarized calculation

#1 Post by chelman » Thu May 12, 2011 12:47 am

Hi
I'm performing spin-orbit calculation to obtain the MAE (magnetic anisotropy energy). The manual states that it is possible to calculate the spin orbit from a converged spin polarized calculation, using the CHGCAR. So therefore I calculated for different directions. For each different direction the initial CHGCAR was spin polarized. But, for reasons that are long to explain in this mail, I did the same calculation but the initial CHGCAR was the previously calculated with spin orbit. To clarify, for the first direction (could be SAXIS=1 0 0) the initial CHGCAR was the obtained in spin polarized calculation. For the second direction (let's said SAXIS = sqr(1/2) sqr(1/2) 0) the initial CHGCAR was the one obtained for the first direction. For the third direction the initial CHGCAR was the one obtained for the second direction, and so on different directions.
What I obtain is:
in-plane angle Energy(1)(eV) Energy(2)(eV)
0 -192.96158 -192.96158
30 -192.96173 -192.96175
45 -192.96189 -192.96193
60 -192.96204 -192.96212
90 -192.96219 -192.96232
120 -192.96204 -192.96225
180 -192.96158 -192.96198
(1) initial CHGCAR from converged spin polarized calculation
(2) initial CHGCAR from the converged for the previous angle (except 0? was from spin pol. calc.)

So three main question:
1) Could somebody explain me why there is an energy difference between columns (1) and (2)?
2) Most important, Since the system has C_2 symmetry, at 180? the energy should be the same at 0?,so what happen in column (2), the energy at 0? should be the same at 180??
3) Could somebody provide me with any reference about the way spin orbit is implemented on non-collinear calculation in VASP?

Thanks in advance
Last edited by chelman on Thu May 12, 2011 12:47 am, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Spin orbit using initial CHGCAR file from previoulsy converged spin-polarized calculation

#2 Post by admin » Thu May 26, 2011 4:28 pm

how large EDIFF in these calculations?
Last edited by admin on Thu May 26, 2011 4:28 pm, edited 1 time in total.

chelman
Newbie
Newbie
Posts: 13
Joined: Fri Mar 11, 2011 4:58 pm
License Nr.: 5-793
Location: Buenos Aires, Argentina

Spin orbit using initial CHGCAR file from previoulsy converged spin-polarized calculation

#3 Post by chelman » Mon May 30, 2011 6:13 pm

The EDIFF is 0.1E-04 for both cases.
Thanks
Last edited by chelman on Mon May 30, 2011 6:13 pm, edited 1 time in total.

Post Reply