Model

Please see our paper for details about the magnetic model implemented by flatspin:

Jensen et al. “flatspin: A Large-Scale Artificial Spin Ice Simulator”

A short summary is provided below, as well as some additional practical information.

Magnets as dipoles

As illustrated in the above figure, magnets are modelled as magnetic dipoles with binary magnetization, i.e., macro spins:

\[s_i \in \{-1, +1\}\]

Each magnetic dipole has an associated position \(\vec{r}_i\) and rotation \(\theta_i\). Using reduced units, the magnetization vector of a single magnet can be expressed as

\[\vec{m}_i=s_i \vec{\hat{m}}_i\]

where \(\vec{\hat{m}}_i\) is the unit vector along \(\vec{m}_i\).

Magnetic Fields and Temperature

External fields and temperature are modeled as a combination of effective magnetic fields. The total field, \(\vec{h}_i\), affecting each magnet \(i\) is the sum of three components:

\[\vec{h}_i = \vec{h}_{dip}^{(i)} + \vec{h}_{ext}^{(i)} + \vec{h}_{th}^{(i)}\]

where \(\vec{h}_\text{dip}^{(i)}\) is the local magnetic field from neighboring magnets (magnetic dipole-dipole interactions), \(\vec{h}_\text{ext}^{(i)}\) is a global or local external field, and \(\vec{h}_\text{th}^{(i)}\) is a random magnetic field representing thermal fluctuations in each magnetic element.

Magnetic dipole-dipole interactions

Spins are coupled through dipole-dipole interactions. Each spin, \(i\), is subject to a magnetic field from all neighboring spins, \(j\neq i\), given by

\[\vec{h}_\text{dip}^{(i)} = \alpha \sum_{j \ne i}\frac{3\vec{r}_{ij}(\vec{m}_j \cdot \vec{r}_{ij})}{|\vec{r}_{ij}|^5} - \frac{\vec{m}_j}{|\vec{r}_{ij}|^3},\]

where \(\vec{r}_{ij}=\vec{p}_i-\vec{p}_j\) is the distance vector from spin \(i\) to \(j\) and \(\alpha\) scales the dipolar coupling strength between spins. The coupling strength \(\alpha\) is given by \(\alpha = \frac{\mu_0 M}{4\pi a^3}\), where \(a\) is the lattice spacing, \(M\) is the absolute magnetic moment of a single magnet, and \(\mu_0\) is the vacuum permeability. The distance \(\vec{r}_{ij}\) is given in reduced units of the lattice spacing. The size of the neighborhood is user-configurable and defined in units of the lattice spacing.

External field

Applying an external magnetic field is the primary mechanism for altering the state of an ASI in a controlled manner. The external field can either be set locally on a per-spin basis, \(\vec{h}_\text{ext}^{(i)}\), globally for the entire system, \(\vec{h}_\text{ext}\), or as a spatial vector field, \(\vec{h}_\text{ext}(\vec{r})\).

Time-dependent external fields are supported, i.e., \(\vec{h}_\text{ext}\) is a discrete time series of either local, global or spatial fields. A variety of time-dependent external fields are provided, including sinusoidal, sawtooth and rotational fields. More complex field-protocols can be generated, e.g., for annealing purposes or probing dynamic response.

Thermal field

Thermal fluctuations are modeled as an additional local field, \(\vec{h}_\text{th}^{(i)}\), applied to each magnet individually. Two orthogonal components of the field are independently drawn from the Normal distribution \(\mathcal{N}(0,\sigma_\text{th}^2)\). The simulated temperature \(T\) is closely related to the value of \(\sigma_\text{th}\), where a large \(\sigma_\text{th}\) corresponds to higher temperatures.

Switching

Magnetization reversal, or switching, may take place when a magnet is subjected to a magnetic field or high temperature. If the field is sufficiently strong (stronger than some critical field) and directed so that the projection onto \(\vec{m}_i\) is in the opposite direction to \(\vec{m}_i\), the magnetization \(\vec{m}_i\) will switch direction.

The critical field strength is referred to as the coercive field \(\vec{h}_\text{c}\), which depends on the angle between the applied field \(\vec{h}_i\) and \(\vec{m}_i\). The external field can be decomposed into two components, \(\vec{h}_\parallel\) and \(\vec{h}_\perp\), corresponding to the field component parallel and perpendicular to the easy axis, respectively. We denote the coercive field strength along the hard axis as \(h_k\).

flatspin uses a generalized Stoner-Wohlfarth (SW) switching model to allow an asymmetry between easy and hard axes. The critical field is described by the equation:

\[\left(\frac{h_{\parallel}}{b h_k}\right)^{2/\gamma} + \left(\frac{h_{\perp}}{c h_k}\right)^{2/\beta} = 1\]

\(b\), \(c\), \(\beta\) and \(\gamma\) are parameters which adjust the shape of the switching astroid: \(b\) and \(c\) define the height and width, respectively, while \(\beta\) and \(\gamma\) adjust the curvature of the astroid at the easy and hard axis, respectively.

