Skip to content

Commit dc32770

Browse files
author
davidcorteso
committed
Micro DMI: Specified interf DMI convention
1 parent 5ba1b49 commit dc32770

File tree

2 files changed

+9
-4
lines changed

2 files changed

+9
-4
lines changed

fidimag/micro/dmi.py

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -24,9 +24,10 @@ class DMI(Energy):
2424
r_ij is the vector from the i-th mesh site to the j-th
2525
neighbour
2626
27-
interfacial :: The energy density of this DMI is
27+
interfacial :: The energy density of this DMI is (according to
28+
Yang's convention; PRL 115, 267210)
2829
29-
w = D * ( L(x)_zx - L(y)_zy )
30+
w = D * ( L(x)_zx + L(y)_zy )
3031
3132
where L are Lifshitz invariants. This DMI is found in
3233
systems with their interface in contact with a heavy

fidimag/micro/lib/dmi.c

Lines changed: 6 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -235,17 +235,21 @@ void dmi_field_interfacial(double *m, double *field, double *energy, double *Ms_
235235
* For the Interfacial DMI, D_{ij} has the structure: D_{ij} = r_{ij} X z
236236
* See Rohart et al. Phys. Rev. B 88, 184422)
237237
*
238-
* but [Yang et al. Phys. Rev. Lett. 115, 267210] uses the opposite sign (?)
238+
* but [Yang et al. Phys. Rev. Lett. 115, 267210] uses the opposite sign.
239+
* We will use Yang's convention for the DMI sign.
239240
*
240-
* so the Dzyaloshinskii vectors are IN plane. This function only works
241+
* The Dzyaloshinskii vectors are IN plane. This function only works
241242
* along the XY plane since we assume there is a non magnetic material below
242243
* with a different SOC which gives the DMI for neighbouring in plane spins
243244
*
244245
* The computation is similar than before, only that we no longer have
245246
* a z component. In the continuum the H_dmi field is:
247+
*
246248
* H_dmi = (2 * D / (Ms a a_z)) * ( x X dM/dy - y X dM/dx )
249+
*
247250
* where "a" is the lattice spacing in the plane and "a_z" the system
248251
* thickness along z
252+
*
249253
* This derivative can be obtained doing the cross product
250254
* as in the atomic model, hence when going through the neighbours
251255
* loop, we can obtain the field doing the following cross products:

0 commit comments

Comments
 (0)