Many-body Green Function Research Articles

A Review on Study Short-Range Correlations and the Isospin Dependence

High-energy electron scattering can be considered a crucial tool for interpreting the function of high-momentum nucleons (and quarks) in nuclei and offers a clean, accurate probing for investigations of hadronic and nuclear organization. The intense repulsion of nuclear contact at small distances is the source of short-range interactions in nuclei and nuclear matter. These connections have been verified in several nuclear tests using hadronic and electroweak probes. Improved knowledge of the underlying structure of these high-momentum and short-range arrangements in nuclei will be crucial for a variety of high-energy observables. The short-range correlations also appear to be associated with changes in the quark distributions in nuclei. Over the past ten years, there has been a significant concentration on both theoretical and experimental efforts aimed at characterizing short-range correlations, namely their isospin dependence. In this review, it will employ a theoretical framework rooted in many-body Green's functions theory. Specifically, review focus on the one-body momentum distributions as a fundamental variable. The objective is to provide a concise overview of the latest research findings in the field.

Graphene with suitable-size nitrogen-doped cavity being an excellent photocatalyst for proton-coupled electron transfer in water splitting

Aiming to attain porous carbon nitride photocatalysts with high specific surface area, abundant active sites, broad absorption range and effective electron-hole separation, we remove carbon rings from graphene and replace all the edge carbon atoms of the cavity with nitrogen atoms. The accurate electronic structures and exciton properties of the proposed materials are revealed by the ab initio many-body Green's function theory. Computational results show that the absorption of the designed materials can cover both the visible and the near-infrared light region. The exciton binding energies of the materials are one order of magnitude lower than those of the typical two-dimensional photocatalyst graphitic carbon nitride. Due to the formation of hydrogen bonds between pyridinic nitrogen atoms and water molecules during water splitting, the key excited-state proton-coupled electron transfer reactions are barrierless. These findings could serve as guidelines for the realization of high-performance metal-free two-dimensional photocatalysts.

BSE@GW Prediction of Charge Transfer Exciton in Molecular Complexes: Assessment of Self-Energy and Exchange-Correlation Dependence.

The Bethe-Salpeter equation using the GW approximation to the self-energy (BSE@GW) is a computationally attractive method for studying electronic excitation from first principles within the many-body Green's function theory framework. We examine its dependence on the underlying exchange-correlation (XC) approximation as well as on the GW approximation for predicting the charge transfer exciton formation at representative type-II interfaces between molecular systems of tetrachloro-1,2-benzoquinone (TCBQ) and acene derivatives. For the XC approximation, we consider several widely used generalized gradient approximation (GGA) and hybrid GGA functionals. For the GW self-energy approximation, we examine the recently proposed renormalized singles approach by Yang and coauthors [J. Phys. Chem. Lett. 2019, 10 (3), 447-452; J. Chem. Theory Comput. 2022, 18, 7570-7585] in addition to other commonly employed approximated GW schemes. We demonstrate a reliable prediction of the charge transfer exciton within the BSE@GW level of theory.

Reference Energies for Valence Ionizations and Satellite Transitions.

Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.

Electronic Couplings and Conversion Dynamics between Localized and Charge Transfer Excitations from Many-Body Green's Functions Theory.

We investigate the determination of electronic coupling between localized excitations (LEs) and charge-transfer (CT) excitations based on many-body Green's functions theory in the GW approximation with the Bethe-Salpeter equation (GW-BSE). Using a small molecule dimer system, we first study the influence of different diabatization methods, as well as different model choices within GW-BSE, such as the self-energy models or different levels of self-consistency, and find that these choices affect the LE-CT couplings only minimally. We then consider a large-scale low-donor morphology formed from rubrene and fullerene and evaluate the LE-CT couplings based on coupled GW-BSE-molecular mechanics calculations. For these disordered systems of bulky molecules, we observe differences in the couplings based on the Edmiston-Ruedenberg diabatization compared to the more approximate Generalize Mulliken-Hush and fragment charge difference diabatization formalisms. In a kinetic model for the conversion between LE and CT states, these differences affect the details of state populations in an intermediate time scale but not the final populations.

Machine Learning Many-Body Green's Functions for Molecular Excitation Spectra.

