skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Stochastic many-body perturbation theory for anharmonic molecular vibrations

Abstract

A new quantum Monte Carlo (QMC) method for anharmonic vibrational zero-point energies and transition frequencies is developed, which combines the diagrammatic vibrational many-body perturbation theory based on the Dyson equation with Monte Carlo integration. The infinite sums of the diagrammatic and thus size-consistent first- and second-order anharmonic corrections to the energy and self-energy are expressed as sums of a few m- or 2m-dimensional integrals of wave functions and a potential energy surface (PES) (m is the vibrational degrees of freedom). Each of these integrals is computed as the integrand (including the value of the PES) divided by the value of a judiciously chosen weight function evaluated on demand at geometries distributed randomly but according to the weight function via the Metropolis algorithm. In this way, the method completely avoids cumbersome evaluation and storage of high-order force constants necessary in the original formulation of the vibrational perturbation theory; it furthermore allows even higher-order force constants essentially up to an infinite order to be taken into account in a scalable, memory-efficient algorithm. The diagrammatic contributions to the frequency-dependent self-energies that are stochastically evaluated at discrete frequencies can be reliably interpolated, allowing the self-consistent solutions to the Dyson equation to be obtained. Thismore » method, therefore, can compute directly and stochastically the transition frequencies of fundamentals and overtones as well as their relative intensities as pole strengths, without fixed-node errors that plague some QMC. It is shown that, for an identical PES, the new method reproduces the correct deterministic values of the energies and frequencies within a few cm{sup −1} and pole strengths within a few thousandths. With the values of a PES evaluated on the fly at random geometries, the new method captures a noticeably greater proportion of anharmonic effects.« less

Authors:
 [1];  [1];  [2]
  1. Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801 (United States)
  2. (Japan)
Publication Date:
OSTI Identifier:
22419820
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Chemical Physics; Journal Volume: 141; Journal Issue: 8; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ALGORITHMS; DEGREES OF FREEDOM; EVALUATION; FREQUENCY DEPENDENCE; MATHEMATICAL SOLUTIONS; MONTE CARLO METHOD; PERTURBATION THEORY; POTENTIAL ENERGY; SELF-ENERGY

Citation Formats

Hermes, Matthew R., Hirata, So, E-mail: sohirata@illinois.edu, and CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012. Stochastic many-body perturbation theory for anharmonic molecular vibrations. United States: N. p., 2014. Web. doi:10.1063/1.4892614.
Hermes, Matthew R., Hirata, So, E-mail: sohirata@illinois.edu, & CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012. Stochastic many-body perturbation theory for anharmonic molecular vibrations. United States. doi:10.1063/1.4892614.
Hermes, Matthew R., Hirata, So, E-mail: sohirata@illinois.edu, and CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012. Thu . "Stochastic many-body perturbation theory for anharmonic molecular vibrations". United States. doi:10.1063/1.4892614.
@article{osti_22419820,
title = {Stochastic many-body perturbation theory for anharmonic molecular vibrations},
author = {Hermes, Matthew R. and Hirata, So, E-mail: sohirata@illinois.edu and CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012},
abstractNote = {A new quantum Monte Carlo (QMC) method for anharmonic vibrational zero-point energies and transition frequencies is developed, which combines the diagrammatic vibrational many-body perturbation theory based on the Dyson equation with Monte Carlo integration. The infinite sums of the diagrammatic and thus size-consistent first- and second-order anharmonic corrections to the energy and self-energy are expressed as sums of a few m- or 2m-dimensional integrals of wave functions and a potential energy surface (PES) (m is the vibrational degrees of freedom). Each of these integrals is computed as the integrand (including the value of the PES) divided by the value of a judiciously chosen weight function evaluated on demand at geometries distributed randomly but according to the weight function via the Metropolis algorithm. In this way, the method completely avoids cumbersome evaluation and storage of high-order force constants necessary in the original formulation of the vibrational perturbation theory; it furthermore allows even higher-order force constants essentially up to an infinite order to be taken into account in a scalable, memory-efficient algorithm. The diagrammatic contributions to the frequency-dependent self-energies that are stochastically evaluated at discrete frequencies can be reliably interpolated, allowing the self-consistent solutions to the Dyson equation to be obtained. This method, therefore, can compute directly and stochastically the transition frequencies of fundamentals and overtones as well as their relative intensities as pole strengths, without fixed-node errors that plague some QMC. It is shown that, for an identical PES, the new method reproduces the correct deterministic values of the energies and frequencies within a few cm{sup −1} and pole strengths within a few thousandths. With the values of a PES evaluated on the fly at random geometries, the new method captures a noticeably greater proportion of anharmonic effects.},
doi = {10.1063/1.4892614},
journal = {Journal of Chemical Physics},
number = 8,
volume = 141,
place = {United States},
year = {Thu Aug 28 00:00:00 EDT 2014},
month = {Thu Aug 28 00:00:00 EDT 2014}
}
  • No abstract prepared.
  • Cited by 1
  • Diagrammatic many-body perturbation theory is used to calculate the electronic energy and the static electric dipole polarizability of the hydrogen molecule in its ground state. An amply extended discrete basis set of Gaussian orbitals is employed to minimize basis set errors and single-electron states are generated by the Hartree--Fock V/sup N/ potential. The correlation energy is evaluated through third order and with some higher-order corrections included by denominator shifts to recover about 95% of the total correlation energy. Dipole polarizabilities are calculated through second order in electron correlation with an accuracy of approx.2%. Also the energy-denominator decoupling theorem is explicitlymore » proved by invoking combinatorial analysis to implement extensive denominator shifts. Considering the values obtained, some comments are given on the application of partial summation techniques to molecular problems.« less
  • The spin-rotational Hamiltonian parameters A{sub ||} and A{sub perpendicular} for the BaF molecule are calculated using four-component relativistic spinors at the second-order many-body perturbation theory (MBPT) level via the Z-vector technique. The second-order MBPT is applied to assess the accuracy of the computed hyperfine-structure constants before studying the problem with the state-of-the-artcoupled cluster with single and double excitations (CCSD) method which is highly accurate but computationally more expensive than MBPT. The hyperfine-structure constants A and A{sub d} resulted from these calculations agree favorably well with experimental findings and with other correlated calculations. The convergence behavior of A and A{sub d}more » with respect to the number of active orbitals used in the perturbative calculations suggests that our estimated A and A{sub d} values should be accurate.« less
  • A formalism for energy-dependent many-body perturbation theory (MBPT), previously indicated in our recent review articles [Lindgren et al., Phys. Rep. 389, 161 (2004); Can. J. Phys. 83, 183 (2005)], is developed in more detail. The formalism allows for an arbitrary mixture of energy-dependent (retarded) and energy-independent (instantaneous) interactions and hence for a merger of quantum-electrodynamics (QED) and standard (relativistic) MBPT. This combined many-body-QED procedure is particularly important for light elements, such as light heliumlike ions, where electron correlation is pronounced. It can also be quite significant in the medium-heavy mass range, as recently discussed by Fritzsche et al. [J. Phys.more » B 38, S707 (2005)], with the consequence that the effects might be significant also in analyzing the data of experiments with highly charged ions. A numerical procedure is described, and some preliminary results are given for heliumlike ions with a single retarded photon. This represent the first numerical evaluation of the combined many-body-QED effect on an atomic system. It is found that for heliumlike neon the effect of one retarded photon with correlation represents more than 99% of the nonradiative effects beyond energy-independent MBPT. The new procedure also allows for the inclusion of radiative effects (self-energy and vacuum polarization) in a more systematic fashion than has previously been possible.« less