非平衡格林函数模型：单粒子negf方程概要

《非平衡格林函数模型：单粒子negf方程概要》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、非平衡格林函数模型：单粒子NEGF方程概要MagnusPaulsson微纳米技术系，NanoDTU，丹麦技术大学2006年1月4日摘要非平衡格林函数方法常用于计算纳米尺寸的电导器件（包括分子和半导体）外加偏压后的电流及电荷密度。该方法主要用于处理弹道输运，但扩展后也可处理存在非弹性散射的输运问题。本文将尽可能清晰的导出NEGF方法中计算电流及电荷密度矩阵的几个重要方程，并加以解释。1、引论非平衡和林函数方法常用于计算纳米尺寸的电导器件（包括分子和半导体）外加偏压后的电流及电荷密度。关于分子电子学的一般性理论可参看文献1，有关半导体纳米器件的相关理论可参看文献2。本文的目的是

2、使读者对于单电子格林函数以及由此得到的电流、电荷密度矩阵表达式有一个直观的理解。本文在内容上并不求全责备，只作为文献1-4的补充。2、格林函数分立的薛定鄂方程(1)我们把系统的哈密顿量以及波函数按三个子空间分割：器件以及两个接触电极(2)其中描述器件与接触电极之间的相互作用。在这里我们假定不同的接触电极之间是相互独立的，亦即，不同的接触间没有相互作用项。我们将格林函数定义为其他文献中的定义可能（也的确有）符号跟我们相反：(3)2.1为什么要计算格林函数？如果在薛定鄂方程中加入一个恒定微扰项，则格林函数正反映了系统对这种微扰的响应。具体来看，薛定鄂方程：(4)对于微扰的响应：

3、(5)(6)为什么我们想要知道针对这种微扰的响应呢？因为通常说来，计算格林函数都比直接处理本征值问题要容易（由下一节的阐述可看出）特别是对无限长系统，更妙的是，体系的大多数性质（对于单粒子系统是所有性质）都可以通过格林函数来得到。举个例子，我们来算算接触电极的波函数。器件中的波函数作为已知条件，根据方程(2)的第三行，有：(7)(8)(9)其中就是接触2孤立时的格林函数。值得注意的是，因为我们的系统无限长，所以我们会从格林函数的定义式中解出两个解当电子的能量与接触电极的能带匹配时，会有两个解，分别对应接触中的入射波和出射波。我们把它们分别称作超前解和推迟解在实际操作中，我们

4、通过在能量上加入一个虚部，继而将虚部对0求极限来得到两个不同的解。极限得到的就是推迟解，得到的就是超前解。我们可以通过傅利叶变换将格林函数变化到时域来检验这一结论。，分别对应接触中的入射波和出射波。标注：我们将用G来表示推迟格林函数，来表示超前格林函数（偶尔可能会使用）。大写的G表示整体系统的格林函数，以及它的子空间，比如。小写的G用来表示孤立子系统的格林函数，比如。还有一点需要注意的是，在方程(9)中，如果使用推迟格林函数,得到的波函数对应接触电极中的一个出射波；如果使用超前格林函数，得到就是入射波。2.2自能计算格林函数的一大原因是，它比求解薛定鄂方程更简单。另外还有一

5、点：想知道器件的格林函数并不需要求解整个系统的格林函数G。从格林函数的定义式可得：(10)选取第二列对应的三个方程：(11)(12)(13)通过方程(11)和(13)可解出：(14)(15)代入方程(12)中：(16)从中易得：(17)其中和即我们所说的自能。粗略的说来，接触对于器件的影响可以看作在器件自身的哈密顿量上加入了一个附加的自能，于是我们计算器件的格林函数时直接采用有效哈密顿量即可。但是，我们仍需保持清醒，这种概念只在计算格林函数的时候才适用。有效哈密顿量的本征值和本征矢并没有直接的物理意义。对于“正常”的接触，电极表面的格林函数和通常利用周期性条件来计算得出。该

6、方法在文献3的附录B和文献2的第三部分中有详细介绍。2.3谱函数格林函数的另一个重要用途是由它来得到谱函数：(18)它给出了态密度和薛定鄂方程的所有解！为了说明这一点，我们首先指出的是，对于任意一个微扰，薛定鄂方程的解有两个（和）：(19)两个解分别由超前格林函数和推迟格林函数给出：(20)(21)两个解的差正是薛定鄂方程的解：(22)这也就是说，任意矢量，左乘谱函数A之后得到的，都是薛定鄂方程的一个解。接下来证明通过谱函数可以得到薛定鄂方程的所有解，稍微有些麻烦，我们需要将格林函数按薛定鄂方程的本征矢进行展开：(23)式中的代表虚部的一个小量（参看注释3），是H的一组正交

7、完备的本征矢，分别对应着本征能量。我们也将谱函数按这一组基展开：(24)(25)(26)令，我们得到：(27)方程(27)说明谱函数可给出薛定鄂方程的所有解。3对入射波的响应对于非平衡问题，不同的接触电极可看作化学势不同的电子库，它们向器件内注入电子，填充能级，对应来自各个电极的入射波。因此，我们希望搞清器件如何响应这些入射波。考虑电极1孤立时的情形。此时薛定鄂方程的解对应入射波在电极边界被完全反射的情形。我们把这些解标记为，其中1为电极编号，n为量子数（电极中存在多种波模）。我们通过孤立电极的谱函数可以把这些波函