We present a machine learning (ML) framework for predicting Green's functions of molecular systems, from which photoemission spectra and quasiparticle energies at quantum many-body level can be obtained. Kernel ridge regression is adopted to predict self-energy matrix elements on compact imaginary frequency grids from static and dynamical mean-field electronic features, which gives direct access to real-frequency many-body Green's functions through analytic continuation and Dyson's equation. Feature and self-energy matrices are represented in a symmetry-adapted intrinsic atomic orbital plus projected atomic orbital basis to enforce rotational invariance. We demonstrate good transferability and high data efficiency of the proposed ML method across molecular sizes and chemical species by showing accurate predictions of density of states (DOS) and quasiparticle energies at the level of many-body perturbation theory (GW) or full configuration interaction. For the ML model trained on 48 out of 1995 molecules randomly sampled from the QM7 and QM9 data sets, we report the mean absolute errors of ML-predicted highest occupied and lowest unoccupied molecular orbital energies to be 0.13 and 0.10 eV, respectively, compared to GW@PBE0. We further showcase the capability of this method by applying the same ML model to predict DOS for significantly larger organic molecules with up to 44 heavy atoms.

Dynamical effects on photoluminescence spectra from first principles: A many-body Green's function approach

Dynamical effects on photoluminescence spectra from first principles: A many-body Green's function approach

Robust relativistic many-body Green's function based approaches for assessing core ionized and excited states.

A two-component contour deformation (CD) based GW method that employs frequency sampling to drastically reduce the computational effort when assessing quasiparticle states far away from the Fermi level is outlined. Compared to the canonical CD-GW method, computational scaling is reduced by an order of magnitude without sacrificing accuracy. This allows for an efficient calculation of core ionization energies. The improved computational efficiency is used to provide benchmarks for core ionized states, comparing the performance of 15 density functional approximations as Kohn-Sham starting points for GW calculations on a set of 65 core ionization energies of 32 small molecules. Contrary to valence states, GW calculations on core states prefer functionals with only a moderate amount of Hartree-Fock exchange. Moreover, modern abinitio local hybrid functionals are also shown to provide excellent generalized Kohn-Sham references for core GW calculations. Furthermore, the core-valence separated Bethe-Salpeter equation (CVS-BSE) is outlined. CVS-BSE is a convenient tool to probe core excited states. The latter is tested on a set of 40 core excitations of eight small inorganic molecules. Results from the CVS-BSE method for excitation energies and the corresponding absorption cross sections are found to be in excellent agreement with those of reference damped response BSE calculations.

Many-Body Green's Function Theory for Electronic Excitations in Complex Chemical Systems.

The GW method and the Bethe-Salpeter equation (BSE) have exhibited excellent performance in computing charged and neutral electronic excitations in materials of various dimensions in the past decades. Extensive benchmark studies have demonstrated that their precision can reach the level of high-level ab initio wave function approaches with a much lower computational cost. GW and BSE outperform the density functional theory because of both the more accurate electronic structures in GW and the capability to treat local, charge-transfer, and Rydberg excitations on the same footing in BSE. Presently, they are the only available first-principles approaches that can study electronic excitations at the surface and in the interior of periodic systems and at the interface between periodic systems with reliable accuracy. In this Perspective, a brief overview of GW and BSE and their applications in complex chemical systems is provided, with the goal of boosting their broader utilization in chemistry.

Computing the Many-Body Green's Function with Adaptive Variational Quantum Dynamics.

We present a method to compute the many-body real-time Green's function using an adaptive variational quantum dynamics simulation approach. The real-time Green's function involves the time evolution of a quantum state with one additional electron with respect to the ground state wave function that is first expressed as a linear-linear combination of state vectors. The real-time evolution and the Green's function are obtained by combining the dynamics of the individual state vectors in a linear combination. The use of the adaptive protocol enables us to generate compact ansatzes on-the-fly while running the simulation. In order to improve the convergence of spectral features, Padé approximants are applied to obtain the Fourier transform of the Green's function. We demonstrate the evaluation of the Green's function on an IBM Q quantum computer. As a part of our error mitigation strategy, we develop a resolution-enhancing method that we successfully apply on the noisy data from the real-quantum hardware.

Exploring Bethe-Salpeter Excited-State Dipoles: The Challenging Case of Increasingly Long Push-Pull Oligomers.

The change of molecular dipole moment induced by photon absorption is key to interpret the measured optical spectra. Except for compact molecules, time-dependent density functional theory (TD-DFT) remains the only theory allowing to quickly predict excited-state dipoles (μES), albeit with a strong dependency on the selected exchange-correlation functional. This Letter presents the first assessment of the performances of the many-body Green's function Bethe-Salpeter equation (BSE) formalism for the evaluation of the μES. We explore increasingly long push-pull oligomers as they present an excited-state nature evolving with system size. This work shows that BSE's μES do present the same evolution with oligomeric length as their CC2 and CCSD counterparts, with a dependency on the starting exchange-correlation functional that is strongly decreased as compared to TD-DFT. This Letter demonstrates that BSE is a valuable alternative to TD-DFT for properties related to the excited-state density and not only for transition energies and oscillator strengths.

Helical Organic and Inorganic Polymers.

