Zeeman Field Tutorial

Note: this page has not been carefully worked out, especially the description of the Stoner parameter may not be accurate.

You can add to the crystal potential an external magnetic field inside each augemetation sphere, in the ASA program lm and full-potential program lmf. Use token BFIELD to tell lm or lmf that a field is to be read. The calculation must be spin polarized, e.g.

HAM NSPIN=2 BFIELD=2
...

BFIELD=1 allows the field to point in any direction; BFIELD=2 restricts the field to the z axis. The former case generates a noncollinear hamiltonian; the latter is spin diagonal. (The hamiltonian may also be noncollinear for another reason, e.g. spin-orbit coupling is turned on.) Caution: lmf does not as yet generate noncollinear output densities. Self-consistent calculations with noncollinear B fields are not possible; also the total energy is not trustworthy if B is not along the z axis.

By selecting BFIELD nonzero, you are telling lm or lmf that a (site-dependent) Zeeman field is to be added. In this case you must also supply a site-dependent external field, which you do through the additional file bfield.ext. If this file is omitted, the program will not execute. bfield.ext contains one row for each site; each row contains three numbers corresponding to x,y,z components of B. (The file must always contain three components Bx.ByBz even if BFIELD=2.) The GdN test case (part of the test suite) is an instructive example. In the top-level directory, run this test:

fp/test/test.fp gdn

This test combines the LDA+U hamiltonian with spin-orbit coupling. Bands are calculated for that case; then a magnetic field of 0.25 eV=0.037/2 Ry is added to the Gd site and the bands are recalculated. GdN has two atoms and file bfield.gdn reads:

0 0 .037/2
0 0 0

The test generates files bnds.bfield.gdn and bnds.nobfield.gdn To analyze the bands you need to use a graphics package. If you have the FPLOT packaged installed, you can draw the picture shown below comparing the two cases, by invoking the following commands after running the test:

In the figure, blue bands are LDA+U results. Because spin-orbit coupling is included both minority and majority spin bands appear. You can see majority Gd f states near -6.5 eV, and minority f states near +5.5 eV. Gd d states are spin-split by the exchange-correlation field originating from spin splitting of the Gd f states; majority and minority d states are found near 2.5 eV and 3.5 eV, respectively, at the Γ point.

Red bands use the same potential as the LDA+U case, except that a Bz field of magnitude 0.5 eV was added to the Gd site. Bz causes the the following shifts:

The majority f shifts down by ~0.22 eV, and the minority f states shift up by ~0.27 eV. States that are purely atomic like and are little affected by the kinetic energy will be shifted by the Zeeman field, ±B/2. This is the case for Gd f orbitals.

The majority d states at Γ shift down by ~0.16 eV; the minority d states shift up by ~0.20 eV.

The As p states, in the (-4,0) eV range, are only slightly shifted.

Longitudinal susceptibility

The magnetic susceptibility is the system’s magnetic response to an external magnetic field:

In linear response theory, χ−1\chi^{-1}χ​−1​​ can be written as a sum of a noninteracting part and an interacting part

χ−1=χ0−1−I\chi^{-1} = \chi_{0}^{-1} - Iχ​−1​​=χ​0​−1​​−I

The kernel III is often associated with the “Stoner parameter.” It is typically about 1 eV in 3d transition metals. The noninteracting part χ0−1\chi_0^{-1}χ​0​−1​​ is connected with the induced magnetization of the noninteracting system, i.e.

δB=χ0−1δM0\delta B = \chi_{0}^{-1} \delta M_0δB=χ​0​−1​​δM​0​​

where δM0\delta M_0δM​0​​ is the magnetization induced by δB\delta BδB without taking into account the internal potential shifts that δM\delta MδM cause. In the density-functional context, δM0\delta M_0δM​0​​ can be obtained from an initially self-consistent potential V0V_0V​0​​, which is generated in the absence of an external field. δM0\delta M_0δM​0​​ is then the change in MMM for the potential V0+δBV_0 + \delta BV​0​​+δB. In practice this is obtained by evaluating the change in MMM from the change in output density of a 1-shot calculation with potential V0+δBV_0 + \delta BV​0​​+δB.

You can use lmf to apply a field and calculate E_LDA and δ_M as a function of B. By comparing self-consistent and 1-shot calculations you can in principle also obtain I. The example below does this in an approximate way for Fe, using input file fp/test/ctrl.felz. Note: a proper calculation of I would require a supercell large enough that the interactions between external B-fields would be sufficiently small. (I is a local function of r in density functional theory; thus I is a site-diagonal matrix.) Here we only consider a single atom cell. The interactions are not small, and I is underestimated. The results below were generated by this command:

As the top figure shows, the total energy E has a minimum near, but not exactly at B=0. This shows that there are some slight errors in the LDA calculation (in an exact self-consistent LDA calculation E should be minimum at exactly B=0). E is smooth in the interval (−0.01, 0.01) Ry, as the figure shows. The behavior is somewhat unsmooth near B=−0.01 Ry, and its discontinuity or near discontinuity is evident near B=0.01 Ry. That means numerical differentiation ∂E/∂B is a bit problematic, particularly outside the (−0.01,0.01) interval. Nevertheless we can perform it. In the bottom figure, M−M(B=0) and ∂E/∂B are plotted respectively as a function of B as dashed lines+circles, and a solid line. The two curves track fairly well, though ∂E/∂B is noisy; nor is M(B) particularly smooth.

We can obtain χ from the slope ∂E/∂B. There is some uncertainty in the calculation since M is not particularly smooth (note in particular that (∂E/∂B)+ and (∂E/∂B)− are a little different).

We can also obtain both χ0 by obtaining δM0 from a 1-shot calculation. The table below shows M0 and M for two values of B: