1
|
Loos PF, Romaniello P. Static and dynamic Bethe-Salpeter equations in the T-matrix approximation. J Chem Phys 2022; 156:164101. [PMID: 35490009 DOI: 10.1063/5.0088364] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
While the well-established GW approximation corresponds to a resummation of the direct ring diagrams and is particularly well suited for weakly correlated systems, the T-matrix approximation does sum ladder diagrams up to infinity and is supposedly more appropriate in the presence of strong correlation. Here, we derive and implement, for the first time, the static and dynamic Bethe-Salpeter equations when one considers T-matrix quasiparticle energies and a T-matrix-based kernel. The performance of the static scheme and its perturbative dynamical correction are assessed by computing the neutral excited states of molecular systems. A comparison with more conventional schemes as well as other wave function methods is also reported. Our results suggest that the T-matrix-based formalism performs best in few-electron systems where the electron density remains low.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Pina Romaniello
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
2
|
Everhart LM, Derteano JA, Bates JE. Tension between predicting accurate ground state correlation energies and excitation energies from adiabatic approximations in TDDFT. J Chem Phys 2022; 156:084116. [DOI: 10.1063/5.0080382] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The connection between the adiabatic excitation energy of time-dependent density functional theory and the ground state correlation energy from the adiabatic connection fluctuation–dissipation theorem (ACFDT) is explored in the limiting case of one excited state. An exact expression is derived for any adiabatic Hartree-exchange–correlation kernel that connects the excitation energy and the potential contribution to correlation. The resulting formula is applied to the asymmetric Hubbard dimer, a system where this limit is exact. Results from a hierarchy of approximations to the kernel, including the random phase approximation (RPA) with and without exchange and the adiabatically exact (AE) approximation, are compared to the exact ones. At full coupling, the numerical results indicate a tension between predicting an accurate excitation energy and an accurate potential contribution to correlation. The AE approximation is capable of making accurate predictions of both quantities, but only in parts of the parameter space that classify as weakly correlated, while RPA tends to be unable to accurately predict these properties simultaneously anywhere. For a strongly correlated dimer, the AE approximation greatly overestimates the excitation energy yet continues to yield an accurate ground state correlation energy due to its accurate prediction of the adiabatic connection integrand. If similar trends hold for real systems, the development of correlation kernels will be important for applications of the ACFDT in systems with large potential contributions to correlation.
Collapse
Affiliation(s)
- Lucas M. Everhart
- Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28607, USA
| | - Julio A. Derteano
- Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28607, USA
| | - Jefferson E. Bates
- Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28607, USA
| |
Collapse
|
3
|
Di Sabatino S, Loos PF, Romaniello P. Scrutinizing GW-Based Methods Using the Hubbard Dimer. Front Chem 2021; 9:751054. [PMID: 34778206 PMCID: PMC8586429 DOI: 10.3389/fchem.2021.751054] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2021] [Accepted: 09/28/2021] [Indexed: 12/03/2022] Open
Abstract
Using the simple (symmetric) Hubbard dimer, we analyze some important features of the GW approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot GW method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard GW approximation and find, in particular, that 1) some neutral excitation energies become complex when the electron-electron interaction U increases, which can be traced back to the approximate nature of the GW quasiparticle energies; 2) the BSE formalism yields accurate correlation energies over a wide range of U when the trace (or plasmon) formula is employed; 3) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); 4) the trace formula has the correct behavior for weak (i.e., small U) interaction, unlike the ACFDT expression.
Collapse
Affiliation(s)
- S. Di Sabatino
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS and ETSF, Toulouse, France
| | - P.-F. Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| | - P. Romaniello
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS and ETSF, Toulouse, France
| |
Collapse
|
4
|
Marie A, Burton HGA, Loos PF. Perturbation theory in the complex plane: exceptional points and where to find them. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:283001. [PMID: 33601362 DOI: 10.1088/1361-648x/abe795] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2020] [Accepted: 02/18/2021] [Indexed: 05/24/2023]
Abstract
We explore the non-Hermitian extension of quantum chemistry in the complex plane and its link with perturbation theory. We observe that the physics of a quantum system is intimately connected to the position of complex-valued energy singularities, known as exceptional points. After presenting the fundamental concepts of non-Hermitian quantum chemistry in the complex plane, including the mean-field Hartree-Fock approximation and Rayleigh-Schrödinger perturbation theory, we provide a historical overview of the various research activities that have been performed on the physics of singularities. In particular, we highlight seminal work on the convergence behaviour of perturbative series obtained within Møller-Plesset perturbation theory, and its links with quantum phase transitions. We also discuss several resummation techniques (such as Padé and quadratic approximants) that can improve the overall accuracy of the Møller-Plesset perturbative series in both convergent and divergent cases. Each of these points is illustrated using the Hubbard dimer at half filling, which proves to be a versatile model for understanding the subtlety of analytically-continued perturbation theory in the complex plane.
Collapse
Affiliation(s)
- Antoine Marie
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France
| | - Hugh G A Burton
- Physical and Theoretical Chemical Laboratory, Department of Chemistry, University of Oxford, Oxford, OX1 3QZ, United Kingdom
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France
| |
Collapse
|
5
|
Monino E, Loos PF. Spin-Conserved and Spin-Flip Optical Excitations from the Bethe-Salpeter Equation Formalism. J Chem Theory Comput 2021; 17:2852-2867. [PMID: 33724811 PMCID: PMC8154368 DOI: 10.1021/acs.jctc.1c00074] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
![]()
Like adiabatic time-dependent
density-functional theory (TD-DFT),
the Bethe–Salpeter equation (BSE) formalism of many-body perturbation
theory, in its static approximation, is “blind” to double
(and higher) excitations, which are ubiquitous, for example, in conjugated
molecules like polyenes. Here, we apply the spin-flip ansatz (which considers the lowest triplet state as the reference configuration
instead of the singlet ground state) to the BSE formalism in order
to access, in particular, double excitations. The present scheme is
based on a spin-unrestricted version of the GW approximation
employed to compute the charged excitations and screened Coulomb potential
required for the BSE calculations. Dynamical corrections to the static
BSE optical excitations are taken into account via an unrestricted
generalization of our recently developed (renormalized) perturbative
treatment. The performance of the present spin-flip BSE formalism
is illustrated by computing excited-state energies of the beryllium
atom, the hydrogen molecule at various bond lengths, and cyclobutadiene
in its rectangular and square-planar geometries.
Collapse
Affiliation(s)
- Enzo Monino
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Route de Narbonne, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Route de Narbonne, 31062 Toulouse, France
| |
Collapse
|
6
|
Abstract
We discuss the physical properties and accuracy of three distinct dynamical (i.e., frequency-dependent) kernels for the computation of optical excitations within linear response theory: (i) an a priori built kernel inspired by the dressed time-dependent density-functional theory kernel proposed by Maitra et al. [J. Chem. Phys. 120, 5932 (2004)], (ii) the dynamical kernel stemming from the Bethe-Salpeter equation (BSE) formalism derived originally by Strinati [Riv. Nuovo Cimento 11, 1-86 (1988)], and (iii) the second-order BSE kernel derived by Zhang et al. [J. Chem. Phys. 139, 154109 (2013)]. The principal take-home message of the present paper is that dynamical kernels can provide, thanks to their frequency-dependent nature, additional excitations that can be associated with higher-order excitations (such as the infamous double excitations), an unappreciated feature of dynamical quantities. We also analyze, for each kernel, the appearance of spurious excitations originating from the approximate nature of the kernels, as first evidenced by Romaniello et al. [J. Chem. Phys. 130, 044108 (2009)]. Using a simple two-level model, prototypical examples of valence, charge-transfer, and Rydberg excited states are considered.
Collapse
Affiliation(s)
- Juliette Authier
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
7
|
Loos PF, Blase X. Dynamical correction to the Bethe-Salpeter equation beyond the plasmon-pole approximation. J Chem Phys 2020; 153:114120. [PMID: 32962392 DOI: 10.1063/5.0023168] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Bethe-Salpeter equation (BSE) formalism is a computationally affordable method for the calculation of accurate optical excitation energies in molecular systems. Similar to the ubiquitous adiabatic approximation of time-dependent density-functional theory, the static approximation, which substitutes a dynamical (i.e., frequency-dependent) kernel by its static limit, is usually enforced in most implementations of the BSE formalism. Here, going beyond the static approximation, we compute the dynamical correction of the electron-hole screening for molecular excitation energies, thanks to a renormalized first-order perturbative correction to the static BSE excitation energies. The present dynamical correction goes beyond the plasmon-pole approximation as the dynamical screening of the Coulomb interaction is computed exactly within the random-phase approximation. Our calculations are benchmarked against high-level (coupled-cluster) calculations, allowing one to assess the clear improvement brought by the dynamical correction for both singlet and triplet optical transitions.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Institut NEEL, F-38042 Grenoble, France
| |
Collapse
|
8
|
Blase X, Duchemin I, Jacquemin D, Loos PF. The Bethe-Salpeter Equation Formalism: From Physics to Chemistry. J Phys Chem Lett 2020; 11:7371-7382. [PMID: 32787315 DOI: 10.1021/acs.jpclett.0c01875] [Citation(s) in RCA: 70] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The Bethe-Salpeter equation (BSE) formalism is steadily asserting itself as a new efficient and accurate tool in the ensemble of computational methods available to chemists in order to predict optical excitations in molecular systems. In particular, the combination of the so-called GW approximation, giving access to reliable ionization energies and electron affinities, and the BSE formalism, able to model UV/vis spectra, has shown to provide accurate singlet excitation energies with a typical error of 0.1-0.3 eV. With a similar computational cost as time-dependent density-functional theory (TD-DFT), BSE is able to provide an accuracy on par with the most accurate global and range-separated hybrid functionals without the unsettling choice of the exchange-correlation functional, resolving further known issues (e.g., charge-transfer excitations). In this Perspective, we provide a historical overview of BSE, with a particular focus on its condensed-matter roots. We also propose a critical review of its strengths and weaknesses in different chemical situations.
Collapse
Affiliation(s)
- Xavier Blase
- Université Grenoble Alpes, CNRS, Institut NEEL, F-38042 Grenoble, France
| | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L Sim, 38054 Grenoble, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
9
|
Litvinova E, Schuck P, Wibowo H. Nuclear response at zero and finite temperature. EPJ WEB OF CONFERENCES 2019. [DOI: 10.1051/epjconf/201922301033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present some recent developments on the nuclear many-body problem, such as the treatment of high-order correlations and finite temperature in the description of in-medium two-nucleon propagators. In this work we discuss two-time propagators of the particle-hole type, which describe the response of finite nuclei to external probes without nucleon transfer. The general theory is formulated in terms of the equation of motion method for these propagators with the only input from the bare nucleon-nucleon interaction. The numerical implementation was performed on the basis of the effective mason-nucleon Lagrangian in order to study the energy-dependent kernels of different complexity. The finite-temperature extension of the theory with ph ⊗ phonon configurations is applied to a study of the multipole response of medium-mass nuclei.
Collapse
|