---
title: Orbital Edelstein Susceptibility
url: https://www.emergentmind.com/topics/orbital-edelstein-susceptibility
type: topic
---

# Orbital Edelstein Susceptibility

Searching arXiv for recent and foundational papers on orbital Edelstein susceptibility to ground the article in published work.
Orbital Edelstein susceptibility is the linear-response coefficient that quantifies how an applied electric field, charge current, or superconducting phase gradient generates a nonequilibrium orbital magnetization or orbital angular momentum in a system with broken inversion symmetry. In its most common normal-state form, it is written as
$$
M_i^{\mathrm{orb}}(\omega)=\sum_j \chi_{ij}^{\mathrm{orb}}(\omega) E_j(\omega),
$$
while in superconductors it is often defined through the phase gradient,
$$
\langle \hat{L}_\mu\rangle=-\frac{1}{2}\chi_{\mu\nu}q_\nu.
$$
The literature shows that the quantity is not tied to a single microscopic mechanism or unit convention: it can be formulated for local atomic projections, unit-cell orbital moments, magnetization densities, or supercurrent-induced orbital polarization, and it has been evaluated with Kubo, semiclassical Boltzmann, modern orbital magnetization, and Bogoliubov–de Gennes response frameworks [1905.08279; 2107.07476; 2504.07665].

## 1. Definition and constitutive relations

The orbital Edelstein effect is the generation of a non-equilibrium orbital polarization by an applied electric field in noncentrosymmetric or gyrotropic materials. In the antiferromagnetic formulation of the generalized Rashba–Edelstein effect, the orbital susceptibility tensor is the coefficient in
$$
M_i^{\mathrm{orb}}(\omega)=\sum_j \chi_{ij}^{\mathrm{orb}}(\omega)E_j(\omega),
$$
with the full induced magnetic polarization given by $\delta M=\mu_B\delta(L+2S)$, so that orbital and spin susceptibilities can be discussed separately as $\chi^L_{ij}$ and $\chi^S_{ij}$ [1905.08279].

The same object is written in several equivalent but context-dependent ways. Carbon nanotube work defines the observable as the non-equilibrium orbital magnetic moment per unit cell along the tube axis, $m_{L_z}^{uc}=\chi_z^{L_z}E_z$, and frequently reports the susceptibility per unit tube length, $\chi_z^{L_z}\cdot a/T$, in units of $\mu_B\cdot\mathrm{m}/\mathrm{V}$ [2504.07665]. Oxide-interface studies write the induced moment per 2D unit cell as $(A_0/A)m=[\chi^s+\chi^l]E\equiv \chi E$, with $\chi^l$ the orbital Edelstein efficiency tensor [2006.14958]. In bilayer Rashba systems, the same distinction appears between a “per-field” susceptibility, $M_i^{\mathrm{orb}}=\alpha^{\mathrm{orb}}_{ij}E_j$, and a “per-current” susceptibility, $M_i^{\mathrm{orb}}=\chi^{\mathrm{orb}}_{ij}j_j$, related by $\chi^{\mathrm{orb}}=\alpha^{\mathrm{orb}}\sigma^{-1}$ [2307.02872].

A central point in the literature is that orbital Edelstein susceptibility is not reducible to spin Edelstein susceptibility with a different prefactor. In noncentrosymmetric antiferromagnets, the dominant staggered orbital response remains finite and essentially unchanged when spin–orbit coupling is turned off, while the spin susceptibility vanishes with spin–orbit coupling [1905.08279]. In Cu\(_2\)WSe\(_4\), the orbital Edelstein effect is described as current-induced orbital magnetization in a noncentrosymmetric conductor, with constitutive relation $M_i=\alpha^{\mathrm{orb}}_{ij}E_j$ and units of $A/V$ for $\alpha^{\mathrm{orb}}_{ij}$ [2512.16387]. This diversity of definitions reflects different observables, not a disagreement about the underlying response.

## 2. Microscopic formulations

