Spin-relaxation time in materials with broken inversion symmetry and large spin-orbit coupling

TL;DR

This paper investigates spin-relaxation times in materials with strong spin-orbit coupling and broken inversion symmetry, using combined Monte Carlo and diagrammatic methods to reveal new regimes and analogies with NMR phenomena.

Contribution

It introduces a combined computational approach that accurately models spin relaxation, identifying novel regimes and providing a comprehensive phase space analysis.

Findings

Quantitative agreement between Monte Carlo and diagrammatic methods.

Identification of non-exponential spin relaxation regimes.

Discovery of regimes where spin relaxation equals momentum relaxation.

Abstract

We study the spin-relaxation time in materials where a large spin-orbit coupling (SOC) is present which breaks the spatial inversion symmetry. Such a spin-orbit coupling is realized in zincblende structures and heterostructures with a transversal electric field and the spin relaxation is usually described by the so-called D'yakonov-Perel' (DP) mechanism. We combine a Monte Carlo method and diagrammatic calculation based approaches in our study; the former tracks the time evolution of electron spins in a quasiparticle dynamics simulation in the presence of the built-in spin-orbit magnetic fields and the latter builds on the spin-diffusion propagator by Burkov and Balents [Burkov and Balents Phys. Rev. B. 69, 245312 (2004).]. Remarkably, we find a parameter free quantitative agreement between the two approaches and it also returns the conventional result of the DP mechanism in the appropriate limit. We discuss the full phase space of spin relaxation as a function of SOC strength, its distribution, and the magnitude of the momentum relaxation rate. This allows us to identify two novel spin-relaxation regimes; where spin relaxation is strongly non-exponential and the spin relaxation equals the momentum relaxation. A compelling analogy between the spin-relaxation theory and the NMR motional narrowing is highlighted.

Tables

Table 1. Table 1: Summary of the relaxation regimes which are encountered for different values of the momentum scattering rate Γ Γ \Gamma , average SOC energy ⟨ ℒ ⟩ delimited-⟨⟩ ℒ \left<\mathcal{L}\right> and its spread Δ ​ ℒ Δ ℒ \Delta\mathcal{L} .

Condition	Relaxation type
	exponential (DP regime), ${\Gamma }_{\text{s}}=\alpha {ℒ}^2/\Gamma$
	oscillatory+exponential decay, ${\Gamma }_{\text{s}}=\alpha \Gamma$
	non-exponential, $T_2^{\ast }=\hslash /\Delta ℒ$, $T_{1,2}=\hslash /\Gamma$

Equations

$${\Gamma }_{\text{s}}=\alpha \frac{⟨{\left|\mathcal{L}\right|}^2⟩}{\Gamma }$$

$${\Gamma }_{\text{s}}=\Delta \mathcal{L}$$

$$H_0=\frac{{\hslash }^2{\mathbf{k}}^2}{2m}+\frac{\mathcal{L}}{k_{\text{F}}}(s_x{k}_y-s_y{k}_x)$$

$$s_z(t)\propto {\mu }_{0}M(t)={\int }_{-\infty }^t\chi \left(t-t^{\prime }\right)B\left(t^{\prime }\right)\text{d}{t}^{\prime }$$

$$D_{zz}(\omega )=\frac{(-i\Gamma +\mathcal{L}-\hslash \omega )(i\Gamma +\mathcal{L}+\hslash \omega )}{-{\hslash }^2{\omega }^2-i\Gamma \hslash \omega +{\mathcal{L}}^2}.$$

$$\hslash {\omega }_{1,2}=\frac{-i\Gamma \pm \sqrt{4{\mathcal{L}}^2-{\Gamma }^2}}{2}.$$

$$\frac{1}{{\tau }_{\text{s}}}=-\text{Im}{\omega }_{1,2}.$$

$$(\hslash \omega {)}^3+2i\Gamma (\hslash \omega {)}^2-\left({\Gamma }^2+{\mathcal{L}}^2\right)\hslash \omega -i\frac{{\mathcal{L}}^2\Gamma }{2}=0.$$

$$\mathbf{\Omega }(\mathbf{k})=\frac{\mathcal{L}}{\hslash k_{\text{F}}^3}\left[k_x\left(k_y^2-k_z^2\right),k_y\left(k_z^2-k_x^2\right),k_z\left(k_x^2-k_y^2\right)\right].$$

$$\omega \cdot \text{Re}\mathcal{S}(\omega )\propto \text{Im}\chi (\omega ).$$

$$H_0=\frac{{\hslash }^2{\mathbf{k}}^2}{2m}+\frac{{\mathcal{L}}_{\mathrm{R}}}{k_{\mathrm{F}}}(s_x{k}_y-s_y{k}_x)+\frac{{\mathcal{L}}_{\mathrm{D}}}{k_{\mathrm{F}}}(s_y{k}_y-s_x{k}_x).$$

$$\displaystyle\chi_{xx}(\omega)=\lim_{\eta\rightarrow 0+}\frac{m}{2\pi\hbar^{2}}\Bigg{(}1+$$

$$\displaystyle\frac{\left(\omega+\mathrm{i}\eta\right)^{2}}{\sqrt{\left(\frac{\mathcal{L}_{\mathrm{R}}+\mathcal{L}_{\mathrm{D}}}{\hbar}\right)^{2}-\left(\omega+\mathrm{i}\eta\right)^{2}}\sqrt{\left(\frac{\mathcal{L}_{\mathrm{R}}-\mathcal{L}_{\mathrm{D}}}{\hbar}\right)^{2}-\left(\omega+\mathrm{i}\eta\right)^{2}}}\Bigg{)}$$

