Dear all,
may I ask if there is a way to compute the occupation matrix, which VASP prints out for LDAUPRINT = 2, in a different frame. What I meant is that VASP computes this matrix for the real spherical harmonics oriented along the global axis. However, if I want to know this projection on a different axis how should I go about it? In principle, I could rotate the global frame via the lattice vectors, which however is quite inconvenient if I have multiple distorted octahedra with different axis.
So is there a way to specify a lokal axis for this projection,
or
does anybody has a reference how to rotate the occupation matrix VASP computes by hand?
thanks for any reply.
occupation matrix rotation
Moderators: Global Moderator, Moderator
-
andreas.singraber
- Global Moderator
- Posts: 236
- Joined: Mon Apr 26, 2021 7:40 am
Re: occupation matrix rotation
Hello!
Sorry for the late reply.. unfortunately there is no feature built into VASP to allow for this rotation. I am not sure if it is relevant for your work but maybe you find this information useful:
https://www.chemistry.tcd.ie/staff/peop ... thodology/
There is also a patch for VASP 5.4.4 available on Github here:
https://github.com/WatsonGroupTCD/Occup ... ol-in-VASP
Best,
Andreas
Sorry for the late reply.. unfortunately there is no feature built into VASP to allow for this rotation. I am not sure if it is relevant for your work but maybe you find this information useful:
https://www.chemistry.tcd.ie/staff/peop ... thodology/
There is also a patch for VASP 5.4.4 available on Github here:
https://github.com/WatsonGroupTCD/Occup ... ol-in-VASP
Best,
Andreas