Edelstein 效应：物理机制、对称性与微观理论的推导

该学习笔记以 V. M. Edelstein, Spin Polarization of Conduction Electrons Induced by Electric Current in Two-Dimensional Asymmetric Electron Systems, Solid State Communications 73, 233–235 (1990) 为主要参考。

Edelstein 效应是什么

Edelstein 效应（Edelstein effect, EE），也称为电流诱导自旋极化（Current-Induced Spin Polarization, CISP），是指在具有适当自旋轨道耦合与晶体对称性的导电体系中，施加电场产生电流的同时，可以诱导一个非平衡自旋极化。其基本响应关系可以写成

$\boxed{\delta S_i=\chi_{ij}^{\mathrm{EE}}E_j}.$

其中：$\delta S_i$ 是电场诱导的自旋角动量密度，$E_j$ 是外加电场，$\chi_{ij}^{\mathrm{EE}}$ 是 Edelstein 响应张量。也可以定义电流诱导的响应：

$\delta S_i=\kappa_{ij}j_j.$

在各向同性的二维 Rashba 体系中，这种关系具有特别简单的形式：

$\boxed{\delta\mathbf S=\chi_R\,\hat{\mathbf z}\times\mathbf E.}$

即电场沿 $x$ 方向时，诱导的自旋极化沿 $y$ 方向。具体正负号由 Rashba 耦合常数的符号及其定义决定。

原始论文的出发点为：二维电子气上下两个方向不等价时，哈密顿量允许出现线性动量的自旋轨道耦合。这个耦合使电子的自旋方向与动量相关。平衡时不同动量的自旋彼此抵消；外加电场改变动量空间的占据分布后，这种抵消不再完全，于是产生净自旋极化。

从对称性理解 Edelstein 效应

为什么普通二维电子气不能产生这种效应？

考虑不包含自旋轨道耦合的二维自由电子模型：

$H_0(\mathbf k)=\frac{\hbar^2k^2}{2m}-\mu.$

其中 $\mathbf k=(k_x,k_y)$，$m$ 为电子有效质量，$\mu$ 为化学势。这个哈密顿量与自旋无关，每个动量态都有两个自旋简并的本征态。如果施加电场，电子会形成非平衡分布，但是由于向上自旋和向下自旋的电子受到完全相同的作用，$\delta n_\uparrow=\delta n_\downarrow$，因此净自旋极化为 $\delta\mathbf S=0$。

更一般地，我们可以直接利用空间反演对称性来判断。空间反演操作 $\mathcal P$ 下：

$\mathcal P:\quad \mathbf E\rightarrow-\mathbf E,\qquad \mathbf S\rightarrow\mathbf S.$

这里电场是极矢量，而自旋是轴矢量。对于线性响应 $S_i=\chi_{ij}^{\mathrm{EE}}E_j$，如果体系具有反演对称性，反演后的物理定律必须保持不变：

$S_i=-\chi_{ij}^{\mathrm{EE}}E_j.$

与原来的关系比较，得到

$\boxed{\chi_{ij}^{\mathrm{EE}}=0.}$

因此，具有空间反演对称性的均匀体系，其均匀线性 Edelstein 响应受到对称性禁止。需要注意，这只是必要的对称性条件。破坏反演对称性并不保证一定存在非零响应；还需要考虑其他晶体对称性、自旋轨道耦合及费米面结构。

为什么时间反演对称性不禁止 Edelstein 效应？

时间反演操作下，

$\mathcal T:\quad \mathbf S\rightarrow-\mathbf S,\qquad \mathbf E\rightarrow\mathbf E,\qquad \mathbf j\rightarrow-\mathbf j.$

如果只考虑形式 $\mathbf S=\chi^{\mathrm{EE}}\mathbf E$，似乎会得到自旋为 TR-odd、电场为 TR-even，因此响应应该在时间反演对称体系中消失。

但是，这一推理不适用于存在输运弛豫的直流非平衡稳态。原因在于，Edelstein 效应的自旋极化不是平衡态自旋磁化，而是由电场持续驱动、散射维持的非平衡自旋积累。从物理上看：时间反演对称性要求平衡时的净自旋极化为零，电场驱动的非平衡分布不再具有平衡态的时间反演对称性，散射弛豫使体系建立具有确定电流方向的稳态。因此，即使哈密顿量具有时间反演对称性，也可以产生 Edelstein 效应。

事实上，文章的原始模型并不包含磁序或外部磁场，本身就是时间反演对称的，这一点对于区分平衡磁电效应与 Edelstein 输运效应非常重要。

镜面对称性决定自旋极化方向

二维 Rashba 体系可以具有平面内连续旋转对称性和垂直镜面对称性。其结构反演不对称方向由极矢量 $\hat{\mathbf z}$ 表示。如果只考虑平面内电场，并要求响应满足该结构的对称性，允许的张量形式为

$\chi_{ij}^{\mathrm{EE}}=\chi_R\epsilon_{izj},\qquad i,j=x,y.$

写成矩阵：

$\begin{pmatrix}S_x\\S_y\end{pmatrix} =\chi_R\begin{pmatrix}0&-1\\1&0\end{pmatrix} \begin{pmatrix}E_x\\E_y\end{pmatrix}.$