$$|\psi (0)⟩=|\uparrow ⟩=\frac{v_2^{(-)}|+⟩-v_2^{(+)}|-⟩}{v_2^{(-)}{v}_1^{(+)}+v_2^{(+)}{v}_1^{(-)}},$$

$$|\psi (t)⟩=\frac{v_2^{(-)}{\mathrm{e}}^{-i{E}_{+}t/\hslash }|+⟩-v_2^{(+)}{\mathrm{e}}^{-i{E}_{-}t/\hslash }|-⟩}{v_2^{(-)}{v}_1^{(+)}+v_2^{(+)}{v}_1^{(-)}}$$

$$S_z(t,\mathbf{k})$$

$$S_z(t)={\int }_{F.S}d\mathbf{k}{S}_z(t,\mathbf{k})$$

$$S_z(t,\mathbf{k})=\cos (\Omega (\mathbf{k})t).$$

$$I(\omega )=\Re \frac{2\mathrm{i}\omega +4{\Gamma }_{\text{c}}}{\left({\Omega }^2-{\omega }^2\right)+2\mathrm{i}\omega {\Gamma }_{\text{c}}}.$$

$$I(\omega )=\Re \left(\frac{A}{\omega -{\omega }_1}+\frac{B}{\omega -{\omega }_2}\right),$$

$${\dot{n}}_1$$

$${\dot{n}}_2$$

$${\partial }_t(n_1-n_2)$$

Full text

Introduction

It is an intriguing possibility to employ the electron spins as information carriers (known as spintronics FabianRMP ). This prospect has revived the experimental and theoretical studies of spin-relaxation in semiconductors and metals. It is the spin-relaxation time which characterizes how rapidly a non-equilibrium spin population decays and it therefore determines whether a material is suitable for spintronics purposes. The theory of spin-relaxation differs for materials with and without spatial inversion symmetry FabianRMP ; WuReview : the Elliott-Yafet (EY) theory Elliott ; yafet1963g describes the former while the so-called D’yakonov-Perel’ (DP) mechanism describes the spin-relaxation for the latter case DyakonovPerelSPSS1972 . We note that the common physical picture to unify the two approaches was developed in Ref. Boross2013, .

The DP theory describes the dominant spin-relaxation mechanism in large band-gap III-V semiconductors (e.g. GaAs) with the zincblende structure (the so-called bulk inversion asymmetry) and for semiconductor heterostructures with an applied transversal electric field (the so-called structure inversion asymmetry). These two cases are known as Dresselhaus or Bychkov-Rashba spin-orbit coupling (SOC), respectively.

The DP spin-relaxation mechanism turned out to be particularly relevant for novel, spintronics candidate materials with two-dimensional structure, such as e.g. mono or bi-layer graphene GeimRMP ; HuertasPRL2009 and transition metal dichalcogenide monolayers Heinz_TMDC ; Wu_TMDC ; Ozyilmaz_TMDC_DP .

The spin-relaxation mechanism, which dominates in materials without inversion symmetry, is depicted schematically in Fig. 1a. The inversion symmetry breaking SOC splits the spin-up/down states at the Fermi level and acts as a $\bm{k}$ dependent magnetic field or SOC field. When the electrons are treated in the quasi-particle approximation, they are assumed to move around in the material and suffer momentum scattering on a timescale of $\tau$. The electron spins precess around the axis of the SOC field with a corresponding Larmor frequency, $\Omega(\bm{k})$. The rigorous derivation of the DP result DyakonovPerelSPSS1972 ; FabianRMP gives that the spin-relaxation time, $\tau_{\text{s}}$ is inversely proportional to $\tau$ if the $\left<\Omega\right>\cdot\tau\ll 1$ holds (here $\left<\Omega\right>$ is an average value of the Larmor frequencies). When rewritten for the quasiparticle momentum scattering rate: $\Gamma=\hbar/\tau$ and the spin scattering rate: $\Gamma_{\text{s}}=\hbar/\tau_{\text{s}}$, the DP result reads:

[TABLE]

where $\alpha$ is a band structure dependent parameter around unity and $\left<\mathcal{L}\right>$ is an average value of the SOC induced splitting and is related to the SOC fields by $\left<\mathcal{L}\right>=\hbar\left<\Omega\right>$.

Rather than a rigorous derivation, we give an illustration of the DP result in Fig. 1b. In the clean limit (i.e. $\Gamma=0$), the dynamic spin-susceptibility, Im$\chi(\omega)$, is an odd function around the origin which is given solely by the distribution of the SOC fields and it describes pure dephasing without relaxation. Im$\chi(\omega)$ could be observed in a clean semiconductor by e.g. an electron spin resonance experiment. The conventional DP theory applies in the dirty limit, i.e. when $\Gamma\gg\hbar\left<\Omega\right>$, which is depicted in Fig. 1b. (dotted curve is the quasi-particle spectral function). This is in fact the textbook situation of motional narrowing, which is well known in NMR spectroscopy AbragamBook and it formally leads to the DP result. In the spintronics literature, the DP regime is also known as motional slowing of the spin-relaxation.

