Spin-helix Larmor mode

Shahrzad Karimi, Carsten A. Ullrich, Irene D 'Amico, Florent Perez

► To cite this version:

Shahrzad Karimi, Carsten A. Ullrich, Irene D 'Amico, Florent Perez. Spin-helix Larmor mode. Scientific Reports, 2018, 8, pp.3470. 10.1038/s41598-018-21818-8. hal-01727534

HAL Id: hal-01727534

https://hal.sorbonne-universite.fr/hal-01727534

Submitted on 9 Mar 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Spin-helix Larmor mode

Shahrzad Karimi<fnref n="1" />, Carsten A. Ullrich<fnref n="1" />, Irene D'Amico² & Florent Perez³

A two-dimensional electron gas (2DEG) with equal-strength Rashba and Dresselhaus spin-orbit coupling sustains persistent helical spin-wave states, which have remarkably long lifetimes. In the presence of an in-plane magnetic field, there exist single-particle excitations that have the character of propagating helical spin waves. For magnon-like collective excitations, the spin-helix texture reemerges as a robust feature, giving rise to a decoupling of spin-orbit and electronic many-body effects. We prove that the resulting spin-flip wave dispersion is the same as in a magnetized 2DEG without spin-orbit coupling, apart from a shift by the spin-helix wave vector. The precessional mode about the persistent spin-helix state is shown to have an energy given by the bare Zeeman splitting, in analogy with Larmor's theorem. We also discuss ways to observe the spin-helix Larmor mode experimentally.

Received: 29 September 2017

Accepted: 1 February 2018

Published online: 22 February 2018

Spin-flip waves are collective excitations of magnetic systems<fnref n="1" />,²: rather than flipping individual magnetic moments, which causes a large exchange energy penalty, the periodic reversal of magnetic moments extends as a precessional wave over the entire system, which is energetically favorable. There has been recent interest in spin waves in ferromagnetic thin films as an information carrier, which constitutes the basis for magnonics³,⁴. Spin waves exist in systems with localized and itinerant magnetic moments. In the latter case, the precession of the interacting spins, the charge motion, and the spin-orbit coupling (SOC) due to inversion asymmetry are all interrelated and lead to novel phenomena. For example, chiral spin waves have been observed in asymmetric monolayers of iron⁵,⁶ and on the surface of topological insulators⁷, helical spin waves have been predicted in a two-dimensional electron gas (2DEG) subject to Rashba SOC⁸, and twisted spin waves have been predicted and observed in magnetized 2DEGs⁹.

The fundamental and practical aspects of spin waves in the presence of SOC have drawn interest recently in the context of spintronics¹⁰⁻¹³. SOC provides the conversion of charge based information into the spin wave¹⁴,¹⁵. However, the presence of both SOC and Coulomb interaction still poses interesting challenges, especially in the dynamical regime.

In this paper, we will study spin waves in a 2DEG in the presence of in-plane magnetization and SOC. This system exhibits a rich interplay between Coulomb many-body effects, Rashba and Dresselhaus SOC, applied magnetic field, and electron density, which we have studied earlier⁹,¹⁶⁻¹⁸. What we found is that the spin waves are modified by the SOC in a subtle manner: the spin waves get a boost of their group velocity whose magnitude and orientation depends on the crystallographic propagation direction in the quantum well plane. This interesting behavior of the spin waves can be understood via a transformation into a spin-orbit twisted reference frame; however, in general this only holds to lowest order in SOC⁹.

In this paper, we consider a very special case in which exact results can be proved to all orders in SOC, namely, the case of a persistent spin helix¹⁹⁻²⁷. The spin helix arises in a 2DEG in which the strengths of the Rashba and Dresselhaus fields, $\alpha$ and $\beta$, are equal (in this paper we only include contributions to the Dresselhaus effect which are linear in the wavevector, i.e., $\beta = \beta_1$; the cubic $\beta_3$ coupling coefficient²³,²⁷ is ignored). We here consider a 2DEG embedded in a zincblende quantum well grown along the [001] direction: SU(2) symmetry is then partially restored, and a helical spin texture can be sustained along the [110] direction. This result is robust against spin-independent disorder scattering and Coulomb interactions²⁰. The main experimental signature of this state is that spin packet excitations are protected from decoherence, leading to extraordinarily long lifetimes²²,²⁵.

If a magnetic field is applied in the plane of a 2DEG with $\alpha = \beta$, spin-packet excitations can sustain long-lived precessional motion²³. Furthermore, some interesting magnetoelectronic effects can occur in 2DEGs²⁸⁻³⁰ or quantum wires³¹,³² which reflect the special condition $\alpha = \beta$. However, to our knowledge the spin-wave dynamics under these circumstances, which involves Coulomb interactions between the electrons at finite magnetic field, has not been explicitly addressed before.

Our treatment of spin waves is based both on time-dependent density-functional theory (TDDFT) in the linear-response regime³³ and on an equation-of-motion approach featuring the full many-body Hamiltonian⁹,³⁴.

<footnote n="1">¹Department of Physics and Astronomy, University of Missouri, Columbia, MO, 65211, USA. ²Department of Physics, University of York, York, YO10 5DD, United Kingdom. ³Institut des Nanosciences de Paris, CNRS/Université Paris VI, Paris, 75005, France. Correspondence and requests for materials should be addressed to C.A.U. (email: ullrichc@missouri.edu)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Reference frames $\mathcal{R}'$ (black) and $\mathcal{R}$ (red) for the electronic states in a 2DEG with SOC and in-plane magnetic field **B**. The striped pattern along the $\varphi = 45^{\circ}$ direction indicates the persistent spin helix state of the 2DEG, with wave vector **Q**, which forms in the absence of **B** if the Rashba and Dresselhaus coupling strengths are equal.</figcaption></figure>

We derive the exact form of the spin-wave dispersion for systems with a spin-helix texture, and find that it is obtained from the dispersion without SOC by a simple wave vector shift. The spin-wave stiffness $S_{sw}$ remains unchanged.

The main result of this paper is that we identify an exact dynamical state of the 2DEG with $\alpha = \beta$, which can be characterized as a collective precession of the 2DEG about the spin-helix state. The precession occurs at the bare Zeeman frequency, and we therefore refer to it as the spin-helix Larmor mode. This mode will be characterized by its long lifetime, and we will discuss ways in which it could be experimentally observed.

## Results

**2DEG with spin-orbit coupling in a magnetic field: the spin helix.** We consider the electronic ground state of a 2DEG in an n-doped zincblende quantum well in the presence of an in-plane magnetic field and Rashba and Dresselhaus SOC. Since the magnetic field is parallel to the 2DEG, it only acts on the spin and there is no Landau level quantization (as long as the magnetic length $l_B = \sqrt{\hbar/eB}$ is not significantly smaller than the well width). We will use the effective-mass approximation and work in units where $\hbar = e^* = m^* = 1$, where $e^*$ and $m^*$ are the effective charge and mass, respectively.

