Nonlinear Magnetoelectric Edelstein Effect

• The nonlinear magnetoelectric Edelstein effect is a second-order response where electric fields induce spin or orbital magnetization through mechanisms like Berry connection polarizability and orbital Rashba textures.
• It encompasses multiple response forms including quadratic electric-field effects, optical rectification, and bilinear electric–magnetic couplings, each governed by symmetry constraints and spin–orbit interactions.
• Experimental platforms such as few-layer WTe2 and chiral CoSi validate the effect by detecting second-harmonic signals and orbital magnetization, offering insights into designing advanced spintronic devices.

The nonlinear magnetoelectric Edelstein effect denotes a family of nonequilibrium responses in which electric driving generates spin or orbital magnetization beyond the conventional linear Edelstein regime. In current usage, the term covers at least three closely related structures: a quadratic electric-field response $M_i=\chi_{ijk}E_jE_k$, an optical rectification response $M_i\propto E(\omega)E(-\omega)$, and a bilinear electric–magnetic response $\delta s^\alpha\propto E_\beta B_\gamma$ generated by the combined action of electric and magnetic fields (Xu et al., 2020, Ye et al., 2024, Jia et al., 31 Jul 2025). Across these variants, the central problem is the same: how broken inversion, spin–orbit coupling, hidden or explicit magnetic order, orbital texture, and dissipation convert charge-sector driving into nonequilibrium magnetic polarization.

1. Terminology, historical development, and scope

The linear Edelstein effect is usually written as $M_i=\alpha_{ij}E_j$ or, equivalently, as current-induced spin polarization in a noncentrosymmetric conductor. A foundational nonlinear extension was developed for a clean Rashba two-dimensional electron gas, where the response was treated nonperturbatively in the electric field and controlled by the parameter $\gamma=eEL_s/E_F$; that work showed an adiabatic regime at small $\gamma$ and a strongly nonadiabatic suppression of spin polarization at large $\gamma$ (Vignale et al., 2015). Subsequent literature broadened the notion of nonlinear Edelstein physics to include optically generated static magnetization, orbital as well as spin magnetization, and transport regimes in which the leading nonvanishing response is second order in the field (Xu et al., 2020).

A useful taxonomy is:

Response class	Constitutive form	Typical interpretation
Linear Edelstein effect	$M_i=\alpha_{ij}E_j$	Current-induced spin or orbital magnetization
Nonlinear Edelstein effect	$M_i=\chi_{ijk}E_jE_k$	Second-order electric or optical generation of magnetization
Nonlinear magnetoelectric Edelstein effect	$$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$	Bilinear $M_i\propto E(\omega)E(-\omega)$0-$M_i\propto E(\omega)E(-\omega)$1 spin magnetization

The terminology is not fully uniform. Some works reserve “Edelstein” for spin density and treat orbital responses separately, whereas others compute total magnetization and explicitly decompose it into spin and orbital parts. Some papers use “nonlinear Edelstein effect” for purely electric $M_i\propto E(\omega)E(-\omega)$2 responses, while a later formulation uses “nonlinear magnetoelectric Edelstein effect” for the mixed $M_i\propto E(\omega)E(-\omega)$3 response with Zeeman coupling (Baek et al., 2023, Xu et al., 17 Jan 2025, Jia et al., 31 Jul 2025). This terminological spread reflects genuine physical variety rather than simple naming inconsistency.

A second historical thread comes from optics. “Light-Induced Static Magnetization: Nonlinear Edelstein Effect” formulated a dc magnetization induced by optical fields through

$M_i\propto E(\omega)E(-\omega)$4

and emphasized that both spin and orbital magnetization can be induced under linearly or circularly polarized light (Xu et al., 2020). Later work connected this optical formulation to first-principles calculations in chiral semimetals, correlated Hubbard systems, centrosymmetric metals, and nonlinear spin–charge conversion at magnetic interfaces (Xu et al., 17 Jan 2025, Ōiké et al., 2024, Baek et al., 2023, Trama et al., 2024).

2. Symmetry structure and response tensors

Symmetry is the organizing principle of the subject. In the broad electric-quadratic form,

$M_i\propto E(\omega)E(-\omega)$5

the output $M_i\propto E(\omega)E(-\omega)$6 is an axial vector, while $M_i\propto E(\omega)E(-\omega)$7 is inversion even. This is the basic reason why a nonlinear Edelstein response can exist in centrosymmetric systems even when the linear Edelstein tensor vanishes. In centrosymmetric metals, the second-order orbital or spin response is therefore symmetry-allowed provided rotational symmetries such as $M_i\propto E(\omega)E(-\omega)$8 or $M_i\propto E(\omega)E(-\omega)$9 do not enforce cancellation (Baek et al., 2023).

