---
title: Madelung–Rashba Equations Overview
url: https://www.emergentmind.com/topics/madelung-rashba-equations
type: topic
---

# Madelung–Rashba Equations Overview

The Madelung–Rashba equations constitute a unified hydrodynamic framework for quantum systems with spin–orbit coupling (SOC), specifically incorporating the Rashba interaction. These equations extend the classical Madelung hydrodynamics of quantum mechanics to the spinorial case, yielding a coupled system for the particle density, hydrodynamic velocity, and spin vector. They arise from the Pauli equation with SOC via either phase–space kinetic expansions or variational/Hamiltonian reduction, and encode the quantum–geometric, semiclassical, and uniquely quantum correlation effects generated by spin–orbit interactions in two-dimensional electron systems and related platforms [2007.15947][2601.10698]. Notably, the Madelung–Rashba system captures the interplay of drift–diffusion, quantum pressure, semiclassical forces, geometric torques, and quantum spin–orbit correlations stemming from the Rashba term, and forms the basis for both analytic understanding and variationally exact numerical simulation.

## 1. Mathematical Formulation of Madelung–Rashba Equations

The starting point is the Pauli equation with spin–orbit interaction:
\[
H_{\mathrm{Rashba}} = -\frac{\hbar^2}{2m^*}\Delta \,\sigma_0 - i\hbar\alpha_R\,\nabla^\perp \cdot \vec{\sigma}
\]
where $\sigma_0$ is the identity, $\vec{\sigma}$ the vector of Pauli matrices, $m^*$ the effective mass, and $\alpha_R$ the Rashba coupling constant.

Applying the Madelung transform to the Pauli spinor $\Psi(\mathbf{x}, t)$,
\[
\Psi(\mathbf{x}, t) = \sqrt{D(\mathbf{x}, t)}\, e^{i S(\mathbf{x}, t)/\hbar} \chi(\mathbf{x}, t) ,
\]
with $D = \Psi^\dagger \Psi$ the probability density and $\chi$ the normalized local spinor, the hydrodynamic fields are:

- $D(\mathbf{x}, t)$: scalar particle (orbital) density
- $\mathbf{v}(\mathbf{x}, t) = \frac{1}{m}(\nabla S + \mathbf{A}) + \langle \chi, \widehat{\mathbf{X}} \chi \rangle$: velocity field with $\mathbf{A}$ as the Berry connection and $\widehat{\mathbf{X}}$ the SOC vector operator
- $\mathbf{s} = \chi^\dagger \vec{\sigma} \chi/2$: local spin vector

The full Madelung–Rashba hydrodynamic system for planar SOC is [2601.10698]:

\[
\begin{align*}
&\text{Continuity:} && \partial_t D + \nabla \cdot (D \mathbf{v}) = 0 \\
&\text{Euler:} && mD (\partial_t + \mathbf{v} \cdot \nabla) \mathbf{v} = - D \nabla V - D \nabla V_Q - m\,\mathrm{Tr}[\nabla \widehat{\mathbf{X}}\cdot J] - \partial_j \left( \frac{\hbar^2}{2m^2} D \,\mathrm{Tr}(\partial_j \hat{\rho} \nabla \hat{\rho}) \right) \\
&\text{Spin evolution:} && (\partial_t + \mathbf{v} \cdot \nabla) \mathbf{s} = \mathbf{T}_{\mathrm{sc}} + \mathbf{T}_{\mathrm{QGT}} + \mathbf{T}_{\mathrm{MCO}}
\end{align*}
\]
with $V_Q = - \frac{\hbar^2}{2m}\frac{\Delta \sqrt{D}}{\sqrt{D}}$ the Madelung–Bohm quantum potential, $\mathbf{T}_{\mathrm{sc}}$ the semiclassical SOC torque, $\mathbf{T}_{\mathrm{QGT}}$ the quantum–metric torque, and $\mathbf{T}_{\mathrm{MCO}}$ a new Mead current operator (MCO) torque.

