
2
Here, we investigate the SOT-driven self-oscillations
in a trilayer system consisting of a thin-film NCAFM
with a kagome structure sandwiched between two met-
als. The external electric field is applied perpendicu-
lar to the thin-film plane (see Fig. 1a). Surprisingly, we
find that the dynamics of the self-oscillations strongly de-
pend on the chirality set by the relativistic Dzyaloshin-
skii–Moriya interaction (DMI) of the system. Despite the
large in-plane and out-of-plane magnetic anisotropies, we
show that one of the two chiral structures hosts gap-
less self-oscillations that are highly tunable via intrinsic
SOTs. In contrast, the structure of opposite chirality
has gapped oscillations. Notably, the gapless oscillations
enable voltage-controlled NCAFM nano-oscillators with
exceptional bandwidths, where the frequency is tunable
from 0 Hz to the THz regime via the applied DC electric
field. Our results thus demonstrate that the NCAFMs
offer distinct chiral magnetic properties that are particu-
larly attractive for bridging the gap between technologies
operating in the microwave and infrared regions.
The material systems we consider are thin-film kagome
AFMs, where the mirror symmetry of the kagome lattice
plane is broken. These systems are described by the point
group D6[44]. Important candidate materials include
Mn3X (X= Ga, Ge, Sn), which in isolation are charac-
terized by the point group D6h[45], sandwiched between
two different metals. The broken spatial inversion sym-
metry of the system has two significant consequences: 1)
it leads to a magnetoelectric effect, and 2) it induces a
DMI. The main effect of the DMI is that it determines the
chirality of the ground state (see Fig 1b-c). The magneto-
electric effect refers to the out-of-equilibrium spin density
produced by electric fields [46], which in magnetic sys-
tems yields an SOT [47–50]. Below, we start by deriving
the magnetoelectric effect of NCAFMs with D6symme-
try from symmetry arguments [51]. Then, based on the
symmetry analysis, we phenomenologically add the cou-
pling terms between the spin system and electric field in
a microscopic model, which is used as starting point for
deriving an effective action and dissipation functional of a
uniform NCAFM. Further, the effective theory is applied
to investigate the voltage-controlled self-oscillations.
In linear response, the out-of-equilibrium spin density
sproduced by the electric field Eis given by [46]
si=ηij Ej.(1)
Here, ηij is a second-rank axial tensor, which satisfies the
following symmetry relationships [48, 49]
ηij =|G|Gii′Gjj′ηi′j′,(2)
dictated by the generators Gof the system’s point group.
|G|represents the determinant of the symmetry oper-
ation G. Throughout, we apply Einstein’s summation
convention for repeated indices. For kagome AFMs de-
scribed by the point group D6, the symmetry relations in
Eq. (2) imply that ηij is diagonal and parameterized by
two independent parameters [52]: ηxx =ηyy ≡η⊥and
ηzz ≡ηz. Here, the xand yaxes span the kagome plane,
whereas the z-axis is perpendicular to the lattice plane
(Fig. 1a). Consequently, the out-of-equilibrium spin den-
sity produced by the electric field can be written as
sx
sy
sz
=
η⊥0 0
0η⊥0
0 0 ηz
Ex
Ey
Ez
.(3)
Interestingly, we see that the electric field in kagome
AFMs can polarize the spin density along any axis (also
the out-of-plane axis z). This is different from most thin-
film systems, which usually are characterized by Dressel-
haus or Rashba SOC where the electric field only gener-
ates spin densities polarized along an in-plane axis of the
thin-film magnet [47–49]. In what follows, we investigate
how the spin density (3) couples to the NCAFM.
The kagome AFM is modeled by the spin Hamiltonian
H=He+Ha+HD+HE.(4)
Here, He=JP⟨ι˜ι⟩Sι·S˜ιdescribes the isotropic
exchange interaction (J > 0) between the neighbor-
ing lattice sites ⟨ι˜ι⟩, whereas Ha=Pι[Kz(Sι·ˆ
z)2−
K(Sι·ˆ
nι)2] represents the easy axes (K > 0) and easy
plane (Kz>0) anisotropy energies. The unit vector
ˆ
nιdenotes the in-plane easy axis at lattice site ι. The
kagome AFM consists of three spin sublattices with in-
plane easy axes ˆ
n1= [0,1,0], ˆ
n2= [√3/2,−1/2,0], and
ˆ
n3= [−√3/2,−1/2,0], respectively (Fig. 1b-c). HD=
P⟨ι˜ι⟩Dι˜ι·(Sι×S˜ι) is the DMI where Dι˜ι=Dzˆ
z[53].
HE=−PιgrSι·ηEexpresses the reactive coupling to
the electric field, where gris the coupling strength.
The ground state of the spin Hamiltonian (4) depends
on the ratio Dz/K. If Dz/K < 1/4√3, the spins are
aligned parallel or anti-parallel to the in-plane easy axes,
i.e., Sι=±ˆ
nι(see Fig. 1b). We will refer to these
two ground states as (+)-chiral. On the other hand,
if Dz/K > 1/4√3, the spins attain a configuration of
opposite chirality, which we will refer to as having (−)-
chirality (Fig. 1c). The (−)-chiral configuration is related
to (+)-chiral structure by a reflection about the xz-plane.
The dynamics of the spin system is described by the ac-
tion S=PιℏRdtA(Sι)·˙
Sι−RdtHand the dissipation
functional G=PιℏRdt[(αG/2) ˙
S2
ι+gd˙
Sι·(ηE×Sι)] [54–
58]. Here, ˙
Sι≡∂tSι,Ais defined via ∇×A(Sι) =
Sι/S,αGis the Gilbert damping parameter, and the
term proportional to gdcharacterizes the dissipative cou-
pling to the current-induced spin density. To derive an
effective description of the dynamics, it is convenient to
express the three sublattice spins as [54]
Sι(t) = SR(t) [ˆ
nι+aL(t)]
∥ˆ
nι+aL(t)∥, ι ∈ {1,2,3}.(5)