Potential energy curve for isomerization of N2H2 and C2H4 using the improved virtual orbital multireference Møller–Plesset perturbation theory

Accurate Calculation of Interatomic Forces with Neural Networks Based on a Generative Transformer Architecture

Using neural networks to express electronic wave functions represents a new paradigm for solving the Schrödinger equation in quantum chemistry. For practical quantum chemistry simulations, one needs to know not only energies of molecules, but also accurate forces acting on constituent atoms. In this work, we achieve the accurate calculation of interatomic forces on QiankunNet, a platform that combines transformer-based deep neural networks with efficient batched autoregressive sampling. Our approach permits the application of the Hellmann-Feynman theorem to force calculations without introducing corrective Pulay terms. The results show that the calculated interatomic forces are in close agreement with those derived from the full configuration interaction method, irrespective of whether the system is a simple molecule or a strongly correlated electron system like a linear hydrogen chain. Furthermore, the calculated interatomic forces are employed for atomic relaxation in the torsional rotation process of ethylene, and the energy barrier obtained from the scanned potential energy surface is in excellent agreement with the experiment. Our work contributes to the application of artificial intelligence to broader quantum chemistry simulations, such as modeling challenging chemical transformations where electron correlations are difficult to describe.

Quantum Algorithm for Full Configuration Interaction Calculations without Controlled Time Evolutions

A quantum phase estimation algorithm allows us to perform full configuration interaction (full-CI) calculations on quantum computers with polynomial costs against the system size under study, but it requires quantum simulation of the time evolution of the wave function conditional on an ancillary qubit, which makes the algorithm implementation on real quantum devices difficult. Here, we discuss an application of the Bayesian phase difference estimation algorithm that is free from controlled time evolution operations to the full-CI calculations.

Single-Root Multireference Brillouin-Wigner Perturbative Approach to Excitation Energies

The state-specific Brillouin-Wigner multireference perturbation theory [which employs Jeziorski-Monkhorst parametrization of the wave function] using improved virtual orbitals, denoted as IVO-BWMRPT, is applied to calculate excitation energies (EEs) for methylene, ethylene, trimethylenemethane, and benzyne systems exhibiting various degrees of diradical character. In IVO-BWMRPT, all of the parameters appearing in the wave function ansatz are optimized for a specific electronic state. For these systems, the IVO-BWMRPT method provides EEs that are in close agreement with the benchmark results and experiments, where available, indicating that the method does not introduce imbalance in the target-specific treatment of closed- and open-shell states involved. The good performance of the present methodology is primarily related to structural compactness of the formalism. Overall, present findings are encouraging for both further development of the approach and chemical applications on the energy differences of strongly correlated systems.

Hartree-Fock on a superconducting qubit quantum computer

The simulation of fermionic systems is among the most anticipated applications of quantum computing. We performed several quantum simulations of chemistry with up to one dozen qubits, including modeling the isomerization mechanism of diazene. We also demonstrated error-mitigation strategies based on N-representability that dramatically improve the effective fidelity of our experiments. Our parameterized ansatz circuits realized the Givens rotation approach to noninteracting fermion evolution, which we variationally optimized to prepare the Hartree-Fock wave function. This ubiquitous algorithmic primitive is classically tractable to simulate yet still generates highly entangled states over the computational basis, which allowed us to assess the performance of our hardware and establish a foundation for scaling up correlated quantum chemistry simulations.

Theoretical investigation of the isomerization pathways of diazenes: torsion vs. inversion