Figure 1 defines two reference frames. The primed frame $\mathcal{R}'$ is fixed with respect to the quantum well: the $x'$-, $y'$- and $z'$-axes point along the crystallographic [100], [010], and [001] directions, respectively; the 2DEG is in the $x' - y'$ plane. The Rashba and Dresselhaus SOC fields will be introduced in $\mathcal{R}'$, but it will be convenient for the discussion of spin waves to work in a coordinate system $\mathcal{R}$ which is oriented such that its x and z axes lie in the quantum well plane, and the z-axis points along the in-plane magnetic field **B**. As shown in Fig. 1, the x-axis is at an angle $\varphi$ with respect to the $x'$-axis.

**Single-particle states.** Without magnetic field or SOC, the lowest electronic conduction subband in the quantum well is spin-degenerate. For simplicity, we will treat the electronic states as purely two-dimensional; however, the main results in this paper will not change qualitatively if one takes the finite well width into account. The magnetic field lifts the degeneracy and splits the lowest subband into two, which we shall denote by the index $p = \pm 1$. In the reference frame $\mathcal{R}'$, the associated single-particle states can be written as $\Phi'_{pk}(\mathbf{r}') = e^{ik \cdot \mathbf{r}} \Psi'_{pk}$, where $\mathbf{r}' = (x', y')$, $k = (k_{x'}, k_{y'})$, and $\Psi'_{pk}$ is a two-component spinor of the form

$$ \Psi'_{pk} = \begin{pmatrix} \psi'_{pk\uparrow} \\ \psi'_{pk\downarrow} \end{pmatrix}. \qquad (1) $$

Here, “spin-up” and “spin-down” ($\uparrow$ and $\downarrow$) refer to the spin quantization axis $z'$.

The states $\Psi'_{pk}$ are obtained from the following Kohn-Sham single-particle equation:

$$ \left[ \frac{k^2}{2} \hat{\sigma}_0 + h_{x'} \hat{\sigma}_{x'} + h_{y'} \hat{\sigma}_{y'} \right] \Psi'_{pk} = E_{pk} \Psi'_{pk}, \qquad (2) $$

where $\hat{\sigma}_{0,x',y',z'}$ are the usual Pauli matrices. The off-diagonal parts in Eq. (2) involve

$$ h_{x'} = - \frac{Z + Z_{xc}}{2} \sin \varphi + \alpha k_{y'} + \beta k_{x'} \qquad (3) $$

$$ h_{y'} = \frac{Z + Z_{xc}}{2} \cos \varphi - \alpha k_{x'} - \beta k_{y'}. \qquad (4) $$

Here, $Z = g^* \mu_B B$ is the bare Zeeman energy ($\mu_B$ is the Bohr magneton, and $g^*$ is the effective g-factor). The presence of Coulomb many-body effects in the interacting 2DEG gives rise to the Zeeman exchange-correlation (xc) energy $Z_{xc}$, which we discuss below. $\alpha$ and $\beta$ are the usual Rashba and Dresselhaus linear coupling parameters;
we ignore contributions to the Dresselhaus effect that are cubic in the wavevector, since these tend to be much smaller than the linear contributions$^{23,27}$.

It is convenient to change the reference frame for the spin, and go over to reference system $\mathcal{R}$, whose z-axis is along the magnetic field direction. We introduce two in-plane vectors, $\mathbf{q}_0$ and $\mathbf{q}_1$, given in $\mathcal{R}'$ by

$$ \mathbf{q}_0 = (\alpha \cos \varphi + \beta \sin \varphi) \hat{e}_{x'} + (\alpha \sin \varphi + \beta \cos \varphi) \hat{e}_{y'} \quad (5) $$

$$ \mathbf{q}_1 = (\beta \cos \varphi - \alpha \sin \varphi) \hat{e}_{x'} + (\alpha \cos \varphi - \beta \sin \varphi) \hat{e}_{y'} \quad (6) $$

With this, Eq. (2) transforms into

$$ \left[ \frac{k^2}{2} \hat{\sigma}_0 + \left( \frac{Z + Z_{xc}}{2} - \mathbf{k} \cdot \mathbf{q}_0 \right) \hat{\sigma}_z + \mathbf{k} \cdot \mathbf{q}_1 \hat{\sigma}_x \right] \Psi_{pk} = E_{pk} \Psi_{pk}, \quad (7) $$

where the scalar products $\mathbf{k} \cdot \mathbf{q}_0$ and $\mathbf{k} \cdot \mathbf{q}_1$ remain invariant under $\mathcal{R}' \to \mathcal{R}$. $\Psi_{pk}$ is now a two-component spinor whose spatial coordinates and spin quantization axes are defined with respect to $\mathcal{R}$.

Let us now discuss the xc contribution. The in-plane magnetic field causes the 2DEG to become uniformly magnetized. The xc energy per particle of a homogeneous 2DEG$^{35}$, $e_{xc}(n, \zeta)$, can be written as a functional of the density $n$ and the spin polarization $\zeta$, where $n = n_↑ + n_↓$ and $\zeta = (n_↑ - n_↓)/n$ ($↑$ and $↓$ are now defined with respect to $\hat{e}_z$). The Zeeman xc energy is then given by

$$ Z_{xc} = 2 \frac{\partial e_{xc}}{\partial \zeta}. \quad (8) $$

The renormalized Zeeman energy$^{36}$ can now be defined as $Z^* = Z + Z_{xc}$.

The general solution of Eq. (7) has been considered elsewhere$^{37}$; instead, we concentrate here on the special case $\alpha = \beta$ and $\varphi = \pi/4$ or $5\pi/4$. Under these circumstances, $\mathbf{q}_1 = 0$ and Eq. (7) simplifies considerably:

$$ \left[ \frac{k^2}{2} \hat{\sigma}_0 + \frac{Z^* - \mathbf{k} \cdot \mathbf{Q}}{2} \hat{\sigma}_z \right] \Psi_{pk} = E_{pk} \Psi_{pk}, \quad (9) $$

where

$$ \mathbf{Q} = \pm 4\alpha \hat{e}_x, \quad \varphi = \begin{cases} \pi/4 \\ 5\pi/4 \end{cases}. \quad (10) $$

To keep the discussion a bit simpler, we will limit ourselves to the case $\varphi = \pi/4$ in the following, so $\mathbf{Q}=4\alpha\hat{e}_x$ (the $\varphi = 5\pi/4$ case is essentially the same, just in the opposite direction).

The solution of Eq. (9) is straightforward. We obtain

$$ \Psi_{+,k} = \begin{pmatrix} 1 \\ 0 \end{pmatrix}, \quad \Psi_{-,k} = \begin{pmatrix} 0 \\ 1 \end{pmatrix}, \quad (11) $$

where

$$ E_{\pm, k} = \frac{k^2}{2} \pm \frac{Z^* - \mathbf{k} \cdot \mathbf{Q}}{2}. \quad (12) $$

The single-particle energies (12) have the important property

$$ E_{+,k+Q} = E_{-,k} + Z^*, \quad (13) $$

which is illustrated in Fig. 2b, using the parameters $Z^* = 0.0381$ and $\alpha = 0.05$. This large value of $\alpha$ (typical experimental values of $\alpha$ are about an order of magnitude smaller) was chosen for clarity of presentation.

**B=0: spin-helical single-particle eigenstates.** In the absence of external magnetic fields, i.e., for $Z^* = 0$ (see Fig. 2a), the degeneracy of the two energy branches $E_{+,k+Q}$ and $E_{-,k}$ gives rise to a persistent spin-helix state of the 2DEG, as illustrated by the stripe-like pattern in Fig. 1. This is easy to see: Due to the degeneracy, linear combinations of the eigenstates $\Psi_{-,k}$ and $\Psi_{+,k+Q}$ are also solutions of Eq. (9). The single-particle wave functions for wave vector $\mathbf{k}$ can therefore be written as

$$ \Phi_{\mathbf{k}}^{\pm}(\mathbf{r}) = a\begin{pmatrix} 1 \\ 0 \end{pmatrix} e^{i(\mathbf{k}+\mathbf{Q})\cdot\mathbf{r}} \pm b\begin{pmatrix} 0 \\ 1 \end{pmatrix} e^{i\mathbf{k}\cdot\mathbf{r}}, \quad (14) $$

where $|a|^2 + |b|^2 = 1$. From the associated spin-density matrix it is straightforward to determine the magnetization in the $\mathcal{R}$ frame. We obtain $m_z = |a|^2 - |b|^2$, and since the macroscopic magnetization must vanish (i.e., $m_z = 0$), this implies $|a| = |b|$. Writing $a, b = e^{i\phi_{a,b}/\sqrt{2}}$ and defining $\delta_{ab} = \phi_a - \phi_b$, we get

$$ m_x = \pm \cos(\mathbf{Q} \cdot \mathbf{r} + \delta_{ab}), \quad m_y = \mp \sin(\mathbf{Q} \cdot \mathbf{r} + \delta_{ab}). \quad (15) $$
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Single-particle energies $E_{-k}$ and $E_{+,k}$ for $\alpha = 0.05$ and $k$ along the [110] direction, see Eq. (13). (a) No magnetic field ($Z^* = 0$). Linear combinations of $E_{-k}$ and $E_{+,k+Q}$ states have spin helix texture, but these cancel out if summed over all occupied states below $E_F$. A persistent spin helix appears if a quasiparticle is injected at the Fermi surface, as shown. (b) Finite magnetic field ($Z^* = 0.0381$). Single-particle excitations across the Fermi energy with momentum transfer $Q$ (thick green arrow) give rise to propagating spin helices.</figcaption></figure>

Accordingly, the stripes in Fig. 1 indicate a periodic rotation of the electronic spin in and out of the plane, with a wave vector $Q = 4\alpha\hat{e}_x$ oriented at a 45° angle with respect to the x'-axis (i.e., the [110] direction). This is the spin helix pattern$^{20,27}$.

The magnetizations $m_x$ and $m_y$ associated with the states $\Phi_k^+$ and $\Phi_k^-$ cancel out; therefore, adding up all occupied spin-helix states below the Fermi energy $E_F$ gives zero. This means that the ground state of the N-electron system has no spin texture. The persistent spin helix pattern can be observed if additional quasiparticles are injected at the Fermi level, as shown in Fig. 2a. Such states will have a very long lifetime$^{20,22,23}$.

**B ≠ 0: spin-helical single-particle excitations.** For a finite magnetic field ($Z^* > 0$), the degeneracy of the two energy branches $E_{+,k+Q}$ and $E_{-k}$ is lifted, as shown in Fig. 2b. As a consequence, the spin-helix pattern (15) is not a property of the ground state anymore: instead, the spin helix becomes a nonequilibrium feature.

To see this, consider a single-particle excitation across the Fermi surface, with wave vector transfer $Q$, as illustrated by the thick green arrow in Fig. 2b. First-order perturbation theory tells us that the time-dependent wave function has the form

$$ \Phi_{k \to k+Q}(\mathbf{r}, t) = \gamma \begin{pmatrix} 1 \\ 0 \end{pmatrix} e^{i(k+Q) \cdot \mathbf{r}} e^{-i(E_{-k}+Z^*)t} + \begin{pmatrix} 0 \\ 1 \end{pmatrix} e^{i\mathbf{k} \cdot \mathbf{r}} e^{-iE_{-k}t}, \quad (16) $$

where $|\gamma| \ll 1$ and we made use of Eq. (13). In the $\mathcal{R}$ frame, the $x$ and $y$ components of the associated magnetization are given by

$$ m_x(t) = 2|\gamma| \cos(\mathbf{Q} \cdot \mathbf{r} - Z^*t + \phi_\gamma) \quad (17) $$

$$ m_y(t) = -2|\gamma| \sin(\mathbf{Q} \cdot \mathbf{r} - Z^*t + \phi_\gamma), \quad (18) $$

where $\phi_\gamma$ is the phase of $\gamma$. This defines a forward propagating spin helix (i.e., a spin-flip wave), with amplitude $2|\gamma|$, wave vector $Q$, and group velocity $Z^*Q/Q^2$. Single-particle excitations of this type, which are long-lived due to the property (13), were experimentally observed by Walser et al.$^{23,38}$ using circularly polarized optical pump pulses and imaging via time-resolved Kerr rotation microscopy. A theoretical study using a drift-diffusion model was recently carried out by Ferreira et al.$^{39}$.

So far, our discussion has been for single-particle excitations, i.e., we did not consider collective excitations. In the following Section, we will present a very special case of a collective mode, which we call the spin-helix Larmor mode, which can be viewed as a coherent superposition of the left- and right-propagating single-particle spin helices considered above. As we will see, this gives rise to a collective, standing precessional wave with wave vector $Q$ which is undamped. Propagating collective spin-flip waves and their wave vector dispersions will be considered in Methods.

## Spin-helix Larmor mode

We consider a 2DEG in the presence of a uniform magnetic field $B = B\hat{e}_z$, in the reference frame $\mathcal{R}$ of Fig. 1. The many-body Hamiltonian without SOC is

$$ \hat{H}_0 = \sum_i \left( \frac{\hat{\mathbf{p}}_i^2}{2} + \frac{1}{2} Z \hat{\sigma}_{z,i} \right) + \frac{1}{2} \sum_{ij} \frac{1}{|\mathbf{r}_i - \mathbf{r}_j|} \quad (19) $$

We consider the spin-wave operator$^{9,34,40,41}$
$$ \hat{S}_{\mathbf{q}}^{+} = \frac{1}{2} \sum_{i} \hat{\sigma}_{+,i} e^{i\mathbf{q} \cdot \mathbf{r}_{i}} \qquad (20) $$

($\hat{\sigma}_{+} = \hat{\sigma}_{x} + i\hat{\sigma}_{y}$), whose Heisenberg equation of motion, in the absence of SOC, is

$$ i \frac{d}{dt} \hat{S}_{\mathbf{q}}^{+} = [\hat{S}_{\mathbf{q}}^{+}, \hat{H}_{0}] = -\omega_{L} \hat{S}_{\mathbf{q}}^{+} - \mathbf{q} \cdot \mathbf{J}_{\mathbf{q}}^{+}. \qquad (21) $$

The second term on the right-hand side of Eq. (21) arises from the commutator of $\hat{S}_{\mathbf{q}}^{+}$ with the kinetic part of $\hat{H}_{0}$, where

$$ \mathbf{J}_{\mathbf{q}}^{+} = \frac{1}{2} \sum_{i} \hat{\sigma}_{+,i} e^{i\mathbf{q} \cdot \mathbf{r}_{i}} (\mathbf{p}_{i} + \mathbf{q}/2) \qquad (22) $$

is the transverse spin-current operator at wave vector $\mathbf{q}$. The Larmor frequency $\omega_{L}$ is equal to the bare Zeeman energy, $\omega_{L} = Z$. For small values of $\mathbf{q}$, the equation of motion (21) can be written as

$$ i \frac{d}{dt} \hat{S}_{\mathbf{q}}^{+} = -\omega_{sw,0}(\mathbf{q}) \hat{S}_{\mathbf{q}}^{+}. \qquad (23) $$

Here, $\omega_{sw,0}(\mathbf{q})$ implicitly contains the coupling between the collective spin waves and the single-particle spin-current dynamics$^{34}$. The real part, $\Re\omega_{sw,0}(\mathbf{q}) = \omega_{L} + S_{sw}q^2/2$, is the spin-wave dispersion, where $S_{sw}$ is the spin-wave stiffness (explicit expressions and discussion of $S_{sw}$ are provided in the Methods section). The imaginary part of $\omega_{sw,0}(\mathbf{q})$ accounts for the damping of the spin wave due to electron-electron interactions$^{42}$.

Now let us include SOC. The SOC Hamiltonian for the spin-helix case is given by [see Eq. (9)]

$$ \hat{H}_{\text{SOC}} = -\frac{1}{2} \sum_{i} \mathbf{Q} \cdot \mathbf{p}_{i} \hat{\sigma}_{z,i} \qquad (24) $$

and the total Hamiltonian of the system is $\hat{H} = \hat{H}_{0} + \hat{H}_{\text{SOC}}$. The equation of motion for $\hat{S}_{\mathbf{q}}^{+}$ in the presence of SOC is given by

$$ i \frac{d}{dt} \hat{S}_{\mathbf{q}}^{+} = [\hat{S}_{\mathbf{q}}^{+}, \hat{H}_{0} + \hat{H}_{\text{SOC}}] = - \omega_{L} \hat{S}_{\mathbf{q}}^{+} - (\mathbf{q} - \mathbf{Q}) \cdot \mathbf{J}_{\mathbf{q}}^{+}. \qquad (25) $$

We now show that the SOC contribution can be transformed away. We introduce the SU(2) unitary transformation $\hat{U} = \exp[-i\sum_{i}\mathbf{Q} \cdot \mathbf{r}_{i}\hat{\sigma}_{z,i}/2]$, which leads to

$$ \hat{U} \hat{S}_{\mathbf{q}}^{+} \hat{U}^{\dagger} = \hat{S}_{\mathbf{q}-\mathbf{Q}}^{+}, \quad \hat{U} \hat{J}_{\mathbf{q}}^{+} \hat{U}^{\dagger} = \hat{J}_{\mathbf{q}-\mathbf{Q}}^{+}. \qquad (26) $$

In other words, $\hat{U}$ causes a boost of the wave vector argument of the spinwave and the spin-current operators, $\mathbf{q} \rightarrow \mathbf{q} - \mathbf{Q}$, and transforms the momentum operator of the $i$th electron into $\hat{U}\mathbf{p}_{i}\hat{U}^{\dagger} = \mathbf{p}_{i} + \mathbf{Q}\hat{\sigma}_{z,i}/2$. On the other hand, $\hat{U}$ leaves the Coulomb and the Zeeman parts of $\hat{H}_{0}$ unchanged. Thus, we obtain

$$ \hat{U} \hat{H} \hat{U}^{\dagger} = \hat{H}_{0} - \frac{Q^2}{8} \hat{\sigma}_{0}. \qquad (27) $$

Since $Q^2\hat{\sigma}_{0}/8$ commutes with $\hat{S}_{\mathbf{q}-\mathbf{Q}}^{+}$, the transformed equation of motion (25) becomes

$$ i \frac{d}{dt} \hat{S}_{\mathbf{q}}^{+} = \hat{U}^{\dagger} [\hat{S}_{\mathbf{q}-\mathbf{Q}}^{+}, \hat{H}_{0}] \hat{U}, \qquad (28) $$

where

$$ [\hat{S}_{\mathbf{q}-\mathbf{Q}}^{+}, \hat{H}_{0}] = -\omega_{L} \hat{S}_{\mathbf{q}-\mathbf{Q}}^{+} - (\mathbf{q} - \mathbf{Q}) \cdot \mathbf{J}_{\mathbf{q}-\mathbf{Q}}^{+}. \qquad (29) $$

For small values of $|\mathbf{q} - \mathbf{Q}|$, this becomes $[\hat{S}_{\mathbf{q}-\mathbf{Q}}^{+}, \hat{H}_{0}] = -\omega_{sw,0}(\mathbf{q} - \mathbf{Q})\hat{S}_{\mathbf{q}-\mathbf{Q}}^{+}$. Substituting this into Eq. (28), we find that the spin waves of the system with $\alpha = \beta$ are those of the system without SOC (governed only by $\hat{H}_{0}$), but where the wave vector is shifted:

$$ \omega_{sw}(\mathbf{q}) = \omega_{sw,0}(\mathbf{q} - \mathbf{Q}). \qquad (30) $$

Let us now consider the important special case $\mathbf{q} = \mathbf{Q}$. It is known$^{43}$ that, in the absence of SOC, the spin-polarized electron system carries out a collective precessional motion with $\omega_{sw,0} = \omega_{L}$. The Larmor frequency $\omega_{L}$ is equal to the bare Zeeman energy, $\omega_{L} = Z$. The Larmor mode has infinite lifetime (zero line width), because it can be represented as a superposition of two exact eigenstates of the system Hamiltonian $\hat{H}_{0}$: the many-body ground state $|0\rangle_{0}$ (the subscript 0 indicates absence of SOC), with ground-state energy $E_{0}$, and the many-body eigenstate $\hat{S}_{0}^{+}|0\rangle_{0}$, with energy $E_{0} + Z$. An important feature of the Larmor's mode is that it does not
<figure><img src="image_5.png" /><figcaption>**Figure 3.** (a) Spin-wave dispersion $\omega_{sw,0}(q)$ (line) and single-particle spin-flip continuum (shaded area) without SOC. (b) Spin-wave dispersion $\omega_{sw}(q)$ and single-particle spin-flip continuum, plotted along [110], for $\alpha = \beta = 0.003$ ($Q = 0.012$). The 2DEG parameters are $r_s = 2$, $\zeta = -0.0762$, $Z = 0.02$, and $Z^* = 0.0381$ (all values are in atomic units). The inset shows the position of Larmor's mode in the wave vector plane ($q_x, q_y$).</figcaption></figure>

carry any spin current in the plane. This is obvious from the equation of motion (21) in the absence of SOC: for the Larmor's mode which occurs at $q = 0$, no current is induced by the homogenous precession.

Coming back to the case with SOC with the full Hamiltonian $\hat{H}$, we can now formulate the spin-helix Larmor theorem. If the spin wave has wave vector $Q$, commensurate with the spin-helix texture, all Coulomb many-body contributions drop out and the frequency is given by the bare Zeeman energy:

$$ \omega_{sw}(Q) = \omega_{sw,0}(0) = Z, \qquad (31) $$

which follows directly from Eq. (30). In the presence of SOC, the Larmor's mode occurs at $q = Q$ and is not a homogenous mode anymore; however, one still has the property that no spin-current is driven by the precession, since the total current term in Eqs (25) and (29) disappears when $q = Q$: this happens because SOC induces spin currents opposite to the spin currents induced by the motion. The spin wave then has vanishing group velocity, $\nabla_q \omega_{sw}(q)|_{q=Q} = 0$, which means that it is a standing wave. All electronic spins precess about their local orientation, given by the spin helix configuration, with the Larmor frequency $\omega_L = Z$. The spin-helix Larmor mode is a superposition of two many-body eigenstates of $\hat{H}$: the ground state $|0\rangle$ and the state $\hat{S}_Q^+|0\rangle$.

We illustrate the spin-wave dispersions with and without SOC in Fig. 3. The left panel shows $\omega_{sw,0}(q)$ (which is independent of the direction of $q$), and the right panel shows $\omega_{sw}(q)$ for $q$ parallel to $Q$, i.e., along the [110] direction. For this particular case, $\omega_{sw}(q)$ is simply obtained by a horizontal shift by $Q$ of $\omega_{sw,0}(q)$ (likewise for the particle-hole continua). The spin-wave dispersions plotted in Fig. 3 are obtained from the numerical solution of Eq. (42), see below. The small-wave vector expansion (44) is very close to the exact result.

## Experimental Schemes

Experimental observation of the spin-helix Larmor mode should be possible in specially designed doped magnetic semiconductor quantum well samples where the $\alpha = \beta$ condition is met. The spin-flip waves under an in-plane magnetic field can be detected using inelastic light scattering, similar to our earlier work$^{9,16-18}$. The Larmor-type precessional mode about the spin-helix state should then be recognizable by a significant narrowing of the linewidth.

Time-resolved Kerr rotation spectroscopy is another suitable technique to detect collective spin excitations in a 2DEG$^{44,45}$. In the experimental setup by Walser et al.$^{23}$ however, only the spin-helical single-particle excitations were observed; at the magnetic field strength of 1 T considered in the experiment, the collective modes were too short-lived to be seen since their separation from the spin-flip continuum was too small and the mode frequency was of the same order as the linewidth.

We also propose a device design which would allow one to excite Larmor's mode optically and probe it electronically. In the absence of SOC, the Larmor's mode is usually probed in a paramagnetic resonance setup, by a microwave (mw) magnetic field $b_{mw} = b \cos(\omega t) \hat{e}_y$ applied in the plane perpendicular to the quantizing field $B = B \hat{e}_z$. The corresponding coupling Hamiltonian is $\hat{H}_{mw} = -\frac{i}{2}g_e \mu_B b \cos(\omega t)(\hat{S}_0^+ - \hat{S}_0^-)$. The unitary transformation $\hat{U}$ maps the situation with SOC to the situation without SOC; hence, the coupling Hamiltonian becomes $\hat{U}\hat{H}_{mw}\hat{U}^\dagger = -\frac{i}{2}g_e \mu_B b \cos(\omega t)(\hat{S}_Q^+ - \hat{S}_Q^-)$. This corresponds to a coupling with a standing helical magnetic field in the plane $(x, \hat{y})$. Such a magnetic field can be generated by a device like the one sketched in Fig. 4. The idea is to deposit metal stripes on top of the sample, separated by a distance $\pi/Q$ (for typical values of $\alpha$, of order ~1 meV Å, this corresponds to a few $\mu m$). The stripes are aligned parallel to the quantizing magnetic field $B = B \hat{e}_z$, perpendicular to the [010] direction, and the spacing of the stripes is commensurate with the standing-wave spin-helix Larmor mode. The metal stripes follow the principle of photo-conductive antennas used in time-domain THz experiments$^{46}$. When an infrared optical pulse hits the sample, the photo-generated carriers in the substrate create a short conductive channel closing the biased circuit of the stripes. A current pulse with frequencies in the THz range will propagate through the stripes in alternating direction. This will, by induction, generate concentric oscillating magnetic fields compatible with the spin-helix pattern as depicted in the right panel of Fig. 4. The magnetic field exerts torques on the spin-polarized 2DEG underneath. If the current pulse spectrum contains the right frequency, $\omega = \omega_L$, this will trigger a helical standing spin wave which will persist after the end of the pulse.
<figure><img src="image_6.png" /><figcaption>Figure 4. Proposed experimental design for the direct excitation of the spin-helix Larmor mode. Left: photo-conductive antenna on top of the sample, converting an infrared optical pulse into a short current pulse. Right: close-up view of the metal stripes on top of the spin-polarized 2DEG. The currents (blue arrows) are in alternating directions in neighboring stripes; the induced magnetic fields (pink circles) trigger a standing spin wave in the 2DEG.</figcaption></figure>

Detection of the spin-helix Larmor mode, and measurement of its lifetime, should then be possible via the currents induced in the metal stripes from the stray magnetic fields associated with the standing spin wave. Such a detection technique is common in ferromagnetic magnonics$^{3,4}$.

In our semiconductor test-bed system$^{9}$, we deal with lower spin densities and higher frequencies: we expect stray fields of the order of nT, and rapidly varying, at a rate of about 50–100 GHz, with a lifetime in the 100 ps to ns range. The measurements will therefore push the limits of present-day electronics, but they should be feasible. Moreover, the physics and concepts presented here are applicable to emerging 2D conducting systems with higher doping and strong spin-orbit effects. The above estimate for the lifetime of the spin-helix Larmor mode is based on spin-wave linewidths of the order of 50 µeV which were measured$^{9}$ for quantum wells with α ≠ β; the corresponding dephasing time (165 ps for a linewidth of 50 µeV) can be viewed as a lower threshold to the lifetime of the Larmor mode, since it includes electronic many-body contributions which drop out if α = β. The lifetime of the spin-helix Larmor mode is mainly determined by cubic Dresselhaus contributions and disorder.

## Discussion

In this paper, we have considered the spin dynamics in a 2DEG in the presence of SOC, under the very special condition where the Rashba and Dresselhaus coupling strengths are equal (α = β) and where an in-plane magnetic field is applied perpendicular to the [110] direction. Without this magnetic field, the system sustains persistent spin-helix states which have been widely studied in the literature. The magnetic field lifts the degeneracy that leads to the persistent spin-helix states; it instead leads to single-particle excitations that have the form of propagating spin helices.

The presence of Coulomb interactions causes these single-particle excitations to combine and form collective spin waves, which are robust against any decoherence caused by SOC. We have found that for the 2DEG with α = β the spin-wave dispersion is the same as for the system without SOC, apart for a rigid wave vector shift by Q (the spin-helix wave vector). The case of q = Q thus produces the special scenario which we have termed the spin-helix Larmor mode, where all many-body effects vanish and the precession frequency is given by the bare Zeeman energy (divided by ħ). This is a new and exact result for electronic many-body systems, which opens up new ways of manipulating and driving electronic spins by optical means.

## Methods

We now discuss the spin-wave dispersions in the quantum well system considered above. We have seen that in the case α = β the single-particle states (11) are pure up and down spinors; therefore, the longitudinal and transverse spin response channels are decoupled. The transverse spin-density response equation reads

$$ \vec{n}_T(\mathbf{q}, \omega) = \chi_T(\mathbf{q}, \omega)[\vec{v}_T(\mathbf{q}, \omega) + \int_T^{xc}(\mathbf{q}, \omega)\vec{n}_T(\mathbf{q}, \omega)], \quad (32) $$

where $\vec{n}_T = (n_{↑↓}, n_{↓↑})$ is the transverse spin-density response, $\vec{v}_T = (v_{↑↓}, v_{↓↑})$ is an external perturbation, and the noninteracting transverse spin response function and the transverse xc kernel are diagonal 2 × 2 matrices in the frame $\mathcal{R}^{41,47}$:

$$ \begin{aligned} \chi_T(\mathbf{q}, \omega) &= \begin{pmatrix} \chi_{↑↓,↑↓}(\mathbf{q}, \omega) & 0 \\ 0 & \chi_{↓↑,↓↓}(\mathbf{q}, \omega) \end{pmatrix}, \\ \int_T^{xc}(\mathbf{q}, \omega) &= \begin{pmatrix} f_{↑↓,↑↓}^{xc}(\mathbf{q}, \omega) & 0 \\ 0 & f_{↓↑,↓↑}^{xc}(\mathbf{q}, \omega) \end{pmatrix}. \end{aligned} \quad (33) $$

The general form of the individual elements of the noninteracting spin-density-matrix response function for a 2DEG is$^{33}$
$$
\begin{equation}
\begin{aligned}
\chi_{\sigma\sigma',\tau\tau}(\mathbf{q}, \omega) = & -\sum_{pp'}^{\pm 1} \int \frac{d^2\mathbf{k}}{(2\pi)^2} \frac{f(E_{p\mathbf{k}})}{\omega - E_{p\mathbf{k}} + E_{p'\mathbf{k}-\mathbf{q}} + i\eta} \\
& \times [\delta_{p,+1}\delta_{\sigma\uparrow} + \delta_{p,-1}\delta_{\sigma\downarrow}][\delta_{p',+1}\delta_{\sigma'\uparrow} + \delta_{p',-1}\delta_{\sigma'\downarrow}] \\
& \times \delta_{\sigma\tau}\delta_{\sigma'\tau'} + \sum_{pp'}^{\pm 1} \int \frac{d^2\mathbf{k}}{(2\pi)^2} \frac{f(E_{p\mathbf{k}})}{\omega + E_{p\mathbf{k}} - E_{p'\mathbf{k}+\mathbf{q}} + i\eta} \\
& \times [\delta_{p',+1}\delta_{\sigma\uparrow} + \delta_{p',-1}\delta_{\sigma\downarrow}] \\
& \times [\delta_{p,+1}\delta_{\sigma'\uparrow} + \delta_{p,-1}\delta_{\sigma'\downarrow}]\delta_{\sigma\tau}\delta_{\sigma'\tau'}
\end{aligned}
\tag{34}
\end{equation}
$$

where $\sigma, \sigma', \tau, \tau'$ are spin indices (↑ or ↓), and $\eta$ is a positive infinitesimal (since we will be considering spin waves outside the particle-hole continuum, we can drop $\eta$). The Fermi function is given by $f(E_{p\mathbf{k}}) = \theta(E_F - E_{p\mathbf{k}})$, where $E_F$ is the Fermi energy of a paramagnetic 2DEG in the presence of SOC. It can be shown that $E_F = E_F^0 - \alpha^2 - \beta^2$, where $E_F^0 = \pi n$ is the Fermi energy of a 2DEG without SOC, and with 2D electronic density $n$.

We recast Eq. (12) as $E_{\pm,\mathbf{k}} = \pm Z^*/2 + \frac{1}{2}|\mathbf{k} \mp \mathbf{Q}|^2 - 2\alpha^2$. With a change of the integration variable, $\mathbf{k} \to \mathbf{k} \pm \mathbf{Q}/2$, the noninteracting spin-flip response functions become

$$
\begin{equation}
\begin{aligned}
\chi_{\uparrow\downarrow,\uparrow\downarrow}(\mathbf{q}, \omega) = & -\int \frac{d^2\mathbf{k}}{(2\pi)^2} \frac{f_0\left(\frac{k^2}{2} + \frac{Z^*}{2}\right)}{\omega - Z^* - \mathbf{k} \cdot (\mathbf{q} - \mathbf{Q}) + |\mathbf{q} - \mathbf{Q}|^2/2} \\
& + \int \frac{d^2\mathbf{k}}{(2\pi)^2} \frac{f_0\left(\frac{k^2}{2} - \frac{Z^*}{2}\right)}{\omega - Z^* - \mathbf{k} \cdot (\mathbf{q} - \mathbf{Q}) - |\mathbf{q} - \mathbf{Q}|^2/2}
\end{aligned}
\tag{35}
\end{equation}
$$

$$
\begin{equation}
\begin{aligned}
\chi_{\downarrow\uparrow,\downarrow\uparrow}(\mathbf{q}, \omega) = & \int \frac{d^2\mathbf{k}}{(2\pi)^2} \frac{f_0\left(\frac{k^2}{2} + \frac{Z^*}{2}\right)}{\omega + Z^* - \mathbf{k} \cdot (\mathbf{q} + \mathbf{Q}) - |\mathbf{q} + \mathbf{Q}|^2/2} \\
& - \int \frac{d^2\mathbf{k}}{(2\pi)^2} \frac{f_0\left(\frac{k^2}{2} - \frac{Z^*}{2}\right)}{\omega + Z^* - \mathbf{k} \cdot (\mathbf{q} + \mathbf{Q}) + |\mathbf{q} + \mathbf{Q}|^2/2},
\end{aligned}
\tag{36}
\end{equation}
$$

where the Fermi function $f_0$ indicates that $E_F$ has been replaced by $E_F^0$.

Looking at Eqs (35) and (36) we immediately see that the spin-flip response functions of the system with SOC can be expressed in terms of the corresponding functions without SOC (denoted by the superscript 0):

$$
\chi_{\uparrow\downarrow,\uparrow\downarrow}(\mathbf{q}, \omega) = \chi_{\uparrow\downarrow,\downarrow\downarrow}^{0}(\mathbf{q} - \mathbf{Q}, \omega) \quad (37)
$$

$$
\chi_{\downarrow\uparrow,\downarrow\uparrow}(\mathbf{q}, \omega) = \chi_{\downarrow\uparrow,\downarrow\downarrow}^{0}(\mathbf{q} + \mathbf{Q}, \omega). \quad (38)
$$

This simple result only holds for the special case $\alpha = \beta$. In TDDFT, the spin-flip linear-response xc kernel of the 2DEG is given in the adiabatic local-density approximation (ALDA) by$^{33}$

$$
f_{\uparrow\downarrow,\uparrow\downarrow}^{xc} = f_{\downarrow\uparrow,\downarrow\uparrow}^{xc} \equiv K_{xc} = \frac{2}{n\zeta} \frac{\partial e_{xc}}{\partial \zeta}. \qquad (39)
$$

The spin-flip wave dispersion now follows from Eq. (32) by setting $\vec{v}_T$ to zero and finding the eigenmodes, which leads to

$$
[K_{xc}\chi_{\uparrow\downarrow,\uparrow\downarrow}^{0}(\mathbf{q} - \mathbf{Q}, \omega_{sw}) - 1][K_{xc}\chi_{\downarrow\uparrow,\downarrow\uparrow}^{0}(\mathbf{q} + \mathbf{Q}, \omega_{sw}) - 1] = 0. \qquad (40)
$$

Upon closer inspection of Eqs (35) and (36), we find

$$
\chi_{\downarrow\uparrow,\downarrow\uparrow}^{0}(\mathbf{q} + \mathbf{Q}, \omega_{sw}) = \chi_{\uparrow\downarrow,\uparrow\downarrow}^{0}(\mathbf{q} - \mathbf{Q}, -\omega_{sw}). \qquad (41)
$$

We are only interested in positive frequencies, so the spin-flip wave dispersion is obtained from Eq. (40) as the implicit solution $\omega_{sw}(\mathbf{q})$ of

$$
\chi_{\uparrow\downarrow,\uparrow\downarrow}^{0}(\mathbf{q} - \mathbf{Q}, \omega_{sw}) = K_{xc}^{-1}. \qquad (42)
$$

The spin-wave dispersion can be analytically determined for small wave vectors. To second order, we find

$$
\omega_{sw}(\mathbf{q}) = Z^* \left( 1 + \frac{K_{xc}}{2\pi} \right) + \frac{E_F^0 |\mathbf{q} - \mathbf{Q}|^2}{Z^*} \left( 1 + \frac{2\pi}{K_{xc}} \right). \qquad (43)
$$
<figure><img src="image_7.png" /><figcaption>Figure 5. Spin-wave stiffness of the 2DEG, obtained with the ALDA, for various values of the spin polarization $\zeta$ between 0.05 and 0.95.</figcaption></figure>

This can be rewritten as

$$ \omega_{sw}(\mathbf{q}) = Z + \frac{1}{2}S_{sw}|\mathbf{q} - \mathbf{Q}|^2, \qquad (44) $$

which corresponds to spin waves propagating with group velocity $v_g(\mathbf{q}) = S_{sw}(\mathbf{q} - \mathbf{Q})$, where the spin-wave stiffness of the 2DEG is

$$ S_{sw} = \frac{2E_F^0}{Z^*} \left( 1 + \frac{2\pi}{K_{xc}} \right). \qquad (45) $$

Since $K_{xc} = Z_{xc}/(n\zeta)$ and $\zeta = -Z^*/(2\pi n)$, this can also be written as

$$ S_{sw} = - \frac{2E_F^0}{Z^*} \left( \frac{Z}{Z^* - Z} \right). \qquad (46) $$

We plot the ALDA spin-wave stiffness for various values of the spin polarization $\zeta$ and as a function of the 2D Wigner-Seitz radius $r_s$ in Fig. 5. The stiffness has negative values for all practically relevant values of $r_s$, and crosses over to $S_{sw} > 0$ for very low densities (at $r_s = 26.96$ for $\zeta = 0.05$ and at $r_s = 24.49$ for $\zeta = 0.95$). In practice, one is interested in studying collective modes in the typical metallic range of $r_s$ between about 2 and 6; this is the intermediate coupling regime where the Coulomb energy (which scales as $r_s^{-1}$) is comparable to the kinetic energy (which scales as $r_s^{-2}$), but which can still be regarded as weakly correlated. As an example, in the quantum well system considered in ref. 37 the spin-wave stiffness was $S_{sw} = -27.6$, for $r_s = 2.2$ and $\zeta = 0.053$.

**Data availability.** All data generated or analysed during this study are included in this article.

## References

1. Stancil, D. D. & Prabhakar, A. *Spin waves: Theory and applications*. (Springer, Berlin, 2010).

2. Zakeri, K. Elementary spin excitations in ultrathin itinerant magnets. *Phys. Rep.* **545**, 47–93 (2014).

3. Kruglyak, V. V., Demokritov, S. O. & Grundler, D. Magnonics. *J. Phys. D: Appl. Phys.* **43**, 264001–1–14 (2010).

4. Chumak, A. V., Vasyuchka, V. I., Serga, A. A. & Hillebrands, B. Magnon spintronics. *Nature Phys.* **11**, 453–61 (2015).

5. Zakeri, K., Zhang, Y., Chuang, T.-H. & Kirschner, J. Magnon lifetimes on the Fe(110) surface: The role of spin-orbit coupling. *Phys. Rev. Lett.* **108**, 197205–1–5 (2012).

6. Di, K. et al. Direct observation of the Dzyaloshinskii-Moriya interaction in a Pt/Co/Ni film. *Phys. Rev. Lett.* **114**, 047201–1–5 (2015).

7. Kung, H.-H. et al. Chiral spin mode on the surface of a topological insulator. *Phys. Rev. Lett.* **119**, 136802–1–6 (2017).

8. Ashrafi, A. & Maslov, D. L. Chiral spin waves in fermi liquids with spin-orbit coupling. *Phys. Rev. Lett.* **109**, 227201–1–5 (2012).

9. Perez, F. et al. Spin-orbit twisted spin waves: group velocity control. *Phys. Rev. Lett.* **117**, 137204–1–5 (2016).

10. Winkler, R. *Spin-orbit coupling effects in two-dimensional electron and hole systems*, vol. 191 of Springer Tracts in Modern Physics. (Springer, Berlin, 2003).

11. Žutić, I., Fabian, J. & Das Sarma, S. Spintronics: Fundamentals and applications. *Rev. Mod. Phys.* **76**, 323–410 (2004).

12. Wu, M. W., Jiang, J. H. & Weng, M. Q. Spin dynamics in semiconductors. *Phys. Rep.* **493**, 61–236 (2010).

13. Bercioux, D. & Lucignano, P. Quantum transport in Rashba spin-orbit materials: a review. *Rep. Prog. Phys.* **78**, 106001–1–31 (2015).

14. Liu, T. & Vignale, G. Electric control of spin currents and spin-wave logic. *Phys. Rev. Lett.* **106**, 247203–1–4 (2011).

15. Kajiwara, Y. et al. Transmission of electrical signals by spin-wave interconversion in a magnetic insulator. *Nature* **464**, 262–266 (2010).

16. Baboux, F. et al. Coulomb-driven organization and enhancement of spin-orbit fields in collective spin excitations. *Phys. Rev. B* **87**, 121303–1–5 (2013).

17. Baboux, F., Perez, F., Ullrich, C. A., Karczewski, G. & Wojtowicz, T. Electron density magnification of the collective spin-orbit field in quantum wells. *Phys. Rev. B* **92**, 125307–1–6 (2015).

18. Baboux, F., Perez, F., Ullrich, C. A., Karczewski, G. & Wojtowicz, T. Spin-orbit stiffness of the spin-polarized electron gas. *Phys. Stat. Solidi RLL* **10**, 315–9 (2016).

19. Schliemann, J., Egues, J. C. & Loss, D. Nonballistic spin-field-effect transistor. *Phys. Rev. Lett.* **90**, 146801–1–4 (2003).
20. Bernevig, B. A., Orenstein, J. & Zhang, S.-C. Exact SU(2) symmetry and persistent spin helix in a spin-orbit coupled system. *Phys. Rev. Lett.* **97**, 236601–1–4 (2006).

21. Liu, M.-H., Chen, K.-W., Chen, S.-H. & Chang, C.-R. Persistent spin helix in Rashba-Dresselhaus two-dimensional electron systems. *Phys. Rev. B* **74**, 235322–1–6 (2006).

22. Koralek, J. D. et al. Emergence of the persistent spin helix in semiconductor quantum wells. *Nature* **458**, 610–3 (2009).

23. Walser, M. P., Reichl, C., Wegscheider, W. & Salis, G. Direct mapping of the formation of a persistent spin helix. *Nature Phys.* **8**, 757–62 (2012).

24. Schönhuber, C. et al. Inelastic light-scattering from spin-density excitations in the regime of the persistent spin helix in a GaAs-AlGaAs quantum well. *Phys. Rev. B* **89**, 085406–1–6 (2014).

25. Sasaki, A. et al. Direct determination of spin-orbit interaction coefficients and realization of the persistent spin helix symmetry. *Nature Nanotech.* **9**, 703–9 (2014).

26. Kammermeier, M., Wenk, P. & Schliemann, J. Control of spin helix symmetry in semiconductor quantum wells by crystal orientation. *Phys. Rev. Lett.* **117**, 236801–1–5 (2016).

27. Schliemann, J. Persistent spin textures in semiconductor nanostructures. *Rev. Mod. Phys.* **89**, 011001–1–17 (2017).

28. Duckheim, M. & Loss, D. Resonant spin polarization and spin current in a two-dimensional electron gas. *Phys. Rev. B* **75**, 201305–1–4 (2007).

29. Nakhmedov, E. P. & Alekperov, O. Interplay of Rashba- and Dresselhaus spin-orbit interactions in a quasi-two-dimensional electron gas of a finite thickness under in-plane magnetic field. *Eur. Phys. J. B* **85**, 298–1–19 (2012).

30. Wilde, M. A. & Grundler, D. Alternative method for the quantitative determination of Rashba- and Dresselhaus spin-orbit interaction using the magnetization. *New J. Phys.* **15**, 115013–1–18 (2013).

31. Nazmitdinov, R. G., Pichugin, K. N. & Valin-Rodríguez, M. Spin control in semiconductor quantum wires: Rashba and Dresselhaus interaction. *Phys. Rev. B* **79**, 193303–1–4 (2009).

32. Micu, C. & Papp, E. Equal coupling strength approach to spin precession effects concerning Rashba and Dresselhaus spin-orbit interactions in the presence of in-plane magnetic fields. *Superlatt. Microstruct.* **51**, 651–662 (2012).

33. Ullrich, C. A. *Time-dependent density-functional theory: concepts and applications*. (Oxford University Press, Oxford, 2012).

34. Perez, F., Cibert, J., Vladimirova, M. & Scalbert, D. Spin waves in magnetic quantum wells with Coulomb interaction and sd exchange coupling. *Phys. Rev. B* **83**, 075311–1–9 (2011).

35. Attaccalite, C., Moroni, S., Gori-Giorgi, P. & Bachelet, G. B. Correlation energy and spin polarization in the 2D electron gas. *Phys. Rev. Lett.* **88**, 256601–1–4 (2002).

36. Perez, F. et al. From spin flip excitations to the spin susceptibility enhancement of a two-dimensional electron gas. *Phys. Rev. Lett.* **99**, 026403–1–4 (2007).

37. Karimi, S. et al. Spin precession and spin waves in a chiral electron gas: beyond Larmor's theorem. *Phys. Rev. B* **96**, 045301–1–14 (2017).

38. Salis, G., Walser, M. P., Altmann, P., Reichl, C. & Wegscheider, W. Dynamics of a localized spin excitation close to the spin-helix regime. *Phys. Rev. B* **89**, 045304–1–7 (2014).

39. Ferreira, G. J., Hernandez, F. G. G., Altmann, P. & Salis, G. Spin drift and diffusion in one- and two-subband helical systems. *Phys. Rev. B* **95**, 125119–1–12 (2017).

40. Holstein, T. & Primakoff, H. Field dependence of the intrinsic domain magnetization of a ferromagnet. *Phys. Rev.* **58**, 1098–113 (1940).

41. Perez, F. Spin-polarized two-dimensional electron gas embedded in a semimagnetic quantum well: Ground state, spin responses, spin excitations, and Raman spectrum. *Phys. Rev. B* **79**, 045306–1–16 (2009).

42. Hankiewicz, E. M., Vignale, G. & Tserkovnyak, Y. Inhomogeneous Gilbert damping from impurities and electron-electron interactions. *Phys. Rev. B* **78**, 020404–1–4 (2008).

43. Lipparini, E. *Modern many-particle physics*. 2nd edn, (World Scientific, Singapore, 2008).

44. Bao, J. M., Pfeiffer, L. N., West, K. W. & Merlin, R. Ultrafast dynamic control of spin and charge density oscillations in a GaAs quantum well. *Phys. Rev. Lett.* **92**, 236601–1–4 (2004).

45. Barate, P. et al. Collective nature of two-dimensional electron gas spin excitations revealed by exchange interaction with magnetic ions. *Phys. Rev. B* **82**, 075306–1–9 (2010).

46. Dreyhaupt, A., Winnerl, S., Dekorsy, T. & Helm, M. High-intensity terahertz radiation from a microstructured large-area photoconductor. *Appl. Phys. Lett.* **86**, 121114–1–3 (2005).

47. Giuliani, G. F. & Vignale, G. *Quantum Theory of the Electron Liquid*. (Cambridge University Press, Cambridge, 2005).

## Acknowledgements

S.K. and C.A.U. are supported by DOE Grant DE-FG02-05ER46213. F.P. acknowledges support from the Fondation CFM, C'NANO IDF and ANR.

## Author Contributions

C.A.U., I.D'A. and F.P. conceived the study, developed the theory, and wrote the manuscript. S.K. carried out the linear response calculations of the spin-wave dispersions, participated in the discussions and contributed to the review of the manuscript.

## Additional Information

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018