Going beyond the Dipole Approximation for the Neutron Form factor

Usually the magnetic form factor is calculated by the dipole approximation (see appendix J). In order to go beyond the dipole approximation (see [26] chapter 11.6.2,equations 11.141-11.143, mind: there is a small error in equ 11.141 and 11.142 in this book, see below) for the magnetic form factor it is necessary to take into account the expectation values and (in addition to and ). Moreover the ground state wave functions, the -coefficients (see appendix K or [26],table 11.1) and the momentum quantum number for each magnetic ion is required. This information has to be given in the input file in the following format:

... # %SECTION 4% DESCRIPTION OF MAGNETIC UNIT CELL AND LIST OF MAGNETIC ATOMS # # # here follows the description of the magnetic unit cell with respect # to the primitive crystallographic unit cell: # 'nr1', 'nr2', 'nr3' ...the crystallographic unit cell has to be taken # nr1 nr2 and nr3 times along r1 r2 and r3, # respectively to get magnetic unit cell # 'nat' denotes the number of magnetic atoms in magnetic unit cell # # Temperature, External Magnetic Field: Magnetic Unit Cell #! T=1.3 K Ha=0 T Hb= 0 T Hc= 0 T: nr1=2 nr2=2 nr3=2 nat=8 # # # It follows a list of nat lines with to describe the magnetic moment configuration # Notes: # 'atom-filename' means the single ion property filename of this magnetic atom: # -it must contain the Formfactor Coefficients (e.g. see international tables) # Lande factor # Neutron Scattering Length (10^-12 cm) # -it may contain a Debey Waller Factor # 'da' 'db' and 'dc' are not used by the program (unless you enter a line #! use_dadbdc=1) # 'dr1','dr2' and 'dr3' refer to the primitive lattice given below # 'Ma','Mb','Mc' denote the magnetic moment components in Bohr magnetons # in case of non orthogonal lattices instead of Ma Mb Mc the components Mi Mj Mk # have to be given, which refer to an right handed orthogonal coordinate system # defined by j||b, k||(a x b) and i normal to k and j # <Sa> <La> <Sb> <Lb > <Sc> <Lc> (optional) denote the spin and orbital angular momentum components # 'Hxc1' 'Hxc2' 'Hxc3' (optional line, used to go beyond dipole approx for formfactor) # denote the exchange fields in meV # #{atom-file} da[a] db[b] dc[c] dr1[r1] dr2[r2] dr3[r3] <Ma> <Mb> <Mc> [mb] [optional <Sa> <La> <Sb> <Lb > <Sc> <Lc>] #{corresponding exchange fields Hxc [meV]- if passed to mcdiff only these are used for calculation (not the magnetic moments)} {Ce3p.sipf} 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 +0.69590 +0.69620 +0.00000 corresponding exchange fields Hxc [meV]--> +0.58640 +0.58660 +0.00000 {Ce3p.sipf} 0.50000 0.50000 0.50000 0.00000 0.00000 1.00000 -0.69590 -0.69620 +0.00000 corresponding exchange fields Hxc [meV]--> -0.58640 -0.58660 +0.00000 {Ce3p.sipf} 0.00000 1.00000 0.00000 0.00000 1.00000 0.00000 -0.69590 -0.69620 +0.00000 corresponding exchange fields Hxc [meV]--> -0.58640 -0.58660 +0.00000 {Ce3p.sipf} 0.50000 1.50000 0.50000 0.00000 1.00000 1.00000 +0.69590 +0.69620 +0.00000 corresponding exchange fields Hxc [meV]--> +0.58640 +0.58660 +0.00000 {Ce3p.sipf} 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 -0.69590 -0.69620 +0.00000 corresponding exchange fields Hxc [meV]--> -0.58640 -0.58660 +0.00000 {Ce3p.sipf} 1.50000 0.50000 0.50000 1.00000 0.00000 1.00000 +0.69590 +0.69620 +0.00000 corresponding exchange fields Hxc [meV]--> +0.58640 +0.58660 +0.00000 {Ce3p.sipf} 1.00000 1.00000 0.00000 1.00000 1.00000 0.00000 +0.69590 +0.69620 +0.00000 corresponding exchange fields Hxc [meV]--> +0.58640 +0.58660 +0.00000 {Ce3p.sipf} 1.50000 1.50000 0.50000 1.00000 1.00000 1.00000 -0.69590 -0.69620 +0.00000 corresponding exchange fields Hxc [meV]--> -0.58640 -0.58660 +0.00000

Such a file can be easily created by using the program *spins*. Alternatively, if no complete *mcphas* model
calculation is available, the meanfields can be varied such as to
give a good fit to experimental data.

If such information is given the formalism used by *mcdiff* changes in the following way.
(i) If possible, the expectation values of
are evaluated by the module functions *Lcalc, Scalc and mcalc*, respectively. These values are used to calculate the neutron intensities in dipole approximation (Idip).
(ii) For the full intensity calculation instead of equation (31) the magnetic structure factor is calculated according to

(46) |

Here the thermal expectation values
of the scattering operators
have to be evaluated
for each ion . Note that the scattering operators correspond to the
Fourier transform of the atomic magnetisation density:
.
If present in the single ion module, the function *mq* is called to calculate this quantity.

For example, for the single ion module *so1ion* the scattering operator
is calculated as follows (for *ic1ion* see section 14.6):
once the eigenstates of each ion are known in terms of
atomic wave functions (with
) this
evaluation can be performed by means of the following expressions for
the matrix elements of the scattering operator (compare chapter 11.6.2
in [26],equations
11.141-11.143, mind: there is a small error in equation 11.141
and 11.142 in this reference - the imaginary has been misplaced ... the following
equations have been corrected for this mistake):

Here the different symbols have the following meaning:

The coefficients have been tabulated (see appendix K or table 11.1 of [26]) and have to be given in the single ion parameter file(s) by entering the following lines

# coefficients of Z(K') according to Lovesey chapter 11.6.1 page 233 Z1c0=1.63636364 Z1c2=2.95041322 Z3c2=-0.20896503 Z3c4=-0.25329095 Z5c4=0.03820789 Z5c6=0.14258681 Z7c6=-0.00614959

The are the spherical harmonics evaluated for the direction of the scattering vector with respect to the crystal field coordinate system . The factor is defined in terms of -symbols.

Martin Rotter 2017-01-10