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.
论文ID : 2507.06933标题 : Photogalvanic effect in hydrodynamic flows of nonreciprocal electron liquids作者 : E. Kirkinis, L. Bonds, A. Levchenko, A. V. Andreev分类 : cond-mat.mes-hall (凝聚态物理-介观和霍尔效应)发表时间 : 2025年9月1日论文链接 : https://arxiv.org/abs/2507.06933 本文研究了由交流电场驱动的非线性流体动力学电子输运现象。在破缺时间反演对称性的非中心对称导体中,这种液体的非线性流动表现出非互易性,产生与交流电场幅值平方成正比的直流电流I D C I_{DC} I D C 。这是线性光伏效应(PGE)的流体动力学类似物,后者出现在破缺时间反演对称性的非中心对称材料中。I D C I_{DC} I D C 的大小既依赖于电子流体的性质,也依赖于流动的几何形状,可以用两个无量纲参数来表征:非互易数N \mathcal{N} N 和频率相关的振动数R \mathcal{R} R 。由于流体动力学输运的非局域性,在低频交流驱动下,I D C I_{DC} I D C 表现为超广延性。电流的交流分量同样受到非互易性的强烈影响:迟滞的电流-电压依赖关系变得倾斜,这可以用记忆保持时间的非互易性来解释。
流体动力学电子输运的兴起 :在电子动量和能量弛豫的时空尺度远大于动量和能量守恒的电子-电子散射尺度的固体系统中,电子输运可以用流体动力学方法描述。非互易效应的重要性 :在破缺时间反演对称性的非中心对称导体中,电子液体的泊肃叶流表现出强烈的非互易非线性双端电导。实验验证的需求 :近年来在石墨烯器件中已经实验证实了各种频率下的非线性电流整流效应,跨越了扩散-流体动力学过渡区域。该研究探索了一种内在的流体动力学电流整流机制,其中所需的反演对称性和时间反演对称性的破缺来自电子液体本身,而非器件几何结构。这与早期基于空间调制2D系统的流体动力学方法形成区别。
建立了非互易电子液体的广义Navier-Stokes方程 ,描述了时间依赖流动的演化。发现了流体动力学光伏效应 :在单色交流驱动下产生与驱动场幅值平方成正比的直流电流。引入了两个关键无量纲参数 :非互易数N = N U 0 \mathcal{N} = NU_0 N = N U 0 振动数R = ω d 2 / ν \mathcal{R} = \omega d^2/\nu R = ω d 2 / ν 揭示了超广延性 :在低频驱动下,直流电流I D C ∝ d 5 I_{DC} \propto d^5 I D C ∝ d 5 ,这是流体动力学输运非局域性的结果。发现了记忆效应 :非互易性影响电流的振荡部分,导致电流-电压迟滞依赖关系的倾斜。研究在霍尔棒几何结构中,由振荡电场驱动的非互易电子液体的非线性流体动力学输运,特别关注产生的直流电流和记忆效应。
电子系统动量密度p = ρ v \mathbf{p} = \rho\mathbf{v} p = ρ v 的演化方程为牛顿第二定律形式:
∂ p ∂ t = − ∇ ⋅ Π ^ + e n E \frac{\partial \mathbf{p}}{\partial t} = -\nabla \cdot \hat{\Pi} + en\mathbf{E} ∂ t ∂ p = − ∇ ⋅ Π ^ + e n E
其中动量通量张量密度为:
Π ^ = P δ i j + ρ v i v j − σ i j \hat{\Pi} = P\delta_{ij} + \rho v_i v_j - \sigma_{ij} Π ^ = P δ ij + ρ v i v j − σ ij
非互易修正δ σ i j \delta\sigma_{ij} δ σ ij 仅在时间反演对称性和反演对称性都破缺时出现:
δ σ i j = η N i j k l m v k ∂ l v m \delta\sigma_{ij} = \eta N_{ijklm} v_k \partial_l v_m δ σ ij = η N ijk l m v k ∂ l v m
对于霍尔棒几何,简化为:
δ σ x y = η N u ∂ y u \delta\sigma_{xy} = \eta N u \partial_y u δ σ x y = η N u ∂ y u
引入特征尺度U 0 = e n E d 2 / η U_0 = enEd^2/\eta U 0 = e n E d 2 / η 和I 0 = 2 3 ( e n ) 2 E d 3 / η I_0 = \frac{2}{3}(en)^2Ed^3/\eta I 0 = 3 2 ( e n ) 2 E d 3 / η ,得到无量纲动量方程:
R ∂ u ∂ τ = cos τ + ∂ ∂ y [ ( 1 + N u ) ∂ u ∂ y ] \mathcal{R}\frac{\partial u}{\partial \tau} = \cos \tau + \frac{\partial}{\partial y}[(1 + \mathcal{N}u)\frac{\partial u}{\partial y}] R ∂ τ ∂ u = cos τ + ∂ y ∂ [( 1 + N u ) ∂ y ∂ u ]
微扰展开求解 :将速度场表示为:
u ( y , t ) = ℜ { u 0 A C ( y ) e − i ω t + N [ u 1 D C ( y ) + u 1 A C ( y ) e − 2 i ω t ] } + O ( N 2 ) 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) u ( y , t ) = ℜ { u 0 A C ( y ) e − iω t + N [ u 1 D C ( y ) + u 1 A C ( y ) e − 2 iω t ]} + O ( N 2 ) 结构因子的引入 :定义结构因子S ± ( x ) = sinh x ± sin x cosh x + cos x S_\pm(x) = \frac{\sinh x \pm \sin x}{\cosh x + \cos x} S ± ( x ) = c o s h x + c o s x s i n h x ± s i n x 来表征频率依赖性。多尺度分析 :通过引入穿透深度δ = 2 ν / ω \delta = \sqrt{2\nu/\omega} δ = 2 ν / ω 和复波数k = ( 1 + i ) / δ k = (1+i)/\delta k = ( 1 + i ) / δ 来处理振荡流动。霍尔棒宽度 :d ∼ 1 d \sim 1 d ∼ 1 μm运动粘度 :ν ∼ 0.1 \nu \sim 0.1 ν ∼ 0.1 m²/s驱动频率 :太赫兹量级(R ∼ 10 R \sim 10 R ∼ 10 )非互易参数 :矢量型N = α B N = \alpha B N = α B ,张量型N = ∣ r ∣ cos ( 3 θ ) N = |r|\cos(3\theta) N = ∣ r ∣ cos ( 3 θ ) 通过数值求解时间依赖的Stokes方程(12)来验证微扰结果,初始条件为u ( y , 0 ) = 1 10 cos ( π y / 2 ) u(y,0) = \frac{1}{10}\cos(\pi y/2) u ( y , 0 ) = 10 1 cos ( π y /2 ) 。
I D C = I 0 3 N 8 [ − 2 R 2 + 2 R 5 / 2 S + ( 2 R ) ] 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] I D C = I 0 8 3 N [ − R 2 2 + R 5/2 2 S + ( 2 R ) ]
低频极限 (R ≪ 1 \mathcal{R} \ll 1 R ≪ 1 ):I D C ≈ − I 0 N / 10 I_{DC} \approx -I_0\mathcal{N}/10 I D C ≈ − I 0 N /10 高频极限 (R ≫ 1 \mathcal{R} \gg 1 R ≫ 1 ):I D C ≈ − 3 I 0 N / ( 4 R 2 ) I_{DC} \approx -3I_0\mathcal{N}/(4\mathcal{R}^2) I D C ≈ − 3 I 0 N / ( 4 R 2 ) 在低频下,直流电流具有超广延性:
I D C ≈ − N ( e n ) 3 E 2 d 5 15 η 2 I_{DC} \approx -\mathcal{N}\frac{(en)^3E^2d^5}{15\eta^2} I D C ≈ − N 15 η 2 ( e n ) 3 E 2 d 5
图3显示了微扰理论结果与完整数值解的良好一致性,验证了理论分析的正确性。
定义记忆时间尺度τ M ∼ d 2 / ( 2 ν ) \tau_M \sim d^2/(2\nu) τ M ∼ d 2 / ( 2 ν ) 在R = 2 \mathcal{R} = 2 R = 2 时迟滞回路面积达到最大 非互易性导致迟滞回路的倾斜和面积增加 Gurzhi开创性工作建立了固体中的流体动力学理论 在石墨烯、强关联电子系统等中得到实验验证 本文扩展到非互易系统的非线性响应 传统光伏效应出现在非中心对称材料中 本文发现的是流体动力学类似物 与器件几何诱导的整流效应不同,这里是内在效应 非互易电子液体在交流驱动下产生直流电流 ,其大小与驱动场幅值平方成正比。直流电流表现出独特的频率依赖性 ,在低频下具有超广延性(∝ d 5 \propto d^5 ∝ d 5 )。非互易性显著影响记忆效应 ,导致电流-电压迟滞关系的倾斜。弱无序的影响可以通过现象学摩擦系数来描述 ,不改定性行为。微扰理论的适用范围 :要求N ≪ 1 \mathcal{N} \ll 1 N ≪ 1 ,虽然数值结果显示在N ∼ 1 \mathcal{N} \sim 1 N ∼ 1 时仍然有效。几何限制 :主要考虑了霍尔棒几何,其他几何结构的效应尚未完全探索。温度效应 :理论分析主要在零温极限,有限温度效应需要进一步研究。实验实现 :在石墨烯等二维材料中验证预测的效应几何优化 :探索不同几何结构对非互易效应的影响量子效应 :考虑量子修正对流体动力学描述的影响理论创新性强 :首次将非互易性引入流体动力学电子输运的交流响应研究数学处理严谨 :微扰展开、结构因子等数学工具使用得当物理图像清晰 :超广延性、记忆效应等概念具有明确的物理意义数值验证充分 :理论预测得到数值模拟的很好验证实验验证缺乏 :主要是理论工作,缺乏直接的实验验证参数估计粗略 :对实际材料参数的估计相对粗糙适用条件限制 :流体动力学近似的适用条件在实际材料中可能难以满足学术价值 :为理解非互易系统的输运性质提供了新的理论框架应用前景 :可能在太赫兹器件、整流器等方面有应用潜力启发意义 :为相关领域的研究提供了新的研究方向高迁移率二维材料 :如石墨烯、过渡金属二硫化物等强关联电子系统 :铜氧化物高温超导体等拓扑材料 :具有内在非互易性的拓扑绝缘体表面态论文引用了60多篇相关文献,涵盖了流体动力学电子输运、非互易效应、光伏效应等多个领域的重要工作,为研究提供了坚实的理论基础。