因此 $S_x=-\chi_RE_y,\quad S_y=\chi_RE_x$。这意味着 Edelstein 自旋极化方向与电场垂直。但这并不是所有 Edelstein 体系的普遍结论。对于低对称晶体，响应张量还可能具有对角元，从而产生与电场平行的自旋极化。

Rashba 自旋轨道耦合模型

考虑电子被限制在 $xy$ 平面内，界面法向为 $\hat{\mathbf z}$。如果上下方向不等价，那么垂直于二维电子气的有效电势梯度不再为零：$\partial_zV\neq0$。从非相对论极限的自旋轨道耦合可以得到如下形式的作用项：

$H_{\mathrm{SOC}}\propto\boldsymbol\sigma\cdot(\nabla V\times\mathbf p).$

当 $\nabla V=(\partial_zV)\hat{\mathbf z}$，并且 $\mathbf p=\hbar(k_x,k_y,0)$，有

$\nabla V\times\mathbf p=\hbar(\partial_zV)(-k_y,k_x,0).$

因此允许的线性自旋轨道耦合为 $H_R=\alpha(k_y\sigma_x-k_x\sigma_y)$，其中 $\alpha$ 吸收了相应的系数与符号，其量纲是能量乘长度。选择这一符号约定后，总哈密顿量写为

${H(\mathbf k)=\left(\frac{\hbar^2k^2}{2m}-\mu\right)\sigma_0+\alpha(k_y\sigma_x-k_x\sigma_y).}\label{eq:rashba_H}$

令 $\varepsilon_0(k)=\hbar^2k^2/(2m)-\mu$，定义 $\mathbf d(\mathbf k)=\alpha(k_y,-k_x,0)$。于是

$H(\mathbf k)=\varepsilon_0(k)\sigma_0+\mathbf d(\mathbf k)\cdot\boldsymbol\sigma.\label{eq:dvector_H}$

这在形式上与电子自旋处于动量依赖的有效磁场中类似。但应注意，$\mathbf d(\mathbf k)$ 并不是真实的外部磁场，其在时间反演下会随动量一起变换。利用 Pauli 矩阵性质：

$(\mathbf d\cdot\boldsymbol\sigma)^2=|\mathbf d|^2\sigma_0,$

可知 $|\mathbf d(\mathbf k)|=|\alpha|k$。选择分支指标 $\lambda=\pm1$，将能谱写为

$\boxed{\varepsilon_\lambda(k)=\frac{\hbar^2k^2}{2m}-\mu+\lambda\alpha k.}\label{eq:dispersion}$

为了明确各分支的标记，后面取 $\alpha>0$。当 $\alpha<0$ 时，结果可以通过保留 $\alpha$ 的符号推广。由于自旋轨道耦合，原来两重简并的抛物线能带发生劈裂。不过在 $k=0$ 处，$\varepsilon_+(0)=\varepsilon_-(0)=-\mu$，仍然存在 Kramers 简并。令 $k_x=k\cos\varphi,\quad k_y=k\sin\varphi$，则

$\mathbf d(\mathbf k)=\alpha k(\sin\varphi,-\cos\varphi,0).$

对两条 Rashba 分支，自旋的 Pauli 期望值满足

$\langle\boldsymbol\sigma\rangle_\lambda=\lambda\frac{\mathbf d}{|\mathbf d|}.$

因此

${\langle\boldsymbol\sigma\rangle_\lambda=\lambda(\sin\varphi,-\cos\varphi,0).}\label{eq:spin_texture}$

各个分量为

$\langle\sigma_x\rangle_\lambda=\lambda\frac{k_y}{k},\quad\langle\sigma_y\rangle_\lambda=-\lambda\frac{k_x}{k},\quad\langle\sigma_z\rangle_\lambda=0.$

可见：电子自旋始终位于二维平面内，自旋与电子动量垂直，两个 Rashba 分支在相同动量处具有相反的自旋方向。这就是 Rashba 自旋—动量锁定（spin-momentum locking）。

为什么平衡态没有净自旋？

平衡时电子遵循 Fermi–Dirac 分布

$f_0(\varepsilon)=\frac{1}{e^{\beta\varepsilon}+1}.$

定义单位面积的自旋角动量密度

$S_i=\frac{\hbar}{2}\sum_\lambda\int\frac{d^2k}{(2\pi)^2}\langle\sigma_i\rangle_\lambda f_0(\varepsilon_\lambda).\label{eq:spin_density}$

由于能谱只依赖 $k$，对 $S_x$ 的角度积分含有 $\int_0^{2\pi}\sin\varphi\,d\varphi=0$，而对 $S_y$ 的角度积分含有 $\int_0^{2\pi}\cos\varphi\,d\varphi=0$。同时 $S_z=0$，因此

$\boxed{\mathbf S^{(0)}=0.}$

从时间反演角度，也可以看到：

$\langle\boldsymbol\sigma\rangle_\lambda(-\mathbf k)=-\langle\boldsymbol\sigma\rangle_\lambda(\mathbf k),\qquad f_0[\varepsilon_\lambda(-\mathbf k)]=f_0[\varepsilon_\lambda(\mathbf k)].$

所以 $\mathbf k$ 和 $-\mathbf k$ 的贡献完全抵消。Edelstein 效应的关键就在于：施加电场后，第二个关系不再适用于非平衡分布。

