格林函数与输运
格林函数及其应用课件
有限差分法
01
有限差分法是将微分方程或积分 方程转化为差分方程，然后求解 差分方程得到格林函数的数值解 。
02
有限差分法适用于求解偏微分方 程，特别是对于具有周期性或对 称性的问题，有限差分法可以大 大简化计算过程。
有限元法
有限元法是将微分方程或积分方程转化为有限元方程，然后求解有限元方程得到 格林函数的数值解。
对于某些领域，需要高精度的格林函数来保证计 算的准确性。
未来格林函数研究的方向与展望
算法优化
寻求更高效、稳定的算法来计算格林函数。
多领域交叉
加强与其他领域的合作，拓展格林函数的应用范围。
数值稳定性
研究如何提高格林函数计算的数值稳定性。
感谢观看
THANKS
量子力学散射问题的格林函数计算
总结词
介绍了量子力学散射问题中格林函数的 计算方法，以及其在散射理论中的应用 。
VS
详细描述
在量子力学中，格林函数用于描述粒子在 相互作用下的运动行为。通过计算格林函 数，可以研究粒子在散射过程中的能量和 动量变化，进一步理解物质的微观结构和 相互作用机制。
流体动力学波动问题的格林函数计算
工程学
在电路分析、控制理论和信号 处理等领域有广泛应用。
生物学
用于研究神经网络的传播和扩 散过程。
金融学
用于描述资产价格波动和风险 评估。
当前格林函数计算中存在的问题与挑战
高维问题
随着问题维度的增加，格林函数的计算变得极为 复杂。
不适定性
在实际应用中，格林函数的求解可能存在数值不 稳定性。
精度要求
有限元法适用于求解复杂的偏微分方程，特别是对于具有复杂边界条件的问题， 有限元法可以更好地处理边界条件。
格林函数方法)
第七章 格林函数方法
第一节 前言
从五十年代开始，量子场论中的格林函数方法被用于研究统计物理学中的 问题。到六十年代后期，格林函数理论在固体物理等多个领域得到了进一步的拓 展，被认为是一种强有力的数学工具[1]。例如，对许多准粒子问题，只需知道相 互作用过程中少数粒子的初态与末态间的跃迁振幅（相应的格林函数），就能得 到体系的一些特征，而对于固体物理中的很多问题，只有对应于费米能量附近的 系统格林函数与我们要研究的性质有关。这样，格林函数方法就成为研究系统性 质的直接有效的方法。 但是在很多的实际问题中，如一些较复杂的有限尺寸量子系统，要得出其格 林函数的解析表达式是很困难的，因此必须要通过数值计算来解决。格点格林函 数方法[2-17]是通过把系统分离成一些格点，然后通过计算这些格点及格点间的格 林函数，进而得出整个体系的格林函数的一种有效数值计算方法。它与其他的一 些数值方法如有限元法[18]、转移矩阵法[19,20]、散射矩阵法[21]、模式匹配法[22]等相 比较，格点格林函数方法能够很方便的处理磁场和无序（掺杂）等问题。在系统 的某个区域加入磁场时，只需要考虑一个 Peierls 相位因子。当系统的自由度很 大时， 用一般的格点格林函数方法求解系统的格林函数就对应一个很大维数的矩 阵计算。虽然计算机技术飞速发展，但是计算机的容量仍然制约着我们所能直接 处理的矩阵的维数。在这种情况下，迭代技术已经被越来越广泛地应用于处理这 一类问题。 递归格林函数方法也在这种要求下得到了很大的发展。 Lee、 Fisher[2,3] 和 MacKinnon[7]等作了开创性的工作，然后人们又发展了各种递归格林函数方法 来处理一些具体的结构或边界条件下的尺寸效应和多终端效应。 如 Soles[5,6]等用 递归格林函数方法计算了有 T-型突起的量子线的电子输运性质， Ando[9]则考虑了 在磁场调制下的量子点接触的电导。 在量子物理中，格林函数常常被定义为 v v v v v [ E − H (r )]G (r , r ' ; E ) = δ (r − r ' ) 其中 E 是复变量，H 是一个厄米的含时算符。 如果 E − H 的本征值是非零的，我们可以写出格林函数的等价定义式： 1 G= (2-2) E−H 如果 H 的本征函数ψ n 是正交完备的，且 λ n 是其相应的本征值，则
数学物理方法 12 格林函数法
（14.2.4）
(r r0 ) 代表三维空间变量的 函数，在直角坐标系中其形式为
(r r0 ) ( x x0 ) ( y y0 ) ( z z0 )
（14.2.4）式中 函数前取负号是为了以后构建格林函数方便 格林函数的物理意义【2】：在物体内部（T 内） r0 处放置一个单位点电荷，而该物体的界面保持电位为零, 那么 该点电荷在物体内产生的电势分布，就是定解问题(14.2.4)的 解――格林函数．由此可以进一步理解通常人们为什么称格林 函数为点源函数．
（14.3.8）
上式正是我们所熟知的静电场的电位表达式
二、二维轴对称情形
用单位长的圆柱体来代替球．积分在单位长的圆柱体内进行，即
G(r ,0)dV (r )dV
T T
因为
(r )dV 1
T
G(r ,0)dV G(r,0)dV
令积分常数为0，得到
G(r ,0)
1 1 G(r ,0) ln c 2π r
1 1 ln 2π r
因此二维轴对称情形的格林函数为
G(r , r0 ) 1 1 ln 2π | r r0 |
（14.3.9）
将（14.3.9）代入式（14.3.1）得到二维无界区域的解为
u (r ) 1 1 f (r0 )ln dS0 S0 2π | r r0 |
点场的影响的总和．两项积分中的格林函数相同．这说明 泊松方程的格林函数是点源在一定的边界条件下所产生的 场． 对于拉普拉斯方程 第一边值问题的解为 第三边值问题的解为
f (r0 ) 0
u (r ) (r0 )
G (r , r0 ) ]dS0 n 0
格林函数法
为第三边值问题的积分表示式
物理意义：右边第一个积分表示区域T中分布的源在r 点产生的场的总和；第二个积分代表边界上的状况对 r点场的影响的总和；两项积分中的格林函数相同。 说明泊松方程的格林函数是点源在一定的边界条件下 所产生的场。
对于拉普拉斯方程，f(r0)=0，因此可得拉普拉斯 方程第一边值问题的解
因此，我们可设想一个等效的点电荷，它位 于球外M1处，且在球面产生的电势与球内点电荷 在球面产生的电势相反。由物理学知识可知，该 设想的点电荷必位于OM0处的延长线上，如图所 示，并记：
OM r, OM0 r0
在∑ε 上的解，该解表示位于球心r=r0处的电量为ε0的 点电荷在半径为ε的球面上产生的电势，根据电磁学 知识，该电势为：
1
G(r, r0 ) 4
因此我们可得∑ε面上的积分
Ò
u(r)
G n
G
u(r) n
dS
Ò
u(r
)
n
(1
4
)
1
4
u(r) n
dS
Ò
u(
r
)
n
(1
4
)
1
4
u(r n
)
2d
(r r0 ) (x x0) ( y y0) (z z0)
格林函数的物理意义：在物体内部（T内）处放置 一个单位点电荷（或热源），而该物体的界面保持 电位为零（或温度为零）, 那么该点电荷（或该点 热源）在物体内产生的电势分布（或稳定温度分 布），就是上述定解问题的解――格林函数。
格林函数互易定理： 格林函数代表r0处的点源在r处 所产生的影响，系统不变，则该影响等同于将移至r 处的该点源在r0处产生影响。故格林函数遵守如下 的互易定理：
多体量子力学中的格林函数方法
多体量子力学中的格林函数方法多体量子力学是研究多粒子体系中粒子之间相互作用的力学理论。
在这个理论框架下，我们需要处理多个粒子的波函数，同时考虑它们之间的相互作用。
为了解决这个问题，物理学家们提出了多种方法，其中一种重要的方法就是格林函数方法。
格林函数方法最早由德国物理学家赫尔曼·哈库斯（Hermann Hankel）于1859年提出，后来由多位物理学家进一步发展和推广。
格林函数可以用来描述量子态的演化和性质，是求解多体问题的有力工具。
在多体量子力学中，格林函数是描述粒子行为的函数。
它可以用来计算不同时间和位置下粒子的性质，比如粒子的动量、位置和电荷等。
格林函数的形式由一般的波函数演化方程决定。
它可以被分为两个部分：单粒子格林函数和相互作用格林函数。
单粒子格林函数描述了一个单粒子在外势场下的行为。
它可以被定义为粒子在某个时刻从一个位置传播到另一个位置的概率幅。
通过计算单粒子格林函数，可以得到粒子的一些重要性质，比如能谱和态密度等。
相互作用格林函数描述了多个粒子之间的相互作用。
在多体问题中，粒子之间的相互作用是一个非常重要的因素。
通过计算相互作用格林函数，可以探究粒子之间的相互作用强度和方式。
相互作用格林函数的求解可以通过一系列的近似方法，比如平均场理论、扰动方法和重整化群等。
格林函数方法在各个领域都有广泛的应用。
在凝聚态物理中，格林函数方法可以用来研究电子系统和其他凝聚态物理体系的性质。
通过计算格林函数，可以得到电子的输运性质、激发态和自能等重要信息。
格林函数方法在量子化学、固体物理、统计物理和粒子物理等领域也都有着重要的应用。
虽然求解格林函数的问题是一个复杂的任务，但是近年来，在计算机科学和数值方法的发展下，越来越多的精确和高效的方法被提出。
比如，基于数值求解的格林函数方法、基于图像处理的格林函数方法、基于机器学习的格林函数方法等。
这些方法为求解多体问题提供了新的思路和工具。
总结起来，格林函数方法是解决多体量子力学问题的一种重要方法。
数学物理方法格林函数
演化问题的格林函数
演化问题的格林函数也可以用冲量定理法得到 问题 等价问题
Gt a 2 G 0 G |x 0 G |x L 0 G | t 0 ( x )
Gtt a 2 G 0 G |x 0 G |x L 0 G |t 0 0 G | t t 0 ( x )
演化问题的基本解
无界输运问题的求解
2 ut a u xx f ( x, t ) u |t 0 0
f ( x, t ) d d f ( , ) ( x ) (t )
0
t
2 Gt a G ( x ) (t ) G |t 0 0
2 ( x ) t exp 2 4a (t ) u d d f ( , ) 2a ( t ) 0
u( x, t ) d d f ( , )G( x, ; t, )
0
t
( x ) 2 exp 2 4 a ( t ) G 2a ( t )
应用（求解数学物理方程的格林函数法）
稳定问题的基本解
稳定问题的基本解可以利用静电场类比法得到 原问题 方程
u f ( r )
点源问题
G ( r r ' )
点电荷电场
V q (r r ' ) / 0
解
u
f (r ' )d ' 1 q G V 4 | r r ' | 4 | r r ' | 4 0 | r r ' |
格林函数
第十四章 格林函数--偏微分方程解的积分表示解偏微分方程主要有两种方法：A: 数理方法中的分离变量法：正交的多项式或无穷级数解，但需要齐次边界条件。
B: 理论物理中的Green 函数方法：既是简单的有理形式解，又允许任意的边界条件！ 1，Green 函数（GF ）的意义：物理上：点源产生的场（函数）在时空中的分布。
特别是它在空间是源函数；在时空是传播函数。
(See below)数学上: 具有点源的偏微分方程在齐次边界条件或者无界区域、初值条件下的解。
2，GF 的分类：边界值GF ：(,')G r r 即源函数；初始值GF ：(,;',')G r t r t 即传播函数。
3，Green 函数的性质：1）对称性：(,')(',)G r r G r r =，它与定解问题相关，即与厄米性相关。
(See 4 below)2）时间传播函数没有对称性：(,;',')(',';,)G r t r t G r t r t ≠.（因果律引起） 3）存在的必要条件：设方程2()(,')(')G r r r r λδ∇+=--，若λ是对应齐次方程的本征值，即2ϕλϕ∇=- 和附加齐次边界条件，则(,')G r r 不存在。
这是因为既有点源：(')r r δ-矛盾于又无流：|0.n G ∂∑∂= 本征值问题存在，但是没有激发，物理上自相矛盾！平面波(),ik xat Ae 球面波1()ik xat Ar e -和柱面波1/2()ik at A e ρρ-均是LaplaceEquation 的解，但不是Possion Equation 的解。
球、柱面波分别来自于1x时（散射问题）渐近行为：(1,2)[(1/2)/2](1,2)[(1)/2]21(),().i x m i x l ml xxH x e h x e πππ±-+±-+4，Green 函数的边值条件：选取边值条件具有人为性，但要求简单并保证算子的厄米性。
第十二章 格林函数法
故得到
( x ) G ( x x ) ( x )d G ( x x ) ( x ) 0 G ( x x ) ( x ) ds n n S
V
这就是用Green函数求解静电问题的一种形式解。 讨论几点：
12
该式左边第二项为 1 1 ( x) ( x x )d ( x )
0
V
0
得到
1 1 G ( x x ) ( x )d ( x ) 0 0 V G ( x x ) ( x ) G ( x x ) ( x ) ds n n S
2 0
2
2、镜像法能解的情况：在求解区域内没有自由电荷， 或者只有有限几个点电荷，并且区域边界或介质界面 规则(电场能用等效电荷代替)+边界条件。
3
二、 Green函数法能解的情况
能用Green定理求解静电边值问题的情况: 给定区域V内电荷分布 (x ) 和区域V的边界面S 上各点的电势 φs 或电势法向导数
1 2 G ( x , x ) ( x x) 0 G ( x , x) 1 G ( x , x) 0, 或 S n 0S S
所在的位置， x 代表观察点，在（3）式和（4）式中，
(5)
7
五、Green公式和边值问题的解
G 0 在一个单位电荷在空间所激发的电势。因此 n S 即代表单位电荷在边界上所激发的电场，由Gauss定 理知道
1 G( x x )ds n 0 S n G( x x )ds 0 S G( x x ) 0 n S
非平衡格林函数方法
非平衡格林函数方法
非平衡格林函数方法是一种量子力学计算方法，用于研究非平衡态下的电子结构和输运性质。
它通过求解非平衡格林函数来描述系统的电子态和输运性质。
非平衡格林函数是描述非平衡态下的电子密度矩阵和电子自能的重要工具。
在非平衡态下，电子系统中存在着电子的注入和抽出，因此电子系统的密度矩阵和自能不再是平衡态下的对角化态。
非平衡格林函数方法通过求解非平衡态下的格林函数，可以得到体系的电子密度矩阵和自能，从而研究体系的输运性质。
非平衡格林函数方法可以用于研究各种体系的输运性质，如半导体器件、分子器件、纳米结构等。
该方法的优点在于可以考虑电子-电子相互作用和电子-声子相互作用等非平衡效应，可以得到更为准确的结果。
非平衡格林函数方法的实现需要使用一系列数学工具，如Keldysh路径积分、费曼图等。
这些工具的使用使得非平衡格林函数方法的计算复杂度较高，但是在研究非平衡态下的电子输运性质时，该方法是一种非常有效的计算工具。
总之，非平衡格林函数方法是一种重要的量子力学计算方法，可以用
于研究非平衡态下的电子结构和输运性质。
在未来的研究中，非平衡格林函数方法将继续发挥重要作用，推动纳米电子学、分子电子学等领域的发展。
格林函数在输运中的应用
I(t)的傅立叶变换
∞ dEdE ' 2e R I β (ω ) = Re{∫ ∑ [Gmn ( E + ω , E ')Σ T<,β nm ( E ', E ) −∞ (2π ) 2 h mn < + Gmn ( E + ω , E ')Σ Ta, β nm ( E ', E )]}
∞ dEdE ' 2e I β (ω ) = Re{∫ Tr[G R ( E + ω , E ')Σ T<, β ( E ', E ) −∞ (2π ) 2 h
G R , A = [ E − ε 0 − Σ R , A ]−1
Σ R , A, < = Σ L R , A, < + Σ R R , A, < 注意这里的 Σ β R , A,< 不同于上边的 Σ T , β R , A,<
相互作用 自能
只有在没有相互作用时它们才相等。通常 Σ R , A,< = Σ T R , A,< + Σ I R , A,<
§20.3 二端输运问题
IL
L
C
R
无相互作用 的自由电子
−i t dt ε (t ) Σ (t , t ') = miθ (±t m t ')∑ V (t )Vα k ,n (t ') exp ∫ 1 α k 1 t' k 0 < ∗ −i t dt ε (t ) Σ T ,mnα (t , t ') = i ∑ f (ε α k )Vm,α k (t )Vα k ,n (t ') exp ∫ 1 α k 1 t' k
什么是格林函数(Green's function)
一般地，点源作用产生的场就是格林函数。
在地震学中，格林函数是单位集中脉冲力产生的场，可以是位移，速度或加速度等，一般指位移场。
集中意味着力只作用于空间中一点，脉冲指力只作用于时间中某一时刻。
在地震学中，应特别注意：1) 集中脉冲型单力产生的位移场是格林函数；2) 一对单力组成的力偶产生的位移场是格林函数空间导数；3) 断层剪切位错所产生的位移场，等效于双力偶所产生的位移场，也等效于单力+单力偶所产生的位移场。
（见《定量地震学》等效体力章节，即3.2节）。
注：单力偶就是一般意义上的力偶，代表一对单力组成的力偶；双力偶是指两个单力偶的组合。
1 什么是格林函数对线性算子 L ，在点源 \delta 作用下的输出（或响应）就是格林函数G，即： LG=\delta 。
不同线性算子对应不同物理问题，也就对应不同性质的方程，如拉普拉斯方程，泊松方程，亥姆霍兹方程，波动方程等，这些方程都对应着各自不同的格林函数（见第二部分Wikipedia汇总）。
如，对声波波动问题，线性算子为 L=\frac{\partial^2}{\partial t^2}-c^2 \nabla^2 .格林函数妙处在于若已知格林函数与源分布（包括时间上与空间上），则可通过格林函数与源的卷积求得在此源作用下系统的输出（或响应）。
郭敦仁先生曾讲：“从物理上看，一个数理方程表示一种特定的场和产生这种场的源之间的关系（如热传导方程表示温度场和热源的关系），而格林函数则代表了一个点源所产生的场。
知道了一个点源的场，就可以用叠加的方法算出任意源的场。
”推导：已知： L\varphi=Q ,其中 L 是线性算子，Q 为源分布， \varphi 为待求输出。
利用卷积的性质，可得： \varphi=\varphi *\delta=\varphi * (LG)=(L\varphi) * G=Q*G .（注：卷积的实质就是把所有源的作用都通过积分叠加起来）因此，问题的关键就是求格林函数。
非平衡格林函数在电子输运问题中的应用
摘要 ： 简单讨 论 了非平衡格 林函数理 论 ， 并利 用非 平衡格 林函数理 论推 导了电子相 互作用体 系的电流佘 式 ， 获
得 了一些一般 性的结果 。 中图分 类号 ： ４ １７ Ｏ ７ ． ． 文献标识码 ： Ａ 关键 词 ： 非平衡格 林函数 ； 电流公 式 ； 电子 输运 ； 态密度 函数
０ 引言
现代纳米技术和材料制备的发展， 得 电子 使 器件越做越小。在这些器件 中， 量子效应和 电荷 量子化效应越来越 明显。这时 ， 电子器件处于非
１ 非平衡格林 函数
在非平衡态 ， 可以像定义平衡 格林 函数一样 定 义 一 种 编 时 格 林 函 数 一闭 路 编 时 格 林 函 数 （ｏｔ ｒ ｏ ｅ） ｃｎｏ — ｒ ｒ 或非平 衡格林 函数 。满足类似 ｕ ｄ 的戴森方程 ， 在非平 衡状态下的哈密顿量中的某
第 ６期
刘 志勇 等 ： 非平 衡 格林 函数 在 电子输 运 问题 中的应 用
・３ ７ Ｉ．
Ｃ ￡ ￡）＝ Ｊ ｄ［ ￡ ￡ （ ，， （ ，， １１ ｔ （ ， ￡ ￡）＋ Ａ １） １１
（１ｔＢ （ ｔ） ｔ，） 。ｔ ］ 数 和推迟 格林 函数 。 同样 有 （） ３ 对非 平 衡 格林 函数 ， 们 还需 要 超前 格 林 函 人
２． ｏ ｅ ｅｏ ｈ ｓｃ Ｃ ￣ ｇ ｆＰ ｙ ｉｓ＆ Ｃ ｍｍｕ ｉａ ｏ ｅ ｔ ｎｃ ，ｉｎ ｘ ｏｍａ ｉｅ ｉ ，ｉｎ ｘ ｎ ｈ ｎ ３ ０ ２Ｐ ｏ ｎｃ ｔ ｎＥｌｃｒｉ ｓＪａ ｇ ｉＮ ｒ ｌＵｎｖｒ ｔ Ｊａ ｇｉＮａ ｅ ａ ｇ３ ０ ２ ＲＣ； ｉ ｏ ｓ ｙ
Ｋｅ ｙ ｗｏｒ ｓ： ｎｅ ｕ ｌｂ ｉ ｍ ｅ ｎ ｆ ｎ ｔ ｎ， ｒｅ ｏ ｍｕ ａ， ｅ ｔｏ ｒ ｎ ｐ ｒ ， ｎ ｉｙ ｆ ｎ ｔ ｎ ｄ Ｎｏ ｑ ｉｉ ｒｕ Ｇｒ ｅ ｕ ｃｉ ｏ Ｃｕ ｒ ｎｔｆｒ ｌ Ｅｌｃ ｒ ｎ ｔａ ｓ ｏ ｔ Ｄｅ ｓｔ ｕ ｃｉ ｏ
数学物理方程第四章_格林函数
1 ⎧ ⎪∆G (r , r0 ) = − δ (r − r0 ) ε ⎨ ⎪G Γ = 0 ⎩
(4.3.7) (4.3.8)
以 G (r , r0 ) 乘式 (4.3.5), u (r ) 乘式 (4.3.7), 二式相减后在 Ω 上对 r 积分 ,以 dr 表示 r 点处的体积微元,有
∫
Ω
(G∆u − u∆G )dr = −
第 4 章 格林函数
在这一章里,我们介绍数学物理方程中另外一种常用的方法—格林函数法.从物理上看, 一个数学物理方程是表示一种特定的“场”和产生这种场的“源”之间的关系.例如,热传导 方程表示温度场和热源之间的关系,泊松方程表示静电场和电荷分布的关系,等等.这样,当源 被分解成很多点源的叠加时,如果能设法知道点源产生的场,利用叠加原理,我们可以求出同 样边界条件下任意源的场,这种求解数学物理方程的方法就叫格林函数法.而点源产生的场就 叫做格林函数. 4.1
⎧0, T ( x) = ⎨ ⎩∞,
x≠0 x=0
且
∫Байду номын сангаас
所以有
+∞
−∞
cρT ( x)dx = Q
T ( x) =
Q δ ( x) cρ
通过以上两个例题,我们对 δ ( x) 有了进一步的认识.如果将坐标平移 x0 ,即集中量 出现在点 x = x 0 处,则有
δ ( x − x0 ) = ⎨
且
⎧0, ⎩∞,
∫
= ∫ (u∆v)dΩ + ∫ gradu ⋅ gradvdΩ
Ω Ω
=∫u
Γ
∂v dS ∂n
或表示为
∫
Ω
(u∆v)dΩ = ∫ u
Γ
∂v dS − ∫ gradu ⋅ gradvdΩ Ω ∂n
量子输运格林函数方法
G(t1, t2 ) = −i
TC [cλ (t1)cλ† (t1)]
≡
⎛ ⎜ ⎝
G++ G−+
G+− G−−
⎞ ⎟ ⎠
（2.2）
其中TC 是回路 C 上的复编时算符，它的作用是将回路 C 上的时间较早的算符排在
右边。+(−) 表示的是时间回路 C 的上（下）支。式（2.1）中 G++ (t1, t2 ) 和 G−− (t1, t2 )
Tp [ A(t1 ) , B (t2 )] ≡ θ (t2 − t1 ) A(t1 ) B (t2 ) ±θ (t1 − t2 ) B (t2 ) A(t1 )
+(−) 分别对应玻色子和费米子。式（2.1）中 G+− (t1, t2 ) 和 G−+ (t1, t2 ) 分别是小于
和大于格林函数： G+− ≡ G< (t1,t2 ) = i cλ† (t2 )cλ (t1)
2.1.1 非平衡格林函数的定义
非平衡格林函数是一种处理非平衡问题的有效方法，它是由平衡态格林函数
推广得到的。首先我们先简单介绍下平衡态理论。平衡态理论中所用到的时间是
定义在实时间轴上的。若系统的哈密顿量可以写为 H = H0 + H ' ，其中 H ' 是相互作
用部分，可以看成是微扰。在薛定谔绘景中 H 0 的基态
∞
∫ Cr (t1, t2 ) = dt[ Ar (t1,t)Br (t, t2 )] −∞
∞
∫ Ca (t1, t2 ) = dt[ Aa (t1, t)Ba (t,t2 )] −∞
（2.9）
同理可以得到常用到的 Lengreth 公式：如果 C(τ ,τ1) = A(τ ,τ1)B(τ1,τ ) ，则有
非平衡格林函数
非平衡格林函数非平衡格林函数（Non-equilibriumGreen'sfunctions,NEGF）是描述非平衡态下系统行为的重要工具。
它是格林函数的一种推广，广泛应用于凝聚态物理、纳米电子学、光电子学等领域。
本文将从NEGF 的基本概念、历史发展、理论框架、应用研究等方面进行介绍和分析。
一、基本概念NEGF是一种描述量子系统非平衡态下的行为的理论工具。
它是格林函数理论的一种推广，用于描述系统中的电荷、能量、自旋等自由度在时间和空间上的演化。
NEGF理论可以用来计算非平衡态下的输运性质，如电导率、热导率等，也可以用于描述非平衡态下的光学性质，如吸收谱、发射谱等。
NEGF理论的核心是非平衡态下的格林函数。
格林函数是描述量子系统中的相互作用效应的数学工具，它反映了系统中某个自由度的激发情况对其他自由度的影响。
在平衡态下，格林函数可以用来描述系统的激发态密度、热力学性质等。
在非平衡态下，格林函数则可以用来描述系统中的输运性质。
二、历史发展NEGF理论的历史可以追溯到20世纪50年代。
当时，人们开始研究电子在晶体中的输运性质，发现传统的电子输运理论无法解释一些实验现象，如局域化、能级移动等。
为了解决这些问题，人们开始研究非平衡态下的电子输运理论。
1960年代初，Kadanoff和Baym等人提出了非平衡态下的格林函数理论，为后来的NEGF理论的发展奠定了基础。
NEGF理论在20世纪80年代得到了快速发展。
当时，人们开始研究纳米电子学、光电子学等领域，需要描述非平衡态下的输运性质。
NEGF理论的优越性质得到了广泛认可，并被应用于多个领域。
目前，NEGF理论已经成为描述非平衡态下的量子系统行为的重要工具。
三、理论框架NEGF理论的核心是非平衡态下的格林函数。
在NEGF理论中，系统的哈密顿量可以表示为H=H0+V其中H0是自由哈密顿量，V是相互作用哈密顿量。
系统的演化可以用密度矩阵来描述。
在NEGF理论中，密度矩阵可以表示为ρ(t)=ρ0+δρ(t)其中ρ0是平衡态下的密度矩阵，δρ(t)是非平衡态下的扰动。
数学物理方法第12章-格林函数
∫∫ ϕ (r )G (r , r )dS . α
0 0 0 Σ
1
12.2
电像法求格林函数 第一边值问 题格林函数
∆v(r , r0 ) = δ (r − r0 )
v(r , r0 ) Σ = 0
⇒
v(r , r0 ) = G (r , r0 )
r r0
导体球内有一个点电荷 ，导体接 地。求球内电势。 电荷的存在，在导体上感应了电荷。 球内的电势为自由电荷和感应电荷电势之和。 将感应电荷的电势由一 “电像电荷”的电势表示
∂G ( r , r0 ) dS . ∂n
第一边值问 题格林函数
u ( r0 ) = ∫∫∫ G ( r , r0 ) f ( r ) dV + ∫∫ ϕ ( r )
T Σ
第三边值问题
[α
∂u + β u ] Σ = ϕ (Σ ) ∂n
∆v(r , r0 ) = δ (r − r0 )
[α ∂v + βv] Σ = 0 ∂n
1 a + 4π r − r0 r0
1 a 4π r − 2 r0 r0
2
1 1 = r − r0 r 2 − 2rr0 cos θ + r02
在球面上
∂ ∂n
Σ
=
∂ ∂r
r =a
[
∂ 1 1 2r − 2r cos θ ]Σ =− ∂n r − r0 2 (r 2 − 2rr0 cos θ + r02 )3 / 2
[
∂ ∂ G (r , r0 )] Σ = [− G (r , r0 )] z =0 ∂n ∂z
Σ
∂u ( r ) ∂v ( r , r0 ) − u (r ) ]dS . ∂n ∂n
格林函数
§2.4 格林函数法 解的积分公式在第七章至第十一章中主要介绍用分离变数法求解各类定解问题，本章将介绍另一种常用的方法——格林函数方法。
格林函数，又称点源影响函数，是数学物理中的一个重要概念。
格林函数代表一个点源在一定的边界条件和（或）初始条件下所产生的场。
知道了点源的场，就可以用迭加的方法计算出任意源所产生的场。
一、 泊松方程的格林函数法为了得到以格林函数表示的泊松方程解的积分表示式，需要用到格林公式，为此，我们首先介绍格林公式。
设u （r ）和v （r ）在区域 T 及其边界 ∑ 上具有连续一阶导数，而在 T 中具有连续二阶导数，应用矢量分析的高斯定理将曲面积分 化成体积积分.)(⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰∇⋅∇+∆=∇⋅∇=⋅∇∑TTTvdV u vdV u dV v u S d v u ϖ（12-1-1）这叫作第一格林公式。
同理，又有.⎰⎰⎰⎰⎰⎰⎰⎰∇⋅∇+∆=⋅∇∑TTvdV u udV v S d u v ϖ（12-1-2）（12-1-1）与（12-1-2）两式相减，得 亦即.)(⎰⎰⎰⎰⎰∆-∆=⎪⎭⎫ ⎝⎛∂∂-∂∂∑T dV u v v u dS n u v n vu（12-1-3）n ∂∂表示沿边界 ∑ 的外法向求导数。
（12-1-3）叫作第二格林公式。
现在讨论带有一定边界条件的泊松方程的求解问题。
泊松方程是)( ),(T r r f u ∈=∆ϖϖ（12-1-4）第一、第二、第三类边界条件可统一地表为),( M u n u ϕβα=⎥⎦⎤⎢⎣⎡+∂∂∑（12-1-5）其中 ϕ（M ）是区域边界 ∑ 上的给定函数。
α＝0，β ≠0为第一类边界条件，α ≠0，β＝0是第二类边界条件，α、β 都不等于零是第三类边界条件。
泊松方程与第一类边界条件构成的定解问题叫作第一边值问题或狄里希利问题，与第二类边界条件构成的定解问题叫作第二边值问题或诺依曼问题，与第三类边界条件构成的定解问题叫作第三边值问题。
量子输运格林函数方法
0
0
间趋于 +∞ 时，系统仍然回到初始时刻的基态，而且只相差一个相位，用公式可
以表示为[1]
∞ = S (+∞, −∞) = eiL
0
0
（2.1）
图 2.1 Contour C
在非平衡状态下，系统并不能保证其基态 在经过 S (+∞, −∞) 作用后不变。 0
人们通过将时间轴扩展到复平面上（如图 2.1 所示），引入了时间回路的概念。这 样系统就可以从 t0 = −∞ 出发沿着 t 轴演化到 t1' = +∞ （上支），然后从 t1' = +∞ 沿着 t 轴演化回到 t0 = −∞ （下支）。这样系统通过时间演化又回到了最初的基态，与平 衡态很类似，所以在这种情形下，在平衡态格林函数基础上发展起来的各种理论 仍然可以方便的使用。此时， S 算符的形式变为 SC = S (−∞, +∞)S (+∞, −∞) 。引入 回路 C 上的非平衡格林函数[2]：
t2
图 2.3 图 2.2 中的回路 C 变形为两个时间回路
其中回路 C1 在时间轴上支，回路 C2 在时间轴的下支。则式子（2.4）变为
∫ C(t1, t2 ) = dtA(t1, t)B(t, t2 ) C
∫ ∫ = dt[ A(t1+ , t)B< (t, t2− ) + dtA< (t1+ , t)B(t, t2− )]
Meir和Wingreen推导出了相互作用区域与理想电极相连时的电流公式 [6,7]。随 后人们沿着Meir的思路和步骤对电流公式进行了推广。孙庆丰等人给出了量子点 多电极（可以是正常电极也可以是超导电极）体系的电流表达式[8]。最近，李 保文等人将非平衡格林函数推广到铁磁电极-正常金属-超导电极构成的异质结 中，并且得到了Landauer-Büttiker型的电流普遍公式[9]。利用他们得出的这个公 式，我们可以用同一套理论来研究自旋相关的电流和Andreev反射电流等输运问
格林函数与输运
《多粒子物理学》读书报告：格林函数与输运内容提要：1概述；2单粒子性质的格林函数表述；3用格林函数推导迁移率中1-α项1概述 1． 1金属中电子输运特性对于金属*m e τμ-=， μσ0en -=，τ是输运驰豫时间，它的物理意义是处在某动量本征态的电子的平均寿命，即0=t 时一个处于某动量本征态的电子在τ=t 时完全失去了对其原有动量的记忆。
输运驰豫时间包括各种相互作用的贡献主要有杂质散射﹑电子-声子相互作用﹑电子-电子相互作用等等：∑=--ii 11ττ即输运驰豫时间由各种机构中i τ最小的决定。
绝对零度时，纯金属晶体中电子不受散射，具有无穷大电导。
T >0时实际金属的电阻是由电子受到杂质和晶格振动的散射引起的。
在室温时，典型金属的电阻率约为10-8Ω.m ，随着温度降低到室温以下，电阻近似线性地减小（图1,see, p.131 in Ref.[1]），在低温时水平地达到一定值。
低温时的电阻率与试样的纯度密切相关，对于高纯度的退火单晶体，约可以达到室温电阻率的10-4倍。
不纯试样中的附加电阻在整个温度范围内近似地与温度无关。
这个事实叫做马赛厄司定则（Mathiessen rule ，又翻译为马提生定则（1862））。
这个附加电阻是由于杂质引起的电子散射，在低温下它构成电阻的主要部分。
杂质散射电阻与温度无关的事实暗示出可动电子的浓度与温度无关，这与半导体中电子浓度与温度呈指数函数关系大不一样。
声子散射电阻依赖于温度，在高温时可变得很大。
这两部分电阻具有可加性，因此可分别处理。
上述金属中的杂质不含磁性杂质。
磁性杂质的散射将导致低温下电阻值的对数上升，称为近藤(Kondo)效应。
1． 2半导体输运特性半导体中的散射仍可分为电离杂质和晶格振动的散射两大类。
晶格振动的散射又分为声学波和光学波散射两种。
声学波通过两种方式散射电子：引起密度变化从而产生形变势（声学声子形变势散射）；在没有反演中心的极性晶体中引起压电极化（压电散射,长声学波明显）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《多粒子物理学》读书报告：格林函数与输运内容提要：1概述；2单粒子性质的格林函数表述；3用格林函数推导迁移率中1-α项1概述 1． 1金属中电子输运特性对于金属*m e τμ-=， μσ0en -=，τ是输运驰豫时间，它的物理意义是处在某动量本征态的电子的平均寿命，即0=t 时一个处于某动量本征态的电子在τ=t 时完全失去了对其原有动量的记忆。
输运驰豫时间包括各种相互作用的贡献主要有杂质散射﹑电子-声子相互作用﹑电子-电子相互作用等等：∑=--ii 11ττ即输运驰豫时间由各种机构中i τ最小的决定。
绝对零度时，纯金属晶体中电子不受散射，具有无穷大电导。
T >0时实际金属的电阻是由电子受到杂质和晶格振动的散射引起的。
在室温时，典型金属的电阻率约为10-8Ω.m ，随着温度降低到室温以下，电阻近似线性地减小（图1,see, p.131 in Ref.[1]），在低温时水平地达到一定值。
低温时的电阻率与试样的纯度密切相关，对于高纯度的退火单晶体，约可以达到室温电阻率的10-4倍。
不纯试样中的附加电阻在整个温度范围内近似地与温度无关。
这个事实叫做马赛厄司定则（Mathiessen rule ，又翻译为马提生定则（1862））。
这个附加电阻是由于杂质引起的电子散射，在低温下它构成电阻的主要部分。
杂质散射电阻与温度无关的事实暗示出可动电子的浓度与温度无关，这与半导体中电子浓度与温度呈指数函数关系大不一样。
声子散射电阻依赖于温度，在高温时可变得很大。
这两部分电阻具有可加性，因此可分别处理。
上述金属中的杂质不含磁性杂质。
磁性杂质的散射将导致低温下电阻值的对数上升，称为近藤(Kondo)效应。
1． 2半导体输运特性半导体中的散射仍可分为电离杂质和晶格振动的散射两大类。
晶格振动的散射又分为声学波和光学波散射两种。
声学波通过两种方式散射电子：引起密度变化从而产生形变势（声学声子形变势散射）；在没有反演中心的极性晶体中引起压电极化（压电散射,长声学波明显）。
光学波也通过两种方式散射电子：二种不等价原子之间的相对移动所引起的形变势（光学波形变势散射）；在极性晶体中伴随光学波的极化所产生的微扰势（极性光学波散射）。
后者只有纵光学波（LO 声子）才产生，横光学波不产生。
对于离子性晶体，这时电子与LO 声子形成极化子（详见后文）。
各种散射机制中，电离杂质散射﹑声学声子形变势散射﹑压电散射中驰豫时间与电子能量的关系可统一地写成[2]r E 0ττ=对于不同的散射机构0τ和r 有不同的值。
电离杂质散射r =3/2；声学声子形变势散射r =-1/2；压电散射r =1/2。
光学声子形变势散射形式比较复杂，这里略。
对于极性光学声子散射不能有驰豫时间的定义。
对于不同的散射机构，i τ随温度的变化关系不同。
电离杂质散射I τ～2/3T ；声学声子形变势散射as τ～2/3-T ；压电散射PZ τ～2/1-T 。
光学声子形变势散射形式比较复杂，这里略。
在低温下I τ值最小，所以低温下主要散射机构是电离杂质散射。
典型半导体中各种散射机构下迁移率温度的变化关系如图2（p.135 in Ref.2）所示。
一般地0en -=σμ中电子浓度0n 可通过霍尔效应精确测定，0n σ叫做霍尔迁移率，它与我们常规定义的电导率迁移率相差一个常数因子。
图3(Fig.7.9 in Ref.3)是CdTe 霍尔迁移率的温度依赖关系。
1． 3极性光学波散射长期以来电导问题是用Boltzmann 方程处理的[4]。
1958年Edward 首先将格林函数方法应用于输运问题。
Kadanoff 和Baym 用格林函数方法证明，在金属中只有k F l >>1时Boltzmann 方程才是正确的，这里k F 是费米波矢，l 是平均自由程。
现在格林函数已经用于推导很多不同系统的输运性质，包括Boltzmann 方程不适用的情况。
格林函数方法的优点是，用它可以推导出输运系数的准确表达式，然后在各种条件下作近似计算。
严格的计算必须同时考虑晶格振动和杂质对电子的散射。
单独考虑这两种散射时计算方法非常相近。
下面只考虑晶格振动（LO 声子）的散射，不考虑杂质散射，电子哈密顿为Frohlich 极化子的哈密顿[3])(12/100+-+++++∑+∑+∑=q q p qpq p q qq p pp p a a c c qM a a c c H νωε (7.2.1)2/12/3020)2()(4B m M ωπα =)11()2(02/102εεωα-=∞ B m eBp m p 22=ε 其中0ω为LO 声子的频率。
我们讨论弱耦合(1<<α)的情况，所以极化子尺寸很大，是所谓的大极化子。
1． 4极化子迁移率理论有许多关于极化子迁移率的理论[3]，例如 （1）通过求解Boltzmann 方程(BE)；（2）通过计算电流-电流关联函数（见后文）；><⎰-=)0().(1)(0u u i j j T e d i ττνωπτβωτ式中u 是(x,y,z )之一，当系统各向同性时><>=<)0().(31)0().(j j T j j T u u ττττ所以 ><⎰-=)0().(31)(0j j T e d i i ττνωπτβωτ。
令δωωi i +→则可由电流-电流关联函数得到推迟的电流-电流关联函数。
]})({Im[lim 0ωωπσωret →-=(久保公式)0→ω说明得到的是直流电导。
在久保公式中，为了求直流电导需要先计算交流电导，然后取极限0→ω。
若开始就从直流电场出发则计算要麻烦一些。
久保的上述理论又叫线性响应理论。
（3）通过计算力-力关联函数；><⎰-=)0().(31)(0F F T e d i R i ττωτβωτ令δωωi i +→则可由力-力关联函数得到推迟的力-力关联函数。
]})({Im[lim 10202ωωρωret R n e →-=力-力关联函数的严格推导是由Mahan 给出的。
（4）通过求解量子Boltzmann 方程(QBE)。
QBE 与BE 的差别是：BE 是关于分布函数 ),,(t r v f 的微分方程；而 QBE 是关于Wigner 分布函数 ),,,(t r k f ω的微分方程。
这些理论各自不同，但在弱耦合 (1<<α)和低温（10>>βω）极限下结论是一致的，这些理论都预言0000)1(2lim 0μωατμβωα≡--=-=→→e m em e B B T (7.2.2)00021ωατN =110-=βωe N(7.2.2)对于检验理论是有用的，但在同实验比较方面无能为力，原因是（1）在温度2/10≈βω范围，迁移率的计算只考虑了光学声子散射。
这是一个低温理论，我们需要能在更高温下计算迁移率的理论（2）我们感兴趣的大多数材料的极化子耦合常数都在中间耦合（31≤≤α）的范畴内。
所以我们需要一个在中间耦合下适用的计算迁移率的理论。
在这方面最成功的是费曼路径积分方法。
对于中间耦合和高温的情况，格林函数并不擅长，所以我们下面的计算仍限于低温弱耦合的情况。
我们通常采用的方法是1964年Langreth 和Kadanoff 用过的方法（1）利用格林函数从久保公式出发推导(7.2.2)式。
迁移率可写成幂级数形式++++=-221001ααααμa a a a(7.2.3)得到(7.2.2)即得到上面级数的首项（1-α项）；（2）在上面基础上获得量级为0α的所有修正项。
定义0002N m e B ωαμ-=则L-K 结果为)(61120ααμμO +-= (7.2.4)这个结果同用下面方法，即由*me τμ-=，然后将*,m τ分别按α展开：)(120αττO +=，α611*+=B m m )(2αO +得到的结果精确符合。
我们将推导L-K 公式中的第一项（即1-α项）。
2 单粒子性质 (see, Section 7.2 of Ref.[3])我们用格林函数来表述单粒子性质。
因为是在低温单声子情况下，所以用零温格林函数。
首先看电子自能。
自能算符∑是一个非局域的且与能量k E 有关的非厄密算符。
由于∑非厄密，k E 一般为复数，∑实部代表多体效应引起的能级移动，而∑虚部为粒子在该状态寿命的倒数。
电子自能的实部为),()(sin )()],(Re[202/101)1(2/12/30ααεωωεεωαωN O p p p p ++-∑-=- 在零温上式退化为)],,(3121[)],(Re[22)1(0p p p O p ωεεωεωωαω+∑-+-=(7.2.6)根据自能我们可写出有效质量*m 和重整化系数（或重整化因子）Z 的表达式)(211)211()1(211000αααεωO Z p p +-=+=∂∑∂-=--== αααεω611211311)1()(00*-≈++=∂∑∂+===p p B z m m (7.2.7)我们用到的另一个量是寿命τ，电子寿命定义为)]},(Im[2){()(1∑-=p E p p Z p τ (7.2.8)1)1()(-=∂∑∂-=p E pp Z ωε 这里Z(p)是重整化系数。
(7.2.8)式中我们计算的是p E =ω点的自能（虚部），而不是0→ω点的自能（虚部）。
E p 是粒子的基态能量，它可从方程∑+=)],(Re[p p p E p E ε中自恰计算求得。
这个方程近似为∑+=)1()],(Re[ret p p p E p E ε。
作为一级近似，E p 也可由下法得出：将E p 展开为幂级数形式)(242*20p O p mE E p ++= 显然p p E E 00lim →==∑)1(0)],0(Re[ret E 。
根据参考文献[3]第六章算出的一阶电子自能∑--=)1(2/1002/300)()],0(Re[retE E ωεω即2/1002/300)(E E --=ωεω从中展开得到020)(ααωO E +-=所以)(242*20p O p mE p ++-= αω。
由(7.2.8)式我们看到Z(p)与寿命的定义有关。
下面我们看Z(p)的物理意义。
由谱函数的定义∑+∑--∑-=-=2)][Im()],(Re {)],(Im[2)],(Im[2),(ωεωωωp p E p G p A p ret (7.2.9)在∑→0Im 的极限下，谱函数变成一个δ函数（表示能量守恒）和重整化系数的乘积)(2)()]},(Re[{2),(lim 0Im p p p Z p p A εωδπωεωπδω-∑⨯=--=∑→推迟格林函数定义为),(2)(),(ωπωθωp A e i d t t p G ti ret -∞∞-⎰= (7.2.10)根据格林函数的衰减我们可定义弛豫时间。