量子输运理论:散射矩阵理论

  1. 4月前

    Phantom_Ghost

    1楼 10月21日 数学版主, 物理版主, 优秀回答者
    4小时前Phantom_Ghost 重新编辑

    - 2014 -

    粒子数密度算符
    一次量子化
    \[
    \rho(x)=\psi^\dagger(x)\psi(x)=\int\psi^\dagger(x')\delta(x-x')\psi(x')dx'
    \]
    二次量子化形式为
    \begin{align*}
    \hat{\rho}(x)&=\int\hat{\psi}^\dagger(x')\delta(x-x')\hat\psi(x')dx'\\
    &=\sum_{\alpha\beta} \rho_{\alpha\beta} b_\alpha^\dagger b_\beta\\
    \rho_{\alpha\beta}(x)&=\int\phi^*_\alpha(x')\delta(x-x')\phi_\beta(x') dx'\\
    &=\phi^*_\alpha(x)\phi_\beta(x)
    \end{align*}
    总粒子数算符为
    \[
    \hat{N}=\int\hat\rho(x) dx=\sum_\alpha b_\alpha^\dagger b_\alpha
    \]
    流密度算符
    在经典力学中 $\mathbf{j}=\rho e\mathbf{v}$,在量子力学中为
    \[
    \hat{j}=\frac{e}{2}\int\hat{\psi}^\dagger(x')\left[v(x)\delta(x-x')+\delta(x-x')v(x')\right]\hat\psi(x') dx'
    \]
    正则量子化中速度算符为 $\hat{v}=-\frac{i\hbar}{m}\nabla$。于是流密度算符表达为
    \[
    \hat{j}=\frac{\hbar e}{2mi}\left[\hat{\psi}^\dagger(\nabla\hat\psi)-(\nabla\hat{\psi}^\dagger)\hat\psi\right]
    \]
    量子统计里面有密度矩阵算符$\hat{\rho}=\frac{1}{Z}e^{-\beta(\hat{H}-\mu\hat{N})}$,配分函数 $Z=\text{Tr}\;e^{-\beta(\hat{H}-\mu\hat{N})}$。于是某算符$\hat{A}$的统计平均值为$\langle\hat{A}\rangle=\text{Tr}[\hat\rho\hat{A}]$,于是有
    \begin{align*}
    \langle b_\alpha^\dagger b_\alpha\rangle&=\text{Tr}\left[\frac{1}{Z}e^{-\beta(\hat{H}-\mu\hat{N})}b_\alpha^\dagger b_\alpha\right]\\
    \hat{H}&=\sum_\alpha \varepsilon_\alpha b_\alpha^\dagger b_\alpha\;,\;\hat{N}=\sum_\alpha b_\alpha^\dagger b_\alpha
    \end{align*}
    可得到Fermi分布
    \begin{align*}
    f_\alpha&=\langle b_\alpha^\dagger b_\alpha\rangle\\
    &=\frac{1}{e^{\beta(\varepsilon_\alpha-\mu)}+1}
    \end{align*}
    在线性相应理论里面,有著名的Kubo公式。假设受微扰系统为 $H=H_0+H_I$,零温密度矩阵算符为$\hat\rho=\sum_i |i\rangle\langle i|$,密度矩阵的运动方程为(注意不同于一般算符的Heisenberg方程 $i\hbar\dot{A}=[A,H]$)
    \[
    i\hbar\frac{d\hat{\rho}}{dt}=[H,\hat\rho]
    \]
    假设微扰在$t=0$时刻引入的 $H_I(t)=H_I(t)\theta(t)$,有限温度当$t<0$,密度矩阵为$\hat{\rho}_0=\frac{1}{Z}e^{-\beta(H_0-\mu N)}$ 不含时。当$t\geq 0$,$\hat\rho(t)$可使用微扰法求解。设$\hat\rho(t)$展开为
    \[
    \hat\rho(t)=\hat{\rho}_0(t)+\hat{\rho}_1(t)+\hat{\rho}_2(t)+...
    \]
    运动方程也展开为
    \[
    i\frac{d}{dt}\left\{\hat{\rho}_0(t)+\hat{\rho}_1(t)+\hat{\rho}_2(t)+...\right\}=[H_0+H_I,\hat{\rho}_0(t)+\hat{\rho}_1(t)+\hat{\rho}_2(t)+...]
    \]
    比较两边,零阶项为 $i\frac{d\rho_0(t)}{dt}=[H_0,\rho_0(t)]=0$
    一阶项为 $i\frac{d\rho_1(t)}{dt}=[H_I,\rho_0(t)]+[H_0,\rho_1(t)]$
    一般解为:
    \begin{align*}
    \hat\rho_0(t)&=\hat\rho_0=\frac{1}{Z}e^{-\beta(H_0-\mu H)}\\
    \hat\rho_1(t)&=-i\int_0^t dt'\;e^{-iH_0(t-t')}[H_I(t'),\rho_0]e^{iH_0(t-t')}
    \end{align*}
    因此一级近似解为 $\hat\rho(t)=\hat\rho_0(t)+\hat\rho_1(t)$
    由于算符$A$的统计平均定义为 $\langle A\rangle=\text{Tr}[\hat\rho A]$ ,因此有
    \[
    \langle A\rangle_t=\text{Tr}[\hat{\rho}_0 A(t)]+\text{Tr}[\hat{\rho}_1(t) A(t)]\Rightarrow \langle A\rangle_t=\langle A\rangle_0+\text{Tr}[\hat{\rho}_1(t) A(t)]
    \]
    于是对于$\hat{\rho}_1(t)$有
    \begin{align*}
    \text{Tr}[\hat{\rho}_1(t) A(t)]&=-i\int_0^t dt'\;\text{Tr}\left\{e^{-iH_0(t-t')}[H_I(t'),\hat{\rho}_0]e^{iH_0(t-t')} A(t)\right\}\\
    &=-i\int_0^t dt'\;\text{Tr}\left\{\hat{\rho}_0[A^{H_0}(t),H_I^{H_0}(t')]\right\}\\
    &=-i\int_0^t dt'\langle[A^{H_0}(t),H_I^{H_0}(t')]\rangle
    \\\\
    A^{H_0}(t)&=e^{iH_0t}A(t)e^{-iH_0 t}\\
    H_I^{H_0}(t')&=e^{iH_0t'}H_I(t')e^{-iH_0 t'}
    \end{align*}
    最后便得到Kubo公式
    \[
    \langle A\rangle_t=\langle A\rangle_0-i\int_0^t dt'\langle[A^{H_0}(t),H_I^{H_0}(t')]\rangle
    \]
    流密度响应
    流密度对外磁场扰动的响应,哈密顿量为
    \[
    \hat{H}=\int \hat{\psi}^\dagger(r,t)\frac{1}{2m}\left[P(r)+\frac{q}{c}A(r,t)\right]^2\hat{\psi}(r,t) dr
    \]
    在哈密顿量里保持$A$的线性项并忽略$A^2$,有
    \begin{align*}
    \hat{H}&=\hat{H}_0+\hat{H}_I\\
    \hat{H}_0&=\frac{1}{2m}\int \hat{\psi}^\dagger(r,t) P^2(r) \hat{\psi}(r,t) dr\\
    \hat{H}_I&=\frac{q}{c}\int A(r,t)\cdot J(r,t) dr\\
    J(r,t)&=-\frac{i\hbar}{2m}\left[\hat{\psi}^\dagger(r,t)(\nabla\hat{\psi}(r,t) )-(\nabla \hat{\psi}^\dagger(r,t))\hat{\psi}(r,t)\right]
    \end{align*}
    根据Kubo定理,流$J(r)$平均为:
    \[
    \langle J(r,t)\rangle_t=\langle J(r)\rangle_0-i\int_0^t dt'\int\langle[J(r,t),J(r',t')]\rangle\cdot A(r',t') dr'
    \]
    $Q(r,t;r',t')=-i\langle[J(r,t),J(r',t')]\rangle$ 为流-流关联函数。

    粒子数密度响应
    粒子数密度对外电场扰动响应,其哈密顿量为
    \begin{align*}
    \hat{H}&=\int\hat{\psi}^\dagger(r,t)\left(\frac{1}{2m}P^2+V_\text{ext}(r,t)\right)\hat{\psi}(r,t) dr\\
    &=\hat{H}_0+\hat{H}_I\\
    \hat{H}_0&=\frac{1}{2m}\int\hat{\psi}^\dagger(r,t) P^2\hat{\psi}(r,t) dr\\
    \hat{H}_I&=\int\hat{\psi}^\dagger(r,t)V_\text{ext}(r,t)\hat{\psi}(r,t) dr=\int\hat{\rho}(r,t)V_\text{ext}(r,t) dr\\
    \hat{\rho}(r,t)&=\hat{\psi}^\dagger(r,t)\hat{\psi}(r,t)
    \end{align*}
    由Kubo定理,粒子数密度平均为
    \[
    \langle\hat{\rho}(r,t)\rangle_t=\langle\hat{\rho}(t)\rangle_0-i\int_0^t dt'\int\langle[\hat{\rho}(r,t),\hat{\rho}(r',t')]\rangle V_\text{ext}(r',t') dr'
    \]
    定义 $\Pi(r,t;r',t')=-i\langle[\hat{\rho}(r,t),\hat{\rho}(r',t')]\rangle$ 为密度-密度关联函数,或者密度响应函数(Lindahard函数)。在单粒子近似下,场算符可写为
    \[
    \hat\psi(r,t)=\sum_i\varphi_i(r)\hat{a}_i(t)\;,\;\hat{a}_i(t)=\hat{a}_i e^{-i\varepsilon_i t}\;(\hbar=1)
    \]
    代入密度响应函数中得到
    \begin{align*}
    \Pi(r,t;r',t')&=-i\langle[\hat{\rho}(r,t),\hat{\rho}(r',t')]\rangle\\
    &=-i\langle[\hat{\psi}^\dagger(r,t)\hat\psi(r,t),\hat{\psi}^\dagger(r',t')\hat\psi(r',t')]\rangle\\
    &=-i\sum_{ijkl}\varphi^*_i(r)\varphi_j(r)\varphi^*_k(r')\varphi_l(r')e^{i\varepsilon_i t-i\varepsilon_j t+i\varepsilon_k t'-i\varepsilon_l t'}\\
    &(\langle a_i^\dagger a_j a_k^\dagger a_l\rangle-\langle a_k^\dagger a_l a_i^\dagger a_j\rangle)
    \end{align*}
    由于 $a_j a_k^\dagger+a_k^\dagger a_j=\delta_{kj}$,及 $\langle a_k^\dagger a_j\rangle=f_k\delta_{kj}$,其中$f_k$为Fermi分布。便得到
    \[
    \Pi(r,t;r',t')=-i\sum_{kl}\varphi^*_l(r)\varphi_k(r)\varphi^*_k(r')\varphi_l(r')e^{i\varepsilon_l(t-t')-i\varepsilon_k(t-t')}(f_l-f_k)
    \]
    若满足时间平移对称则可作Fourier变换到频率空间
    \begin{align*}
    \Pi(r,r';\omega)&=\int \Pi(r,t;r',t')e^{i\omega(t-t')}d(t-t')\\
    &=\sum_{kl}\varphi^*_l(r)\varphi_k(r)\varphi^*_k(r')\varphi_l(r')\frac{f_k-f_l}{\omega+(\varepsilon_l-\varepsilon_k)+i0^+}\\
    -i\lim_{\eta\to 0^+}&\int_{\infty}^0 e^{i((\varepsilon_l-\varepsilon_k)x}e^{i\omega x}e^{-\eta x}dx=-\frac{1}{\omega+(\varepsilon_l-\varepsilon_k)+i0^+}\\
    &(0<t'<t\;,\;t\to\infty)
    \end{align*}

    Büttiker提出散射矩阵理论(PRB,Vol.46,12485(1992))及由Datta等发展(PRB,Vol.53,No.24(1996)),讨论低频、弱非线性介观体系中量子输运问题。纳米电子器件在此情况下要保持电流守恒以及规范不变性则须考虑电子间Coulomb作用,可在平均场意义下约化成有效势(需自洽求解),有效势和外电场(导致非平衡态)改变了波函数,波函数也反过来影响有效势。散射矩阵理论可以简明地给出自洽计算以讨论量子输运问题。

    研究一个多端口介观导体,端口以$\alpha$标记,每个端口有多种模式。现在先对模型做一些理想处理:假设各端口间存在相干透射,则介观导体须足够小;每个端口与处于热平衡的热库连接,暂时假定端口均处于相同温度$T$,且每个端口载流子化学势各为$\mu_\alpha$,平衡时化学势为$\mu$,弱外场(各端口电压$V_\alpha$)激励下有$\mu_\alpha=\mu+e V_\alpha$;在介观导体上加静磁场,各端口间透射振幅可用散射矩阵$S$描述。$S$是幺正矩阵,加磁场时左右入射体系不对称 $S_{\alpha\beta}(-B)=S_{\beta\alpha}(B)$,$S(E_F,U)$与Fermi能以及有效势都有关系从而依赖于端口电势。散射矩阵理论在平均场意义下考虑电子间相互作用,形式上仍是单粒子理论,引入自洽平均场在于既要考虑相互作用又要确保流守恒以及规范不变性。
    现在来具体计算端口$\alpha$的电流,可由电流守恒再推出流经器件的电流。设端口上电子的色散关系为 $E_{\alpha m}(k)=E_{\alpha m}(0)+\frac{\hbar^2 k_{\alpha m}^2}{2m}$,散射态为 $\chi_{\alpha m}^\pm(x,y)=e^{\pm ik_{\alpha m}x}\phi_{\alpha m}(y)$,这里由于是散射态故不需要归一化,$\phi_{\alpha m}(y)$为$y$轴方向驻波(量子限域),$e^{\pm ik_{\alpha m}x}$为$x$轴方向平面波(行波)。通常端口内电子波函数一般形式写为:
    \[
    \psi(r,t)=\sum_n [a_n(t)\chi_n^+(r) +b_n(t)\chi_n^-(r)]
    \]
    因此在端口内散射态的场算符为
    \begin{align*}
    \hat{\psi}(r,t)&=\sum_n\int\frac{dE_{\alpha m}}{2\pi\hbar}\sqrt{v_{\alpha m}}\left\{\chi^+_{\alpha m}(r)\hat{a}_{\alpha m}(E_{\alpha m})+\chi^-_{\alpha m}(r)\hat{b}_{\alpha m}(E_{\alpha m})\right\}e^{-i\frac{E_{\alpha m}}{\hbar}t}\\
    \hat{\psi}^\dagger(r,t)&=\sum_n\int\frac{dE_{\alpha m}}{2\pi\hbar}\sqrt{v_{\alpha m}}\left\{\chi^{+ *}_{\alpha m}(r)\hat{a}^\dagger_{\alpha m}(E_{\alpha m})+\chi^{- *}_{\alpha m}(r)\hat{b}^\dagger_{\alpha m}(E_{\alpha m})\right\}e^{i\frac{E_{\alpha m}}{\hbar}t}
    \end{align*}
    $v_{\alpha m}=\frac{\hbar k_{\alpha m}}{m}$为速度,作为入射流归一化常数。二次量子化下电流密度算符为
    \[
    \hat{J}(r,t)=-i\frac{e\hbar}{2m}\left\{\hat{\psi}^\dagger(\nabla\hat{\psi})-(\nabla\hat{\psi}^\dagger)\hat{\psi}\right\}
    \]
    则通过$\alpha$端口电流算符对应着电流密度在$x_\alpha$上沿截面的积分 $I_\alpha(t)=\int \hat{J}(r,t)\left|\right._{x_\alpha}\;dy_{\alpha}$。注意只有散射态对电流有贡献(尽管电流 密度分布可依赖于局域态),于是可把散射态场算符代入电流表达式:
    \[
    \frac{\partial\hat{\psi}^\dagger}{\partial x}\hat{\psi}=\sum_{mn}\sqrt{v_{\alpha m} v_{\alpha n}}e^{-i(E_{\alpha n}-E_{\alpha m})t/\hbar}\left(\frac{d\chi_{\alpha m}^{+ *}}{dx}\hat{a}_{\alpha m}^\dagger+\frac{d\chi_{\alpha m}^{- *}}{dx}\hat{b}_{\alpha m}^\dagger\right)\left(\chi_{\alpha n}^+\hat{a}_{\alpha n}+\chi_{\alpha n}^-\hat{b}_{\alpha n}\right)
    \]
    乘积展开得到四个项
    \begin{align*}
    \frac{d\chi_{\alpha m}^{+ *}}{dx}\chi_{\alpha n}^+\hat{a}_{\alpha m}^\dagger\hat{a}_{\alpha n}\;&,\;\frac{d\chi_{\alpha m}^{+ *}}{dx}\chi_{\alpha n}^-\hat{a}_{\alpha m}^\dagger\hat{b}_{\alpha n}\\
    \frac{d\chi_{\alpha m}^{- *}}{dx}\chi_{\alpha n}^+\hat{b}_{\alpha m}^\dagger\hat{a}_{\alpha n}\;&,\;\frac{d\chi_{\alpha m}^{- *}}{dx}\chi_{\alpha n}^-\hat{b}_{\alpha m}^\dagger\hat{b}_{\alpha n}
    \end{align*}
    同理
    \[
    \hat{\psi}^\dagger\frac{\partial\hat{\psi}}{\partial x}=\sum_{mn}\sqrt{v_{\alpha m} v_{\alpha n}}e^{-i(E_{\alpha n}-E_{\alpha m})t/\hbar}\left(\chi_{\alpha m}^{+*}\hat{a}_{\alpha m}^\dagger+\chi_{\alpha m}^{-*}\hat{b}_{\alpha m}^\dagger\right)\left(\frac{d\chi_{\alpha n}^{+}}{dx}\hat{a}_{\alpha n}+\frac{d\chi_{\alpha n}^{-}}{dx}\hat{b}_{\alpha n}\right)
    \]
    展开也得四项
    \begin{align*}
    \chi_{\alpha m}^{+*}\frac{d\chi_{\alpha n}^{+}}{dx}\hat{a}_{\alpha m}^\dagger\hat{a}_{\alpha n}\;&,\;\chi_{\alpha m}^{+*}\frac{d\chi_{\alpha n}^{-}}{dx}\hat{a}_{\alpha m}^\dagger\hat{b}_{\alpha n}\\
    \chi_{\alpha n}^{+*}\frac{d\chi_{\alpha n}^{-}}{dx}\hat{b}_{\alpha m}^\dagger\hat{a}_{\alpha n}\;&,\;\chi_{\alpha m}^{-*}\frac{d\chi_{\alpha n}^{-}}{dx}\hat{b}_{\alpha m}^\dagger\hat{b}_{\alpha n}
    \end{align*}
    由于 $\hat{J}\sim\hat{\psi}^\dagger(\nabla\hat{\psi})-(\nabla\hat{\psi}^\dagger)\hat{\psi}$,由乘积项差再作横向积分可得
    \[
    I_{mn}^{\sigma_1\sigma_2}(E,E')=\frac{e\hbar}{2mi}\int dy_\alpha\left[\chi_{\alpha m}^{\sigma_1 *}(E)\frac{d\chi_{\alpha n}^{\sigma_2}(E')}{dx}-\frac{d\chi_{\alpha m}^{\sigma_1 *}(E)}{dx}\chi_{\alpha n}^{\sigma_2 *}(E')\right]
    \]
    设$E=k^2,E'=q^2$(原子单位),求出微分后积分化为 $I_{mn}^{\sigma_1\sigma_2}(E,E')=\frac{e\hbar}{2m}\delta_{mn}(\sigma_2 q+\sigma_1 k) e^{-i(\sigma_1 k-\sigma_2 q)x}$,当$E\simeq E',k\simeq q$,得到 $I_{mn}^{\sigma_1\sigma_2}(E,E')=\frac{e\hbar}{m}\delta_{mn}\delta_{\sigma_1\sigma_2}\sigma_1 k$,再利用速度定义 $v_n(E)=\frac{\hbar k}{m}$,得到 $I_{mn}^{\sigma_1\sigma_2}(E,E')=\delta_{mn}\delta_{\sigma_1\sigma_2}\sigma_1 e v_n(E)$。若假定是低频 $E'=E+\hbar\omega$,则$E\simeq E'$,上述积分可近似为
    \[
    I_{mn}^{\sigma_1\sigma_2}(E,E')=\delta_{mn}\delta_{\sigma_1\sigma_2}\sigma_1 e \sqrt{v_{\alpha m}(E)v_{\alpha n}(E')}
    \]
    注意到所有$a^\dagger b,b^\dagger a$项涉及的积分是$\sigma_1\neq\sigma_2$,故这些项为零。最后电流算符表达为
    \[
    \hat{I}_\alpha(t)=\frac{e}{(2\pi\hbar)^2}\sum_m\iint dE dE'\;e^{i(E-E')t/\hbar}\left[\hat{a}_{\alpha m}^\dagger(E)\hat{a}_{\alpha m}(E')-\hat{b}_{\alpha m}^\dagger(E)\hat{b}_{\alpha m}(E')\right]
    \]
    $\hat{a}_{\alpha m}^\dagger\hat{a}_{\alpha m}$为入射流的粒子数密度,$\hat{b}_{\alpha m}^\dagger\hat{b}_{\alpha m}$则为出射流的粒子数密度。净电流由通过端口的电子数确定。
    对电流作Fourier变换
    \begin{align*}
    \hat{I}_\alpha(\omega)&=\int dt\;\hat{I}_\alpha(t)e^{i\omega t}\\
    &=\frac{1}{2\pi}\frac{e}{\hbar^2}\int dE dE'\sum_m\left[\hat{a}_{\alpha m}^\dagger(E)\hat{a}_{\alpha m}(E')-\hat{b}_{\alpha m}^\dagger(E)\hat{b}_{\alpha m}(E')\right]\int e^{i(E-E')t/\hbar}e^{i\omega t}dt\\
    \hat{I}_\alpha(\omega)&=\frac{e}{2\pi\hbar}\sum_m\iint dE dE'\;\delta(E-E'+\hbar\omega)\left[\hat{a}_{\alpha m}^\dagger(E)\hat{a}_{\alpha m}(E')-\hat{b}_{\alpha m}^\dagger(E)\hat{b}_{\alpha m}(E')\right]\\
    &=\frac{e}{h}\sum_m\int dE\left[\hat{a}_{\alpha m}^\dagger(E)\hat{a}_{\alpha m}(E+\hbar\omega)-\hat{b}_{\alpha m}^\dagger(E)\hat{b}_{\alpha m}(E+\hbar\omega)\right]
    \end{align*}
    矩阵表达形式 $\hat{\boldsymbol{a}}_\alpha^\dagger=\{\hat{a}_{\alpha m}^\dagger\}$,求和变为矩阵求迹,电流公式有简洁的表达式
    \[
    \hat{I}_\alpha(\omega)=\frac{e}{h}\int dE\;\text{Tr}\left[\hat{\boldsymbol{a}}_\alpha^\dagger(E)\hat{\boldsymbol{a}}_\alpha(E+\hbar\omega)-\hat{\boldsymbol{b}}_\alpha^\dagger(E)\hat{\boldsymbol{b}}_\alpha(E+\hbar\omega)\right]
    \]
    关于低频情况,能量$E$接近Fermi能,若$k\sim q$,则要求$\hbar\omega\ll E_F$,对通常金属$\omega\sim 10^{10}Hz$也满。由于波函数一般为$\psi\sim a\chi^++b\chi^-$,其中$a$为入射波振幅,$b$为出射波振幅。根据散射矩阵的定义 $b=S a$,考虑二次量子化场算符表达下,$\hat\psi=\hat{a}\chi^++\hat{b}\chi^-$,故也有$\hat{\boldsymbol{b}}=\boldsymbol{S}\hat{\boldsymbol{a}}$,得到
    \[
    \hat{I}_\alpha=\frac{e}{h}\int dE\;\text{Tr}\left[\hat{\boldsymbol{a}}_\alpha^\dagger(E)\hat{\boldsymbol{a}}_\alpha(E)-\sum_{\beta\gamma}S_{\alpha\beta}^\dagger\hat{\boldsymbol{a}}_\beta^\dagger(E)S_{\alpha\gamma}\hat{\boldsymbol{a}}_\gamma(E)\right]
    \]
    于是直流情形($\omega=0$)的多端口Büttiker公式便导出为
    \[
    \hat{I}_\alpha(\omega)=\frac{e}{h}\int dE\;\text{Tr}\left[\hat{\boldsymbol{a}}_\alpha^\dagger(E)\hat{\boldsymbol{a}}_\alpha(E)-\hat{\boldsymbol{b}}_\alpha^\dagger(E)\hat{\boldsymbol{b}}_\alpha(E)\right]
    \]
    此时 $\hat{b}_{\alpha m}=\sum_{\beta n}S_{\alpha m,\beta n}\hat{a}_{\beta n}\;,\;\hat{\boldsymbol{b}}_\alpha=S_{\alpha\beta}\hat{\boldsymbol{a}}_\beta$,$\boldsymbol{S}$即为散射矩阵。
    注意到流过端口的电流因为电流守恒的缘故必然流过器件中心区,而上面电流时在端口内计算的,可观测量必须是统计平均值,在端口内做统计平均 $I_\alpha=\langle\hat{I}_\alpha\rangle$。假设散射与自旋无关,因此散射矩阵中不含自旋指标,所以引入自旋简并因子2得到
    \begin{align*}
    I_\alpha&=\frac{2e}{h}\int dE\;\text{Tr}\left[\langle\hat{\boldsymbol{a}}_\alpha^\dagger\hat{\boldsymbol{a}}_\alpha\rangle-\sum_{\beta\gamma}S^\dagger_{\alpha\beta}S_{\alpha\gamma}\langle\hat{\boldsymbol{a}}_\beta^\dagger\hat{\boldsymbol{a}}_\gamma\rangle\right]\\
    &=\frac{2e}{h}\int dE\;\sum_{\beta}\text{Tr}\left[\delta_{\alpha\beta}\delta_{\alpha\gamma}-S^\dagger_{\alpha\beta}S_{\alpha\gamma}\right]\langle\hat{\boldsymbol{a}}_\beta^\dagger\hat{\boldsymbol{a}}_\gamma\rangle\\
    \langle\hat{\boldsymbol{a}}_\beta^\dagger\hat{\boldsymbol{a}}_\gamma\rangle&=f_\beta\delta_{\beta\gamma}\;,\;f_\beta=\frac{1}{e^{\frac{E-\mu_\beta}{k_B T}}+1}\\
    I_\alpha&=\frac{2e}{h}\int dE\;\sum_{\beta\gamma}\text{Tr}\left[\delta_{\alpha\beta}\delta_{\alpha\gamma}-S^\dagger_{\alpha\beta}S_{\alpha\gamma}\right]f_\beta\delta_{\beta\gamma}\\
    &=\frac{2e}{h}\sum_\beta\int dE\;\text{Tr}\left[\delta_{\alpha\beta}-S^\dagger_{\alpha\beta}S_{\alpha\beta}\right]f_\beta
    \end{align*}
    定义矩阵 $A_{\alpha\beta}=\delta_{\alpha\beta}-S^\dagger_{\alpha\beta}S_{\alpha\beta}$,直流Büttiker公式简化为
    \[
    I_\alpha=\frac{2e}{h}\sum_\beta\int dE\;\text{Tr}A_{\alpha\beta}f_\beta
    \]

 

后才能发言