An official website of the United States government
The accurate description of electron correlation and excitation energies remains a fundamental challenge in quantum chemistry. The particle–particle random phase approximation (ppRPA) has emerged as a promising method for capturing a broad range of excited-state properties. However, the implementation of ppRPA has been largely limited to in-house software, restricting its accessibility and usability. In this work, we present LibppRPA, an open-source and lightweight Python library designed for efficient and flexible ppRPA calculations of (1) electronic excitation energy and its associated analytical gradients and (2) the ground state correlation energy and its associated analytical gradients. LibppRPA enables seamless integration with existing quantum chemistry packages, such as PySCF, by utilizing occupation numbers, molecular orbital coefficients, and three-center electron repulsion integrals. We implement both direct diagonalization and the iterative Davidson algorithm for solving the ppRPA equations, as well as active-space approximations, allowing users to balance accuracy and computational efficiency. We demonstrate the performance of LibppRPA through benchmark calculations on singlet–triplet gaps, double excitations, charge-transfer excitations, and valence/Rydberg excitations, showcasing its reliability across diverse molecular systems. The library provides a robust platform for studying electronic excitations and offers new opportunities for future developments in electronic structure theory.
more »
« less
Accurate and efficient prediction of double excitation energies using the particle–particle random phase approximation
Double excitations are crucial to understanding numerous chemical, physical, and biological processes, but accurately predicting them remains a challenge. In this work, we explore the particle–particle random phase approximation (ppRPA) as an efficient and accurate approach for computing double excitation energies. We benchmark ppRPA using various exchange-correlation functionals for 21 molecular systems and two point defect systems. Our results show that ppRPA with functionals containing appropriate amounts of exact exchange provides accuracy comparable to high-level wave function methods such as CCSDT and CASPT2, with significantly reduced computational cost. Furthermore, we demonstrate the use of ppRPA starting from an excited (N − 2)-electron state calculated by ΔSCF for the first time, as well as its application to double excitations in bulk periodic systems. These findings suggest that ppRPA is a promising tool for the efficient calculation of double and partial double excitation energies in both molecular and bulk systems.
more »
« less
- Award ID(s):
- 2337991
- PAR ID:
- 10640333
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 162
- Issue:
- 9
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
It is demonstrated that a double hybrid density functional approximation, ωB88PTPSS, that incorporates equipartition of density functional theory and the non-local correlation, however with a meta-generalized gradient approximation correlation functional, as well as with the range-separated exchange of ωB2PLYP, provides accurate excitation energies for conventional systems, as well as correct prescription of negative singlet–triplet gaps for non-conventional systems with inverted gaps, without any necessity for parametric scaling of the same-spin and opposite-spin non-local correlation energies. Examined over “safe” excitations of the QUESTDB set, ωB88PTPSS performs quite well for open-shell systems, correctly and fairly accurately [relative to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) reference] predicts negative gaps for 50 systems with inverted singlet–triplet gaps, and is one of the leading performers for intramolecular charge-transfer excitations and achieves near-second-order approximate coupled cluster (CC2) and second-order algebraic diagrammatic construction quality for the Q1 and Q2 subsets. Subsequently, we tested ωB88PTPSS on two sets of real-life examples from recent computational chemistry literature–the low energy bands of chlorophyll a (Chl a) and a set of thermally activated delayed fluorescence (TADF) systems. For Chl a, ωB88PTPSS qualitatively and quantitatively achieves DLPNO-STEOM-CCSD-level performance and provides excellent agreement with experiment. For TADF systems, ωB88PTPSS agrees quite well with spin-component-scaled CC2 (SCS-CC2) excitation energies, as well as experimental values, for the gaps between the S1 and T1 excited states.more » « less
-
We applied localized orbital scaling correction (LOSC) in Bethe–Salpeter equation (BSE) to predict accurate excitation energies for molecules. LOSC systematically eliminates the delocalization error in the density functional approximation and is capable of approximating quasiparticle (QP) energies with accuracy similar to or better than GW Green’s function approach and with much less computational cost. The QP energies from LOSC, instead of commonly used G 0 W 0 and ev GW, are directly used in BSE. We show that the BSE/LOSC approach greatly outperforms the commonly used BSE/ G 0 W 0 approach for predicting excitations with different characters. For the calculations of Truhlar–Gagliardi test set containing valence, charge transfer, and Rydberg excitations, BSE/LOSC with the Tamm–Dancoff approximation provides a comparable accuracy to time-dependent density functional theory (TDDFT) and BSE/ev GW. For the calculations of Stein CT test set and Rydberg excitations of atoms, BSE/LOSC considerably outperforms both BSE/ G 0 W 0 and TDDFT approaches with a reduced starting point dependence. BSE/LOSC is, thus, a promising and efficient approach to calculate excitation energies for molecular systems.more » « less
-
Abstract It is commonly believed that there are only two types of particle exchange statistics in quantum mechanics, fermions and bosons, with the exception of anyons in two dimensions1–5. In principle, a second exception known as parastatistics, which extends outside two dimensions, has been considered6but was believed to be physically equivalent to fermions and bosons7–9. Here we show that non-trivial parastatistics inequivalent to either fermions or bosons can exist in physical systems. These new types of identical particle obey generalized exclusion principles, leading to exotic free-particle thermodynamics distinct from any system of free fermions and bosons. We formulate our theory by developing a second quantization of paraparticles that naturally includes exactly solvable non-interacting theories and incorporates physical constraints such as locality. We then construct a family of exactly solvable quantum spin models in one and two dimensions, in which free paraparticles emerge as quasiparticle excitations, and their exchange statistics can be physically observed and are notably distinct from fermions and bosons. This demonstrates the possibility of a new type of quasiparticle in condensed matter systems and—more speculatively—the potential for previously unconsidered types of elementary particle.more » « less
-
We present an implementation of relativistic double ionization potential (DIP) equation-of-motion coupled cluster (EOMCC) with up to 4-hole–2-particle (4h2p) excitations that makes use of the molecular mean-field exact two-component (mmfX2C) framework. We apply mmfX2C-DIP-EOMCC to several neutral atoms and diatomic molecules to obtain the ground and first few excited states of the corresponding dication species, and we observe excellent agreement (to within 0.001 eV) between DIPs obtained from mmfX2C- and four-component DIP-EOMCC calculations that include 3-hole–1-particle (3h1p) excitations, with either the Dirac–Coulomb or Dirac–Coulomb–Gaunt Hamiltonians. We also compare double IPs for mmfX2C-DIP-EOMCC calculations with the full Dirac–Coulomb–Breit Hamiltonian to those from experiment. The mmfX2C-DIP-EOMCC with 3h1p excitations leads to errors in absolute double IPs that generally overestimate experimental data for noble gases by 0.1–0.4 eV, whereas the inclusion of 4h2p excitations results in double IPs that are too low by 0.1–0.2 eV, at the large basis set limit.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1063/5.0251418