Diazenes are an important family of organic compounds used widely in synthetic and materials chemistry. These molecules have a planar geometry and exhibit cis-trans isomerization. The simplest of all these molecules - diazene (N2H2) - has been subjected to several experimental and theoretical studies. Two mechanisms have been proposed for the cis-trans isomerization of diazene, which are an in-plane inversion and an out-of-plane torsion. The activation energies for these pathways are similar and the competition between these two mechanisms has been discussed in the literature based on electronic structure theory calculations. Three decades ago, a classical dynamics investigation of diazene isomerization was carried out using a model Hamiltonian and it was indicated that the in-plane inversion is forbidden classically because of a centrifugal barrier and the out-of-plane torsion is the only isomerization pathway. In the present work, we investigated the cis-trans isomerization dynamics of diazene using ab initio classical trajectory simulations at the CASSCF(2,2)/aug-cc-pVDZ level of electronic structure theory. The simulation results confirmed the presence of the aforementioned centrifugal barrier for the inversion and torsion was the only observed pathway. The calculations were repeated for a similar system (difluorodiazene, N2F2) and again the centrifugal barrier prevented the inversion pathway.

Ab Initio Probing of the Ground State of Tetraradicals: Breakdown of Hund's Multiplicity Rule

The electronic structure of organic σ-type polyradical including 2,4,6-tridehydropyridine radical cation (246-TDHP) and three isomers of tetradehydrobenzene (TDHB) have been studied using a computationally robust and cost-effective second-order multireference perturbative model which provides a balanced treatment of nondynamic and dynamic contributions to the electron correlation problem in the ground or excited electronic states which are imperative for predicting structural properties (e.g., ground state multiplicity, energy gaps between high-spin and low-spin states, etc.) of polyradicals. Energy gaps are useful to capture insight into the degree of interaction between the radical sites. An important finding of this study is that the tetraradicals considered here possess singlet ground states, contrary to Hund's rule. Present findings are in close agreement with the available high-level ab initio estimates at attainable cost implying that a perturbative description of the systems is adequate. The impact of N⁺ on the nature of ground state for the 246-TDHP have also analyzed. The singlet-triplet energy gaps for 1245- and 1234-TDHB are smaller than for o-benzyne mainly due to the ring strain. 1235-TDHB is 14.42 and 11.05 kcal/mol lower in energy than 1245- and 1234-isomers, respectively. IVO-SSMRPT predicts ¹A₁-³B₂ and ¹A₁-⁵B₂ gaps of 25.84 and 105.15 kcal/mol, respectively for the 246-TDHP cation.

The Hunt for Elusive Molecules: Insights from Joint Theoretical and Experimental Investigations

Rotational spectroscopy is an invaluable tool to unambiguously determine the molecular structure of a species, and sometimes even to establish its very existence. This article illustrates how experimental and theoretical state-of-the-art tools can be used in tandem to investigate the rotational structure of molecules, with particular emphasis on those that have long remained elusive. The examples of three emblematic species-gauche-butadiene, disilicon carbide, and germanium dicarbide-highlight the close, mutually beneficial interaction between high-level theoretical calculations and sensitive microwave measurements. Prospects to detect other elusive molecules of chemical and astronomical interest are discussed.

Bias-Free Chemically Diverse Test Sets from Machine Learning

Current benchmarking methods in quantum chemistry rely on databases that are built using a chemist's intuition. It is not fully understood how diverse or representative these databases truly are. Multivariate statistical techniques like archetypal analysis and K-means clustering have previously been used to summarize large sets of nanoparticles however molecules are more diverse and not as easily characterized by descriptors. In this work, we compare three sets of descriptors based on the one-, two-, and three-dimensional structure of a molecule. Using data from the NIST Computational Chemistry Comparison and Benchmark Database and machine learning techniques, we demonstrate the functional relationship between these structural descriptors and the electronic energy of molecules. Archetypes and prototypes found with topological or Coulomb matrix descriptors can be used to identify smaller, statistically significant test sets that better capture the diversity of chemical space. We apply this same method to find a diverse subset of organic molecules to demonstrate how the methods can easily be reapplied to individual research projects. Finally, we use our bias-free test sets to assess the performance of density functional theory and quantum Monte Carlo methods.

Four-Component Relativistic State-Specific Multireference Perturbation Theory with a Simplified Treatment of Static Correlation

