计算机科学 ›› 2023, Vol. 50 ›› Issue (6): 81-85.doi: 10.11896/jsjkx.220500252
刘人僪1,2, 徐直前2,3, 商红慧2, 张云泉2
LIU Renyu1,2, XU Zhiqian2,3, SHANG Honghui2, ZHANG Yunquan2
摘要: 针对密度泛函微扰理论中响应密度矩阵的计算问题,提出了一种全新的Sternheimer方程的并行求解方法,即通过共轭梯度算法和矩阵直接分解算法对Sternheimer方程进行求解,并且在第一性原理的分子模拟软件FHI-aims中实现了这两种算法。实验结果表明采用共轭梯度算法和矩阵直接分解算法的计算结果精度较高,相比传统方法的计算结果误差较小,且具有可扩展性,验证了新的Sternheimer方程中线性方程求解的正确性和有效性。
中图分类号:
[1]SHAM L J,KOHN W.One-particle properties of an inhomogeneous interacting electron gas[J].Physical Review,1966,145(2):561. [2]ALZAHRANI A Z.Structural and electronic properties of graphene upon molecular adsorption:DFT comparative analysis[J].Graphene Simulation,2011,1:21-38. [3]CHAKRABORTY D,CHATTARAJ P K.Conceptual density functional theory based electronic structure principles[J].Chemical Science,2021,12(18):6264-6279. [4]BESLEY N A.Modeling of the spectroscopy of core electrons with density functional theory[J].Wiley Interdisciplinary Reviews:Computational Molecular Science,2021,11(6):e1527. [5]GONZE X.First-principles responses of solids to atomic dis-placements and homogeneous electric fields:Implementation of a conjugate-gradient algorithm[J].Physical Review B,1997,55(16):10337. [6]GONZE X,LEE C.Dynamical matrices,Born effective charges,dielectric permittivity tensors,and interatomic force constants from density-functional perturbation theory[J].Physical Review B,1997,55(16):10355. [7]BARONI S,DE GIRONCOLI S,DAL CORSO A,et al.Phonons and related crystal properties from density-functional perturbation theory[J].Reviews of Modern Physics,2001,73(2):515. [8]OCHSENFELD C,HEAD-GORDON M.A reformulation of the coupled perturbed self-consistent field equations entirely within a local atomic orbital density matrix-based scheme[J].Chemical Physics Letters,1997,270(5/6):399-405. [9]FRISCH M,HEAD-GORDON M,POPLE J.Direct analyticSCF second derivatives and electric field properties[J].Chemical Physics,1990,141(2/3):189-196. [10]LIANG W Z,ZHAO Y,HEAD-GORDON M.An efficient approach for self-consistent-field energy and energy second derivatives in the atomic-orbital basis[J].The Journal of Chemical Physics,2005,123(19):194106. [11]SHANG H,CARBOGNO C,RINKE P,et al.Lattice dynamics calculations based on density-functional perturbation theory in real space[J].Computer Physics Communications,2017,215:26-46. [12]STERNHEIMER R M.Electronic polarizabilities of ions from the Hartree-Fock wave functions[J].Physical Review,1954,96(4):951. [13]MOTORNYI O,VAST N,TIMROV I,et al.Electron energyloss spectroscopy of bulk gold with ultrasoft pseudopotentials and the Liouville-Lanczos method[J].Physical Review B,2020,102(3):035156. [14]HOFMANN F,SCHELTER I,KÜMMEL S.Efficiently evaluating the Krieger-Li-Iafrate and common-energy-denominator approximations in the frequency-dependent Sternheimer scheme[J].Physical Review A,2019,99(2):022507. [15]KAMEARI A.Improvement of ICCG convergence for thin elements in magnetic field analyses using the finite-element method[J].IEEE Transactions on Magnetics,2008,44(6):1178-1181. [16]HESTENES M R,STIEFEL E.Methods of conjugate gradients for solving linear systems[M].Washington:NBS,1952. |
|