I noticed the following description given by pyprocar documentation:
OTOH, I also noticed the following description on vaspwiki:
In case of non-collinear calculations, it seems that the forms of the two equations mentioned above do not correspond exactly.
Any tips/comments on the differences and connections between them will be highly appreciated?
Regards,
HZ
About the equation used by non-collinear electronic structure calculation.
Moderators: Global Moderator, Moderator
-
- Full Member
- Posts: 191
- Joined: Tue Oct 13, 2020 11:32 pm
-
- Full Member
- Posts: 212
- Joined: Tue Jan 19, 2021 12:01 am
Re: About the equation used by non-collinear electronic structure calculation.
Hi HZ,
Could you please clarify your question? In particular, can you specify what you mean by "it seems that the forms of the two equations mentioned above do not correspond exactly"?
Cheers,
Marie-Therese
Could you please clarify your question? In particular, can you specify what you mean by "it seems that the forms of the two equations mentioned above do not correspond exactly"?
Cheers,
Marie-Therese
-
- Full Member
- Posts: 191
- Joined: Tue Oct 13, 2020 11:32 pm
Re: About the equation used by non-collinear electronic structure calculation.
I mean whether the two equations in the red box below are equivalent and can be derived from each other:
1. The version used by pyprocar:
2. The version described by vaspwiki:
1. The version used by pyprocar:
2. The version described by vaspwiki:
-
- Full Member
- Posts: 212
- Joined: Tue Jan 19, 2021 12:01 am
Re: About the equation used by non-collinear electronic structure calculation.
Well, the first red box is not an equation.
If you are wondering about the sum over k points, it is required as the projected charge density is given in real space. fnk is the occupation number, to show that the sum is over occupied states.
The second red box are the projected components. To obtain the magnetization as a vector, look at the definition of m^j_alpha,l below. If you were to sum over all k points in the first red box, then it equals to m^i_alpha,l.
Best regards,
Marie-Therese
If you are wondering about the sum over k points, it is required as the projected charge density is given in real space. fnk is the occupation number, to show that the sum is over occupied states.
The second red box are the projected components. To obtain the magnetization as a vector, look at the definition of m^j_alpha,l below. If you were to sum over all k points in the first red box, then it equals to m^i_alpha,l.
Best regards,
Marie-Therese