Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1909.06456

Scalable approximation of Green's function for estimation of anharmonic energy corrections

Prashant Rai, Khachik Sargsyan, Habib N. Najm, So Hirata

Published Sep 13, 2019• 0 citations
arXivHTMLPDF

Abstract

A method based on separated integration to estimate anharmonic corrections to energy and vibration of molecules in a second-order diagrammatic vibrational many-body Green's function formalism has already been presented. A severe bottleneck in extending this approach to bigger molecules is that the storage of the Green's function scales exponentially with the number of atoms in the molecule. In this article, we present a method that overcomes this limitation by approximating the Green's function in the Hierarchical Tucker tensor format. We illustrate that the storage cost is linear in dimension and hence one can obtain accurate representations of the Green's function for a molecule of any size. Application of this method to estimate the second-order correction to energy of molecules illustrates the advantage of this approach.

DOI
10.48550/arxiv.1909.06456

Related papers

Low-rank canonical-tensor decomposition of potential energy surfaces: application to grid-based diagrammatic vibrational Green's function theory

Prashant Rai, Khachik Sargsyan, Habib N. Najm +2 more

A new method is proposed for a fast evaluation of high-dimensional integrals of potential energy surfaces (PES) that arise in many areas of quantum dynamics. It decomposes a PES into a canonical low-rank tensor format, reducing its integral into a relatively short sum of products of low-dimensional integrals. The decomposition is achieved by the alternating least squares (ALS) algorithm, requiring only a small number of single-point energy evaluations. Therefore, it eradicates a force-constant evaluation as the hotspot of many quantum dynamics simulations and also possibly lifts the curse of dimensionality. This general method is applied to the anharmonic vibrational zero-point and transition energy calculations of molecules using the second-order diagrammatic vibrational many-body Green's function (XVH2) theory with a harmonic-approximation reference. In this application, high dimensional PES and Green's functions are both subjected to a low-rank decomposition. Evaluating the molecular integrals over a low-rank PES and Green's functions as sums of low-dimensional integrals using the Gauss–Hermite quadrature, this canonical-tensor-decomposition-based XVH2 (CT-XVH2) achieves an accuracy of 0.1 cm−1 or higher and nearly an order of magnitude speedup as compared with the original algorithm using force constants for water and formaldehyde.

Molecular Physics•Mar 7, 2017•8 citations

Diagrammatic theories for the vibrational many-body problem

Jacob A. Faucheaux

Anharmonic vibrational many-body methods are developed for and applied to small molecules and extended systems in a bound potential energy surface (PES). Diagrammatically size-consistent and basis-set-free vibrational coupled-cluster (XVCC) theory for both zero-point energies and transition frequencies, the latter through the equation-of-motion (EOM) formalism, is defined for an $n$th-order Taylor-series PES. Quantum-field-theoretical tools (the rules of normal-ordered second quantization and Feynman--Goldstone diagrams) for deriving their working equations are established. The equations of XVCC and EOM-XVCC are derived and implemented with the aid of computer algebra. Algorithm optimizations known as strength reduction, intermediate reuse, and factorization are carried out before code generation,producing algorithms with optimal cost scaling. A similarity-transformed equation-of-motion vibrational coupled-cluster (STEOM-XVCC) method is introduced as a one-mode theory with an effective vibrational Hamiltonian, which is similarity transformed twice so that its lower-order operators are dressed with higher-order anharmonic effects. From diagonalization of this doubly similarity-transformed Hamiltonian in the small one-mode excitation space, the method simultaneously computes accurate anharmonic vibrational frequencies of all fundamentals, which have unique significance in vibrational analyses. We establish a diagrammatic method of deriving the working equations of STEOM-XVCC and prove their connectedness and thus size-consistency as well as the exact equality of its frequencies with the corresponding roots of EOM-XVCC. An extended STEOM-XVCC (Ext-STEOM-XVCC) method is defined as an $m$th order configuration interaction method with the doubly similarity-transformed Hamiltonian including up to $m$th-order excitation operators. Because the doubly transformed Hamiltonian is dressed with higher-order anharmonic effects, the frequencies of overtones and combinations obtained are different and superior to the corresponding EOM-XVCC method. We compare and contrast the Ext-STEOM-XVCC method to its electronic counterpart. We apply the previously established second-order size-extensive vibrational many-body perturbation (XVMP2) method to the anharmonic phonon dispersion curves of a model two-mass system and the optical phonons of polyethylene. We obtain accurate results despite the presence of multiple Fermi-resonances in the crystalline systems.