The opposite limit, i.e. when $\Gamma\lesssim\hbar\left<\Omega\right>$ is less well-known and studied although recent experiments on ultra-pure GaAs indicates Oscillatory_DP_exp1 ; Oscillatory_DP_exp2 that this regime is accessible. In addition, the theoretical foundations of this regime were presented in Ref. Oscillatory_DP_theor, . This regime corresponds to semiconductors with a large SOC and low quasiparticle scattering. The $\Gamma\ll\hbar\left<\Omega\right>$ limit is discussed in Ref. FabianRMP, and it is argued that the distribution of the SOC fields causes a rapid dephasing of a spin ensemble orientation on the timescale of $\tau_{\text{s}}=1/\Delta\Omega$, where the latter quantity is the characteristic distribution width or variance of the SOC fields. This process is analogous to the so-called reversible dephasing in NMR spectroscopy with the timescale known as $T_{2}^{*}$ (Ref. SlichterBook, ). A momentum scattering after the dephasing causes memory loss therefore the observable spin-decay time roughly equals the spin-dephasing time. In our notation, the spin-relaxation rate was suggested to read in this limit as FabianRMP :

[TABLE]

where $\Delta\mathcal{L}=\hbar\Delta\Omega$ is the distribution width of the SOC splitting.

While we believe this qualitative picture to be correct, this regime requires a more quantitative description and as we show herein, the simple description breaks down depending on the relative magnitude of $\Gamma$ and $\Delta\mathcal{L}$. Another unresolved issue is whether the motional narrowing description of D’yakonov and Perel’ could be continued for cases when its conditions are not fulfilled. This regime would be relevant for clean semiconductors.

These shortcomings motivated us to study the full phase space of $\left(\Gamma,\left<\mathcal{L}\right>,\Delta\mathcal{L}\right)$ with emphasis on the case of large SOC, i.e. when $\Gamma\ll\mathcal{L}$. We compare two distinct approaches: a Monte Carlo method which is essentially a quasiparticle kinetics based approach using numerical simulation and the diagrammatic technique, accounting for both SOC and impurity scattering. We studied the Bychkov-Rashba Hamiltonian for a two-dimensional electron gas with both methods. The two approaches yield quantitatively identical results for the spin-relaxation time without adjustable parameters. We observe a non single-exponential spin decay with both methods in certain cases. This validates the Monte Carlo approach and lets us present a method which allows calculation of $\tau_{\text{s}}$ (or $\Gamma_{\text{s}}$) for an arbitrary SOC model including e.g. the Dresselhaus SOC. We identify a yet unknown regime: when the distribution of the SOC is sharper than the broadening parameter, i.e. $\Gamma\gg\Delta\mathcal{L}$, then $\Gamma_{\text{s}}\approx\Gamma$ is realized. For the case of $\Gamma\lesssim\left(\left<\mathcal{L}\right>,\Delta\mathcal{L}\right)$ a strongly non-exponential spin decay is observed. All the results can be elegantly visualized by considering the evolution of the dynamic spin-susceptibility from the clean to the dirty limit. We note that we focus on the effect of inversion symmetry breaking SOC fields on spin-relaxation and we do not discuss other mechanisms, such as e.g. the Elliott-Yafet mechanism Elliott ; yafet1963g .

Methods

The Monte Carlo simulations

The time evolution of the electron spin direction in the presence of internal spin-orbit coupling is studied with a Monte Carlo approach which essentially mimics the mathematical description of D’yakonov and Perel’ DyakonovPerelSPSS1972 : an initially polarized electron spin ensemble travels in a solid where the SOC related and momentum ($\bm{k}$) dependent magnetic fields ($\bm{B}(\bm{k})$) are present. The spins precess freely around the SOC fields as classical variables with angular momentum $\bm{\Omega}(\bm{k})$ between two scattering events, where $\bm{k}$ points to the Fermi surface. The scattering, what we model as a stochastic process, induces a new random $\bm{k}$ on the Fermi surface thus precession starts around a new axis. The momentum-relaxation is described by a Poisson process with an expectation value of $1/\tau$. Thus the time between successive scattering events follows an exponential distribution with the same $1/\tau$ parameter. We also assume that the momentum scatters uniformly on the Fermi surface and that the scattering does not affect the spin direction, i.e. the Elliott-Yafet type spin-flip mechanisms Elliott ; yafet1963g are not considered.

After running the simulation on several independent spins, the mean spin component is sampled in uniform time intervals. We keep the momentum relaxation $\tau=1$ constant in the simulations and vary the strength of the spin-orbit coupling and the proper time dependent data are obtained by a rescaling. The spin is measured in units of $\hbar/2$.

The physical picture behind this approach is that i) the system is at $T=0\,\text{K}$, ii) $\Gamma$ appears solely through the momentum scattering time without considering its origin (i.e. impurities, phonons etc) and the $\Gamma$ related state broadening, iii) we consider only $k=k_{\text{F}}$ states, i.e. the SOC induced band splitting is disregarded in the initial state, and that iv) all electrons on the Fermi surface are spin polarized and the rest of the Fermi sea is unpolarized in the initial state. Neglecting the effect of $\Gamma$ on the band structure for the Monte Carlo approach is justified by comparing the result with that of the diagrammatic technique (which correctly accounts for the finite life-time effects) in Ref. BurkovBalents, . Considering $k=k_{\text{F}}$ states only for the band structure without SOC splitting is a common approach in similar spin dynamics BurkovBalents studies and when calculating the SOC related dynamic spin-susceptibility ErlingssonPRB . It essentially corresponds to neglecting corrections on the order of $\mathcal{L}/E_{\text{F}}$, where $E_{\text{F}}$ is the Fermi energy and further details about this approximation are provided in the Supplementary Material. The assumption iv) corresponds to the application of a small magnetic field for the time $t<0$ which causes all spins on the Fermi surface to be polarized and the rest of the Fermi sea to remain unpolarized and the magnitude of the required magnetic field is discussed below.