利用 Boltzmann 方程推导 Edelstein 效应

接下来将从电子分布函数的运动方程出发，不预先假设费米面如何移动，而是推导出电场诱导的非平衡分布，再计算净自旋。在均匀、静态电场 $\mathbf E$ 中，电子的半经典运动满足

$\hbar\dot{\mathbf k}=-e\mathbf E,\label{eq:semiclassical_motion}$

其中 $e>0$ 表示元电荷，电子的实际电荷为 $-e$。电子在第 $\lambda$ 个能带上的群速度为

$\mathbf v_\lambda(\mathbf k)=\frac{1}{\hbar}\nabla_{\mathbf k}\varepsilon_\lambda(\mathbf k).$

将式 $(\ref{eq:dispersion})$ 代入：

$\nabla_{\mathbf k}\varepsilon_\lambda=\left(\frac{\hbar^2k}{m}+\lambda\alpha\right)\hat{\mathbf k}.$

于是

${\mathbf v_\lambda=\left(\frac{\hbar k}{m}+\frac{\lambda\alpha}{\hbar}\right)\hat{\mathbf k}.}\label{eq:band_velocity}$

注意两个能带的速度除了自由电子贡献，还包含 Rashba 自旋轨道耦合产生的贡献。例如 $x$ 方向速度为

$v_{\lambda x}=\left(\frac{\hbar k}{m}+\frac{\lambda\alpha}{\hbar}\right)\cos\varphi.$

定义第 $\lambda$ 个能带上的分布函数为 $f_\lambda(\mathbf r,\mathbf k,t)$，半经典 Boltzmann 方程为

$\frac{\partial f_\lambda}{\partial t}+\dot{\mathbf r}\cdot\nabla_{\mathbf r}f_\lambda+\dot{\mathbf k}\cdot\nabla_{\mathbf k}f_\lambda=I_{\mathrm{coll}}[f_\lambda].\label{eq:boltzmann}$

碰撞积分 $I_{\mathrm{coll}}$ 描述杂质或其他散射机制对分布函数的影响。我们先作如下近似：（1）电场空间均匀且不随时间变化，（2）系统已经达到非平衡稳态，（3）电子分布相对于平衡态的偏离较小，（4）可以用一个有效输运弛豫时间 $\tau$ 描述分布函数的恢复。其中，最后一个条件是弛豫时间近似，而不是一般含自旋轨道耦合体系的严格散射解。

在这些条件下，$\partial_t f_\lambda=0$，$\nabla_{\mathbf r}f_\lambda=0$，碰撞项近似为

$I_{\mathrm{coll}}[f_\lambda]= -\frac{f_\lambda-f_0}{\tau}.$

Boltzmann 方程因此简化为

$\dot{\mathbf k}\cdot\nabla_{\mathbf k}f_\lambda=-\frac{f_\lambda-f_0}{\tau}.$

利用式 $(\ref{eq:semiclassical_motion})$，得到

$-\frac{e}{\hbar}\mathbf E\cdot\nabla_{\mathbf k}f_\lambda=-\frac{f_\lambda-f_0}{\tau}.$

在电场下展开分布函数

由于我们研究线性响应，令 $f_\lambda=f_0+\delta f_\lambda$，其中 $\delta f_\lambda=O(E)$。将其代入：

$-\frac{e}{\hbar}\mathbf E\cdot\nabla_{\mathbf k}(f_0+\delta f_\lambda)=-\frac{\delta f_\lambda}{\tau}.$

其中 $\mathbf E\cdot\nabla_{\mathbf k}\delta f_\lambda$ 属于电场的二阶项，在一次线性响应中忽略。因此

$-\frac{e}{\hbar}\mathbf E\cdot\nabla_{\mathbf k}f_0=-\frac{\delta f_\lambda}{\tau}.$

整理得

$\delta f_\lambda=\frac{e\tau}{\hbar}\mathbf E\cdot\nabla_{\mathbf k}f_0.$

由于平衡分布是能量的函数，

$\nabla_{\mathbf k}f_0=\frac{\partial f_0}{\partial\varepsilon_\lambda}\nabla_{\mathbf k}\varepsilon_\lambda.$

而 $\nabla_{\mathbf k}\varepsilon_\lambda=\hbar\mathbf v_\lambda$，所以

${\delta f_\lambda(\mathbf k)=e\tau(\mathbf E\cdot\mathbf v_\lambda)\frac{\partial f_0}{\partial\varepsilon_\lambda}.}\label{eq:delta_f}$

这是 Edelstein 效应半经典推导中最核心的公式之一。在 $T=0$ 时，$f_0(\varepsilon)=\Theta(-\varepsilon)$，其中我们已经将化学势包含在 $\varepsilon_\lambda$ 中。因此 $\partial f_0/\partial\varepsilon=-\delta(\varepsilon)$。将其代入式 $(\ref{eq:delta_f})$，得到

${\color{blue}\delta f_\lambda=-e\tau(\mathbf E\cdot\mathbf v_\lambda)\delta(\varepsilon_\lambda).}\label{eq:delta_f_T0}$

这说明电场诱导的非平衡占据变化主要集中在费米面上。对于 $\mathbf E=E_x\hat{\mathbf x}$，