Two broad microscopic routes dominate the subject. The first is Kubo-type linear response, where the orbital operator and the current or velocity operator appear explicitly in interband and intraband terms. For noncentrosymmetric antiferromagnets, the susceptibility is evaluated from a Kubo-like expression with
$$
\chi^{B}_{\alpha\beta}(\omega)=\chi^{B,\mathrm{inter}}_{\alpha\beta}(\omega)+\chi^{B,\mathrm{intra}}_{\alpha\beta}(\omega),
$$
where $B=\hat L$ or $\hat S$, $p_j=m_ev_j$, and a single phenomenological broadening $\hbar/\tau\approx 0.41\ \mathrm{eV}$ is used in both channels [1905.08279]. In the diffusive parabolic-band problem with an asymmetric scalar potential, the orbital Edelstein effect is also derived from a Kubo bubble between the current vertex and the orbital magnetic moment operator, yielding a result that scales as a cube of the momentum relaxation time and does not rely on spin–orbit coupling [2407.00516].

The second route is semiclassical Boltzmann transport in the relaxation-time approximation. In the oxide two-dimensional electron gas at AlO\(_x\)/SrTiO\(_3\), the orbital susceptibility is computed as
$$
\chi^l_{ij}=\frac{A_0 e\mu_B}{A\hbar}\sum_k \langle l\rangle_{k,i}\Lambda_{k,j}\delta(\mathcal{E}_k-\mathcal{E}_F),
$$
with $\Lambda_k=\tau_0v_k$ and the orbital moment built from atomic-like orbital angular momentum expectation values in a realistic eight-band \(t_{2g}\) model [2006.14958]. The (111) LaAlO\(_3\)/SrTiO\(_3\) interface uses an analogous Boltzmann formula based on $\langle L_\alpha\rangle_n(\mathbf{k})$ and the shifted distribution function, evaluated at $T=10\ \mathrm{K}$ and $\tau_0=3.4\times10^{-12}\ \mathrm{s}$ [2207.07663]. Chiral CNTs, Cu\(_2\)WSe\(_4\), analytically solvable helical models for chirality-induced spin selectivity, and gated monolayer TMDs all employ related Boltzmann expressions, although the orbital moment entering the transport kernel may come from modern orbital magnetization rather than from an atomic operator [2504.07665; 2512.16387; 2502.04978; 2509.25634].

A further distinction concerns how the orbital moment itself is represented. Some works use the expectation value of an orbital angular momentum operator in a restricted orbital basis, as in oxide \(t_{2g}\) models [2006.14958; 2207.07663]. Others use the gauge-invariant Bloch-band orbital magnetic moment,
$$
m_n^{O}(\mathbf{k})=-\frac{e}{2\hbar}\,\mathrm{Im}\,\langle \nabla_k n|\times\{H_k-\varepsilon_{nk}\}|\nabla_k n\rangle,
$$
or its modern-theory counterpart in terms of interband velocity matrix elements [2512.16387; 2504.07665]. Edge-state studies emphasize that the relevant orbital degree of freedom can be entirely inter-atomic, with the orbital angular momentum per Bloch eigenstate computed from velocity matrix elements and no intra-atomic contribution in the tight-binding basis [2411.00353].

