Search

We introduce \textsc{qcmath}, a user-friendly quantum chemistry software tailored for electronic structure calculations, implemented using the Wolfram Mathematica language and available at \url{https://github.com/LCPQ/qcmath}. This software, designed with accessibility in mind, takes advantage of the symbolic capabilities intrinsic to Mathematica. Its primary goal is to provide a supportive environment for newcomers to the field of quantum chemistry, enabling them to easily conceptualize, develop, and test their own ideas. The functionalities of \textsc{qcmath} encompass a broad spectrum of methods, catering to both ground- and excited-state calculations. We provide a comprehensive overview of these capabilities, complemented by essential theoretical insights. To facilitate ease of use, we offer an exhaustive blueprint of the software's architecture. Furthermore, we provide users with comprehensive guides, addressing both the operational aspects and the more intricate programming facets of \textsc{qcmath}., Comment: 9 pages

In recent years, Green's function methods have garnered considerable interest due to their ability to target both charged and neutral excitations. Among them, the well-established $GW$ approximation provides accurate ionization potentials and electron affinities and can be extended to neutral excitations using the Bethe-Salpeter equation (BSE) formalism. Here, we investigate the connections between various Green's function methods and evaluate their performance for charged and neutral excitations. Comparisons with other widely-known second-order wave function methods are also reported. Additionally, we calculate the singlet-triplet gap of cycl[3,3,3]azine, a model molecular emitter for thermally activated delayed fluorescence, which has the particularity of having an inverted gap thanks to a substantial contribution from the double excitations. We demonstrate that, within the $GW$ approximation, a second-order BSE kernel with dynamical correction is required to predict this distinctive characteristic., Comment: 18 pages, 1 figure

Here, we build on the works of Scuseria (et al.) http://dx.doi.org/10.1063/1.3043729 and Berkelbach https://doi.org/10.1063/1.5032314 to show connections between the Bethe-Salpeter equation (BSE) formalism combined with the $GW$ approximation from many-body perturbation theory and coupled-cluster (CC) theory at the ground- and excited-state levels. In particular, we show how to recast the $GW$ and Bethe-Salpeter equations as non-linear CC-like equations. Similitudes between BSE@$GW$ and the similarity-transformed equation-of-motion CC method introduced by Nooijen are also put forward. The present work allows to easily transfer key developments and general knowledge gathered in CC theory to many-body perturbation theory. In particular, it may provide a path for the computation of ground- and excited-state properties (such as nuclear gradients) within the $GW$ and BSE frameworks., Comment: 6 pages

Cyclobutadiene is a well-known playground for theoretical chemists and is particularly suitable to test ground- and excited-state methods. Indeed, due to its high spatial symmetry, especially at the $D_{4h}$ square geometry but also in the $D_{2h}$ rectangular arrangement, the ground and excited states of cyclobutadiene exhibit multi-configurational characters and single-reference methods, such as adiabatic time-dependent density-functional theory (TD-DFT) or equation-of-motion coupled cluster (EOM-CC), are notoriously known to struggle in such situations. In this work, using a large panel of methods and basis sets, we provide an extensive computational study of the automerization barrier (defined as the difference between the square and rectangular ground-state energies) and the vertical excitation energies at $D_{2h}$ and $D_{4h}$ equilibrium structures. In particular, selected configuration interaction (SCI), multi-reference perturbation theory (CASSCF, CASPT2, and NEVPT2), and coupled-cluster (CCSD, CC3, CCSDT, CC4, and CCSDTQ) calculations are performed. The spin-flip formalism, which is known to provide a qualitatively correct description of these diradical states, is also tested within TD-DFT (combined with numerous exchange-correlation functionals) and the algebraic diagrammatic construction [ADC(2)-s, ADC(2)-x, and ADC(3)] schemes. A theoretical best estimate is defined for the automerization barrier and for each vertical transition energy., Comment: 15 pages, 4 figures (supporting information available)

By recasting the non-linear frequency-dependent $GW$ quasiparticle equation into a linear eigenvalue problem, we explain the appearance of multiple solutions and unphysical discontinuities in various physical quantities computed within the $GW$ approximation. Considering the $GW$ self-energy as an effective Hamiltonian, it is shown that these issues are key signatures of strong correlation in the $(N\pm1)$-electron states and can be directly related to the intruder state problem. A simple and efficient regularization procedure inspired by the similarity renormalization group is proposed to avoid such issues and speed up convergence of partially self-consistent $GW$ calculations., Comment: 8 pages, 5 figures (supplementary material available)

In this contribution we start by proving and generalizing a conjecture that has been established few decades ago, relating the value of the integral of the detachment/attachment density in two pictures - one accounting for transition-induced basis relaxation and one which does not account for such a relaxation. To this end, we combine an extension to Cauchy's interlacing theorem and a brother theorem that we bring in the paper. Using this result together with a corollary to Lidskii-Wielandt theorem allows us to derive an upper bound for the electronic charge that is effectively displaced during the molecular electronic transition from one quantum state to another. This quantity can be regarded as the neat charge that has been transferred during the transition. Our derivations show that this boundary value can be determined from a simple singular value decomposition and two trace operations.

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 \textit{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., Comment: 14 pages, 3 figures and 3 tables