DFT-10复杂体系的O(N)算法.ppt

第10章 复杂体系的O(N)算法 1。引言 2。O(N)算法的物理基础 －量子力学局域性 3。O(N)算法的基本策略 4。DFT框架下的O(N)算法 5。计算流程和主要步骤 1。引言 Order-N算法或O(N)算法的必要性 目前，DFT第一性原理计算方法，如fplapw, fplmto, Car-Parrinello, 从头赝势以及许多量子化学计算方法，对于由大量原子组成的复杂体系已经不能满足需要。 原因是以上传统方法的数值运算工作量（操作数）～Nat3。即体系的原子数增加一倍，必须消耗8倍cpu时间。 研究计算操作数与体系原子数成比例的方法，即O(N)算法对于研究复杂体系十分迫切。 本章着重与分析O(N)算法的物理基础、实现O(N)算法的基本策略以及把O(N)算法纳入DFT框架的方法。 2。O(N)算法的物理基础－量子力学局域性 Kohn的“近视原理(near-sightedness principle)” （Kohn, PRL 76, 3168 （1996）） Kohn证明了如下原理：多电子体系的某部分的物理性质，不因远离它的区域有势的变化而受影响。 量子力学局域性 例：大多数量子力学静态性质有局域性： 分子或固体中的化学键 局域态密度 电荷密度分布 局域磁矩 结合能，。。。 它们都只依赖于几个近邻原子“壳层”内的局域环境。 局域性的描述 一般采用广义Wannier函数（GWF,wi, 非周期系的局域Wannier函数）构造密度矩阵： 3。O(N)算法的基本策略 如何实现O(N)算法？ O(N)算法的基本策略 考虑到处理波函数和密度矩阵的具体要求，已经提出了多种O(N)算法方案： Fermi算符展开方法（FOE） Fermi算符投影方法（FOP） 分治（Divide and conquer, DC）方法 密度矩阵最小化方法（DMM） 轨道最小化方法（OM） 优化基密度矩阵最小化方法（OBDMM） McWeeny净化算法 McWeeny提出一种将接近等幂的密度矩阵 变换为更接近等幂的密度矩阵?的算法： LNV密度矩阵最小化方法 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) 线性标度的HGG方法-1 HGG（Hernández-Gillan-Goringe）方法属于自洽第一性原理方法，并与LNV方法密切相关。 方法的特点： 基态的描述：把DFT中关于总能Etot是KS占有轨道?i 或电子密度的泛函，等价地表述为密度矩阵的泛函。并要求密度矩阵是等幂的。 采用局域化支持函数（support function）?i?和空间有限的变分参数矩阵Li?j?来表示密度矩阵。 用变分法求总能关于支持函数和Li?j? 均为最小。HGG方法采用的是固定电子数而不是固定化学势。计算上更为方便。 线性标度的HGG方法-2 用KS占有轨道定义密度矩阵 线性标度的HG方法-3 净化之后称为等幂的密度矩阵 线性标度的HG方法-4 以下的步骤就是采用变分法，在电子数固定的条件下，求总能关于支持函数和L矩阵为最小，由此得到真正基态能量的上限。 目前的O(N)算法，仅限于基态性质的研究。 4。DFT框架下的O(N)算法 以上基本原理的实际执行，可以在LDA近似下采用赝势法。但是，是在实空间的网格点上计算。 以每一个原子为中心，取半径为Rreg的球作为支持区。每一个支持区包含一定数量v的支持函数，并且各区的v都一样。 实际执行表明，支持函数的总数 ? 0.5Nel（val）。 在原来的方法中，每一个支持函数?i?(r)都用它在该区的网格点rl 上之值?i?(rl)表示。后来发现这一方法在动能计算精度及不同的网格点总能出现不连续性等问题。新方法中采用一种局域函数将?i?展开。 支持函数的表达式 支持函数用局域化的基函数展开，他们称这种基函数为“视点函数（blip function)”。 DFT框架下的O(N)算法-2 主要计算公式： DFT框架下的O(N)算法-3 通过变分法使总能最小，在HGG方法中，采用共轭梯度近似，涉及如下解析式： 以及 矩阵乘积中Hi?j?是支持函