SAXIS

From VASP Wiki
Revision as of 09:55, 21 September 2023 by Huebsch (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

SAXIS = [real array]
Default: SAXIS = (0, 0, 1) 

Description: Set the global spin-quantization axis w.r.t. Cartesian coordinates.


SAXIS specifies the relative orientation of spinor space spanned by the Pauli matrices , , with respect to Cartesian coordinates . The default is , , . The direction of the spin-quantization axis with respect to Cartesian coordinates is set

 SAXIS =   sx sy sz    ! global spin-quantization axis

such that , i.e., points along . The directions of and are a consequence of rotating to point along as described below.

The relative orientation of spinor space with respect to real space becomes important in case spin-orbit coupling is included (LSORBIT=True). All magnetic moments and spinor-like quantities written or read by VASP are given in the basis of the spinor space , , . This includes the MAGMOM tag in the INCAR file, the total and local magnetizations in the OUTCAR and PROCAR file, the spinor-like orbitals in the WAVECAR file, and the magnetization density in the CHGCAR file.

Coordinate system

Fig 1. Euler angles and defined by .

The default orientation is , , . To set , VASP applies two rotations with Euler angles

Here, is the angle between the projection of SAXIS onto the xy plane (sx,sy,0) and the Cartesian vector , and is the angle between the vector SAXIS and the Cartesian vector , see Fig. 1. Search for `Euler angles` in the OUTCAR file to see what VASP uses. For the default , and .

The transformation of a vector given in the basis , , into in Cartesian coordinates and its inverse transformation read

where the rotation matrices are

Mind: Apply the proper basis transformation when comparing vector-like quantities and spinor-like quantities.

For instance, when LORBMOM=True the orbital angular momentum is written to the OUTCAR file in Cartesian coordinates. Thus, when comparing the orbital angular momentum (vector-like quantity) and the magnetization (spinor-like quantity), one has to perform a basis transformation on one of the quantities unless the bases agree (default).

Example

  • In case the bases have the same orientation, i.e., , , (default)
For a single site this implies setting
MAGMOM = mx my mz ! magnetic moment in Cartesian coordinates
SAXIS =  0 0 1   ! default
Fig 2. Example with and Euler angles and .
  • Another good choice is setting to point along the direction of the on-site magnetic moment such that
where is the total on-site magnetic moment.
For a single site, this case implies setting
MAGMOM = 0 0 m   ! magnetic moment along sigma3
SAXIS =  sx sy sz ! direction of sigma3
Thus, there are two methods to rotate the initial magnetization in an arbitrary direction: either by changing the initial magnetic moments MAGMOM or by changing SAXIS. Both methods should, in principle, yield exactly the same energy, but for implementation reasons, the second method might be more precise.
  • In case
SAXIS =  1 1 0   ! alpha=pi/4, beta=pi/2
the spinor space , , will be rotated with respect to real space as shown in Fig. 2.

Related tags and articles

LNONCOLLINEAR, MAGMOM, LSORBIT

Examples that use this tag