Constraining the local magnetic moments: Difference between revisions
No edit summary |
No edit summary |
||
(3 intermediate revisions by 2 users not shown) | |||
Line 1: | Line 1: | ||
Description: Constraining the local magnetic moments on NiO. | Description: Constraining the local magnetic moments on NiO. | ||
To | To orient locally each spin, a penalty energy <math>E_p</math> is added, inversely proportional to the <math>\lambda</math> ([[LAMBDA]]) parameter. Thus, the process has to begin with a low [[LAMBDA]], which is increased step by step. | ||
For additional | For additional information, go to the [[I_CONSTRAINED_M|I_CONSTRAINED_M]] page. | ||
''<u>Exercise :</u>'' Verify the | ''<u>Exercise :</u>'' Verify the <math>E_p=f(1/\lambda)</math> relation by constraining the spin directions. Check the efficiency of the method by constraining only the direction, then the direction and the size of spins. | ||
---- | ---- | ||
Line 49: | Line 49: | ||
{{TAGBL|LDAUU}} = 5.00 0.00 | {{TAGBL|LDAUU}} = 5.00 0.00 | ||
{{TAGBL|LDAUJ}} = 0.00 0.00 | {{TAGBL|LDAUJ}} = 0.00 0.00 | ||
{{TAGBL|LDAUPRINT}} = | {{TAGBL|LDAUPRINT}} = 1 | ||
{{TAGBL|LMAXMIX}} = 4 | {{TAGBL|LMAXMIX}} = 4 | ||
Line 75: | Line 75: | ||
1.5 1.5 1.5 | 1.5 1.5 1.5 | ||
</pre> | </pre> | ||
== Download == | == Download == | ||
[ | [[Media:nio_constr.tgz| nio_constr.tgz]] | ||
---- | ---- | ||
[[Category:Examples]] | [[Category:Examples]] |
Latest revision as of 08:49, 11 April 2023
Description: Constraining the local magnetic moments on NiO.
To orient locally each spin, a penalty energy is added, inversely proportional to the (LAMBDA) parameter. Thus, the process has to begin with a low LAMBDA, which is increased step by step.
For additional information, go to the I_CONSTRAINED_M page.
Exercise : Verify the relation by constraining the spin directions. Check the efficiency of the method by constraining only the direction, then the direction and the size of spins.
NiO GGA+U Constr. SYSTEM = "NiO" Electronic minimization ENCUT = 450 EDIFF = 1E-5 LORBIT = 11 LREAL = .False. ISYM = -1 NELMIN = 6 LSORBIT = .True. GGA_COMPAT = .FALSE. DOS ISMEAR = -5 Magnetism ISPIN = 2 I_CONSTRAINED_M = 1 # direction # I_CONSTRAINED_M = 2 # size and direction M_CONSTR = 2 0 0 0 -2 0 6*0 LAMBDA = 1 RWIGS = 1.30 0.70 Orbital moment LORBMOM = T Mixer AMIX = 0.2 BMIX = 0.00001 AMIX_MAG = 0.8 BMIX_MAG = 0.00001 GGA+U LDAU = .TRUE. LDAUTYPE = 2 LDAUL = 2 -1 LDAUU = 5.00 0.00 LDAUJ = 0.00 0.00 LDAUPRINT = 1 LMAXMIX = 4
k-points 0 gamma 4 4 4 0 0 0
NiO 4.17 1.0 0.5 0.5 0.5 1.0 0.5 0.5 0.5 1.0 2 2 Cartesian 0.0 0.0 0.0 1.0 1.0 1.0 0.5 0.5 0.5 1.5 1.5 1.5