The relativistic multireference (MR) perturbative approach is one of the most successful tools for the description of computationally demanding molecular systems of heavy elements. We present here the ground state dissociation energy surfaces, equilibrium bond lengths, harmonic frequencies, and dissociation energies of Ag₂, Cu₂, Au₂, and I₂ computed using the four-component (4c) relativistic spinors based state-specific MR perturbation theory (SSMRPT) with improved virtual orbital complete active space configuration interaction (IVO-CASCI) functions. The IVO-CASCI method is a simple, robust, useful and lower cost alternative to the complete active space self-consistent field approach for treating quasidegenerate situations. The redeeming features of the resulting method, termed as 4c-IVO-SSMRPT, lies in (i) manifestly size-extensivity, (ii) exemption from intruder problems, (iii) the freedom of convenient multipartitionings of the Hamiltonian, (iv) flexibility of the relaxed and unrelaxed descriptions of the reference coefficients, and (v) manageable cost/accuracy ratio. The present method delivers accurate descriptions of dissociation processes of heavy element systems. Close agreement with reference values has been found for the calculated molecular constants indicating that our 4c-IVOSSMRPT provides a robust and economic protocol for determining the structural properties for the ground state of heavy element molecules with eloquent MR character as it treats correlation and relativity on equal footing.

A simplified ab initio treatment of diradicaloid structures produced from stretching and breaking chemical bonds

With a proper choice of active spaces, the single root perturbation theory employing improved virtual orbitals can flawlessly describe the ground, excited, ionized, and dissociated states having varying degrees of degeneracy at the expense of low computational cost.

Cyclobutadiene automerization and rotation of ethylene: energetics of the barriers by using spin-polarized wave functions

Spin-projected spin polarized Møller-Plesset and spin polarized coupled clusters calculations have been made to estimate the cyclobutadiene automerization, the ethylene torsion barriers in their ground state, and the gap between the singlet and triplet states of ethylene. The results have been obtained optimizing the geometries at MP4 and/or CCSD levels, by an extensive Gaussian basis set. A comparative analysis with more complex calculations, up to MP5 and CCSDTQP, together with others from the literature, have also been made, showing the efficacy of using spin-polarized wave functions as a reference wave function for Møller-Plesset and coupled clusters calculations, in such problems.

Prediction of electronic structure of organic radicaloid anions using efficient, economical multireference gradient approach

The improved virtual orbital-complete active space configuration interaction (IVO-CASCI) method enables an economical and reasonably accurate treatment of static correlation in systems with significant multireference character, even when using a moderate basis set. This IVO-CASCI method supplants the computationally more demanding complete active space self-consistent field (CASSCF) method by producing comparable accuracy with diminished computational effort because the IVO-CASCI approach does not require additional iterations beyond an initial SCF calculation, nor does it encounter convergence difficulties or multiple solutions that may be found in CASSCF calculations. Our IVO-CASCI analytical gradient approach is applied to compute the equilibrium geometry for the ground and lowest excited state(s) of the theoretically very challenging 2,6-pyridyne, 1,2,3-tridehydrobenzene and 1,3,5-tridehydrobenzene anionic systems for which experiments are lacking, accurate quantum calculations are almost completely absent, and commonly used calculations based on single reference configurations fail to provide reasonable results. Hence, the computational complexity provides an excellent test for the efficacy of multireference methods. The present work clearly illustrates that the IVO-CASCI analytical gradient method provides a good description of the complicated electronic quasi-degeneracies during the geometry optimization process for the radicaloid anions. The IVO-CASCI treatment produces almost identical geometries as the CASSCF calculations (performed for this study) at a fraction of the computational labor. Adiabatic energy gaps to low lying excited states likewise emerge from the IVO-CASCI and CASSCF methods as very similar. We also provide harmonic vibrational frequencies to demonstrate the stability of the computed geometries.