The role of time reversal is subtler and depends on context. In optical NLEE, one formulation argues that linearly polarized light can induce magnetization even in nonmagnetic materials because absorption and dissipation render the driven steady state effectively irreversible; in the clean limit $\delta s^\alpha\propto E_\beta B_\gamma$0, the linearly polarized response in a $\delta s^\alpha\propto E_\beta B_\gamma$1-symmetric system vanishes (Xu et al., 2020). In a first-principles study of CoSi, the authors similarly argued that both linear and nonlinear Edelstein responses require time-reversal breaking, but that in nonmagnetic materials such breaking can be supplied by heat and dissipation in the driven state; within the same study, the linear tensor is inversion odd and reverses sign between enantiomers, whereas the nonlinear tensor is inversion even and has the same sign in both enantiomers (Xu et al., 17 Jan 2025). By contrast, the mixed $\delta s^\alpha\propto E_\beta B_\gamma$2-$\delta s^\alpha\propto E_\beta B_\gamma$3 NMEE was introduced precisely because the intrinsic components of ordinary linear and quadratic Edelstein response are generally forbidden in $\delta s^\alpha\propto E_\beta B_\gamma$4-invariant or $\delta s^\alpha\propto E_\beta B_\gamma$5-symmetric systems, while the bilinear magnetoelectric variant can remain intrinsic in $\delta s^\alpha\propto E_\beta B_\gamma$6-invariant but inversion-breaking materials, including insulators (Jia et al., 31 Jul 2025).

Specific point-group constraints are often decisive. In few-layer $\delta s^\alpha\propto E_\beta B_\gamma$7-phase WTe$\delta s^\alpha\propto E_\beta B_\gamma$8, mirror and glide symmetries enforce

$\delta s^\alpha\propto E_\beta B_\gamma$9

so an out-of-plane magnetization $M_i=\alpha_{ij}E_j$0 is allowed only for an in-plane electric field tilted away from the crystalline $M_i=\alpha_{ij}E_j$1 and $M_i=\alpha_{ij}E_j$2 axes (Ye et al., 2024). In magnetic-octupolar antiferromagnets, the nonlinear magnetoelectric tensor is written as

$M_i=\alpha_{ij}E_j$3

with $M_i=\alpha_{ij}E_j$4 a rank-3, $M_i=\alpha_{ij}E_j$5-odd axial tensor transforming like a magnetic octupole; this is why the effect is symmetry-selected in magnetic point groups with lowest-rank octupolar order (Ōiké et al., 2024). At altermagnetic interfaces, the out-of-plane nonlinear Edelstein response

$M_i=\alpha_{ij}E_j$6

is allowed even though the corresponding linear out-of-plane Edelstein coefficient vanishes by $M_i=\alpha_{ij}E_j$7 symmetry (Trama et al., 2024).

3. Microscopic mechanisms

A major theme of the literature is that nonlinear Edelstein physics is often orbital-first rather than spin-first. In few-layer WTe$M_i=\alpha_{ij}E_j$8, the experimentally observed quadratic magnetization is interpreted primarily as a nonlinear orbital Edelstein effect governed by the Berry connection polarizability tensor $M_i=\alpha_{ij}E_j$9. The applied field induces a positional shift $\gamma=eEL_s/E_F$0, which produces a field-induced orbital moment texture $\gamma=eEL_s/E_F$1; under the transport-induced Fermi-surface shift, that texture yields a net orbital magnetization $\gamma=eEL_s/E_F$2, and spin polarization follows secondarily through spin–orbit coupling (Ye et al., 2024).

The same hierarchy appears in centrosymmetric metals. There, a first power of the electric field generates an “electric-field-induced orbital Rashba texture,” summarized as

$\gamma=eEL_s/E_F$3

A second electric-field factor then shifts the nonequilibrium distribution on this emergent orbital-Rashba-textured Fermi surface, producing a nonlinear orbital Edelstein effect; spin density appears only after orbital-to-spin conversion by $\gamma=eEL_s/E_F$4 coupling. In this formulation, NOEE survives even without SOC, whereas NSEE vanishes in the zero-SOC limit (Baek et al., 2023).

A distinct intrinsic route is provided by quantum geometry. In the magnetic-octupole formulation, the intrinsic nonlinear magnetoelectric tensor contains derivatives of the Berry-connection polarizability $\gamma=eEL_s/E_F$5 and the mixed spin–momentum object $\gamma=eEL_s/E_F$6,

$\gamma=eEL_s/E_F$7

In a pyrochlore lattice with all-in/all-out order, the intrinsic part is strongly enhanced near Weyl points, with an asymptotic logarithmic growth $\gamma=eEL_s/E_F$8 as the chemical potential approaches the Weyl node (Ōiké et al., 2024).

Another mechanism is purely dynamical. In odd-wave magnets, the steady-state Edelstein magnetization is not a weak-field linear law but

