In this work, we develop two many-body quantum electrodynamic methods to calculate the ground-state energies of strongly coupled light-matter molecular systems. Specifically, we extend the second-order many-body Green's function theory (GF2) for electronic systems to incorporate electron-boson couplings. We employ two ansätze to treat the bosonic part of the system, namely the coherent-state (CS) and Lang-Firsov (LF) transformed vacuum state. These are combined with the GF2 method to construct two new approaches, which we refer to as CS-GF2 and LF-GF2. We benchmark CS- and LF-GF2 by studying various molecular systems inside an optical cavity. We investigate $\mathrm{H}_2$ and $\mathrm{LiH}$ potential energy surfaces, keto-eneol tautomerization energy barrier, van-der Waals interactions between two $\mathrm{H_2}$ molecules and the torsional potential energy surface of the ethylene molecule, $\mathrm{C_2H_4}$. Both methods provide highly accurate energies, with only modest additional improvement observed in LF-GF2.
In this work, we develop two many-body quantum electrodynamic methods to calculate the ground-state energies of strongly coupled light-matter molecular systems. Specifically, we extend the second-order many-body Green's function theory (GF2) for electronic systems to incorporate electron-boson couplings. We employ two ansätze to treat the bosonic part of the system, namely the coherent-state (CS) and Lang-Firsov (LF) transformed vacuum state. These are combined with the GF2 method to construct two new approaches, which we refer to as CS-GF2 and LF-GF2. We benchmark CS- and LF-GF2 by studying various molecular systems inside an optical cavity. We investigate H2 and LiH potential energy surfaces, keto-eneol tautomerization energy barrier, van-der Waals interactions between two H2 molecules and the torsional potential energy surface of the ethylene molecule, C2H4. Both methods provide highly accurate energies, with only modest additional improvement observed in LF-GF2.
We report an implementation of self-consistent Green's function many-body theory within a second-order approximation (GF2) for application with molecular systems. This is done by iterative solution of the Dyson equation expressed in matrix form in an atomic orbital basis, where the Green's function and self-energy are built on the imaginary frequency and imaginary time domain, respectively, and fast Fourier transform is used to efficiently transform these quantities as needed. We apply this method to several archetypical examples of strong correlation, such as a H32 finite lattice that displays a highly multireference electronic ground state even at equilibrium lattice spacing. In all cases, GF2 gives a physically meaningful description of the metal to insulator transition in these systems, without resorting to spin-symmetry breaking. Our results show that self-consistent Green's function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.
Recent progress in materials synthesis combined with advances in computational science have opened up an opportunity to design new materials with specific properties. Semiconductors, heterostructures, and thin layered films are used for energy harvesting and solar cells. Transition metal oxides and metallic surfaces with deposited magnetic molecules are used for signal conversion, nonvolatile memory modules and spintronics. Due to growing complexity of the new materials the most efficient approach to their design requires combined experimental and theoretical effort. In order to advance our understanding of the complicated physics of the modern materials new theoretical methods based on controlled, reliable, computationally efficient and systematically improvable approximations describing correlation and finite-temperature effects are needed. Density functional theory (DFT) remains a method of choice for materials science calculations due to its low computational cost. However, despite all recent efforts, the majority DFT applications are limited to weakly correlated systems and zero temperature. The lack of systematic improvability, inaccurate treatment of strongly correlated systems will continue to restrict the domain of DFT applications. The multi-reference wave-function methods can provide an accurate description of the strongly correlated systems. Unfortunately, these methods are computationally expensive. In spite of recent improvement the computational complexity will likely to hinder the applications of wave-function based methods to materials science in foreseeable future. Finite-temperature Green's function methods offer several advantages such as systematic improvability, access to the excitation spectrum, thermodynamic properties and straightforward implementation of embedding frameworks. It has been realized long time ago and has been successfully used in condensed matter physics to study model systems. The applications of such methods to realistic systems are impeded due to lack of theoretical approaches providing chemical accuracy and efficient computational algorithms. In this dissertation, these shortcomings are addressed by developing self-energy embedding theory and its extensive comparison to established quantum chemistry methods. New numerical algorithms for efficient application of the finite-temperature Green's function methods to molecular and extended systems were developed. Properties of the second-order Green's function theory (GF2) in the context of quantum chemistry were investigated. Additionally, a new range-separated hybrid functional combining GF2 and DFT was developed and benchmarked.
Achieving materials-specific predictions for large, realistic molecules and solids with strong electron correlations remains a long-standing challenge in quantum chemistry. This dissertation develops scalable ab initio quantum many-body methods for crystalline materials, demonstrating that rigorous optimization of theoretical approximations and numerical algorithms enables predictive computation even in strongly correlated systems. The work tackles three interconnected challenges, namely describing material-specific properties of strongly correlated superconductors, achieving linear scaling with k-points for large-scale periodic calculations, and establishing reliable methods for electron-phonon interactions. The following chapters address each challenge in turn. Chapter 2 demonstrates that ab initio density matrix embedding theory, combined with symmetry-breaking strategies for superconducting order parameters, reproduces experimental trends in cuprate superconductivity, including pressure and layer effects on transition temperatures. By directly solving the electronic Schrödinger equation from material structures, these calculations capture multi-orbital covalency effects that simplified model Hamiltonians neglect, revealing their essential role in the spin fluctuations driving pairing. This predictive capability establishes a foundation for efficient computational screening of superconducting materials. Extending these calculations to larger k-point meshes and broader phase diagram explorations, however, demands more efficient computational algorithms. To overcome this computational bottleneck, Chapter 3 introduces interpolative separable density fitting framework. This development enables thermodynamic limit calculations with up to 1000 k-points, as validated on diamond, carbon dioxide, nickel monoxide, and cuprate systems. Integration with density matrix embedding theory and local natural orbital methods further yields converged ground-state energies from accurate correlated wavefunction methods such as CCSD(T) for crystalline materials in the thermodynamic limit. With computational efficiency established, the remaining challenge lies in extending the theoretical framework beyond purely electronic degrees of freedom. Chapter 4 addresses this need by establishing exponential ansatz wavefunctions as competitive approaches for electron-phonon systems. Unlike previous treatments based on perturbation theory, systematic benchmarking of coupled cluster theory and variational Lang-Firsov methods on the Holstein model demonstrates, for the first time, that these approaches accurately describe polaron formation across coupling regimes when combined with appropriate reference state optimization. These results provide a foundation for incorporating phonon effects into ab initio studies of superconducting materials. In addition to these methodological advances, GPU acceleration and MPI parallelization enhance computational throughput, enabling efficient large-scale calculations. Together, these developments constitute a computational framework for predictive ab initio quantum many-body calculations in crystalline materials, opening new avenues for understanding and designing strongly correlated quantum materials.
We investigate analytically the performance of many-body energy functionals, derived by Klein, Luttinger, and Ward, at different levels of diagrammatic approximations, ranging from second Born, to , to the so-called T matrix, for the calculation of total energies and potential energy surfaces. We benchmark our theoretical results on the extended two-site Hubbard model, which is analytically solvable and for which several exact properties can be calculated. Despite its simplicity, this model displays the physics of strongly correlated electrons: it is prototypical of the H 2 dissociation, a notoriously difficult problem to solve accurately for the majority of mean-field-based approaches. We show that both functionals exhibit good to excellent variational properties, particularly in the case of the Luttinger-Ward one, which is in close agreement with fully self-consistent calculations, and we elucidate the relation between the accuracy of the results and the different input one-body Green's functions. Provided that these are wisely chosen, we show how the Luttinger-Ward functional can be used as a computationally inexpensive alternative to fully self-consistent many-body calculations, without sacrificing the precision of the results obtained. Furthermore, in virtue of this accuracy, we argue that this functional can also be used to rank different many-body approximations at different regimes of electronic correlation, once again bypassing the need for self-consistency. Published by the American Physical Society 2024
We present a scalable single-particle framework to treat electronic correlation in molecules and materials motivated by Green's function theory. We derive a size-extensive Brillouin-Wigner perturbation theory from the single-particle Green's function by introducing the Goldstone self-energy. This new ground state correlation energy, referred to as Quasi-Particle MP2 theory (QPMP2), avoids the characteristic divergences present in both second-order Møller-Plesset perturbation theory and Coupled Cluster Singles and Doubles within the strongly correlated regime. We show that the exact ground state energy and properties of the Hubbard dimer are reproduced by QPMP2 and demonstrate the advantages of the approach for larger Hubbard models where the metal-to-insulator transition is qualitatively reproduced, contrasting with the complete failure of traditional methods. We apply this formalism to characteristic strongly correlated molecular systems and show that QPMP2 provides an efficient, size-consistent regularization of MP2.
The coupled coherent states method has demonstrated itself as an accurate and efficient method of studying the quantum dynamics of various systems. In recent years, its applicability has been extended by incorporating a number of new numerical expansions and modifications to generate a closely related family of methods. In this thesis, two new augmentations are developed to further broaden the scope of problems that are able to be treated. The first of these is a 2-layer extension of coupled coherent states, capable of providing an increased mathematical description of a degree or degrees of freedom within a quantum mechanical system, as well as beneficial numerical and scalability properties. The newly developed method is tested on a model system-bath Hamiltonian consisting of a tunnelling mode governed by an asymmetric double well potential coupled to a harmonic bath. It is found to compare well to previous methods of studying the Hamiltonian, as well as a benchmark calculation on the system conducted in this thesis, and demonstrate the beneficial numerical and scalability properties expected. The second development is to extend coupled coherent states to treat systems of indistinguishable bosons in the second quantisation representation. The method is tested on the same Hamiltonian as the 2-layer coupled coherent states scheme, where the harmonic bath is second quantised as it is comprised of oscillators of the same frequency, so they may be thought of as indistinguishable. Exploiting this symmetry property is found to be extremely advantageous, with remarkable agreement to the benchmark calculation. The method is then tested on a model Hamiltonian consisting of 100 bosons in a shifted harmonic trap, with oscillations in the 1-body density calculated. The results are found to compare favourably with a multiconfigurational time-dependent Hartree for bosons calculation that is equivalent to the Gross-Pitaevskii equation, providing impetus for future studies on systems of Bose-Einstein condensates. The existing ab initio multiple cloning extension of coupled coherent states for nonadiabatic dynamics is also used to study the ultrafast photodissociation of 2-ethylpyrrole. The results are compared to experimental data, and a novel insight into the dissociation mechanism is obtained, with it shown to be composed of a two step process. Firstly, molecules that are able to dissociate immediately over the barrier along the N-H coordinate do so in < 50 fs, and this is followed by a second slower dissociation process from molecules that must sample the potential energy surface before finding a way around the barrier. This is not observed experimentally due to the temporal widths of the laser pulses obscuring the dynamics in the < 50 fs window.