In flatspin, the generalized SW model is used as the switching criteria, i.e., a spin may flip if the left hand side of the above equation is greater than one. Additionaly, the projection of \(\vec{h}_\text{i}\) onto \(\vec{m}_i\) must be in the opposite direction of \(\vec{m}_i\):

\[\vec{h}_i \cdot \vec{m}_i < 0.\]

Imperfections and disorder

Variations in each magnet are modelled as a disorder in the coercive fields. A user-defined parameter, \(k_\text{disorder}\), defines the distribution of coercive fields, i.e., \(h_{k}^{(i)}\) is sampled from a normal distribution \(\mathcal{N}(h_{k}, \sigma)\), where \(\sigma = k_\text{disorder} \cdot h_k\).

Dynamics

flatspin employs deterministic single spin flip dynamics. At each simulation step, we calculate the total magnetic field, \(\vec{h}_i\), acting on each spin. Next, we determine which spins may flip according to the switching criteria above. Finally we flip the spin where \(\vec{h}_i\) is the furthest outside its switching astroid, i.e., where the left hand side of the SW equation is greatest. Ties are broken in a deterministic, arbitrary manner, although with non-zero disorder such occurrences are rare. The above process is repeated until there are no more flippable spins.

This relaxation process is performed with constant external and thermal fields. To advance the simulation, the fields are updated and relaxation is performed again. Hence a simulation run consists of a sequence of field updates and relaxation processes.

Geometries

The particular spatial arrangement of the magnets is referred to as the geometry. flatspin includes the most commonly used ASI geometries: square, kagome, pinwheel and ising.

flatspin can be used to model any two-dimensional ASI comprised of identical elements. New geometries can easily be added by extending the model with a new set of positions \(\vec{r}_i\) and rotations \(\theta_i\).

Parameter and variable listing

The following table lists model parameters with their equivalents in the Python code:

Parameter	Python parameter	Description
\(\alpha\)	alpha	Dipolar coupling strength
\(h_k\)	hc	Mean coercive field (along hard axis)
\(k_\text{disorder}\)	disorder	Disorder in coercive fields
\(b\)	sw_b	Height of switching astroid (easy axis)
\(c\)	sw_c	Width of switching astroid (hard axis)
\(\beta\)	sw_beta	Curvature of switching astroid (easy axis)
\(\gamma\)	sw_gamma	Curvature of switching astroid (hard axis)
\(\sigma_\text{th}\)	thermal_std	Thermal field standard deviation

Additional parameters in the Python code include:

Python parameter	Description
size	Size of ASI (geometry specific)
lattice_spacing	Spacing between each spin (in reduced units)
neighbor_distance	Neighborhood to consider when calculating dipole interactions

The following table lists model variables with their equivalents in the Python code:

Variable	Python name	Description
\(s_i\)	spin[i]	Magnet spin
\(\vec{r}_i\)	pos[i]	Magnet position
\(\theta_i\)	angle[i]	Magnet angle
\(\vec{m}_i\)	vectors[i]	Magnet magnetization
\(h_k^{(i)}\)	hc[i]	Magnet coercive field (along hard axis)
\(\vec{h}_i\)	total_field(i)	Total magnetic field
\(\vec{h}_{dip}^{(i)}\)	dipolar_field(i)	Magnetic field from neighboring magnets
\(\vec{h}_{ext}^{(i)}\)	h_ext[i]	External magnetic field
\(\vec{h}_{th}^{(i)}\)	h_therm[i]	Thermal field

Choosing parameters

Here we give some advice on how to set some of the parameters.

Coupling parameters \(\alpha\) and lattice spacing

Although alpha defines the lattice spacing between magnets, it is also possible to change the lattice spacing via the lattice_spacing parameter. This has the effect of scaling the magnetic coupling as if changing lattice spacing directly. A neat trick is to calculate alpha based on a lattice spacing of 1 nm, set lattice_spacing=L to obtain an effective lattice spacing of L nm. Hence lattice_spacing allows setting the lattice spacing directly, e.g., for sweeps, as opposed to indirectly through alpha.

Switching parameters \(h_k\), \(b\), \(c\), \(\beta\) and \(\gamma\)

These parameters define the switching characteristic of the magnets. They can be estimated from micromagnetic simulations of a single magnet, when subject to a gradually increasing external field at different angles.

Here we provide a list some example magnet sizes corresponding switching parameters:

Magnet shape	Magnet size	\(h_k\)	\(b\)	\(c\)	\(\beta\)	\(\gamma\)
Elliptical	220x80x10 nm	0.2	0.77	1.0	2.2	4.2
Stadium	220x80x25 nm	0.2	0.38	1.0	1.3	3.6
Stadium	220x80x20 nm	0.2	0.41	1.0	1.5	3.9
Stadium	470x170x10 nm	0.14	0.28	1.0	4.8	3.0