$\gamma=eEL_s/E_F$9

This implies linear response for $\gamma$0, a maximum at $\gamma$1, and a strong-field decay $\gamma$2 associated with Wannier–Stark localization (Habel et al., 29 Jun 2026). Here the nonlinearity is generated by Bloch-periodic band motion plus relaxation rather than by a small set of symmetry-allowed tensor coefficients.

Correlation effects supply yet another layer. In a Hubbard-interacting Kane–Mele model solved with DMFT, the dc current-induced nonlinear Edelstein coefficient is enhanced approximately as $\gamma$3 by the real part of the self-energy, while the imaginary part suppresses dissipative optical channels. The result is a clean separation between renormalization-driven enhancement and lifetime-driven suppression (Ōiké et al., 2024).

4. Electronic-structure platforms and model systems

The current literature spans a broad set of electronic environments. Few-layer WTe$\gamma$4 is the clearest experimental platform for a quadratic current-induced out-of-plane magnetization: the response is symmetry-selected, detected electrically with a perpendicular Fe$\gamma$5GeTe$\gamma$6 probe, and interpreted through Berry-connection-driven orbital magnetization (Ye et al., 2024). Chiral CoSi provides a first-principles platform in which both linear and nonlinear Edelstein effects were computed for spin and orbital magnetization. In that study, cubic symmetry leaves only one independent linear tensor element and one independent symmetric and antisymmetric nonlinear optical element,

$\gamma$7

with optically induced moments estimated to reach $\gamma$8 for fields of order $\gamma$9, subject to the perturbative caveat that nonperturbative methods are needed when external fields approach internal atomic fields (Xu et al., 17 Jan 2025).

Centrosymmetric metals constitute a second major class. Tight-binding and first-principles studies of orthorhombic metals such as IrW, MoRh, MoIr, and TiPt, together with (110) films of Pt, Pd, Ir, and Rh, found nonlinear orbital and spin Edelstein densities of roughly $\gamma$0 for $\gamma$1, $\gamma$2, and $\gamma$3. These calculations are notable because they do not require either global or local inversion breaking to generate the orbital precursor (Baek et al., 2023).

Magnetic systems add qualitatively new possibilities. At altermagnetic interfaces, the coexistence of $\gamma$4-wave altermagnetic exchange splitting and interfacial Rashba SOC produces chiral spin textures and an anomalous out-of-plane Edelstein response, both as a mixed $\gamma$5 term and as a purely electric nonlinear term $\gamma$6; a continuum model fitted to RuO$\gamma$7 bilayers was used as the concrete realization (Trama et al., 2024). In magnetic-octupole systems, model calculations were carried out for a $\gamma$8-wave altermagnet and a pyrochlore lattice with all-in/all-out order, where the nonlinear magnetoelectric tensor functions as a probe of octupolar symmetry (Ōiké et al., 2024). In odd-wave magnets, the nonlinear Edelstein response arises without SOC from exchange-driven odd-parity spin textures (Habel et al., 29 Jun 2026).

The bilinear $\gamma$9-$M_i=\alpha_{ij}E_j$0 NMEE has so far been established theoretically through a two-band Dirac model and a tight-binding model on a honeycomb lattice. Its distinctive claim is that the intrinsic $M_i=\alpha_{ij}E_j$1 spin magnetization can remain finite in $M_i=\alpha_{ij}E_j$2-invariant, inversion-breaking systems, including insulators, whereas the extrinsic component is proposed as a detector of Néel-vector reversal in $M_i=\alpha_{ij}E_j$3-symmetric antiferromagnets (Jia et al., 31 Jul 2025).

5. Experimental manifestations and detection

The experimentally established reference point is few-layer WTe$M_i=\alpha_{ij}E_j$4. Under an ac current at frequency $M_i=\alpha_{ij}E_j$5, an out-of-plane magnetization with second-harmonic response at $M_i=\alpha_{ij}E_j$6 was detected using a perpendicular Fe$M_i=\alpha_{ij}E_j$7GeTe$M_i=\alpha_{ij}E_j$8 electrode separated by h-BN. The hysteretic second-harmonic voltage $M_i=\alpha_{ij}E_j$9 scaled quadratically with the applied field, $M_i=\chi_{ijk}E_jE_k$0, and the signal did not reverse sign under current reversal, as expected for a second-order effect. The response was strong for oblique current injection, strongly suppressed for current near the $M_i=\chi_{ijk}E_jE_k$1 axis, and became barely discernible above about $M_i=\chi_{ijk}E_jE_k$2 (Ye et al., 2024).

