A fully unconstrained noncollinear density functional theory (DFT) is supported including the spinorbit coupling (SOC) [8,9,10,11,16]. When the noncollinear DFT is performed, the following option for the keyword 'scf.SpinPolarization' is available.
scf.SpinPolarization NC # OnOffNCWhen the option 'NC' is specified, wave functions are expressed by a two component spinor. An initial spin orientation of each site is given by
<Atoms.SpeciesAndCoordinates # Unit=Ang 1 Mn 0.00000 0.00000 0.00000 8.0 5.0 45.0 0.0 45.0 0.0 1 on 2 O 1.70000 0.00000 0.00000 3.0 3.0 45.0 0.0 45.0 0.0 1 on Atoms.SpeciesAndCoordinates>
1: sequential serial number 2: species name 3: xcoordinate 4: ycoordinate 5: zcoordinate 6: initial occupation for up spin 7: initial occupation for down spin 8: Euler angle, theta, of the magnetic field for spin magnetic moment 9: Euler angle, phi, of the magnetic field for spin magnetic moment Also, the 8th and 9th are used to generate the initial noncollinear spin charge distribution 10: the Euler angle, theta, of the magnetic field for orbital magnetic moment 11: the Euler angle, phi, of the magnetic field for orbital magnetic moment 12: switch for the constraint schemes specified by the keywords 'scf.Constraint.NC.Spin', 'scf.NC.Zeeman.Orbital' and 'scf.NC.Zeeman.Orbital'. '1' means that the constraint is applied, and '0' no constraint. 13: switch for enhancement of orbital polarization in the LDA+U method, 'on' means that the enhancement is made, 'off' no enhancement.The initial Euler angles, and , for orientation of the spin and orbital magnetic moment are given by the 8th and 9th columns, and 10th and 11th columns, respectively. The 12th column is a switch for a constraint scheme that a constraint (penalty or Zeeman) functional to the spin and orbital orientation is added on each site, where '1' means that the constraint functional is added, and '0' means no constraint. For the details of the constraint DFT for the spin orientation, see Sec. 38 'Constraint DFT for noncollinear spin orientation'. The final 13th column is a switch for enhancement of orbital polarization in the LDA+U method, 'on' means that the enhancement is made, 'off' no enhancement. Figure 32 shows the spin orientation in a MnO molecule calculated by the noncollinear DFT. You can follow the calculation using an input file 'Mol_MnO_NC.dat' in the directory 'work'. To visualize the spin orientation in real space, two files are generated:
System.Name.nc.xsf System.Name.ncsden.xsfwhere System.Name means 'System.Name' you specified. Two files 'System.Name.nc.xsf' and 'System.Name.ncsden.xsf' store a projected spin orientation to each atom by Mulliken analysis and the spin orientation on real space grids in a vector file format (XSF) supported by XCrySDen. Both the files can be visualized using 'Display Forces' in XCrySDen as shown in Fig. 32.
The spin moment and Euler angles of each atom, which are calculated by Mulliken analysis, are found in the file 'System.Name.out' as follows:
*********************************************************** *********************************************************** Mulliken populations *********************************************************** *********************************************************** Total spin moment (muB) 4.998503442 Angles (Deg) 44.991211196 0.000000000 Up Down Sum Diff theta phi 1 Mn 9.59803 4.76902 14.36705 4.82901 44.99208 0.00000 2 O 3.40122 3.23173 6.63295 0.16949 44.96650 0.00000Also it should be noted that it is difficult to achieve a self consistent field in the noncollinear DFT more than the collinear DFT calculation, since there are many minima, having almost comparable energy, in the spin orientation space, while the constraint DFT is useful for such a case.
In the noncollinear DFT, the inclusion of spinorbit coupling is supported, while it is not supported for the collinear DFT. See also the Section 'Relativistic effects' for the issue.
