2025-11-18T15:07:13.798459

Photogalvanic effect in hydrodynamic flows of nonreciprocal electron liquids

Kirkinis, Bonds, Levchenko et al.
We study nonlinear hydrodynamic electron transport driven by an AC electric field. In noncentrosymmetric conductors with broken time-reversal (TR) symmetry the nonlinear flow of such liquids is nonreciprocal, giving rise to a DC current $I^{DC}$ that is quadratic in the amplitude of the AC electric field. This is the hydrodynamic analogue of the linear photogalvanic effect (PGE), which arises in bulk noncentrosymmetric materials with broken TR symmetry. The magnitude of $I^{DC}$ depends on both the properties of the electron fluid and the geometry of the flow, and may be characterized by two dimensionless parameters: the nonreciprocity number $\mathcal{N}$, and the frequency-dependent vibrational number $\mathcal{R}$. Due to nonlocality of hydrodynamic transport, at low frequencies of the AC drive, $I^{DC}$ is super-extensive. The AC component of the electric current is likewise strongly affected by nonreciprocity: the hysteretic current-voltage dependence becomes skewed, which can be interpreted in terms of nonreciprocity of the memory retention time.
academic

Photogalvanic effect in hydrodynamic flows of nonreciprocal electron liquids

Basic Information

  • Paper ID: 2507.06933
  • Title: Photogalvanic effect in hydrodynamic flows of nonreciprocal electron liquids
  • Authors: E. Kirkinis, L. Bonds, A. Levchenko, A. V. Andreev
  • Classification: cond-mat.mes-hall (Condensed Matter Physics - Mesoscopic and Hall Effects)
  • Publication Date: September 1, 2025
  • Paper Link: https://arxiv.org/abs/2507.06933

Abstract

This paper investigates nonlinear hydrodynamic electron transport phenomena driven by alternating electric fields. In non-centrosymmetric conductors with broken time-reversal symmetry, the nonlinear flow of such liquids exhibits nonreciprocal behavior, generating a direct current IDCI_{DC} proportional to the square of the alternating field amplitude. This represents a hydrodynamic analog of the linear photogalvanic effect (PGE), which appears in non-centrosymmetric materials with broken time-reversal symmetry. The magnitude of IDCI_{DC} depends on both the properties of the electron fluid and the geometry of the flow, characterized by two dimensionless parameters: the nonreciprocal number N\mathcal{N} and the frequency-dependent oscillation number R\mathcal{R}. Due to the nonlocal nature of hydrodynamic transport, under low-frequency alternating driving, IDCI_{DC} exhibits superextensive behavior. The alternating component of the current is similarly strongly affected by nonreciprocity: the hysteretic current-voltage dependence becomes tilted, which can be explained by the nonreciprocal nature of the memory retention time.

Research Background and Motivation

Problem Background

  1. Rise of hydrodynamic electron transport: In solid-state systems where the spatiotemporal scales of electron momentum and energy relaxation far exceed those of electron-electron scattering that conserves momentum and energy, electron transport can be described using hydrodynamic methods.
  2. Importance of nonreciprocal effects: In non-centrosymmetric conductors with broken time-reversal symmetry, the Poiseuille flow of electron liquids exhibits strong nonreciprocal nonlinear two-terminal conductance.
  3. Need for experimental verification: Recent experiments in graphene devices have confirmed various frequency-dependent nonlinear current rectification effects, spanning the diffusive-hydrodynamic transition region.

Research Motivation

This study explores an intrinsic hydrodynamic current rectification mechanism, where the required breaking of inversion symmetry and time-reversal symmetry originates from the electron liquid itself rather than from device geometry. This distinguishes it from earlier hydrodynamic approaches based on spatially modulated 2D systems.