In superconductors the susceptibility is formulated against a supercurrent source rather than an electric field. For multiorbital, inversion-breaking superconductors with orbital Rashba coupling, the static susceptibility is extracted from Bogoliubov–de Gennes theory and has the band-resolved form
$$
\chi^o_{yx}=-2\sum_k\sum_{\lambda>\lambda'}\left[1-\frac{\epsilon_\lambda\epsilon_{\lambda'}+\Delta^2}{E_\lambda E_{\lambda'}}\right]
\frac{\langle\phi_\lambda|\hat v_x|\phi_{\lambda'}\rangle\langle\phi_{\lambda'}|\hat L_y|\phi_\lambda\rangle}{E_\lambda+E_{\lambda'}},
$$
which makes explicit the role of interband matrix elements and superconducting coherence factors [2408.08151]. Closely related work on non-centrosymmetric superconductors writes the response as $\langle \hat{L}_\mu\rangle=-\frac{1}{2}\chi_{\mu\nu}q_\nu$ and finds a supercurrent-induced orbital magnetization more than one order of magnitude greater than the spin contribution [2107.07476].

## 3. Symmetry, gyrotropy, and what is actually required

Broken inversion symmetry is the recurring symmetry condition, but the literature repeatedly distinguishes inversion breaking from chirality and from spin–orbit coupling. In Cu\(_2\)WSe\(_4\), which has space group \(I\bar{4}2m\) and point group \(D_{2d}\), the orbital Edelstein tensor has the symmetry-allowed form
$$
\alpha^{\mathrm{orb}}(D_{2d})=
\begin{bmatrix}
\alpha & 0 & 0\\
0 & -\alpha & 0\\
0 & 0 & 0
\end{bmatrix},
$$
so that $\alpha_{xx}^{\mathrm{orb}}=-\alpha_{yy}^{\mathrm{orb}}\neq 0$, $\alpha_{zz}^{\mathrm{orb}}=0$, and all off-diagonal components vanish. The same work stresses that chirality is not required for $\alpha_{ij}\neq 0$; gyrotropy is the essential requirement [2512.16387].

Other settings enforce different tensor structures. In the (111) LaAlO\(_3\)/SrTiO\(_3\) interface with point group \(C_{3v}\), the in-plane susceptibility is antisymmetric, \(M\propto \hat z\times E\), so only \(\chi_{XY}\) needs to be evaluated [2207.07663]. Gated monolayer TMDs with \(C_{3v}\) symmetry likewise allow only antisymmetric in-plane off-diagonal components,
$$
\chi_{ij}=
\begin{pmatrix}
0 & \chi_{xy} & 0\\
-\chi_{xy} & 0 & 0\\
0 & 0 & 0
\end{pmatrix},
$$
so that \(M_y=\chi_{xy}E_x\) and \(M_x=-\chi_{xy}E_y\) [2509.25634]. In the analytically solvable chiral-helix model and in chiral CNTs, the response is effectively one-dimensional, and only the longitudinal component along the chiral axis is nonzero [2502.04978; 2504.07665].

In antiferromagnets, symmetry must also be resolved at the sublattice level. For CuMnAs and Mn\(_2\)Au, the orbital susceptibility decomposes naturally into staggered and uniform parts,
$$
\chi^{\mathrm{orb,stag}}_{ij}=\frac{\chi^A_{ij}-\chi^B_{ij}}{2},\qquad
\chi^{\mathrm{orb,uni}}_{ij}=\frac{\chi^A_{ij}+\chi^B_{ij}}{2},
$$
with the staggered in-plane components controlling sublattice torques and the uniform components tilting moments [1905.08279]. By contrast, edge-state work shows that even a fixed bulk Hamiltonian does not determine whether an orbital Edelstein response exists at all: zigzag-type edges exhibit an \(L_z(k_\parallel)\) texture and hence an edge susceptibility, whereas straight edges can give \(L_z(k_\parallel)\equiv0\). The paper explicitly characterizes this as an absence of bulk–boundary correspondence in the accumulation process [2411.00353].

Several common misconceptions are therefore incorrect. Spin–orbit coupling is not universally required; chirality is not universally required; and bulk orbital Berry curvature is not, by itself, a predictor of edge Edelstein accumulation. These conclusions are explicit in the antiferromagnetic, scalar-potential, Cu\(_2\)WSe\(_4\), and edge-state formulations [1905.08279; 2407.00516; 2512.16387; 2411.00353].

## 4. Representative material platforms and magnitudes

The reported susceptibility depends strongly on dimensionality, symmetry class, and choice of observable. Representative results span antiferromagnets, oxide interfaces, superconductors, chiral nanotubes, chiral helices, semiconductors, and topological materials.

| Platform | Representative response statement | Paper |
|---|---|---|
| CuMnAs, Mn\(_2\)Au | Orbital response dominates spin; in CuMnAs with moments along \(c\), \(\chi^L_{xy}(\omega=0)\approx 60\times \chi^S_{xy}(\omega=0)\) | [1905.08279] |
| AlO\(_x\)/SrTiO\(_3\) 2DEG | Orbital Edelstein effect exceeds spin Edelstein effect by more than one order of magnitude | [2006.14958] |
| (111) LaAlO\(_3\)/SrTiO\(_3\) | Orbital Edelstein effect an order of magnitude larger then the spin one | [2207.07663] |
| Non-centrosymmetric superconductors | Supercurrent-induced orbital magnetization more than one order of magnitude greater; up to approximately 25 times and approximately 60 times larger in compared cases | [2107.07476] |
| Metallic chiral CNTs | \(\chi_z^{L_z}\cdot(a/T)\) reaches up to about \(10^{-8}\ \mu_B\cdot\mathrm{m}/\mathrm{V}\) in the low-energy range and \(\approx10^{-7}\ \mu_B\cdot\mathrm{m}/\mathrm{V}\) at higher energies | [2504.07665] |
| Chiral-helix model for Te | \(\chi_0=\pm125.0\times10^{-9}\ \mu_B\,\mathrm{m}/\mathrm{V}\) spinless; \(\approx\pm250.0\times10^{-9}\ \mu_B\,\mathrm{m}/\mathrm{V}\) including spin degeneracy | [2502.04978] |

Several platform-specific features recur. In noncentrosymmetric antiferromagnets, the dominant components are staggered off-diagonal terms, and the induced local magnetization can contain both staggered in-plane components and non-staggered out-of-plane components [1905.08279]. In oxide interfaces, the enhancement is traced to multi-orbital hybridization and unequal orbital moment magnitudes within Rashba-split band pairs, rather than simply to large spin–orbit coupling [2006.14958; 2207.07663]. In bilayer Rashba models, the sign of the orbital response can be enhanced, suppressed, or reversed depending on the relation of the effective Rashba parameters of each layer, and a sign change is related to an interchange of the corresponding layer localization of the states [2307.02872].

Chiral systems furnish especially transparent scaling laws. In metallic chiral CNTs near the Fermi level,
$$
\chi_z^{L_z}(E,\theta,R)\propto R\cdot \sin(\theta/6)\cdot (E-E_F)^2,
$$
the susceptibility is an odd function of chirality angle, and it vanishes for zigzag and armchair tubes [2504.07665]. In the analytically solvable chiral-helix model,
$$
\chi_z^{L_z}(E)=\frac{e^2}{2\pi\sqrt{3}\hbar^2}\tau |t|ca^2\sqrt{1-(E/2t)^2},
$$
so the response is linear in the chirality parameter \(c\), linear in \(\tau\), and maximal near the band center [2502.04978].

Recent first-principles work extends these trends to gated monolayer TMDs, where the gate-induced broken mirror symmetry produces a Rashba-type chiral spin/orbital angular momentum texture and the predicted OEE is described as an order of magnitude larger than values reported in previously studied systems [2509.25634]. In optically driven monolayer jacutingaite, the orbital Edelstein susceptibility instead becomes a probe of a Floquet-engineered topological transition: it vanishes at the onset of the semimetallic regime because the orbital magnetic moment is proportional to the light-controlled mass term [2505.06144].

## 5. Nonlinear, reciprocal, and edge-resolved extensions

The susceptibility concept has expanded beyond linear dc bulk response. In few-layer WTe\(_2\), the reported effect is explicitly nonlinear:
$$
M_i=\chi_{ijk}E_jE_k,
$$
with the second-harmonic component under an ac drive given by
$$
M_i^{(2\omega)}=\frac{1}{2}\beta^{\mathrm{orb}}_{ijk}E_j^\omega E_k^\omega.
$$
The out-of-plane magnetization is explained by current-induced orbital magnetization via the Berry connection polarizability tensors in WTe\(_2\), and the orbital degree of freedom plays the primary role in the observed nonlinear Edelstein effect [2412.02938].

Reciprocity has also moved to the forefront. Nonlocal measurements on an orbital Edelstein system based on oxidized Cu demonstrate that the direct and inverse orbital-charge conversion processes produce identical electric voltages, confirming Onsager reciprocity. At \(d=140\ \mathrm{nm}\) with FM = Co\(_{25}\)Fe\(_{75}\) and \(t_{\mathrm{Cu}}=40\ \mathrm{nm}\),
\(2\Delta R_{\mathrm{DOEE}}\approx +0.22\ \mathrm{m}\Omega\) and
\(2\Delta R_{\mathrm{IOEE}}\approx -0.22\ \mathrm{m}\Omega\),
while the lateral orbital decay length is \(\lambda_{\mathrm{orb}}\approx100\ \mathrm{nm}\) at room temperature and is independent of Cu thickness [2502.11040]. A plausible implication is that orbital Edelstein susceptibility is not only a local constitutive parameter but also part of a nonlocal transport reciprocity structure once orbital accumulation and its conjugate chemical potential are resolved experimentally.

Edge-local susceptibility provides another extension. For finite slabs, the nonequilibrium OAM profile is written as
$$
M_x=M_x^{(0)}+\zeta_x e\tau E_y,
$$
so that the local orbital Edelstein susceptibility can be identified as \(\chi^{\mathrm{orbE}}_{zy}(x)=\zeta_x e\tau\). The effect is controlled by the existence of an \(L_z(k_\parallel)\) texture in edge-localized states and can therefore be present for zigzag or irregular terminations and absent for straight edges of the same bulk model [2411.00353].

Superconducting and driven systems show that the “source” field need not be a dc electric field. In superconductors it is the phase gradient \(q\) or the supercurrent; in Floquet-engineered jacutingaite it is a dc probe response carried by a light-modified band structure; in WTe\(_2\) it is a quadratic response to an ac current. Orbital Edelstein susceptibility is therefore best understood as a family of symmetry-allowed magnetoelectric coefficients rather than a single fixed transport number [2107.07476; 2505.06144; 2412.02938].

## 6. Experimental probes, open problems, and current research directions

The experimental literature and the theory papers converge on a set of favored probes. Magnetometry under applied electric field or current, Kerr or Faraday rotation, XMCD, torque magnetometry, NV-center magnetometry, lateral SQUID magnetoscopy, and nonlocal electrical detection all appear as viable approaches depending on whether the target is bulk magnetization, site-resolved orbital moments, or reciprocal charge–orbital conversion [2512.16387; 2107.07476; 2006.14958; 2207.07663; 2502.11040]. In few-layer WTe\(_2\), the nonlinear orbital Edelstein effect is electrically probed using an Fe\(_3\)GeTe\(_2\) electrode with perpendicular magnetic anisotropy, and the second-harmonic response scales quadratically with the applied current [2412.02938].

Several unresolved issues are explicit in the literature. For antiferromagnets, open questions include how efficiently orbital polarization transfers torque to spins, how to quantify and experimentally detect OREE versus SREE contributions at optical and THz frequencies, and how uniform components affect switching pathways [1905.08279]. In reciprocal nonlocal transport, the experiment constrains the product of interfacial charge-to-orbital conversion and orbital-to-charge detection efficiencies, but does not isolate an absolute \(\chi_{\mathrm{OE},j}\) because the detector calibration is not independently known [2502.11040]. In WTe\(_2\), the experiment establishes quadratic scaling, symmetry, and temperature dependence, but does not convert the second-harmonic voltage into an absolute material \(\chi_{ijk}\) because interfacial orbital-to-spin and spin-to-voltage conversion factors are not known [2412.02938].

A further open problem concerns the appropriate microscopic orbital variable. Some calculations use atomic-like \(L\) operators; others use gauge-invariant Bloch orbital moments; still others emphasize purely inter-site self-rotation or effective orbital moments generated by Berry connection polarizability. This suggests that “orbital Edelstein susceptibility” is not a single microscopic observable but a response class whose concrete realization depends on the effective low-energy description. A plausible implication is that cross-comparison between materials requires care not only with units but also with which orbital quantity is being transported or accumulated [2006.14958; 2512.16387; 2411.00353; 2412.02938].

Taken together, the published results establish a consistent core picture. Orbital Edelstein susceptibility is the response coefficient linking inversion-breaking electronic structure to current- or field-induced orbital polarization. Its dominant tensor components are fixed by crystal or boundary symmetry; its microscopic value depends on orbital texture, interband hybridization, and scattering; it can persist without spin–orbit coupling; and in multiple concrete systems it is reported as larger than the corresponding spin Edelstein response [1905.08279; 2107.07476; 2006.14958; 2502.04978].

Source: https://www.emergentmind.com/topics/orbital-edelstein-susceptibility