$\delta f_\lambda=-e\tau E_x\left(\frac{\hbar k}{m}+\frac{\lambda\alpha}{\hbar}\right)\cos\varphi\,\delta(\varepsilon_\lambda).$

当 $\cos\varphi>0$ 时，分布变化与 $\cos\varphi<0$ 时符号相反。因此动量空间中原来对称的占据发生了不对称变化，正负动量的自旋贡献不再完全抵消。

电场诱导的自旋密度

自旋密度定义为

$S_i=\frac{\hbar}{2}\sum_\lambda\int\frac{d^2k}{(2\pi)^2}\langle\sigma_i\rangle_\lambda f_\lambda.$

将 $f_\lambda=f_0+\delta f_\lambda$ 代入，并减去平衡贡献，得到

$\delta S_i=\frac{\hbar}{2}\sum_\lambda\int\frac{d^2k}{(2\pi)^2}\langle\sigma_i\rangle_\lambda\delta f_\lambda.\label{eq:delta_spin}$

利用式 $(\ref{eq:delta_f})$：

${\delta S_i=\frac{e\tau\hbar}{2}E_j\sum_\lambda\int\frac{d^2k}{(2\pi)^2}\langle\sigma_i\rangle_\lambda v_{\lambda j}\frac{\partial f_0}{\partial\varepsilon_\lambda}.}$

由此读出 Edelstein 响应系数：

${\chi_{ij}^{\mathrm{EE}}=\frac{e\tau\hbar}{2}\sum_\lambda\int\frac{d^2k}{(2\pi)^2}\langle\sigma_i\rangle_\lambda v_{\lambda j}\frac{\partial f_0}{\partial\varepsilon_\lambda}.}\label{eq:ee_chi_boltzmann}$

这是在能带对角、常数弛豫时间近似下的一般表达式，可以推广到具有更复杂能带结构和自旋纹理的体系。它揭示了 Edelstein 效应的三个基本因素：自旋纹理 + 费米面输运速度 + 非平衡占据变化。由前面得到的 Rashba 自旋纹理：

$\langle\sigma_x\rangle_\lambda=\lambda\sin\varphi,\qquad\langle\sigma_y\rangle_\lambda=-\lambda\cos\varphi.$

对于 $E_x$，$v_{\lambda x}\propto\cos\varphi$，因此 $\delta S_x$ 的角度积分含有 $\int_0^{2\pi}\sin\varphi\cos\varphi\,d\varphi=0$，所以 $\delta S_x=0$。而 $\delta S_y$ 的角度积分含有 $\int_0^{2\pi}\cos^2\varphi\,d\varphi=\pi$，因此 $\delta S_y\neq0$ 是允许的。同时 $\delta S_z=0$。从而，在当前模型中：

$\boxed{\mathbf E\parallel x\quad\Longrightarrow\quad\delta\mathbf S\parallel y.}$

这与第二节的对称性分析一致。

完整计算 Edelstein 响应系数

接下来不再停留于一般积分表达式，而是显式完成角向积分和径向积分，求出 $\delta S_y$。将式 $(\ref{eq:spin_texture})$、$(\ref{eq:band_velocity})$、$(\ref{eq:delta_f_T0})$ 代入式 $(\ref{eq:delta_spin})$：

$\begin{aligned} \delta S_y={}&\frac{\hbar}{2}\sum_\lambda\int\frac{d^2k}{(2\pi)^2}(-\lambda\cos\varphi)\left[-e\tau E_x\left(\frac{\hbar k}{m}+\frac{\lambda\alpha}{\hbar}\right)\cos\varphi\,\delta(\varepsilon_\lambda)\right]. \end{aligned}$

整理：

$\delta S_y=\frac{e\tau E_x\hbar}{2}\sum_\lambda\lambda\int\frac{d^2k}{(2\pi)^2}v_\lambda(k)\cos^2\varphi\,\delta(\varepsilon_\lambda),$

其中 $v_\lambda(k)=\hbar k/m+\lambda\alpha/\hbar$。在二维极坐标中 $d^2k=k\,dk\,d\varphi$，所以

$\begin{aligned} \delta S_y={}&\frac{e\tau E_x\hbar}{2}\sum_\lambda\lambda\int_0^\infty k\,dk\int_0^{2\pi}\frac{d\varphi}{(2\pi)^2}v_\lambda(k)\cos^2\varphi\,\delta(\varepsilon_\lambda). \end{aligned}$

因为 $\int_0^{2\pi}\cos^2\varphi\,d\varphi=\pi$，有

${\delta S_y=\frac{e\tau E_x\hbar}{8\pi}\sum_\lambda\lambda\int_0^\infty k\,dk\,v_\lambda(k)\delta(\varepsilon_\lambda).}\label{eq:Sy_radial}$

现在只剩下两个 Rashba 分支的径向积分。费米面条件为 $\varepsilon_\lambda(k)=0$，因此

$\frac{\hbar^2k^2}{2m}+\lambda\alpha k-\mu=0.$

两边乘以 $2m/\hbar^2$：

$k^2+2\lambda\frac{m\alpha}{\hbar^2}k-\frac{2m\mu}{\hbar^2}=0.$

定义 Rashba 动量尺度

$k_R=\frac{m\alpha}{\hbar^2}.\label{eq:kR}$