Typical time dependent magnetization curve for an ensemble of electrons is shown in Fig. 2. In this case, we assumed a two-dimensional electron gas with a Bychkov-Rashba type SOC ($x$ and $y$ denotes coordinates within the plane):

[TABLE]

where the first term is the kinetic energy, $k_{\text{F}}$ is the Fermi wavenumber, $\bm{k}^{2}=k_{x}^{2}+k_{y}^{2}$, $s_{x,y}$ and $k_{x,y}$ are the components of the spin and momentum, respectively. The corresponding SOC related Larmor frequencies read: $\bm{\Omega}(\bm{k})=\frac{\mathcal{L}}{\hbar k_{\text{F}}}\left[k_{y},-k_{x},0\right]$. We also assumed that the spins only on the Fermi surface are fully polarized along $z$ at $t=0$. We considered two cases, a large SOC ($\hbar\Omega=\Gamma$) and moderate SOC ($\hbar\Omega=0.4\Gamma$) and the result is shown in Fig. 2. As we show below, these are archetypes of the different relaxation regimes. A decaying magnetization is observed for both cases with and without an oscillating component. The real part of the Fourier transform of the time dependent $s_{z}(t)$ data, $\text{Re}\mathcal{S}(\omega)$ is also shown in Fig. 2, which depicts better the presence of the oscillation (peaks at $\omega\neq 0$) and a single decay (a nearly Lorentzian peak at $\omega=0$). The details of these calculations are analyzed below by fitting two Lorentzian curves to them, whose position and width give the frequency of the damped oscillations and the damping, respectively. These parameters are compared to the analytic calculations for the same Bychkov-Rashba Hamiltonian performed in Ref. BurkovBalents, . The clean limit, i.e. $\Gamma=0$ would give a $\text{Re}\mathcal{S}(\omega)$ with two Dirac-delta peaks at $\omega=\pm\Omega$ upon neglecting terms of the order $\mathcal{L}/E_{\text{F}}$.

We briefly discuss the relationship between $\mathcal{S}(\omega)$ and the dynamic spin-susceptibility, $\chi(\omega)$. The time dependent net magnetic moment, $s_{z}(t)$ is proportional to the sample magnetization. Assuming linear response theory to apply, the magnetization of the sample, $M$ is related to an external magnetic field, $B$ by:

[TABLE]

where $\mu_{0}$ is the vacuum permeability. Our calculation assumes that the strength of the magnetic field is such that electrons on the Fermi surface are spin-polarized and electrons on the next level are already compensated by spin-degeneracy. The corresponding magnetization is given by: $M\approx\frac{g\mu_{\text{B}}}{V_{\text{c}}}\frac{\Delta k}{k_{\text{F}}}$, where $\Delta k$ denotes the typical distance between allowed momentum space points. A magnetic field of $g\mu_{\text{B}}B_{0}=\frac{\Delta k}{k_{\text{F}}g(E_{\text{F}})}$ (where $g(E_{\text{F}})$ is the density of states at the Fermi energy) is required to sustain such a magnetization at $T=0$. This magnetic field tends to zero in the thermodynamic limit, i.e. it is sufficiently small for the linear response theory to be valid.

Considering that in our model $B(t^{\prime})=B_{0}\Theta(-t^{\prime})$ (where $\Theta(t^{\prime})$ is the Heaviside function), i.e. it is switched off at $t=0$, a straightforward calculation leads to $i\omega\cdot\mathcal{S}(\omega)\propto\chi(\omega)$ the proportionality constant being the infinitesimal magnetic field. We numerically confirm this proportionality for the two-dimensional Bychkov-Rashba and Dresselhaus Hamiltonians without impurities (data and additional discussion is given in the Supplementary Material) i.e. that $\omega\cdot\text{Re}\mathcal{S}(\omega)\propto\text{Im}\chi(\omega)$ holds in the clean limit. We used $\chi(\omega)$ data published in Ref. ErlingssonPRB, for the comparison within the above discussed $\mathcal{O}\left(\mathcal{L}/E_{\text{F}}\right)$ approximation. From Eq. (4), the above relation is expected to remain also valid in the dirty limit, i.e. $\omega\cdot\text{Re}\mathcal{S}(\omega)\propto\text{Im}\chi(\omega)$. We note that since $\text{Im}\chi(\omega)$ is an odd function of frequency, $\text{Re}\mathcal{S}(\omega)$ must be even. In general, $\text{Im}\chi(\omega)$ is known to carry information about the spin dynamics and losses, which explains why we present the numerically obtained $\text{Re}\mathcal{S}(\omega)$ throughout this contribution.

The Monte Carlo method we use can be readily applied for an arbitrary distribution of the SOC fields, e.g. for the three-dimensional Dresselhaus Hamiltonian as we discuss below. The system specific parameters are present through the actual $\bm{\Omega}(\bm{k})$ function. We note that in general the spin-relaxation is anisotropic, i.e. different oscillations and damping are observed for the same SOC distribution depending on the starting polarization direction of the spins.

Results of the diagrammatic technique

Burkov and Balents calculated the spin-relaxation in a two-dimensional electron gas BurkovBalents for a Rashba type SOC for an arbitrary value of the magnetic field. Their calculation turns out to be very general and although they applied it in the D’yakonov-Perel’ regime (i.e. for weak SOC), it can be also used for an arbitrary strength of the SOC. The model Hamiltonian was the same as Eq. (3), i.e. the 2DEG with Rashba SOC. The so-called spin-diffusion propagator, $D(\omega)$ was calculated in Ref. BurkovBalents, with a diagrammatic technique. We use Eqs. 36 and Eq. 38 in Ref. BurkovBalents, to obtain the spin-relaxation times from the spin-diffusion propagators, which reads in the $z$ direction as:

[TABLE]

The real and imaginary parts of the two poles ($\omega_{1,2}$) of Eq. (5) describe the oscillation frequency and the damping of the spin-propagation. The poles read:

[TABLE]

The spin-relaxation time is obtained from the poles as:

[TABLE]

Similar results can be obtained for the diffusion propagator which is perpendicular to the quantization axis ($D_{xy}$) and the poles are the roots of a third order polynomial:

[TABLE]

We emphasize that the result of Burkov and Balents is valid for any value of $\mathcal{L}$ and $\Gamma$ for the studied Hamiltonian. We use this result for comparison with the Monte Carlo simulations on the same Bychkov-Rashba SOC model system. We presume that analytic results could be obtained for other simple models of the SOC using the formalism of Ref. BurkovBalents, but these are beyond our scope.

Results and discussion

Validation of the Monte Carlo approach

As a first step, we validate the above described Monte Carlo approach by comparing the numerical results with that of the analytic calculations. The comparison is shown in Fig. 3. The analytic result in Ref. BurkovBalents, yields two parameters: the real and imaginary parts of poles of the spin-diffusion propagator ($\omega_{1,2}$) which are shown with solid curves in Fig. 3. The earlier represents the frequency of the oscillating component, whereas the latter describes the damping or relaxation and it is the spin-relaxation rate in frequency units. A fit to the numerically obtained $\text{Re}\mathcal{S}(\omega)$ data with several Lorentzian components (including both Kramer-Kronig pairs) yield the position of the Lorentzian as well as its width. This result is shown in Fig. 3. with symbols. We observe a surprisingly good agreement between the two types of data, which as we note is obtained without any adjustable paramaters. Essentially, the two types of calculations consider the same SOC Hamiltonian (the Bychkov-Rashba) and a two-dimensional electron gas and the same approximation (neglect of SOC splitting of the band structure, i.e. zeroth order in $\mathcal{L}/E_{\text{F}}$) but the methods are quite different. We can omit the $\left<...\right>$ notation for the Bychkov-Rashba model as the SOC field consists of two delta functions when the $\mathcal{O}\left(\mathcal{L}/E_{\text{F}}\right)$ corrections are neglected.

A horizontal arrow in Fig. 3. indicates the regime where the conventional DP mechanism is realized and both the numerical and Monte Carlo methods give $\Gamma_{\text{s}}=\mathcal{L}^{2}/\Gamma$ (i.e. herein $\alpha=1$). However, the figure also indicates the presence of two additional, previously unknown spin-relaxation regimes: on the far right, when $\Gamma\lesssim\mathcal{L}$, two peaks are present in with a broadening of $\Gamma_{\text{s}}=\alpha\Gamma$ (where $\alpha$=0.5). Another spin-relaxation regime occurs right below the ”bifurcation point” ($\mathcal{L}=0.5\Gamma$): therein two relaxation components with different broadening are present, i.e. it represents a non single-exponential relaxation. This regime crosses over smoothly to the conventional DP regime where a single exponential is observed: the weight of the component with larger broadening gradually disappears, which is apparent from the scatter in the parameters which are obtained by fitting the Monte Carlo data.

The conventional DP mechanism is often referred to as a motional narrowing effect, which gives rise to the $\Gamma_{\text{s}}\propto\Gamma^{-1}$ behavior. In fact, motional narrowing is the extremal case of a theory known as ”the effect of motion on the spectral lines”, which is well developed in e.g. NMR spectroscopy AbragamBook . We show in the Supplementary Material that a conventional textbook description of the so-called two-site moving nuclei problem (after Ref. AbragamBook, , Chapter X.) gives numerically identical results to the parameters of the above described Bychkov-Rashba model for the full motional range. While this is an interesting analogy, it is a straightforward argument that the two situations, i.e. momentum scattering of an electron between different $\bm{\Omega}(\bm{k})$ and motion of nuclei between different sites with different Larmor frequencies, leads to the same result. A more interesting consequence of this analogy is that the full phase space of $\left(\Gamma,\mathcal{L}\right)$ can be regarded as a motional problem and it is not restricted to the DP regime.

Spin-relaxation for the Dresselhaus spin-orbit coupling

The agreement between the spin-relaxation parameters as obtained from the Monte Carlo and from the analytic calculations validates the use of the numerical method for the two-dimensional electron gas with the Bychkov-Rashba SOC. Although it represents no formal proof, we believe that it justifies the use of the Monte Carlo method to obtain spin-relaxation parameters and eventually $\text{Re}\mathcal{S}(\omega)$ for more complicated distributions of the SOC fields, where analytic calculations are not available.

First, we note that the choice of the 2D electron gas with the Bychkov-Rashba SOC is somewhat exceptional as its $\text{Re}\mathcal{S}(\omega)$ in the clean limit consists of two delta functions (when quantization is along the $z$ axis), i.e. this model has a zero width of the SOC field distribution. Generally, the width of the SOC field distribution ($\Delta\mathcal{L}=\hbar\Delta\Omega$) could be sizeable, i.e. comparable to the average SOC field ($\left<\mathcal{L}\right>=\hbar\left<\Omega\right>$). In fact, the Dresselhaus SOC represents such a case. Its Larmor frequency distribution reads:

[TABLE]

The corresponding $\text{Re}\mathcal{S}(\omega)$ in the clean limit is shown in Fig. 4 obtained with $\mathcal{L}=1$. The SOC field distribution is sizeable and $\Delta\mathcal{L}$ and $\left<\mathcal{L}\right>$ have the same order of magnitude.

Two features are observed for this type of SOC in the clean limit: a Dirac-delta function in $\text{Re}\mathcal{S}(\omega=0)$ and ”beating” in $s_{z}(t)$. While similar in its form, the $\text{Re}\mathcal{S}(\omega)$ function is not identical to the histogram of the $\left|\bm{\Omega}(\bm{k})\right|$ Larmor frequency distribution, which is evident by the presence of the Dirac-delta function in the earlier. This is due to some geometric factors which appear in the calculation of $\text{Re}\mathcal{S}(\omega)$ and is discussed in depth in the Supplementary Material. The beats in the time domain data are the consequence of coherent spin oscillations and its details are specific for the angular distribution of the SOC field.

Fig. 4. also presents the time dependence of $s_{z}$ and $\text{Re}\mathcal{S}(\omega)$ for finite $\Gamma$. The D’yakonov-Perel’ limit is recovered when $\Gamma$ is much larger than $\left<\mathcal{L}\right>$ and $\Delta\mathcal{L}$. However, a novel regime is identified when $\Gamma\lesssim\left(\left<\mathcal{L}\right>,\Delta\mathcal{L}\right)$. Then, an initial rapid dephasing due to the distribution of SOC fields is present in agreement with Ref. FabianRMP, , however the dephasing is not complete (i.e. $s_{z}\neq 0$) and the remaining ensemble $s_{z}$ decays on the timescale of $\tau=\hbar/\Gamma$ only. This feature is due to the presence of the Dirac-delta function in addition to the SOC fields at $\omega\neq 0$. This observation mimics the situation encountered in pulsed NMR spectroscopySlichterBook : therein a rapid dephasing (on a timescale denoted as $T_{2}^{*}$) is caused by local magnetic field inhomogeneities which is not accompanied by a true information loss. It is followed by a true relaxation (a timescale denoted as $T_{2}$or in the absence of an external field $T_{1}=T_{2}$) where the information is inevitably lost. This phenomenon leads to the presence of NMR spin-echo, i.e. the spins can be restored in-phase on a timescale within $T_{2}$ with a suitable external excitation. Our observation predicts that a similar scheme may lead to the observation of spin-echo in semiconductors under the circumstances which correspond to the situation shown in Fig. 4.

We believe that the well-know Dresselhaus Hamiltonian is general enough to properly capture the essential features of spin-relaxation for the entire $\left(\Gamma,\left<\mathcal{L}\right>,\Delta\mathcal{L}\right)$ phase space. It is also the most important Hamiltonian which is relevant for most III-V semiconductors where bulk spatial inversion symmetry breaking occurs. We summarize our qualitative findings in a compact form in Table LABEL:Summary_table., which is briefly repeated herein: i) the conventional DP regime occurs when $\Gamma$ is much larger than the SOC, ii) the spin decay is oscillatory and spin-relaxation time has the same order of magnitude as momentum relaxation time when the SOC is significant but its distribution is sharp: $\Delta\mathcal{L}\ll\Gamma\lesssim\left<\mathcal{L}\right>$, iii) an NMR-like rapid spin dephasing followed by a true spin-relaxation occurs when $\Gamma\lesssim\left(\left<\mathcal{L}\right>,\Delta\mathcal{L}\right)$.

We finally comment on the most important approximation of our approach, i.e. that the SOC is smaller than the kinetic energy. This approximation is valid for most technically relevant semiconductors but the opposite is true for heavy elements such as e.g. Bi, where the SOC becomes the leading energy term. Those materials, however, are characterized by very unconventional spin-dynamical properties whose study was attempted in e.g. Ref. DoraSimonTopIns, .

Conclusions

In conclusion, we studied spin-relaxation in materials without a spatial inversion symmetry considering the full phases space of the momentum scattering, the SOC field strength and its distribution properties. Our main tool was a Monte Carlo method which conceptually follows the D’yakonov-Perel’ relaxation mechanism. The method is validated by a comparison to a diagrammatic approach based calculation of the spin-relaxation parameters. We identify a compelling analogy between spin-relaxation and the textbook description of NMR spectroscopy in the presence of motion. The method, when applied for the Dresselhaus SOC Hamiltonian, allowed to expand our knowledge about the different spin-relaxation regimes.

Acknowledgements

Work supported by the ERC Grant Nr. ERC-259374-Sylo and by the Hungarian National Research, Development and Innovation Office (NKFIH) Grant Nr. K119442. A.K. acknowledges the Bolyai Program of the Hungarian Academy of Sciences.

Author Contributions

