Basic theory

Local fields in MuSR

Muon spin rotation and relaxation spectroscopy is mainly used to probe magnetic materials. We briefly describe here the interactions between the magnetic moments of the muon and the electrons in the host system that produce a local magnetic field at the muon site in magnetically ordered samples.

Dipolar Field

The dipolar field is produced by the magnetic dipolar interaction between the spin polarized electronic orbitals and the muon spin. Even though the interaction is best described with quantum mechanics, for the sake of simplicity, here we approximate the spin polarized electronic orbitals with classical dipoles centered at the nuclei of the magnetic atoms. This approximation is also implicit in the code and works rather well in many cases.

The dipolar field is given by

\[\mathbf{B_{\mathrm{dip}}^\prime} = \frac{\mu_0}{4 \pi} \sum _{i=1} ^N \left( -\frac{\mathbf{m}}{r^3 _i} + \frac{3 (\mathbf{m}_i \cdot \mathbf{r}_i)\mathbf{r}_i }{r^5 _i} \right)\]

where, from a quantum perspective, \(\mathbf{m}_i = -g_i \mu_\mathrm{B} \mathbf{J}_i\) and \(\mathbf{J}_i\) is the total angular momentum of the i-th atom. Finally, the radius \(\mathbf{r}_i\) is the distance between the muon and the i-th atom of N magnetic ions in the sample.

When the above sum is performed in real space, it is customary to select a spherical portion of the sample (smaller than a magnetic domain) centered at the muon site and subdivide \(\mathbf{B_{\mathrm{dip}}}\) in three terms:

\[\mathbf{B_{\mathrm{dip}}^\prime} = \mathbf{B_{\mathrm{dip}}} + \mathbf{B_{\mathrm{Lor}}} + \mathbf{B_{\mathrm{dem}}}\]

The first term originates from the magnetic moments inside the sphere of radius \(R_\mathrm{sphere}\), i.e.:

\[\mathbf{B_{\mathrm{dip}}} = \frac{\mu_0}{4 \pi} \sum _{r_i<R_\mathrm{sphere}} \left( -\frac{\mathbf{m}}{r^3 _i} + \frac{3 (\mathbf{m}_i \cdot \mathbf{r}_i)\mathbf{r}_i }{r^5 _i} \right)\]

The second and the term originate from magnetic moments outside the sphere and are evaluated in the continuum approximation. They are

\[\mathbf{B_{\mathrm{Lor}}} = \frac{\mu_0}{3} \mathbf{M}_{\mathrm{Lor}} = \frac{\mu_0}{3 V_\mathrm{sphere}} \sum _{r_i < R_\mathrm{sphere}} \mathbf{m}_i\]
\[\mathbf{B_{\mathrm{dem}}} = - \mu_0 \mathbf{N} \mathbf{M}_\mathrm{meas}\]

where \(\mathbf{N}\) is the demagnetization tensor and \(\mathbf{M}_\mathrm{meas}\) is the bulk magnetization of the sample.

Note

muesr only estimates \(\mathbf{B}_\mathrm{dip}\) and \(\mathbf{B}_\mathrm{Lor}\). The demagnetisation field depends on both the sample shape and the experiment conditions and it must be evaluated case by case.

Contact Hyperfine field

A distinct contribution to the local magnetic field at the muon site is referred to as Fermi contact hyperfine field. It accounts for the finite probability for the quantum electron with wavefunction \(\psi_s (\mathbf{r})\) to share the classical muon position \(\mathbf{r}_\mu\) and it amounts to

\[\mathbf{B_{\mathrm{cont}}} = \frac{2 \mu_0}{3} \vert \psi_s (\mathbf{r}_\mu) \vert ^2 \mathbf{m}_e ^s\]

In muesr, only a scalar coupling between \(\mathbf{B_{\mathrm{cont}}}\) and \(\mathbf{m}_e\) is allowed, proportional to \(\vert \psi_s (\mathbf{r}_\mu) \vert ^2\).

In principle the quantum nature of the contact hyperfine interaction requires the knowledge of the electronic distribution around the muon site for an accurate description. Each magnetic atom that is a muon neighbor may contribute with a different coupling value, while the current version of muesr allows for just one average value. Furthermore the coupling may produce contributions from one or more neighbor magnetic atoms. They add up differently, depending to the magnetic structure.

At the moment this is only implemented by varying the number of nearest neighbours considered in the above sum. It is very badly approximated by considerig that each magnetic atom within a given radius contributes to the total hyperfine field by an amount inversely proportional to the cube of its distance from the muon. The total is then scaled by the common factor ACont.

[TODO]

Improve the implementation of an effective contact interaction in muesr!!!!

Description of Magnetic Structures

There are two possibilities to describe a magnetic structure: by using the color (Shubnikov) group theory or by defining one (or more) propagation vector(s) and using the Fourier coefficients formalism. muesr opts for the latter, limited to single wavevector (1-k) structures for the time being. A magnetic structure is defined as

\[\mathbf{\mu_{n \nu}} = \sum _{\mathbf{k}} \mathbf{m}_{\nu \mathbf{k}} e ^{- 2 \pi i \mathbf{k} \cdot \mathbf{R}_n}\]

where \(\nu\) runs over the atoms of the unit cell and \(n\) identifies the n-th cell where atomic positions \(\mathbf{R}_{n\nu}\) are obtained according to

\[\mathbf{R}_{n\nu} = \mathbf{R}_{n} + \mathbf{r_\nu}\]

with \(\mathbf{R}_{n} = n_a \mathbf{a} + n_b \mathbf{b} + n_c \mathbf{c}\) and \(\mathbf{r}_\nu = x_\nu \mathbf{a} + y_\nu \mathbf{b} + z_\nu \mathbf{c}\).

The fourier coefficients \(m_{\nu \mathbf{k}}\) are three dimensional complex vectors. They are related to the irreducible representations of the so called “little groups” i.e. the subgroup of the crystallographic space group formed by the operators leaving invariant the propagation vector.

[TODO] Discuss the phase!

As we said muesr can only handle 1-k magnetic structures. However, since local field are linear in the magnetic moment, the results for multiple-k magnetic orders can be obtained by performing multiple simulations for each of the k vectors and Fourier components which describe the system and summing the results.

Implementation details

muesr is a tool to analyze muon sites and local field contributions generated by a known magnetic structure. It is intended to be used in an interactive python environment such as IPython or Jupyter notebooks.

Internally, muesr uses Tesla units and Angstrom for lengths if not specified. Magnetic moments are specified in units of Bohr magnetons.