## 2. Quantum Drift–Diffusion and Kinetic Derivation

In the phase–space (Wigner function) approach [2007.15947], the density-matrix evolution is cast into the Wigner representation and expanded via the Chapman–Enskog method with a BGK-type relaxation toward the quantum maximum entropy equilibrium. The Wigner function $w(x, p, t) = w_0(x, p, t) \sigma_0 + \vec{w}(x, p, t) \cdot \vec{\sigma}$ yields:

- $n_0(x, t) = \int w_0(x, p, t) dp$ (scalar density)
- $S(x, t) = \int \vec{w}(x, p, t) dp$ (spin-density vector)

Projecting the kinetic equation onto these moments produces systems of the form
\[
\partial_t n = T g + \tau T T g - \tau \langle T \delta g/\delta n \rangle \odot T g
\]
with the Wigner-transport operator $T$ encoding both drift and Rashba-driven terms. A crucial feature is the presence of a non-vanishing current at $\mathcal{O}(1)$ in $\tau$ due to SOC, an effect not present in spinless or axis-aligned spin models.

Explicitly, for the charge continuity:
\[
\partial_t n_0(x, t) + \nabla \cdot \mathbf{J}_c(x, t) = 0
\]
with
\[
\mathbf{J}_c = -\tau \left[ n_0 \nabla(V + a_0) + S \cdot \nabla \vec{a} \right] + 2 \tau \alpha \nabla^\perp \cdot (\vec{a} \times S) + \mathcal{O}(\tau^2)
\]
where $a_0$ and $\vec{a}$ are Lagrange multipliers, and the Rashba SOC enters through $\alpha$ as well as direct cross-terms in the fluxes.

The corresponding spin-density equation includes leading-order spin torque $Tg = \frac{1}{2}\vec{a} \times S$ and spin currents with further spin–orbit and quantum geometric structure.

## 3. Quantum–Geometric and Correlation Effects

Variational reduction of the Pauli action exposes two principal forms of quantum–geometric forces induced by SOC [2601.10698]:

1. **Quantum Geometric Tensor (QGT) Forces:** Terms proportional to $\frac{\hbar^2}{4m}\|\nabla \hat{\rho}\|^2$ (with $\hat{\rho}$ the spinor density matrix) generate quantum–metric forces and torques, including diffusive spin precession and corrections to the quantum pressure.

2. **Mead Current Operator (MCO) Forces:** Terms of the form $\frac{\hbar}{2}\mathrm{Tr}(i \widehat{\mathbf{X}}\cdot [\hat{\rho}, \nabla \hat{\rho}])$ yield strictly quantum, non-factorizing spin–orbit correlation forces and novel torques. In the hydrodynamic equations, these manifest as terms with a single derivative of the spin field, distinct from the classical and QGT contributions.

This structure leads to three evolutionary mechanisms for the spin vector:
- Semiclassical SOC torque ($\propto m \mathbf{F} \times \mathbf{s}$)
- Quantum geometric metric torque ($\propto -\frac{\hbar^2}{2mD} \nabla \cdot (D \nabla \mathbf{s})$)
- MCO quantum correlation torque (proportional to $\nabla \times \mathbf{s}$ and directly encoding spin–orbit correlations)

## 4. Structural Properties and Lie–Poisson Formulation

The full Madelung–Rashba system possesses a Lie–Poisson Hamiltonian structure [2601.10698]:
\[
h[\mathbf{M}, D, \tilde{\mathbf{s}}] = \int \left[ \frac{|\mathbf{M}|^2}{2 m D} + \alpha\,\mathbf{M} \cdot (\mathbf{e}_z \times \mathbf{s}) + \frac{\hbar^2}{8m} \frac{|\nabla D|^2}{D} + D\,\mathcal{E}(\tilde{\mathbf{s}}, \nabla \tilde{\mathbf{s}}) \right] d^2 x
\]
for momentum density $\mathbf{M} = D \mathbf{v}$ and spin density $\tilde{\mathbf{s}} = D \mathbf{s}$.