二次方程变成

$k^2+2\lambda k_Rk-\frac{2m\mu}{\hbar^2}=0.$

利用求根公式：

$k=-\lambda k_R\pm\sqrt{k_R^2+\frac{2m\mu}{\hbar^2}}.$

在 $\mu>0$ 的情况下，两条分支各具有一个正的费米波矢。定义

$k_0=\sqrt{k_R^2+\frac{2m\mu}{\hbar^2}},$

则

${k_{F,\lambda}=k_0-\lambda k_R.}\label{eq:kF}$

因此 $k_{F,+}=k_0-k_R$，$k_{F,-}=k_0+k_R$。两个费米圆之间的波矢差为

${k_{F,-}-k_{F,+}=2k_R.}\label{eq:kF_difference}$

注意，这个结果不显含化学势 $\mu$。由式 $(\ref{eq:band_velocity})$：

$v_\lambda(k)=\frac{\hbar k}{m}+\frac{\lambda\alpha}{\hbar}.$

在费米面 $k=k_{F,\lambda}$ 处：

$v_{F,\lambda}=\frac{\hbar(k_0-\lambda k_R)}{m}+\frac{\lambda\alpha}{\hbar}.$

利用 $\hbar k_R/m=\alpha/\hbar$，得到

$v_{F,\lambda}=\frac{\hbar k_0}{m}-\frac{\lambda\alpha}{\hbar}+\frac{\lambda\alpha}{\hbar}.$

所以

${v_{F,+}=v_{F,-}=\frac{\hbar k_0}{m}\equiv v_F.}\label{eq:vF}$

这个结果很有意义：虽然 Rashba 自旋轨道耦合使两个费米圆的半径不同，但在当前抛物线模型且两条分支均被占据时，它们的费米速度大小相同。

对于任意满足 $\varepsilon_\lambda(k_F)=0$ 的简单根，有

$\delta[\varepsilon_\lambda(k)]=\frac{\delta(k-k_{F,\lambda})}{|\partial_k\varepsilon_\lambda(k_{F,\lambda})|}.$

因为 $\partial_k\varepsilon_\lambda=\hbar v_\lambda(k)$，所以

$\delta[\varepsilon_\lambda(k)]=\frac{\delta(k-k_{F,\lambda})}{\hbar v_F}.$

于是

$\begin{aligned} I_\lambda&=\int_0^\infty k\,dk\,v_\lambda(k)\delta[\varepsilon_\lambda(k)]=\int_0^\infty k\,dk\,v_\lambda(k)\frac{\delta(k-k_{F,\lambda})}{\hbar v_F}. \end{aligned}$

利用 Dirac delta 函数积分性质：

$I_\lambda=\frac{k_{F,\lambda}v_{F,\lambda}}{\hbar v_F}.$

因此

${I_\lambda=\frac{k_{F,\lambda}}{\hbar}.}$

代回式 $(\ref{eq:Sy_radial})$：

$\delta S_y=\frac{e\tau E_x}{8\pi}\sum_\lambda\lambda k_{F,\lambda}.$

显式求和：

$\sum_\lambda\lambda k_{F,\lambda}=k_{F,+}-k_{F,-}.$

利用式 $(\ref{eq:kF_difference})$：

$k_{F,+}-k_{F,-}=-2k_R.$

所以

$\delta S_y=-\frac{e\tau E_x}{4\pi}k_R.$

再利用式 $(\ref{eq:kR})$：

${\delta S_y=-\frac{e\tau m\alpha}{4\pi\hbar^2}E_x.}\label{eq:Sy_final}$

这就是在我们选定的哈密顿量符号、自旋角动量归一化以及常数弛豫时间近似下，二维 Rashba 电子气的 Edelstein 响应。由旋转对称性得到一般矢量形式：

${\delta\mathbf S=-\frac{e\tau m\alpha}{4\pi\hbar^2}\hat{\mathbf z}\times\mathbf E.}\label{eq:EE_final}$

因此

$\chi_{yx}^{\mathrm{EE}}=-\frac{e\tau m\alpha}{4\pi\hbar^2},\qquad\chi_{xy}^{\mathrm{EE}}=+\frac{e\tau m\alpha}{4\pi\hbar^2}.$

并且 $\chi_{xx}^{\mathrm{EE}}=\chi_{yy}^{\mathrm{EE}}=0$。上述系数以物理自旋角动量密度 $\mathbf S=\hbar\langle\boldsymbol\sigma\rangle/2$ 定义。如果改用 Pauli 自旋极化 $\langle\boldsymbol\sigma\rangle$，或定义为磁矩密度，必须相应修改归一化因子。

为什么结果与化学势无关？

式 $(\ref{eq:Sy_final})$ 中没有显式出现 $\mu$，这是因为 $k_{F,-}-k_{F,+}=2k_R$ 与 $\mu$ 无关。在计算中，两条费米面的自旋贡献具有相反符号，但因费米圆周长不同，无法完全抵消。更具体地说：

$\delta S_y^{(+)}=\frac{e\tau E_x}{8\pi}k_{F,+},\qquad\delta S_y^{(-)}=-\frac{e\tau E_x}{8\pi}k_{F,-}.$

二者相加：

$\delta S_y=\frac{e\tau E_x}{8\pi}(k_{F,+}-k_{F,-}).$

