... Rotter¹

Independent international Scientific Consultant, email: martin_rotter@mcphase.de

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... Le²

Dept. Physics and Astronomy, Seoul National University, Seoul 151-742, Korea

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... Keller³

Universität Regensburg

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... Hoffmann⁴

University of Oxford, Physics Department, Clarendon Laboratory, Parks Road, Oxford, UK

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... Schedler⁵

Institut für Festkörperhysik, Technische Universität Dresden, D-01062 Dresden, Germany

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... Hoffmann⁶

Forschungszentrum Jülich, D-52425 Jülich, Germany

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... Rotter⁷

Wien, Austria

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... Banks⁸

Max Planck Institute, Stuttgart, Germany

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...üver ⁹

Gymnasium Dresden-Cotta

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...#tex2html_wrap_inline683#¹⁰

Note that these conditions are essential and put a limit to the applicability of the theory, for example in the case of charge transfer excitations from one subsystem to the next.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...cfze). ¹¹

In addition to $\langle \mathbf J_i \rangle$ the module also returns the partition sum $z$ and the magnetic energy $u=\sum_{i=1}^{2J+1} p_{i} \epsilon_{i}$.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.... ¹²

Note that if you use the module cfield, the choice is more unconventional:$\vec a\vert\vert y$, $\vec b\vert\vert z$ and $\vec c\vert\vert x$ Tools for rotating crystal field parameters are described in appendix I.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... edited).¹³

Alternatively, one may make use of the program spins.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...Eel). ¹⁴

$R^{\alpha}R^{\delta}\omega_{\sigma} \epsilon_{\sigma\alpha\beta}\delta_{\beta\... ...\sigma\eta}\delta_{\alpha\delta})= (\vec R . \vec \omega)^2-R^2\omega^2 \neq 0$, if $\vec R$ is not parallel to the rotation axis. Thus for some rotation the elastic energy will depend on rotation angle if transversal springs are introduced.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... using ¹⁵

S. Bluegel, Juelich, private communication

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... first¹⁶

Note that we do not calculate all the terms in the braces in the first line of equation 4. In particular using the central field approximation the first summation produces Hydrogen-like energy levels, called configurations, which are split into terms by the second summation. What we refer to as the Coulomb interaction is only this second summation, and we shall consider the lowest energy configuration only. This configuration corresponds to the outer most electrons in ion, and is labelled $nl^\nu$, for example, $4f^1$ for Ce$^{3+}$.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... properties¹⁷

For the case of $LS$-coupling this corresponds to ignoring the Coulomb and SO interactions, which are in this limit both much larger than the CF. For $jj$-coupling, the Coulomb interaction is treated as small and neglected, but the spin-orbit interaction is considered.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...#tex2html_wrap_inline9331#¹⁸

Also denoted $m_J$.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...racah49-1352¹⁹

Equation 11

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... formulae²⁰

Except for the case of the rotation group in 3 dimensions, SO(3) (whose representations are labelled by the angular momentum quantum numbers, $L$ and $J$), where they are simply the $3j$ symbols

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...allison74²¹

http://cpc.cs.qub.ac.uk/cpc/cgi-bin/showversions.pl/?catid=acry

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...elliot57-509²²

Eqn 16 for SO, and 25-27 for CF, and reproduced in more modern notation in appendix C of [45]

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...racah49-1352²³

Eqns 66 for the conversion of the Racah parameters to Slater integrals, 63 for the $\hat e_0$ operator, 69 for $\hat e_1$, 73-74 for $\hat e_2$, and 78, 28 and 80-87 for $\hat e_3$

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...judd88-1²⁴

Eqn 7-58

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... intervals²⁵

Depending on the Java setup on your computer, a conflict may occur between the McPhase Java programs and other Java programs, such as the Matlab GUI. If you find that no graphics windows open, and have other Java programs running in the background, try to close all other Java programs.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...wybourne65²⁶

Note that table 6-1 and equ 6-7 are not correct in this reference.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...#tex2html_wrap_inline16852#²⁷

Note that these conditions are essential and put a limit to the applicability of the theory, for example in the case of charge transfer excitations from one subsystem to the next [65].

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... susceptibility²⁸

the $-$ on top of $\chi$ indicates matrix notation for $\chi_{\alpha\beta}^{ss'}$

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... investigated²⁹

supplementary material - screen shot movie comparing the speed of traditional Green's function method and DMD

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... cell³⁰

In case of magnetic order. In general this will be the unit cell of the Bravais lattice in section M, which is a superlattice of the crystal lattice.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... elements³¹

using the appropriate expression for the matrix elements of of the scattering operator as given in [29, equ. (11.86) or (11.87b)].

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.