最近在学习 如何以DFT第一性原理计算出的高对称路径上的能量本征值为数据源,进行晶体TB模型参数的拟合,为有效学习,记录之。
本文以Marzari发表在RMP上的文章:Maximally localized Wannier functions: Theory and applications 为主要参考文献,行文顺序大致也和引文一致,加以自己的理解,有误之处欢迎各位网友指正。
一、introduction:(导言很多都是直译,待完全理解后适当更改)
在有N个初基原胞的晶体中,当规定周期性边界条件时,可得到N个平移不变操作组成的N阶平移群G,且在此平移群中,每个元素都自成一个共轭类。因此平移群有N个不可约表示,每个不可约表示的维数都为1,可找到平移群的一个1次实表示(不可约表示的维数):
其中k为周期性边界条件限定的倒空间中N个波矢,用k来区分平移算符体现在波函数中(解本征方程)就得到bloch定理,详细可见李正中老师的《固体理论》p.5。
由bloch定理,我们可以用n,k来作为满足平移不变系统的量子数,其中n为主量子数,由求解有限区域(单个元胞)内的厄密本征值问题带来(无穷多分立解),k为倒空间波矢,由平移不变性带来。
在求解多体问题时,通常考虑单电子近似,我们首选单电子能量本征态的bloch波函数
wannier 函数(WF)提供了这个可能(关于wannier函数可见一般固体物理书籍),在紧束缚近似(TB)中便是利用原子波函数来直接替代wannier函数(实际上由此反推出的类bloch态并严格为系统能量本征态)。bloch表象到wannier表象变换的意义便在于次,牺牲能量上的局域性换取空间上的局域性。
WFs最初用于化学上的分子轨道计算“localized molecular orbitals"(LMOs),空间上的局域性使得它在表征化学键时无比方便。而在基于Kohn-Sham eq.的DFT第一性原理计算的应用上发展缓慢。其中一大困难在于对应同一组bloch态的WFs选择并不唯一,这是由于bloch函数的相位不确定性导致的(即 在bloch态中作任意幺正变换,其表征的物理系统的观测量并不改变(规范自由))。更重要的,高对称K点上的bloch能级简并会造成更大的麻烦。
为确定一套WFs用于计算,Marzari and Vanderbilt (1997)拟定了"maximal-localization"的标准,为WFs的应用做出一大贡献。Marzari和V证明了局部泛函的最小化,对应于每一个Wannier电荷密度在其自身电荷中心的二阶矩分布的和的最小化,这在计算上是易于处理的。
Souza, Marzari, and Vanderbilt (2001)将这个方法进行了推广,按需要建立了由WFs张成的子空间(如占据态所处空间)。
值得一提的是,bloch的berry phase 和 WFs的中心电荷有关。
下面我们介绍一下文章结构:
sec. 2:介绍BFs与WFs之间的变换。讨论他们的规范自由(gauge freedom),讨论构建WFs的两种方法:projection 和 maximal localization。A 'disentangling procedure' for constructing WFs for a nonisolated set of bands (e.g., in metals) is also described.
sec. 3:介绍构建过程中使用的不同的标准以及不用的算法,并讨论'downfolding'和 linear-scaling methonds之间的关系。
sec. 4:简要证明WFs在分析晶体化学键中的重要性,尤其是WFs在于电子极化上的应用。
sec. 5:讨论Wannier表象在描述 轨道磁化, NMR化学转变,轨道磁电响应,以及TI中的运用。
sec. 6:介绍Wannier插值法(可将粗糙的计算结果扩展成更细致的网格)。
sec. 7:讨论WFs在量子输运,半经验势的推导以及强关联系统中的运用。
sec. 8:简要介绍WFs在除了电子结构以外的应用,如声子晶体,cold atoms in optical lattices。
sec. 9:summary and conclusions.