Optical control has become a second experimental frontier. In ferromagnetic Cr$M_i=\chi_{ijk}E_jE_k$3Ge$M_i=\chi_{ijk}E_jE_k$4Te$M_i=\chi_{ijk}E_jE_k$5, time-domain THz emission spectroscopy under 1.2 eV, 160 fs excitation was interpreted as direct evidence for a light-driven nonlinear Edelstein effect producing an Edelstein–Zeeman field. The emitted THz waveform was modeled as magnetic dipole radiation from ultrafast magnetization dynamics, with polarization, fluence, and temperature dependences captured by a weakly anisotropic Heisenberg ferromagnet subject to the optically induced exchange field (Lv et al., 5 Mar 2026).

Several additional detection routes have been proposed. In magnetic-octupole systems, the induced spin density was estimated to lie in the $M_i=\chi_{ijk}E_jE_k$6 range for $M_i=\chi_{ijk}E_jE_k$7, increasing to $M_i=\chi_{ijk}E_jE_k$8 for terahertz fields exceeding $M_i=\chi_{ijk}E_jE_k$9; these values were argued to be accessible to magneto-optical Kerr measurements (Ōiké et al., 2024). In odd-wave magnets, higher-harmonic generation in THz sub-cycle lightwave spectroscopy, read out by time-resolved Kerr or Faraday rotation, was proposed as the signature of magnetic Bloch oscillations and of the associated nonlinear Edelstein magnetization (Habel et al., 29 Jun 2026). At altermagnetic interfaces, the finite-frequency extension of $$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$0 yields both rectification and second-harmonic generation of the out-of-plane spin accumulation (Trama et al., 2024).

The interpretation of experiments requires careful control of competing mechanisms. The WTe$$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$1 study explicitly argued against contact nonlinearity, magnetic proximity, nonlinear Hall accumulation, and spin-Seebeck scenarios by using h-BN spacers, nonmagnetic control electrodes, and angular systematics inconsistent with thermal or edge-accumulation explanations (Ye et al., 2024). Similar caution is likely essential for any nonlinear Edelstein measurement because even-in-current magnetization, rectified Hall signals, nonlinear optical torques, and thermally driven spin signals can coexist.

6. Conceptual issues, related phenomena, and open directions

Two conceptual issues dominate current discussion. The first is the spin–orbital hierarchy. Across WTe$$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$2, centrosymmetric metals, and CoSi, the orbital channel is often primary: orbital magnetization or orbital density is generated directly by field-induced band geometry, while spin follows through SOC (Ye et al., 2024, Baek et al., 2023, Xu et al., 17 Jan 2025). This makes “nonlinear spin Edelstein effect” an incomplete description in many materials. The second issue is the distinction between intrinsic and extrinsic response. Some theories are Fermi-surface and $$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$3-dependent, some are band-geometric and persist in insulators, and some mix both sectors. The mixed $$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$4-$$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$5 NMEE is notable because it was introduced precisely to recover an intrinsic contribution in $$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$6-invariant inversion-breaking systems where intrinsic ordinary Edelstein response is forbidden (Jia et al., 31 Jul 2025).

The relationship to neighboring magnetoelectric phenomena is also active. Optical NLEE overlaps conceptually with the inverse Faraday effect and inverse Cotton–Mouton effect, but it is broader because it includes linearly polarized-light-induced magnetization and does not require preexisting magnetic order in the same way (Xu et al., 2020). In correlated systems, optical spin injection and the inverse Faraday effect emerge as different sectors of the same second-order spin susceptibility (Ōiké et al., 2024). By contrast, superconducting Edelstein physics in gyrotropic superconductors is usually formulated as a linear supercurrent-to-magnetization tensor $$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$7, while the magnetoelectric Andreev effect is odd in supercurrent and probes current-induced non-unitary triplet pairing rather than a quadratic electric response (He et al., 2019, Tkachov, 2016).

A further open direction is geometric and mesoscopic design. A strictly linear study of a V-shaped chain showed that local inversion breaking induced by geometry can generate effective SOC concentrated at an apex, with the geometric contribution entering as a $$\delta s^\alpha=\mu_B E_\beta \bar B_\gamma[\Gamma^{\rm in}_{\alpha\gamma,\beta}+\tau\Gamma^{\rm ext}_{\alpha\gamma,\beta}]$$8-matrix associated with local translational-symmetry breaking (Kanda et al., 18 Nov 2025). This suggests, as a plausible implication, that corners, bends, and other local geometric singularities may become useful control parameters for nonlinear Edelstein coefficients as well, although that extrapolation has not yet been computed explicitly.

At present, the field is converging on a more general picture: nonlinear magnetoelectric Edelstein phenomena are not a single mechanism but a response family. Depending on symmetry and frequency regime, they may originate from Berry connection polarizability, quantum metric, hidden orbital Rashba textures, exchange-generated odd-parity spin structures, coherent interband optical processes, or mixed electric–magnetic couplings. What unifies them is the conversion of electric-sector driving into nonequilibrium magnetization through symmetry-lowered spin–orbital electronic structure.