虽然每一条分支的贡献都依赖 $\mu$，它们之和仅依赖两条费米圆之间的半径差。不过，该化学势无关性不是 Edelstein 效应的普遍性质。它依赖于抛物线色散、线性 Rashba 耦合、常数弛豫时间以及两条 Rashba 分支均被占据等条件。如果化学势降低至只有下 Rashba 分支被占据，或者出现非抛物线色散，该结果一般会改变。

从电场响应转换为电流响应

Edelstein 原始论文的标题强调的是电流诱导自旋极化，而不仅是电场诱导自旋极化。因此，我们还需要建立 $\delta\mathbf S$ 与电流密度 $\mathbf j$ 的关系。电子的电流密度为

$\mathbf j=-e\sum_\lambda\int\frac{d^2k}{(2\pi)^2}\mathbf v_\lambda(\mathbf k)f_\lambda(\mathbf k).\label{eq:current}$

由于平衡态无净电流，$\mathbf j^{(0)}=0$，因此线性电流为

$j_i=-e\sum_\lambda\int\frac{d^2k}{(2\pi)^2}v_{\lambda i}\delta f_\lambda.$

代入式 $(\ref{eq:delta_f})$：

${j_i=e^2\tau E_j\sum_\lambda\int\frac{d^2k}{(2\pi)^2}v_{\lambda i}v_{\lambda j}\left(-\frac{\partial f_0}{\partial\varepsilon_\lambda}\right).}$

于是电导张量为

$\sigma_{ij}=e^2\tau\sum_\lambda\int\frac{d^2k}{(2\pi)^2}v_{\lambda i}v_{\lambda j}\left(-\frac{\partial f_0}{\partial\varepsilon_\lambda}\right).\label{eq:sigma_boltzmann}$

这是常数弛豫时间近似下的能带对角 Drude 电导公式。对于各向同性 Rashba 模型，

$\sigma_{xx}=e^2\tau\sum_\lambda\int\frac{d^2k}{(2\pi)^2}v_\lambda^2(k)\cos^2\varphi\,\delta(\varepsilon_\lambda).$

完成角度积分：

$\sigma_{xx}=\frac{e^2\tau}{4\pi}\sum_\lambda\int_0^\infty k\,dk\,v_\lambda^2(k)\delta(\varepsilon_\lambda).$

径向积分为

$\int_0^\infty k\,dk\,v_\lambda^2(k)\delta(\varepsilon_\lambda)=\frac{k_{F,\lambda}v_F}{\hbar}.$

因此

$\sigma_{xx}=\frac{e^2\tau v_F}{4\pi\hbar}(k_{F,+}+k_{F,-}).$

利用 $k_{F,+}+k_{F,-}=2k_0$，$v_F=\hbar k_0/m$，得到

${\sigma_{xx}=\frac{e^2\tau k_0^2}{2\pi m}.}\label{eq:sigma_final}$

于是

$j_x=\frac{e^2\tau k_0^2}{2\pi m}E_x.$

这里需要注意：式 $(\ref{eq:sigma_final})$ 是我们所采用的简单能带 Boltzmann 近似的结果。若定义真实二维粒子数密度

$n=\frac{k_{F,+}^2+k_{F,-}^2}{4\pi}=\frac{k_0^2+k_R^2}{2\pi},$

则式 $(\ref{eq:sigma_final})$ 一般不严格等于 $ne^2\tau/m$。包含杂质诱导的自旋分布修正和相应顶角修正时，电导表达式需要一致地重新计算，不能直接将不同近似下的结果混用。将式 $(\ref{eq:sigma_final})$ 解出的 $E_x$ 代入式 $(\ref{eq:Sy_final})$：

$E_x=\frac{2\pi m}{e^2\tau k_0^2}j_x.$

于是

$\begin{aligned} \delta S_y&=-\frac{e\tau m\alpha}{4\pi\hbar^2}\frac{2\pi m}{e^2\tau k_0^2}j_x=-\frac{m^2\alpha}{2e\hbar^2k_0^2}j_x. \end{aligned}$

因此

${\delta S_y=-\frac{m^2\alpha}{2e\hbar^2k_0^2}j_x.}\label{eq:EE_current}$

这个结果清楚地展示出：在当前近似下，将 Edelstein 自旋极化表示为电流的函数以后，显式的弛豫时间 $\tau$ 被消去了。其物理原因是 $\delta S_y\propto\tau E_x$，$j_x\propto\tau E_x$。两者的比值不再显式依赖 $\tau$。这也正是原始论文讨论的一个重要现象：电流诱导的自旋极化与动量弛豫时间的关系，可以比电场诱导形式更加简单。原文式 (13)、(14) 采用了不同的自旋归一化、SOC 参数定义及杂质处理方法，因此在逐项比较数值系数前，需要先统一这些约定。

从 Kubo 线性响应理论理解 Edelstein 效应

前面使用 Boltzmann 理论完成了计算，但原始论文采用的是格林函数与顶角修正方法。为了理解两种方法之间的关系，我们接下来建立 Edelstein 响应的 Kubo 表达式。对标准 Rashba 哈密顿量

$H(\mathbf k)=\frac{\hbar^2k^2}{2m}\sigma_0+\alpha(k_y\sigma_x-k_x\sigma_y),$