Core Contributions

  1. Established generalized Navier-Stokes equations for nonreciprocal electron liquids, describing the evolution of time-dependent flows.
  2. Discovered the hydrodynamic photogalvanic effect: generation of a direct current proportional to the square of the driving field amplitude under monochromatic alternating driving.
  3. Introduced two key dimensionless parameters:
    • Nonreciprocal number N=NU0\mathcal{N} = NU_0
    • Oscillation number R=ωd2/ν\mathcal{R} = \omega d^2/\nu
  4. Revealed superextensive behavior: under low-frequency driving, the direct current IDCd5I_{DC} \propto d^5, resulting from the nonlocal nature of hydrodynamic transport.
  5. Discovered memory effects: nonreciprocity affects the oscillating component of the current, leading to tilting of the current-voltage hysteresis dependence.

Methodology Details

Problem Definition

Investigation of nonlinear hydrodynamic transport in nonreciprocal electron liquids driven by oscillating electric fields in Hall bar geometry, with particular focus on the generated direct current and memory effects.

Theoretical Framework

Fundamental Equations

The evolution equation for electron momentum density p=ρv\mathbf{p} = \rho\mathbf{v} follows Newton's second law:

pt=Π^+enE\frac{\partial \mathbf{p}}{\partial t} = -\nabla \cdot \hat{\Pi} + en\mathbf{E}

where the momentum flux tensor density is: Π^=Pδij+ρvivjσij\hat{\Pi} = P\delta_{ij} + \rho v_i v_j - \sigma_{ij}

Nonreciprocal Stress Tensor

The nonreciprocal correction δσij\delta\sigma_{ij} appears only when both time-reversal and inversion symmetries are broken:

δσij=ηNijklmvklvm\delta\sigma_{ij} = \eta N_{ijklm} v_k \partial_l v_m

For Hall bar geometry, this simplifies to: δσxy=ηNuyu\delta\sigma_{xy} = \eta N u \partial_y u

Dimensionless Equations

Introducing characteristic scales U0=enEd2/ηU_0 = enEd^2/\eta and I0=23(en)2Ed3/ηI_0 = \frac{2}{3}(en)^2Ed^3/\eta, the dimensionless momentum equation becomes:

Ruτ=cosτ+y[(1+Nu)uy]\mathcal{R}\frac{\partial u}{\partial \tau} = \cos \tau + \frac{\partial}{\partial y}[(1 + \mathcal{N}u)\frac{\partial u}{\partial y}]

Technical Innovations

  1. Perturbative expansion solution: expressing the velocity field as: u(y,t)={u0AC(y)eiωt+N[u1DC(y)+u1AC(y)e2iωt]}+O(N2)u(y,t) = \Re\{u_0^{AC}(y)e^{-i\omega t} + \mathcal{N}[u_1^{DC}(y) + u_1^{AC}(y)e^{-2i\omega t}]\} + O(\mathcal{N}^2)
  2. Introduction of structure factors: defining structure factors S±(x)=sinhx±sinxcoshx+cosxS_\pm(x) = \frac{\sinh x \pm \sin x}{\cosh x + \cos x} to characterize frequency dependence.
  3. Multi-scale analysis: employing penetration depth δ=2ν/ω\delta = \sqrt{2\nu/\omega} and complex wavenumber k=(1+i)/δk = (1+i)/\delta to handle oscillatory flows.

Experimental Setup

Physical Parameters

  • Hall bar width: d1d \sim 1 μm
  • Kinematic viscosity: ν0.1\nu \sim 0.1 m²/s
  • Driving frequency: terahertz scale (R10R \sim 10)
  • Nonreciprocal parameters: vector-type N=αBN = \alpha B, tensor-type N=rcos(3θ)N = |r|\cos(3\theta)

Numerical Verification

Time-dependent Stokes equation (12) is solved numerically to verify perturbative results, with initial condition u(y,0)=110cos(πy/2)u(y,0) = \frac{1}{10}\cos(\pi y/2).

Experimental Results

Main Results

Direct Current Expression

IDC=I03N8[2R2+2R5/2S+(2R)]I_{DC} = I_0 \frac{3\mathcal{N}}{8}\left[-\frac{2}{\mathcal{R}^2} + \frac{\sqrt{2}}{\mathcal{R}^{5/2}}S_+(\sqrt{2}\mathcal{R})\right]

