Non-collinear magnetism setup
To configure a Fleur calculation incorporating non-collinear magnetism, some parameters have to be set in
the calculationSetup section and further parameters have to be set for each atomGroup in
the atomGroups section. Templates with
default parameters are generated by using the input generator with the -explicit or -noco command line options.
An example for the magnetic input elements in the calculation setup section is:
<magnetism jspins="1" l_noco="F" l_ss="F" lflip="F">
<qss>.0000000000 .0000000000 .0000000000</qss>
<mtNocoParams l_mperp="F" l_mtNocoPot="F" l_relaxSQA="F" mag_mixing_scheme="0" mix_RelaxWeightOffD="1.00000000" l_constrained="F" mix_constr="1.00000000"/>
<sourceFreeMag l_sourceFree="F" l_scaleMag="F" mag_scale="1.00000000"/>
</magnetism>
The following attributes can/have to be modified here:
| Tag | Attribute | Description |
|---|---|---|
| magnetism | jspins | Number of spins to consider, i.e. jspins="1" for non-magnetic, and ="2" for magnetic calculations |
| l_noco | This boolean switch activated non-collinear calculations | |
| l_ss | This boolean switch is used to activate spin-spiral calculations. | |
| lfip | Flip the magnetisation. See section below for details | |
| qss | In case of a spin-spiral calculation you have to give the q-vector of the spiral here | |
| mtNocoParams | In this tag flags controlling the treatment of magnetism in the MT-spheres can be set. | |
| l_mperp | Here the output of the magnetization perpendicular to the chosen axis can be activated. | |
| l_mtNocoPot | Switching this to "T" will allow a fully unconstrained calculation in which the MT spin-offdiagonal potential is calculated and use in the SCF. (Needs l_mperp="T") | |
| l_relaxSQA | Adjust the local spin-quantization axis in the sphere in each iteration to the direction of the local magnetic moment. Needed for some features like LDA+U in a l_mtNocoPot="t" calculation. | |
| mag_mixing_scheme | In case of l_relaxSQA="T" different preconditioners can be tried out to accelerate the convergence of the magnetic direction. Currently values of 1-3 are implemented | |
| mix_RelaxWeightOffD | Weight used to increase the local off-diagonal magnetisation in the mag_mixing_scheme not "0" cases. Should be >1 to have an effect. | |
| l_constrained | Perform a calculation in which an additional constraining field is used to fix the direction of the moments. Typical value "0.5" | |
| mix_constr | Mixing factor used to determine the constraining field self-consistently | |
| sourceFreeMag | Switches to use a source-free correction to the XC-B-Field. (Experimental) |
In addition to these global switches there is the nocoParams tag in each atomGroup to control magnetic calculations.
<nocoParams alpha="Pi/2.0" beta="Pi/4.0"
l_constrained="F" l_mtNocoPot="F" l_relaxSQA="F"/>
The following attributes have to be set here:
| Attribute | Description |
|---|---|
| alpha | The 1st angle that determines the magnetic structure. It is equal to in spherical coordinates. |
| beta | The 2nd angle that determines the magnetic structure. It is equal to (measured from the z axis) in spherical coordinates. |
| Please note, that in a l_relaxSQA="T" calculation, these are only used in constructing the starting charge as the correspong angles are determined from the local magnetic moment | |
| l_constrained | These switches overwrite the corresponding global switches (optional). |
| l_mtNocoPot | These switches overwrite the corresponding global switches (optional). |
| l_relaxSQA | These switches overwrite the corresponding global switches (optional). |
Manipulating the magnetism
If you set the lfip="T" switch in the magnetism tag, you can manipulate the magnetisation in the spheres. To do so you need to specify a modInitDen flag in the atomGroup tag. Please be careful about the workflow as you do not want to do this contineously.
<modInitDen flipSpinTheta="0.0" flipSpinPhi="0.0" flipSpinScale="F"/>
| Attribute | Description |
|---|---|
| flipSpinTheta | Rotate the magnetisation by this angle around the local "y" axis. This is done first. |
| flipSpinPhi | Rotate the magnetisation by this angle around local "z" axis. Done after "theta" rotation. |
| flipSpinScale | If "T" flip the local spin. |
| beta | The 2nd angle that determines the magnetic structure. It is equal to (measured from the z axis) in spherical coordinates. |
| magMom | (Experimental) |