自旋角动量算符是 $\hat S_i=\hbar\sigma_i/2$。速度算符需要由哈密顿量对动量求导：

$\hat v_i=\frac{1}{\hbar}\frac{\partial H}{\partial k_i}.\label{eq:velocity_operator}$

因此

$\boxed{\hat v_x=\frac{\hbar k_x}{m}\sigma_0-\frac{\alpha}{\hbar}\sigma_y,}\label{eq:vx_operator}$ $\boxed{\hat v_y=\frac{\hbar k_y}{m}\sigma_0+\frac{\alpha}{\hbar}\sigma_x.}\label{eq:vy_operator}$

这里特别需要注意：速度算符本身包含自旋矩阵。电流算符为 $\hat j_i=-e\hat v_i$。因此，外加电场不仅改变电子的动量分布，同时也通过自旋轨道耦合影响电流与自旋算符之间的响应关联。在空间均匀的电场下，可以使用随时间变化的矢势：

$\mathbf E(t)=-\frac{\partial\mathbf A(t)}{\partial t}.$

以最小耦合引入矢势：

$\hbar\mathbf k\rightarrow\hbar\mathbf k+e\mathbf A.$

将哈密顿量按 $\mathbf A$ 展开到一阶：

$H(\mathbf k+e\mathbf A/\hbar)\simeq H(\mathbf k)+eA_j\hat v_j.$

由于 $\hat j_j=-e\hat v_j$，得到微扰项