This semidirect product structure ensures
- conservation of total energy,
- explicit identification of Casimir invariants $C_\Phi = \int D\,\Phi(D^{-1}\tilde{\mathbf{s}}) d^2x$ for arbitrary scalar functions $\Phi$,
- the presence of hydrodynamic transport, precessional, and geometric flow properties inherent to the spin–orbit coupled system.

## 5. Specialization to Two-Dimensional Electron Gases and Planar SOC

For planar Rashba coupling $H_R = \alpha(p_y \sigma_x - p_x \sigma_y)$, the hydrodynamic velocity reduces to:
\[
\mathbf{v} = \frac{1}{m}\nabla S + \alpha\, \mathbf{e}_z \times \mathbf{s}
\]
and the system further simplifies to explicit $2$D forms for continuity, Euler, and spin–evolution equations, including:
\[
\partial_t D + \nabla \cdot [ D( (1/m) \nabla S + \alpha \mathbf{e}_z \times \mathbf{s} ) ] = 0
\]
with new terms arising in the spin equation from the MCO, which have no classical analogue and encode the quantum–geometric SOC correlation dynamics.

Drift–diffusion approximations of this system produce coupled equations for the particle and spin densities (including full spin vector, not just projections), suitable for modeling two-dimensional electron gases with Rashba SOC, and generalize previous treatments that were either semiclassical or only included axis-projected spin degrees of freedom [2007.15947].

## 6. Numerical Implementation: Bohmion Particle Scheme

A particle-based (bohmion) discretization enables faithful and Hamiltonian-preserving simulation [2601.10698]. Key steps include:

- Mollification of density and spin densities with a regularizing kernel $K$, e.g. (Gaussian, width $\Delta$)
- Singular ansatz: $D(\mathbf{x}, t) = \sum_a w_a \delta(\mathbf{x} - \mathbf{x}_a (t))$, $\tilde{\mathbf{s}}(\mathbf{x}, t) = \sum_a w_a \mathbf{s}_a (t) \delta(\mathbf{x} - \mathbf{x}_a (t))$
- Insertion into a regularized action yields finite-dimensional Hamiltonian ODEs for particle positions $\mathbf{x}_a$, momenta $\mathbf{p}_a$, and spins $\mathbf{s}_a$
- Algorithmic steps (pseudo-code outlined in [2601.10698]) maintain exact conservation of discrete Hamiltonian and Casimirs

This numerical framework leverages the variational foundation to robustly model spin–orbit coupled hydrodynamics, especially relevant for quantum fluids and twodimensional electron gases.

## 7. Physical Interpretation and Applications

The Madelung–Rashba equations provide a macroscopic, yet fully quantum, description of spin–orbit coupled dynamics in low-dimensional systems. They resolve both the transport and geometric torque effects driving phenomenology such as the spin Hall effect, persistent spin helix, and quantum geometrically-induced spin precession [2601.10698][2007.15947]. The explicit inclusion of quantum–geometric (QGT) and MCO correlation forces distinguishes this framework from drift–diffusion and semiclassical models, allowing the study of strictly quantum correlation phenomena in SOC systems. Analysis of these equations supports theoretical and computational efforts in topological spintronics, quantum transport, and hydrodynamic electron flow in materials with strong SOC.

**References:**  
L. Barletti, P. Holzinger, A. Jüngel, "Quantum drift-diffusion equations for a two-dimensional electron gas with spin-orbit interaction" [2007.15947];  
G. Tronci et al., "Madelung hydrodynamics of spin-orbit coupling: action principles, currents, and correlations" [2601.10698].

Source: https://www.emergentmind.com/topics/madelung-rashba-equations