LS and BD performed the analytic calculations, LS made the numerical studies. AK validated the results. FS outlined the overall composition of the paper and all authors contributed to writing the manuscript.

Additional Information

Competing financial interests: The authors declare no competing financial interests.

Note added

We became aware of a contributionCosacchi , which investigates similar phenomena such as our work, shortly prior to submitting our manuscript.

Appendix A The effect of the $\mathcal{O}\left(\mathcal{L}/E_{\text{F}}\right)$ approximation

It was mentioned in the main text that both the Monte Carlo and the diagrammatic technique in Ref. BurkovBalents, neglects the effect of the SOC on the Fermi surface. A proper calculation should consider that the bands are split due to the SOC and therefore the corresponding $\Omega(\bf{k})$ is $\left|\bf{k}\right|$ dependent due to this effect. In contrast, our approximation neglects this effect and Fig. 5 depicts this approximation for a one-dimensional band dispersion. The figure suggests that corrections to our approximation are on the order of $\mathcal{O}\left(\mathcal{L}/E_{\text{F}}\right)$.

Appendix B Relationship between $\mathcal{S}$ and the dynamic spin-susceptibility

We discussed in the main text that the time decay of a spin-polarized ensemble (described by $s_{z}(t)$) is calculated with a Monte Carlo approach for both the clean and dirty cases. We found that the real part of its Fourier transform, $\text{Re}\mathcal{S}(\omega)$ can be conveniently displayed in order to demonstrate the spin-relaxation properties. We also showed that a relation between $\text{Re}\mathcal{S}(\omega)$ and the dynamic spin-susceptibility, $\chi(\omega)$ holds:

[TABLE]

We present the above relationship on a quantitative agreement which we numerically obtained by comparing our Monte Carlo results on $\mathcal{S}(\omega)$ with analytic calculations of $\chi(\omega)$ for a particular Hamiltonian which is available in the literature (Ref. ErlingssonPRB, ). The result is shown in Fig. 6 with the calculation details as follows.

Erlingsson et. al. (Ref. ErlingssonPRB, ) calculated $\chi(\omega)$ for a two-dimensional electron gas for a Hamiltonian containing both Bychkov-Rashba and Dresselhaus type SOC terms as follows:

[TABLE]

where $\mathcal{L}_{\mathrm{R}}$ and $\mathcal{L}_{\mathrm{D}}$ are the strengths of the Bychkov-Rashba and the Dresselhaus type spin-orbit couplings, respectively.

The dynamic spin-susceptibility was calculated in the absence of momentum relaxation and when $\mathcal{L}_{\mathrm{R}}\sim\mathcal{L}_{\mathrm{D}}\ll E_{\mathrm{F}}$. This limit matches the above discussed $\mathcal{O}\left(\mathcal{L}/E_{\text{F}}\right)$ approximation. Eq. 15 in Ref. ErlingssonPRB, gives the dynamic spin-susceptibility for the $x$ direction as:

[TABLE]

Fig. 6. demonstrates a good agreement between the two kinds of data which supports the statement in the main text concerning the connection between $\mathcal{S}$ and $\chi$.

Appendix C Analytic treatment of the spin dynamics by the time evolution in the clean limit

We give an analytic description of the time evolution of the magnetization of an electron ensemble subjected to internal SOC fields to verify further the Monte Carlo method applied in the main text since this description is equivalent to the numerical method in the clean limit. We consider the SOC in the form $H_{\rm 0}=\frac{1}{2}\mathbf{\Omega}(\mathbf{k})\cdot{\bm{\sigma}}$, where ${\bm{\sigma}}$ is a vector composed by the Pauli matrices, and $\mathbf{\Omega}(\mathbf{k})$ is the $\mathbf{k}$-dependent internal SOC field.

The time evolution of the state of an electron under the SOC is determined by the time evolution operator $U(t)={\rm exp}(-iH_{\rm 0}t/\hbar)$. Supposing that the electron is initially in spin-up state, i.e. its spin is polarized along the $z$ direction,

[TABLE]

its state ket at time $t$ is obtained as

[TABLE]

by applying the time evolution operator. Here, $|\pm\rangle=[v_{1}^{(\pm)},v_{2}^{(\pm)}]$ and $E_{\pm}$ are the eigenkets and eigenenergies of the Hamiltonian $H_{\rm 0}$. Finally, the time development of the $z$ component of the electron spin is obtained as

[TABLE]

where $\Omega=\sqrt{\Omega_{x}^{2}+\Omega_{y}^{2}+\Omega_{z}^{2}}$. The quantity $S_{z}(t)$, i.e. the $z$ component of a spin ensemble, calculated by Monte Carlo method in the main text is obtained as

[TABLE]

within this approach, i.e. by integration over $\mathbf{k}$ on the Fermi surface. Arbitrary $\mathbf{\Omega}({\mathbf{k}})$ can be considered in the Hamiltonian $H_{\rm 0}$ such as the two-dimensional Bychkov-Rashba SOC or the three-dimensional Dresselhaus case discussed in the main text. We note that for a complicated distribution of the SOC fields the $\mathbf{k}$ integration in Eq. (16) might not be performed analytically. In such a case, the integration can be performed by choosing random $\mathbf{k}$ values on the Fermi surface.

When this calculation is performed according to Eq. (LABEL:eq-expSz1) for different model Hamiltonians such as those given in the main text, we obtained numerically identical results (data not shown) as for the Monte Carlo, and the analytic result for the dynamic spin susceptibility given in Eq. (12) is reproduced as well.