Frequency-Dependent Asymptotic Behavior

  • Low-frequency limit (R1\mathcal{R} \ll 1): IDCI0N/10I_{DC} \approx -I_0\mathcal{N}/10
  • High-frequency limit (R1\mathcal{R} \gg 1): IDC3I0N/(4R2)I_{DC} \approx -3I_0\mathcal{N}/(4\mathcal{R}^2)

Superextensive Behavior

At low frequencies, the direct current exhibits superextensive behavior: IDCN(en)3E2d515η2I_{DC} \approx -\mathcal{N}\frac{(en)^3E^2d^5}{15\eta^2}

Numerical Verification

Figure 3 demonstrates good agreement between perturbative theory results and complete numerical solutions, validating the theoretical analysis.

Memory Effects

  • Memory timescale defined as τMd2/(2ν)\tau_M \sim d^2/(2\nu)
  • Hysteresis loop area reaches maximum at R=2\mathcal{R} = 2
  • Nonreciprocity causes tilting and area increase of hysteresis loops

Hydrodynamic Electron Transport

  • Gurzhi's pioneering work established hydrodynamic theory in solids
  • Experimentally verified in graphene, strongly correlated electron systems, etc.
  • This work extends to nonlinear response in nonreciprocal systems

Photogalvanic Effects

  • Traditional photogalvanic effects appear in non-centrosymmetric materials
  • The discovered effect is a hydrodynamic analog
  • Unlike device geometry-induced rectification, this is an intrinsic effect

Conclusions and Discussion

Main Conclusions

  1. Nonreciprocal electron liquids under alternating driving generate direct current, with magnitude proportional to the square of the driving field amplitude.
  2. Direct current exhibits unique frequency dependence, showing superextensive behavior (d5\propto d^5) at low frequencies.
  3. Nonreciprocity significantly affects memory effects, leading to tilting of current-voltage hysteresis relationships.
  4. Effects of weak disorder can be described through phenomenological friction coefficients, without altering qualitative behavior.

Limitations

  1. Applicability range of perturbative theory: requires N1\mathcal{N} \ll 1, though numerical results show validity even at N1\mathcal{N} \sim 1.
  2. Geometric constraints: primarily considers Hall bar geometry; effects in other geometries remain incompletely explored.
  3. Temperature effects: theoretical analysis mainly at zero temperature; finite temperature effects require further investigation.

Future Directions

  1. Experimental realization: verification of predicted effects in 2D materials such as graphene
  2. Geometric optimization: exploration of nonreciprocal effects in different geometric structures
  3. Quantum effects: consideration of quantum corrections to hydrodynamic descriptions

In-Depth Evaluation

Strengths

  1. Strong theoretical innovation: first incorporation of nonreciprocity into alternating response studies of hydrodynamic electron transport
  2. Rigorous mathematical treatment: appropriate use of perturbative expansion, structure factors, and other mathematical tools
  3. Clear physical picture: concepts such as superextensive behavior and memory effects have well-defined physical meanings
  4. Sufficient numerical verification: theoretical predictions are well-validated by numerical simulations

Weaknesses

  1. Lack of experimental verification: primarily theoretical work lacking direct experimental validation
  2. Rough parameter estimates: relatively coarse estimation of actual material parameters
  3. Limited applicability conditions: hydrodynamic approximation conditions may be difficult to satisfy in real materials

Impact

  1. Academic value: provides new theoretical framework for understanding transport properties in nonreciprocal systems
  2. Application prospects: potential applications in terahertz devices, rectifiers, etc.
  3. Inspirational significance: provides new research directions for related fields

Applicable Scenarios

  1. High-mobility 2D materials: such as graphene, transition metal dichalcogenides
  2. Strongly correlated electron systems: cuprate high-temperature superconductors, etc.
  3. Topological materials: topological insulator surface states with intrinsic nonreciprocity

References

The paper cites over 60 relevant references, covering important works in hydrodynamic electron transport, nonreciprocal effects, photogalvanic effects, and other fields, providing a solid theoretical foundation for the research.