${H'(t)=-\hat j_jA_j(t).}\label{eq:perturbation}$

这样就把电场诱导的自旋极化转换成了自旋算符与电流算符之间的线性响应问题。

定义推迟自旋—电流关联函数：

$\boxed{\Pi_{ij}^{R}(t)=-\frac{i}{\hbar}\Theta(t)\left\langle[\hat S_i(t),\hat j_j(0)]\right\rangle.}\label{eq:Pi_retarded}$

这里的算符应理解为总自旋和总电流；对均匀体系还需除以面积，才能得到自旋密度响应。进行 Fourier 变换：

$\Pi_{ij}^{R}(\omega)=\int_{-\infty}^{\infty}dt\,e^{i\omega t}\Pi_{ij}^{R}(t).$

根据 Kubo 公式，电场诱导自旋密度为

$\delta S_i(\omega)=\Pi_{ij}^{R}(\omega)A_j(\omega).$

采用时间依赖 $e^{-i\omega t}$，有 $E_j(\omega)=i\omega A_j(\omega)$。

于是

$\delta S_i(\omega)=\frac{\Pi_{ij}^{R}(\omega)}{i\omega}E_j(\omega).$

因此形式上可以定义

$\chi_{ij}^{\mathrm{EE}}(\omega)=\frac{\Pi_{ij}^{R}(\omega)}{i\omega}.$

更严谨地讲，在规范不变的完整响应理论中，需要正确处理静态矢势及平衡响应的抵消，可以写成

${\chi_{ij}^{\mathrm{EE}}(\omega)=\frac{\Pi_{ij}^{R}(\omega)-\Pi_{ij}^{R}(0)}{i\omega}.}\label{eq:Kubo_EE}$

在使用这一表达式时，频率极限、无序平均以及顶角修正必须采用一致的处理方式。非相互作用电子的 Matsubara 格林函数为

$G(\mathbf k,i\omega_n)=\left[i\hbar\omega_n-H(\mathbf k)\right]^{-1}.$

引入 Rashba 分支的投影算符：

$P_\lambda(\mathbf k)=\frac12\left[\sigma_0+\lambda\hat{\mathbf d}(\mathbf k)\cdot\boldsymbol\sigma\right].\label{eq:projectors}$

其中 $\hat{\mathbf d}=\mathbf d/|\mathbf d|$。格林函数可以写成谱分解形式：

${G(\mathbf k,i\omega_n)=\sum_{\lambda=\pm}\frac{P_\lambda(\mathbf k)}{i\hbar\omega_n-\varepsilon_\lambda(k)}.}\label{eq:greens_projected}$

忽略顶角修正时，自旋—电流关联函数由单个泡图给出：

$\begin{aligned} \Pi_{ij}(i\Omega_m)={}&-\frac1\beta\sum_{\omega_n}\int\frac{d^2k}{(2\pi)^2}\operatorname{Tr}\left[\hat S_iG(\mathbf k,i\omega_n+i\Omega_m)\hat j_jG(\mathbf k,i\omega_n)\right], \end{aligned}\label{eq:kubo_bubble}$

这里默认采用单粒子迹及与式 $(\ref{eq:Pi_retarded})$ 一致的归一化约定。代入式 $(\ref{eq:greens_projected})$，得到

$\begin{aligned} \Pi_{ij}(i\Omega_m)={}&-\frac1\beta\sum_{\omega_n}\sum_{\lambda,\lambda'}\int\frac{d^2k}{(2\pi)^2}\frac{\operatorname{Tr}[\hat S_iP_\lambda\hat j_jP_{\lambda'}]}{(i\hbar\omega_n+i\hbar\Omega_m-\varepsilon_\lambda)(i\hbar\omega_n-\varepsilon_{\lambda'})}. \end{aligned}$

执行 Matsubara 频率求和之后：

${\Pi_{ij}(i\Omega_m)=\sum_{\lambda,\lambda'}\int\frac{d^2k}{(2\pi)^2}\mathcal M_{ij}^{\lambda\lambda'}(\mathbf k)\frac{f_0(\varepsilon_\lambda)-f_0(\varepsilon_{\lambda'})}{i\hbar\Omega_m+\varepsilon_\lambda-\varepsilon_{\lambda'}}.}\label{eq:Pi_band}$

其中矩阵元 $\mathcal M$ 根据上述算符顺序定义，适当交换分支指标后即可写为这一形式。由式 $(\ref{eq:Pi_band})$，可以区分两类过程：

$\lambda=\lambda’$：同一能带内的响应，与费米面附近的占据变化密切相关。
$\lambda\neq\lambda’$：不同 Rashba 分支之间的相干响应，包含带间矩阵元。

这种分解在复杂多能带体系中非常有用，因为 Edelstein 效应不一定完全由简单的费米面移动描述。

Edelstein 效应与时间反演对称性

平衡自旋极化为什么被时间反演禁止?

自旋是时间反演奇的物理量：

$\mathcal T\hat{\mathbf S}\mathcal T^{-1}=-\hat{\mathbf S}.$

如果哈密顿量具有时间反演对称性，

$\mathcal TH(\mathbf k)\mathcal T^{-1}=H(-\mathbf k),$

则热平衡态也满足相应的时间反演不变性。因此平衡自旋期望值必须满足

$\langle\mathbf S\rangle_{\rm eq}=-\langle\mathbf S\rangle_{\rm eq},$

得到

$\boxed{\langle\mathbf S\rangle_{\rm eq}=0.}$

为什么非平衡 Edelstein 自旋极化可以不为零？

关键是电流本身属于时间反演奇的物理量：

$\mathcal T\mathbf j=-\mathbf j.$

Edelstein 效应也可以定义为 $S_i=\kappa_{ij}j_j$。在时间反演下，自旋和电流同时改变符号：

$-S_i=\kappa_{ij}(-j_j).$

因此这种关系并不与时间反演对称的哈密顿量矛盾。不过，必须强调：这只是直观的对称性说明，不能取代对非平衡响应函数的严格分析。有限电流的稳态不是平衡态，它的维持依赖外部驱动和耗散。真正能够判断响应是否存在的，是完整动力学和相应的输运系数。

Edelstein 效应是否一定属于 TR-even 响应？

这里最重要的是先明确对响应系数的定义。如果所谓 TR-even 是指：将某一体系的全部磁性序参量 $M$ 反转以后，比较响应系数 $\chi^{\mathrm{EE}}(M)$ 与 $\chi^{\mathrm{EE}}(-M)$，则可以定义

$\chi_{\mathrm{even}}=\frac{\chi(M)+\chi(-M)}2,\qquad\chi_{\mathrm{odd}}=\frac{\chi(M)-\chi(-M)}2.\label{eq:TR_decomposition}$

它们分别满足

$\chi_{\mathrm{even}}(M)=\chi_{\mathrm{even}}(-M),\qquad\chi_{\mathrm{odd}}(M)=-\chi_{\mathrm{odd}}(-M).$

对于本文研究的普通 Rashba 电子气，不存在磁性序参量，因此传统 Edelstein 效应无需时间反演对称性破缺即可存在。但是，在磁性 Rashba 体系或更复杂的自旋轨道耦合磁体中，可以进一步研究反转磁序时自旋响应系数的不同变化。不能简单根据自旋是 TR-odd、电场是 TR-even，就断言所有电场诱导自旋响应系数必然是 TR-odd。

Edelstein 效应的几个重要物理特点

由式 $(\ref{eq:Sy_final})$ 可知：$\delta S_y\propto\alpha$。因此，在当前模型及近似下：

$\lim_{\alpha\to0}\delta S_y=0.$

如果没有 Rashba 自旋轨道耦合，动量与自旋不再锁定，电场驱动的动量分布不对称不能产生这一自旋极化机制。此外，当 $\alpha$ 反号时，$\delta S_y(\alpha)=-\delta S_y(-\alpha)$，即诱导自旋的方向反转。

在固定电场的条件下：$\delta S_y\propto\tau$。弛豫时间越长，电场引起的动量分布偏离越明显，非平衡自旋极化也越大，但在固定电流下，简单 Drude 近似给出的显式 $\tau$ 依赖会消去。这两个结论并不冲突，因为它们描述的是不同的实验控制条件。还应强调，固定电场下的 $\tau\to\infty$ 极限不能直接用上述稳态公式外推。无散射体系中的直流电场未必能建立同样的稳态分布。

与自旋 Hall 效应的区别

虽然 Edelstein 效应和自旋 Hall 效应都与自旋轨道耦合有关，但它们研究的响应量不同。

性质	Edelstein 效应	自旋 Hall 效应
外部驱动	电场或电流	电场或电流
核心响应	自旋密度 $S_i$	自旋流 $J_i^{s_a}$
典型结果	均匀非平衡自旋极化	横向自旋输运
常见后果	电流诱导自旋积累或自旋力矩	边界自旋积累、自旋注入
一般条件	受到反演及其他点群对称性约束	可以存在于某些中心对称体系