Development of efficient computational methods for many-body perturbation calculations of extended systems

Marco Govoni, Giulia Galli

Katoda, Michio

制度:新 ; 報告番号: 乙2324号 ; 学位の種類:博士 (理学) ; 授与年月日: 2011/7/28 ; 早大学位記番号:新5742

Zhendong Li

Treating electron correlation more accurately and efficiently is at the heart of the development of electronic structure methods. In the present work, we explore the use of stochastic approaches to evaluate high-order electron correlation energies, whose conventional computational scaling is unpleasantly steep, being O(Nn+3) with respect to the system size N and the perturbation order n for the Møller-Plesset (MP) series. To this end, starting from Goldstone’s time-dependent formulation of ab initio many-body perturbation theory (MBPT), we present a reformulation of MBPT, which naturally leads to a Monte Carlo scheme with O(nN2 + n2N + f(n)) scaling at each step, where f(n) is a function of n depending on the specific numerical scheme. Proof-of-concept calculations demonstrate that the proposed quantum Monte Carlo algorithm successfully extends the previous Monte Carlo approaches for MP2 and MP3 to higher orders by overcoming the factorial scaling problem. For the first time, Goldstone’s time-dependent formulation is made useful numerically for electron correlation energies, not only being purely as a theoretical tool.

Annabelle Canestraight, Xiaohe Lei, Khaled Z. Ibrahim +1 more

Calculations of excited states in the Green’s function formalism often invokes the diagonal approximation, in which the quasiparticle states are taken from a mean-field calculation. In this paper, we extend the stochastic approaches applied in the many-body perturbation theory and overcome this limitation for large systems in which we are interested in a small subset of states. We separate the problem into a core subspace whose coupling to the remainder of the system environment is stochastically sampled. This method is exemplified on computing hole injection energies into CO 2 on an extended gold surface with nearly 3000 electrons. We find that in the extended system the size of the problem can be compressed up to 95% using stochastic sampling. This result provides a way forward for self-consistent stochastic methods and determination of Dyson orbitals in large systems.

Ionel Stetcu, BRUCE R. BARRETT, U. van Kolck