•Nov 28, 2017•0 citations

Stochastic Self-Consistent Second-Order Green’s Function Method for Correlation Energies of Large Electronic Systems

Daniel Neuhauser, Roi Baer, Dominika Zgid

The second-order Matsubara Green's function method (GF2) is a robust temperature-dependent quantum chemistry approach, extending beyond the random-phase approximation. However, until now the scope of GF2 applications was quite limited as they require computer resources that rise steeply with system size. In each step of the self-consistent GF2 calculation there are two parts: estimating of the self-energy from the previous step's Green's function and updating the Green's function from the self-energy. The first part formally scales as the fifth power of the system size, while the second has a much gentler cubic scaling. Here, we develop a stochastic approach to GF2 (sGF2), which reduces the fifth power scaling of the first step to merely quadratic, leaving the overall sGF2 scaling as cubic. We apply the method to linear hydrogen chains with up to 1000 electrons, showing that the approach is numerically stable, efficient, and accurate. The stochastic errors are very small, on the order of 0.1% or less of the correlation energy for large systems, with only a moderate computational effort. The first iteration of GF2 is an MP2 calculation that is done in linear scaling; hence we obtain an extremely fast stochastic MP2 (sMP2) method as a byproduct. While here we consider finite systems with large band gaps where at low temperatures effects are negligible, the sGF2 formalism is temperature dependent and general and can be applied to finite or periodic systems with small gaps at finite temperatures.

arXiv:1707.08296•Journal of Chemical Theory and Computation•Sep 29, 2017•54 citations

Spectral properties from Matsubara Green's function approach: Application to molecules

Michael Schüler, Y. Pavlyukh

We present results for many-body perturbation theory for the one-body Green's function at finite temperatures using the Matsubara formalism. Our method relies on the accurate representation of the single-particle states in standard Gaussian basis sets, allowing to efficiently compute, among other observables, quasiparticle energies and Dyson orbitals of atoms and molecules. In particular, we challenge the second-order treatment of the Coulomb interaction by benchmarking its accuracy for a well-established test set of small molecules, which includes also systems where the usual Hartree-Fock treatment encounters difficulties. We discuss different schemes how to extract quasiparticle properties and assess their range of applicability. With an accurate solution and compact representation, our method is an ideal starting point to study electron dynamics in time-resolved experiments by the propagation of the Kadanoff-Baym equations.

arXiv:1710.08660•Physical review. B./Physical review. B•Mar 30, 2018•26 citations

Computing many-body wave functions with guaranteed precision: The first-order Møller-Plesset wave function for the ground state of helium atom

Florian A. Bischoff, Robert J. Harrison, Edward F. Valeev

We present an approach to compute accurate correlation energies for atoms and molecules using an adaptive discontinuous spectral-element multiresolution representation for the two-electron wave function. Because of the exponential storage complexity of the spectral-element representation with the number of dimensions, a brute-force computation of two-electron (six-dimensional) wave functions with high precision was not practical. To overcome the key storage bottlenecks we utilized (1) a low-rank tensor approximation (specifically, the singular value decomposition) to compress the wave function, and (2) explicitly correlated R12-type terms in the wave function to regularize the Coulomb electron-electron singularities of the Hamiltonian. All operations necessary to solve the Schrödinger equation were expressed so that the reconstruction of the full-rank form of the wave function is never necessary. Numerical performance of the method was highlighted by computing the first-order Møller-Plesset wave function of a helium atom. The computed second-order Møller-Plesset energy is precise to ~2 microhartrees, which is at the precision limit of the existing general atomic-orbital-based approaches. Our approach does not assume special geometric symmetries, hence application to molecules is straightforward.

The Journal of Chemical Physics•Sep 10, 2012•45 citations

