Core equations¶
Landau-Lifshitz-Gilbert (LLG) equation¶
The dynamic of magnetic moment \(\vec{\mu}_i\) is governed by LLG equation,
where \(\vec{S}_i\) is the unit vector of magnetic moment,
and \(\vec{\mu}_s = |\vec{\mu}_i|\), the effective field \(\vec{H}_i\) is defined as
Interactions¶
In the level of atomic moments, \(\vec{\mu}_s\), originated from the angular momentum of electrons in atoms, could be employed as the base unit in simulations. There are several typical interactions between magnetic moments. The total Hamiltonian is the summation of them.
Exchange interaction¶
The classical Heisenberg Hamiltonian with the nearest-neighbor exchange interaction,
where the summation is taken only once for each pair, so the effective field is
In the continuum limit the exchange energy could be written,
so the corresponding effective field is
Once we implemented the Heisenberg exchange interaction, the effective field in the continuum case can be computed by the same codes with
Note that we needs the factor of \(\mu_0\) to convert the units from T to A/m.
Anisotropy¶
The Hamiltonian for uniaxial anisotropy with easy axis in x direction is expressed as,
and the corresponding field is
UniaxialAnisotropy¶
The UniaxialAnisotropy energy of the magnetic system is defined as
the the effective field is,
Dipolar interaction¶
The Hamiltonian for dipolar interaction is,
Dzyaloshinskii-Moriya interaction (DMI)¶
DMI is an antisymmetric, anisotropic exchange coupling beteen spins (magnetic moments),
Note that \(\vec{a}\cdot(\vec{b}\times\vec{c})=(\vec{a}\times\vec{b})\cdot\vec{c}\), the effective field can be computed by
For bulk materials \(\vec{D}_{ij} = D \vec{r}_{ij}\) and for interfacial DMI one has \(\vec{D}_{ij} = D \vec{r}_{ij} \times \vec{e}_z\), in both cases the vector \(\vec{D}_{ij}\) such that \(\vec{D}_{ij}=-\vec{D}_{ji}\).
In the continuum limit the bulk DMI energy could be written,
where \(D_a = -D/a^2\) and the effective field is
For the interfacial case, the effective field thus becomes,
Compared with the effective field [PRB 88 184422]
we have \(D_a = D/a^2\), note that there is no negative sign for the interfacial case.
Zeeman energy¶
The zeeman energy is,
Basically, we will follow the above equations to write codes.