Despite being a staple of synthetic plastics and biomolecules, helical polymers are scarcely studied with Gaussian-basis-set ab initio electron-correlated methods on an equal footing with molecules. This article introduces an ab initio second-order many-body Green's function [MBGF(2)] method with nondiagonal, frequency-dependent Dyson self-energy for infinite helical polymers using screw-axis-symmetry-adapted Gaussian-spherical-harmonics basis functions. Together with the Gaussian-basis-set density-functional theory for energies, analytical atomic forces, translational-period force, and helical-angle force, it can compute correlated energy, quasiparticle energy bands, structures, and vibrational frequencies of an infinite helical polymer, which smoothly converge at the corresponding oligomer results. These methods can handle incommensurable structures, which have an infinite translational period and are hard to characterize by any other method, just as efficiently as commensurable structures. We apply them to polyethylene (2/1 helix), polyacetylene (Peierls' system) and polytetrafluoroethylene (13/6 helix) to establish the quantitative accuracy of MBGF(2)/cc-pVDZ in simulating their (angle-resolved) ultraviolet photoelectron spectra and of B3LYP/cc-pVDZ or 6-31G** in reproducing their structures, infrared and Raman band positions, phonon dispersions, and (coherent and incoherent) inelastic neutron scattering spectra. We then predict the same properties for infinitely catenated chains of nitrogen or oxygen and discuss their possible metastable existence under ambient conditions. They include planar zigzag polyazene (N2)x (Peierls' system), 11/3-helical isotactic polyazane (NH)x, 9/4-helical isotactic polyfluoroazane (NF)x, and 7/2-helical polyoxane (O)x as potential high-energy-density materials.

Excess and excited-state dipole moments of real-life dyes: a comparison between wave-function, BSE/GW, and TD-DFT values.

In this work, we assess the accuracy of the Bethe-Salpeter equation (BSE) many-body Green's function formalism, adopting the eigenvalue-self-consistent evGW exchange-correlation kernel, for the calculation of the excited-state (μES) and excess dipole moments (Δμ), the latter ones being the changes of dipole amplitude between the ground and excited states (ES), in organic dyes. We compare the results obtained with wave-function methods [ADC(2), CC2, and CCSD], time-dependent density functional theory (TD-DFT), and BSE/evGW levels of theory. First, we compute the evolution of the dipole moments of the two lowest singlet excited states of 4-(dimethylamino)benzonitrile (DMABN) upon twisting of the amino group. Next, we use a set of 25 dyes having ES characters ranging from locally excited to charge transfer to determine both μES and Δμ. For DMABN our results show that BSE/evGW provides Δμ values closer to the CCSD reference and more consistent trends than TD-DFT. Moreover, a statistical analysis of both Δμ and μES for the set of 25 dyes shows that the BSE/evGW accuracy is comparable or sometimes slightly better than that of TD-M06-2X and TD-CAM-B3LYP, BSE/evGW outperforming TD-DFT in challenging cases (zwitterionic and cyanine transitions). Finally, the starting point dependency of BSE/evGW seems to be larger for Δμ, ES dipoles, and oscillator strengths than for transition energies.

Exact factorization of the many-body Green's function theory of electrons and nuclei

We combine the recently developed many-body Green's function theory for electrons and nuclei with the exact factorization of the wave function. The existing Born-Oppenheimer Green's functions are shown to be special cases of our exact approach. We consider the limitations of the laboratory frame formulation of the Green's function theory and discuss why the body-fixed frame formulation is needed in order to go beyond the Born-Oppenheimer theory. We give exact forms of the electronic and nuclear Green's functions written in terms of the exact factorized states, providing a systematic approach beyond the Born-Oppenheimer approximation. The lowest order approximation to the exact electronic Green's function is found to be an expected value of the Born-Oppenheimer electronic Green's function with respect to the nuclear density.

Training biases in machine learning for the analytic continuation of quantum many-body Green's functions

We address the problem of analytic continuation of imaginary-frequency Green's functions, which is crucial in many-body physics, using machine learning based on a multi-level residual neural network. We specifically address potential biases that can be introduced due to the use of artificially created spectral functions that are employed to train the neural network. We also implement an uncertainty estimation of the predicted spectral function, based on Monte Carlo dropout, which allows to identify frequency regions where the prediction might not be accurate, and we study the effect of noise, in particular also for situations where the noise level during training is different from that in the actual data. Our analysis demonstrates that this method can indeed achieve a high quality of prediction, comparable or better than the widely used maximum entropy method, but that further improvement is currently limited by the lack of true data that can be used for training. We also benchmark our approach by applying it to the case of SrVO$_3$, where an accurate spectral function has been obtained from dynamical mean-field theory using a solver that works directly on the real frequency axis.

One-particle Green's functions from the quantum equation of motion algorithm

Many-body Green's functions encode all the properties and excitations of interacting electrons. While these are challenging to be evaluated accurately on a classical computer, recent efforts have been directed towards finding quantum algorithms that may provide a quantum advantage for this task, exploiting architectures that will become available in the near future. In this work we introduce a novel near-term quantum algorithm for computing one-particle Green's functions via their Lehmann representation. The method is based on a generalization of the quantum equation of motion algorithm that gives access to the charged excitations of the system. We demonstrate the validity of the present proposal by computing the Green's function of a two-site Fermi-Hubbard model on a IBM quantum processor.

The charge-transfer states and excitation energy transfers of halogen-free organic molecules from first-principles many-body Green’s function theory

The organic solar cells based on halogen-free components, have been the new favorites to develop green and renewable energy. PBDB-T and its derivatives are considered the superior electron donors to construct the solar cells. Although there are plenty of researches about them, the charge-transfer mechanisms and excitation energy transfers of relative organic solar cells are still unclear, the developments of photovoltaic devices are restricted consequently. In this work, we calculate the electronic structures and excited-state properties of PBDB-T, PBT1-C, PBT1-O and PBT1-S donors in the gas phase from the many-body Green's function theory. With BTP-IC and BTP-IS as the acceptors, we consider the Förster, Dexter, and overlap electronic couplings to compute the excitation energy transfers of the dimers. The ionization energies and excited-state energies of the four donors calculated by GW + BSE are in good agreement with experiments, and they are sensitive to the functionals in the computation. We find two charge transfer schemes. The thienyl of PBDB-T molecule makes its charge-transfer state at the lowest energy, and the total electronic coupling of PBDB-T based dimer is the strongest. The Dexter, and overlap types electronic couplings are significant to study the excitation energy transfer of organic heterojunctions. We provide a theoretical guide in the design and synthesis of higher-performance halogen-free donors.

Identifying topological superconductivity in two-dimensional transition-metal dichalcogenides

We study the superconducting pairing instabilities and gap functions for prototypical two-dimensional (2D) transition-metal dichalcogenides (TMDCs) WS$_2$, MoTe$_2$, and MoS$_2$ in the 2H phase under both hole and electron doping at 10 K. Our first-principles quantum many-body Green's function approach allows us to treat the full $d$ and $p$ manifold of orbitals with strong spin-orbit coupling, yielding pairing predictions with material specific detail. The resulting gap functions exhibit a variety of mixed-parity superconducting states, including $s$, $p$, $d$, $f$, $d\pm id$, and $p\pm ip$ pairing modes. In particular, we predict 3% and 4% hole-doped WS$_2$ to be a chiral $p\pm ip$ topological superconductor. For 1% hole-doped MoS$_2$, we find a competition between three doubly degenerate chiral and non-chiral instabilities. Overall, the relative pairing strengths are found to follow the Fermi surface topology, due to nesting between the Fermi surface sheets. Finally, we discuss our predictions in relation to available experimental data and classify the topology of the predicted superconducting pairing symmetries.

Ana_cont: Python package for analytic continuation

We present the Python package ana_cont for the analytic continuation of fermionic and bosonic many-body Green's functions by means of either the Padé approximants or the maximum entropy method. The determination of hyperparameters and the implementation are described in detail. The code is publicly available on GitHub, where also documentation and learning resources are provided. Program summaryProgram Title: ana_contCPC Library link to program files:https://doi.org/10.17632/vphzykvvf2.1Developer's repository link:https://github.com/josefkaufmann/ana_contLicensing provisions: MITProgramming language:PythonExternal dependencies:Python(≥3.6), numpy, scipy, matplotlib, h5py, PyQt5, CythonSupplementary material: Test case files, tutorials, and instructionsNature of problem: Analytic continuation of correlation functions from Matsubara frequencies/imaginary time to real frequencies.Solution method: Padé interpolation, maximum entropy methodAdditional comments including restrictions and unusual features: The most important features can be accessed through the graphical user interface. For more flexibility, it is recommended to use the code as a library and write problem specific scripts.

Unified Approach to Polarons and Phonon-Induced Band Structure Renormalization.

Abinitio calculations of the phonon-induced band structure renormalization are currently based on the perturbative Allen-Heine theory and its many-body generalizations. These approaches are unsuitable to describe materials where electrons form localized polarons. Here, we develop a self-consistent, many-body Green's function theory of band structure renormalization that incorporates localization and self-trapping. We show that the present approach reduces to the Allen-Heine theory in the weak-coupling limit, and to total energy calculations of self-trapped polarons in the strong-coupling limit. To demonstrate this methodology, we reproduce the path-integral results of Feynman and diagrammatic MonteCarlo calculations for the Fröhlich model at all couplings, and we calculate the zero point renormalization of the band gap of an ionic insulator including polaronic effects.