Efficient computation of the second-Born self-energy using tensor-contraction operations

Riku Tuovinen, Fabio Covito, Michael A. Sentef

In the nonequilibrium Green’s function approach, the approximation of the correlation self-energy at the second-Born level is of particular interest, since it allows for a maximal speed-up in computational scaling when used together with the generalized Kadanoff-Baym ansatz for the Green’s function. The present day numerical time-propagation algorithms for the Green’s function are able to tackle first principles simulations of atoms and molecules, but they are limited to relatively small systems due to unfavorable scaling of self-energy diagrams with respect to the basis size. We propose an efficient computation of the self-energy diagrams by using tensor-contraction operations to transform the internal summations into functions of external low-level linear algebra libraries. We discuss the achieved computational speed-up in transient electron dynamics in selected molecular systems. ACKNO

•Nov 7, 2019•6 citations
A method based on separated integration to estimate anharmonic corrections to energy and vibration of molecules in a second-order diagrammatic vibrational many-body Green's function formalism has already been presented. A severe bottleneck in extending this approach to bigger molecules is that the storage of the Green's function scales exponentially with the number of atoms in the molecule. In this article, we present a method that overcomes this limitation by approximating the Green's function in the Hierarchical Tucker tensor format. We illustrate that the storage cost is linear in dimension and hence one can obtain accurate representations of the Green's function for a molecule of any size. Application of this method to estimate the second-order correction to energy of molecules illustrates the advantage of this approach.
A new method is proposed for a fast evaluation of high-dimensional integrals of potential energy surfaces (PES) that arise in many areas of quantum dynamics. It decomposes a PES into a canonical low-rank tensor format, reducing its integral into a relatively short sum of products of low-dimensional integrals. The decomposition is achieved by the alternating least squares (ALS) algorithm, requiring only a small number of single-point energy evaluations. Therefore, it eradicates a force-constant evaluation as the hotspot of many quantum dynamics simulations and also possibly lifts the curse of dimensionality. This general method is applied to the anharmonic vibrational zero-point and transition energy calculations of molecules using the second-order diagrammatic vibrational many-body Green's function (XVH2) theory with a harmonic-approximation reference. In this application, high dimensional PES and Green's functions are both subjected to a low-rank decomposition. Evaluating the molecular integrals over a low-rank PES and Green's functions as sums of low-dimensional integrals using the Gauss–Hermite quadrature, this canonical-tensor-decomposition-based XVH2 (CT-XVH2) achieves an accuracy of 0.1 cm−1 or higher and nearly an order of magnitude speedup as compared with the original algorithm using force constants for water and formaldehyde.
Anharmonic vibrational many-body methods are developed for and applied to small molecules and extended systems in a bound potential energy surface (PES). Diagrammatically size-consistent and basis-set-free vibrational coupled-cluster (XVCC) theory for both zero-point energies and transition frequencies, the latter through the equation-of-motion (EOM) formalism, is defined for an nnnth-order Taylor-series PES. Quantum-field-theoretical tools (the rules of normal-ordered second quantization and Feynman--Goldstone diagrams) for deriving their working equations are established. The equations of XVCC and EOM-XVCC are derived and implemented with the aid of computer algebra. Algorithm optimizations known as strength reduction, intermediate reuse, and factorization are carried out before code generation,producing algorithms with optimal cost scaling. A similarity-transformed equation-of-motion vibrational coupled-cluster (STEOM-XVCC) method is introduced as a one-mode theory with an effective vibrational Hamiltonian, which is similarity transformed twice so that its lower-order operators are dressed with higher-order anharmonic effects. From diagonalization of this doubly similarity-transformed Hamiltonian in the small one-mode excitation space, the method simultaneously computes accurate anharmonic vibrational frequencies of all fundamentals, which have unique significance in vibrational analyses. We establish a diagrammatic method of deriving the working equations of STEOM-XVCC and prove their connectedness and thus size-consistency as well as the exact equality of its frequencies with the corresponding roots of EOM-XVCC. An extended STEOM-XVCC (Ext-STEOM-XVCC) method is defined as an mmmth order configuration interaction method with the doubly similarity-transformed Hamiltonian including up to mmmth-order excitation operators. Because the doubly transformed Hamiltonian is dressed with higher-order anharmonic effects, the frequencies of overtones and combinations obtained are different and superior to the corresponding EOM-XVCC method. We compare and contrast the Ext-STEOM-XVCC method to its electronic counterpart. We apply the previously established second-order size-extensive vibrational many-body perturbation (XVMP2) method to the anharmonic phonon dispersion curves of a model two-mass system and the optical phonons of polyethylene. We obtain accurate results despite the presence of multiple Fermi-resonances in the crystalline systems.
The second-order Matsubara Green's function method (GF2) is a robust temperature-dependent quantum chemistry approach, extending beyond the random-phase approximation. However, until now the scope of GF2 applications was quite limited as they require computer resources that rise steeply with system size. In each step of the self-consistent GF2 calculation there are two parts: estimating of the self-energy from the previous step's Green's function and updating the Green's function from the self-energy. The first part formally scales as the fifth power of the system size, while the second has a much gentler cubic scaling. Here, we develop a stochastic approach to GF2 (sGF2), which reduces the fifth power scaling of the first step to merely quadratic, leaving the overall sGF2 scaling as cubic. We apply the method to linear hydrogen chains with up to 1000 electrons, showing that the approach is numerically stable, efficient, and accurate. The stochastic errors are very small, on the order of 0.1% or less of the correlation energy for large systems, with only a moderate computational effort. The first iteration of GF2 is an MP2 calculation that is done in linear scaling; hence we obtain an extremely fast stochastic MP2 (sMP2) method as a byproduct. While here we consider finite systems with large band gaps where at low temperatures effects are negligible, the sGF2 formalism is temperature dependent and general and can be applied to finite or periodic systems with small gaps at finite temperatures.
We present results for many-body perturbation theory for the one-body Green's function at finite temperatures using the Matsubara formalism. Our method relies on the accurate representation of the single-particle states in standard Gaussian basis sets, allowing to efficiently compute, among other observables, quasiparticle energies and Dyson orbitals of atoms and molecules. In particular, we challenge the second-order treatment of the Coulomb interaction by benchmarking its accuracy for a well-established test set of small molecules, which includes also systems where the usual Hartree-Fock treatment encounters difficulties. We discuss different schemes how to extract quasiparticle properties and assess their range of applicability. With an accurate solution and compact representation, our method is an ideal starting point to study electron dynamics in time-resolved experiments by the propagation of the Kadanoff-Baym equations.
We present an approach to compute accurate correlation energies for atoms and molecules using an adaptive discontinuous spectral-element multiresolution representation for the two-electron wave function. Because of the exponential storage complexity of the spectral-element representation with the number of dimensions, a brute-force computation of two-electron (six-dimensional) wave functions with high precision was not practical. To overcome the key storage bottlenecks we utilized (1) a low-rank tensor approximation (specifically, the singular value decomposition) to compress the wave function, and (2) explicitly correlated R12-type terms in the wave function to regularize the Coulomb electron-electron singularities of the Hamiltonian. All operations necessary to solve the Schrödinger equation were expressed so that the reconstruction of the full-rank form of the wave function is never necessary. Numerical performance of the method was highlighted by computing the first-order Møller-Plesset wave function of a helium atom. The computed second-order Møller-Plesset energy is precise to ~2 microhartrees, which is at the precision limit of the existing general atomic-orbital-based approaches. Our approach does not assume special geometric symmetries, hence application to molecules is straightforward.
In the nonequilibrium Green’s function approach, the approximation of the correlation self-energy at the second-Born level is of particular interest, since it allows for a maximal speed-up in computational scaling when used together with the generalized Kadanoff-Baym ansatz for the Green’s function. The present day numerical time-propagation algorithms for the Green’s function are able to tackle first principles simulations of atoms and molecules, but they are limited to relatively small systems due to unfavorable scaling of self-energy diagrams with respect to the basis size. We propose an efficient computation of the self-energy diagrams by using tensor-contraction operations to transform the internal summations into functions of external low-level linear algebra libraries. We discuss the achieved computational speed-up in transient electron dynamics in selected molecular systems. ACKNO