DFT复杂体系的ON算法.pptxVIP

会计学1DFT复杂体系的ON算法 21。引言Order-N算法或O(N)算法的必要性 目前，DFT第一性原理计算方法，如fplapw, fplmto, Car-Parrinello, 从头赝势以及许多量子化学计算方法，对于由大量原子组成的复杂体系已经不能满足需要。 原因是以上传统方法的数值运算工作量（操作数）～Nat3。即体系的原子数增加一倍，必须消耗8倍cpu时间。研究计算操作数与体系原子数成比例的方法，即O(N)算法对于研究复杂体系十分迫切。本章着重与分析O(N)算法的物理基础、实现O(N)算法的基本策略以及把O(N)算法纳入DFT框架的方法。第1页/共27页 32。O(N)算法的物理基础－量子力学局域性Kohn的“近视原理(near-sightedness principle)” （Kohn, PRL 76, 3168 （1996）） Kohn证明了如下原理：多电子体系的某部分的物理性质，不因远离它的区域有势的变化而受影响。?r’Δv(r’)r考虑一个量子多粒子系，在r处的静态物理性质为F，它依赖于r周围线度为?的体积内的坐标，?为de Broglie波长量级。Kohn证明，F对于r’处势的变化Δv(r’)是不敏感的。所以，r处的势保持不变，但比?更远的区域会变。Δv(r)=0第2页/共27页 4量子力学局域性例：大多数量子力学静态性质有局域性：分子或固体中的化学键局域态密度电荷密度分布局域磁矩结合能，。。。 它们都只依赖于几个近邻原子“壳层”内的局域环境。2。局域性的描述 主要方案：采用局域化的Wannier函数和密度矩阵方法。 Wannier函数的衰减行为： 有带隙的绝缘体（1D, 3D, 无序，团簇，缺陷和表面），都有 指数衰减行为第3页/共27页 5局域性的描述一般采用广义Wannier函数（GWF,wi, 非周期系的局域Wannier函数）构造密度矩阵：N是体系每个自旋的电子数。因为wi是局域化的，?将按|r-r’|衰减。对于绝缘体和金属， |r-r’|都表现出指数衰减率。T?0时，衰减甚至更迅速。在DFT下，核心问题是使?成为一个投影算符，其作用是把它投射到占有态空间。这在数学上等价于要求? 必须是等幂的 （Idempotent）,即要求其本征值在（0,1）区间。如何把一个接近等幂的密度矩阵 变为等幂矩阵?将在下面介绍。(10.1)第4页/共27页 63。O(N)算法的基本策略如何实现O(N)算法？vivi’根据Kohn近视原理，把体积为V的体系分成N个子体积vi(i-1..N)vi~?3. 在vi处取体积为vi’的区域，它包括vi和一个缓冲区，然后解出每一个vi’的性质。如果vi’vi，那么vi内的性质是相当精确的。由于计算每一个vi的工作量完全独立于体系的大小，只要知道vi’内的资料即可。整个体系的大小~vi的数目N, 于是得到线性标度算法。V=N*vi第5页/共27页 7O(N)算法的基本策略考虑到处理波函数和密度矩阵的具体要求，已经提出了多种O(N)算法方案：Fermi算符展开方法（FOE）Fermi算符投影方法（FOP）分治（Divide and conquer, DC）方法密度矩阵最小化方法（DMM）轨道最小化方法（OM）优化基密度矩阵最小化方法（OBDMM）采用Chebyshev多项式将DM展开杨伟涛教授与DFT密切结合第6页/共27页 8McWeeny净化算法McWeeny提出一种将接近等幂的密度矩阵 变换为更接近等幂的密度矩阵?的算法：用?和 分别表示?和 的本征值，这两个本征值的关系是可见所以，这种映射迭代将驱使本征值趋于0或1，由此得到符合等幂要求的?。(10.2)(10.3)第7页/共27页 9LNV密度矩阵最小化方法Li, Nunes, Vanderbilt (LNV)提出DMM方法，文献上常称LNV方法。它所采用的净化方法有完全不同的方式。其线性标度是通过对密度矩阵的空间截断得到的。 Ref. Li, Nunes, Vanderbilt : Phys. Rev. B47, 10891 (1993)LNV方法已经在TB方法的框架下得到广泛应用。采用化学势固定使总能最小的方法。(后来发现，固定化学势方法在实际计算上并不是最方便的)。 Ref. Nunes, Vanderbilt : Phys. Rev. B50, 17611 (1994)第8页/共27页 10线性标度的HGG方法-1HGG（Hernández-Gillan-Goringe）方法属于自洽第一性原理方法，并与LNV方法密切相关。方法的特点：基态的描述：把DFT中关于总能Etot是KS占有轨道?i