From expression (LABEL:eq-expSz1) it is obvious that in a two-dimensional case, i.e. when $\Omega_{z}=0$:

[TABLE]

By taking the Bychkov-Rashba SOC, $\mathbf{\Omega}(\mathbf{k})=\frac{{\cal L}}{\hbar k_{\rm F}}[k_{x},k_{y},0]$, $\Omega({\mathbf{k}})$ becomes $\mathbf{k}$ independent as $\Omega({\mathbf{k}})={\cal L}/\hbar=\Omega$, which means a single oscillating component in $S_{z}(t)$ as it is shown by dashed curve in Fig.2 of the main text, and two Dirac-delta peaks at $\pm\Omega$ in the real part of the Fourier transform $\mathcal{S}(\omega)$. In three-dimensional cases, there is always a $t$-independent non-zero term in $S_{z}(t)$ coming from the last term in Eq. (LABEL:eq-expSz1), which explains the finite $S_{z}$ value in Fig.4 of the main text with $\Gamma=0$. This non-zero and time independent term corresponds to a Dirac-delta function centered on $\omega=0$ in $\mathcal{S}(\omega)$.

The origin of this effect is further supported by a geometric consideration which is depicted in Fig. 7. The presence of the Dirac-delta peak for $\mathcal{S}(\omega=0)$ is a generic feature and its absence for the two-dimensional electron gas and the Bychkov-Rashba SOC is an exeption. For the latter, when the spins are aligned perpendicular to the 2D plane, all SOC fields are in the plane, i.e. the precession of the spins around the built in $\mathbf{\Omega}(\mathbf{k})$ results in a zero-averaged net magnetization. However, for a general distribution of the SOC fields and the corresponding $\mathbf{\Omega}(\mathbf{k})$ vectors, the precession of the spins retains a finite positive $s_{z}$ component as Fig. 7. depicts. A straightforward geometric consideration shows that the $\left<s_{z}\right>$, i.e. the Dirac delta function strength is given by $\Omega_{z}^{2}/\Omega^{2}$ ($\Omega_{z}$ and $\Omega$ are the $z$ component and the magnitude of the $\mathbf{\Omega}(\mathbf{k})$ vector, respectively) for a particular $\mathbf{\Omega}(\mathbf{k})$ component. Similarly, we obtain that the amplitude of the oscillation goes as $1-\Omega_{z}^{2}/\Omega^{2}$, in full agreement with Eq. (LABEL:eq-expSz1).

In Fig. 8., we present the comparison between the $\text{Re}\mathcal{S}(\omega)$ and the histogram of the internal Larmor frequency distribution for the three-dimensional Dresselhaus model. The latter data integrated for the positive frequencies gives the average of the SOC fields which coincides with the frequency value where it is peaked and its weighted integral, i.e the integral of the histogram values multiplied by the related frequency, gives unity since all weights are summed up in this way.

Appendix D Analogy between the motional narrowing and the spin-relaxation for the 2D BR modell

Abragam AbragamBook considered the so-called two-site NMR motional narrowing problem: a nuclei is allowed to jump with the transition rate $\Gamma_{\text{c}}=1/\tau_{\text{c}}$ between two sites with different local Larmor frequencies: $\pm\Omega$ around a central Larmor frequency (defined as zero in this case). The resulting NMR lineshape is shown in Fig. 9 for a fixed $\Omega=1$ and different values of the jumping frequency, $\Gamma_{\text{c}}$. The analogy between the spin-relaxation and the motional narrowing is clear: the $\pm\Omega$ local Larmor frequencies correspond to the built-in Zeeman field distribution of the spin-relaxation problem and the jumping frequency ($\Gamma_{\text{c}}$) of the motional narrowing problem corresponds to the $\Gamma$ momentum relaxation rate (besides a factor 2 which is discussed below). The analogy can be quantified for the simplest case as follows.

[TABLE]

The denominator of Eq. (18) has poles at $\omega_{1,2}=\mathrm{i}\left(\Gamma_{\text{c}}\pm\sqrt{\Gamma_{\text{c}}^{2}-\Omega^{2}}\right)$, i.e. Eq. (18) can be rewritten as:

[TABLE]

where $A=-\mathrm{i}\frac{\Gamma_{\text{c}}+\Delta}{\Delta}$ and $B=\mathrm{i}\frac{\Gamma_{\text{c}}-\Delta}{\Delta}$. Herein, we introduced $\Delta=\sqrt{\Gamma_{\text{c}}^{2}-\Omega^{2}}$. Evaluation of Eqs. (18) and (19) yields the curves shown in Fig. 9.

Our definition of $\Gamma$ differs from that of $\Gamma_{\text{c}}$ in Abragam’s work in a factor 2 as discussed herein, as $\Gamma$ corresponds to the momentum relaxation rate and $\Gamma_{c}$ corresponds to the transition rate between the two states. The rate equations for the two state’s populations:

[TABLE]

We can see that $\frac{1}{\tau_{\text{m}}}=\frac{2}{\tau_{c}}$, so $\Gamma=2\Gamma_{c}$. With this change, the real and imaginary values of the above defined $\omega_{1,2}$ are shown in Fig. 10. The real part of the roots describe the position of the two peaks and the imaginary parts describe the linewidths in agreement with Fig. 10. Remarkably, this figure is identical to the spin-relaxation problem for the 2D Rashba model in the main paper with a straightforward identification of the correspondence of the two parameters.