Ground- and Excited-State Geometry Optimization of Small Organic Molecules with Quantum Monte Carlo

Optimization Stability in Excited-State-Specific Variational Monte Carlo

We investigate the issue of optimization stability in variance-based state-specific variational Monte Carlo, discussing the roles of the objective function, the complexity of wave function ansatz, the amount of sampling effort, and the choice of minimization algorithm. Using a small cyanine dye molecule as a test case, we systematically perform minimizations using variants of the linear method as both a standalone algorithm and in a hybrid combination with accelerated descent. We demonstrate that adaptive step control is crucial for maintaining the linear method's stability when optimizing complicated wave functions and that the hybrid method enjoys both greater stability and minimization performance.

A structural optimization algorithm with stochastic forces and stresses

We propose an algorithm for optimizations in which the gradients contain stochastic noise. This arises, for example, in structural optimizations when computations of forces and stresses rely on methods involving Monte Carlo sampling, such as quantum Monte Carlo or neural network states, or are performed on quantum devices that have intrinsic noise. Our proposed algorithm is based on the combination of two ingredients: an update rule derived from the steepest-descent method, and a staged scheduling of the targeted statistical error and step size, with position averaging. We compare it with commonly applied algorithms, including some of the latest machine learning optimization methods, and show that the algorithm consistently performs efficiently and robustly under realistic conditions. Applying this algorithm, we achieve full-degree optimizations in solids using ab initio many-body computations, by auxiliary-field quantum Monte Carlo with plane waves and pseudopotentials. A potential metastable structure in Si is discovered using density-functional calculations with synthetic noisy forces.

Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry

The accurate prediction of reaction mechanisms in heterogeneous (surface) catalysis is one of the central challenges in computational chemistry. Quantum Monte Carlo methods─Diffusion Monte Carlo (DMC) in particular─are being recognized as higher-accuracy, albeit more computationally expensive, alternatives to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the broader adoption of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that it is possible to significantly reduce this computational cost by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like adsorption (binding) energies and mapping potential energy surfaces. We analyze the cancellation and convergence of many-body finite-size errors in two well-known adsorbate/slab systems, H₂O/LiH(001) and CO/Pt(111). Based on this analysis, we identify strategies for obtaining binding energies in the thermodynamic limit that optimally utilize error cancellation to balance accuracy and computational efficiency. Using one such strategy, we then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for a wide range of density functionals. Our accurate and inexpensive DMC calculations are found to unambiguously recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to use this DMC method to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.

Near-Exact Nuclear Gradients of Complete Active Space Self-Consistent Field Wave Functions

In this paper, we study the nuclear gradients of heat bath configuration interaction self-consistent field (HCISCF) wave functions and use them to optimize molecular geometries for various molecules.We show that the HCISCF nuclear gradients are fairly insensitive to the size of the "selected" variational space, which allows us to reduce the computational cost without introducing significant error.The ability of HCISCF to treat larger active spaces combined with the flexibility for users to control the computational cost makes the method very attractive for studying strongly correlated systems which require a larger active space than possible with complete active space self-consistent field (CASSCF).Finally, we study the realistic catalyst, Fe(PDI), and highlight some of the challenges this system poses for density functional theory (DFT).We demonstrate how HCISCF can clarify the energetic stability of geometries obtained from DFT when the results are strongly dependent on the functional.We also use the HCISCF gradients to optimize geometries for this species and study the adiabatic singlet-triplet gap. During geometry optimization, we find that multiple near-degenerate local minima exist on the triplet potential energy surface.

Excited-State Geometry Optimization of Small Molecules with Many-Body Green's Functions Theory

We present a benchmark study of gas phase geometry optimizations in the excited states of carbon monoxide, acetone, acrolein, and methylenecyclopropene using many-body Green’s functions theory within the GW approximation and the Bethe–Salpeter equation (BSE) employing numerical gradients. We scrutinize the influence of several typical approximations in the GW-BSE framework; we used one-shot G₀W₀ or eigenvalue self-consistent evGW, employing a fully analytic approach or plasmon-pole model for the frequency dependence of the electron self-energy, or performing the BSE step within the Tamm–Dancoff approximation. The obtained geometries are compared to reference results from multireference perturbation theory (CASPT2), variational Monte Carlo (VMC) method, second-order approximate coupled cluster (CC2) method, and time-dependent density-functional theory (TDDFT). We find overall a good agreement of the structural parameters optimized with the GW-BSE calculations with CASPT2, with an average relative error of around 1% for the G₀W₀ and 1.5% for the evGW variants based on a PBE0 ground state, respectively, while the other approximations have negligible influence. The relative errors are also smaller than those for CC2 and TDDFT with different functionals and only larger than VMC, indicating that the GW-BSE method does not only yield excitation energies but also geometries in good agreement with established higher-order wave function methods.

Accurate spectroscopic properties by diffusion quantum Monte Carlo calculations

The capability of Diffusion Quantum Monte Carlo (DMC) to produce high quality potential energy curve (PEC) was evaluated. H₂⁺, HeH⁺ and LiH PECs were built by all-electron fixed-node DMC calculations. Trial wave functions were obtained from Hartree-Fock (HF) (H₂⁺), MCSCF and CI (HeH⁺ and LiH) calculations multiplied by Jastrow factor. The quality of these generated PECs was analyzed throughout equilibrium distance, dissociation energy, vibrational energies and rovibrational spectroscopic constants (ω_e, ω_ex_e, ω_ey_e, α_e, γ_e and B_e). The Discrete Variable Representation (DVR) and the Dunham approaches were used to determine the rovibrational spectroscopic constants. The PECs and the aforementioned properties were also obtained by the following methods: MCSCF/aug-cc-pV5Z (LiH), CCSD(T)/aug-cc-pV5Z (HeH⁺ and LiH) and HF (H₂⁺ and HeH⁺) levels. The results of these DMC computations, specially the DMC-DVR procedure, are the most accurate among others DMC calculations available in the literature for these systems. They suggest that DMC can be used to achieve accurate PECs to produce spectroscopic properties with the same level of accuracy of theoretical benchmarks and experimental data of the literature.

Rational design of photo-cleavable alkoxyamines for polymerization and synthesis

Theoretical calculations have been performed in order to investigate the impact of different substitution patterns on predicted photoreactivity of alkoxyamines fused to an anthraquinone chromophore. Amino and hydroxy groups (similar to those which have been previously synthesized) are introduced and their effect on excited state energies and charge transfer is assessed. Analogous to formally oxidized alkoxyamines, the charge-separated nNπ* state can undergo mesolytic cleavage or bimolecular or SN2 reactions with nucleophiles, according to the substitution patterns and other reagents present. While homolytic cleavage is in principle promoted by triplet ππ* states, the accessible ππ* triplet states in this system are centered on the chromophore and unreactive. We show that the reactive nNπ* state, which bears a negative charge, is stabilized by hydroxy substitution while amino substitution will destabilize it. After mesolysis to a carbon centred radical, the nitroxide radical re-forms; however, when carbocations are produced the remaining open-shell singlet is stable and unable to undergo coupling with the carbocation.

The Quest for Highly Accurate Excitation Energies: A Computational Perspective

We provide an overview of the successive steps that made it possible to obtain increasingly accurate excitation energies with computational chemistry tools, eventually leading to chemically accurate vertical transition energies for small- and medium-size molecules. First, we describe the evolution of ab initio methods employed to define benchmark values, with the original Roos CASPT2 method, then the CC3 method as in the renowned Thiel set, and more recently the resurgence of selected configuration interaction methods. The latter method has been able to deliver consistently, for both single and double excitations, highly accurate excitation energies for small molecules, as well as medium-size molecules with compact basis sets. Second, we describe how these high-level methods and the creation of representative benchmark sets of excitation energies have allowed the fair and accurate assessment of the performance of computationally lighter methods. We conclude by discussing possible future theoretical and technological developments in the field.

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules

Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.

Singlet-triplet gaps in diradicals obtained with diffusion quantum Monte Carlo using a Slater-Jastrow trial wavefunction with a minimum number of determinants

Diradicals are essential species in a wide range of chemical processes, whereas the computational study of their electronic structure often remains a challenge due to near-degeneracy of the frontier molecular orbitals. The fixed-node diffusion quantum Monte Carlo (FN-DMC) method is employed to calculate adiabatic energy gaps of some typical diradicals with the Slater-Jastrow trial wavefunction. The antisymmetrized part of the trial wavefunction is taken to be a linear combination of a minimum number of determinants using RB3LYP orbitals from the closed-shell singlet state or ROB3LYP orbitals from the triplet state. Our results show that using the two-determinant-Jastrow trial wavefunction is necessary to achieve reliable energy differences between closed-shell singlet states. The energy of the triplet state with MS = 1 is calculated to be lower than that with MS = 0 with FN-DMC even using trial wavefunctions with spin-pure states as their antisymmetrized parts and this difference is reduced with better orbitals. This indicates that the fixed-node error is smaller for the triplet state with MS = 1. Adiabatic energy gaps obtained from the present FN-DMC calculations are in reasonable agreement with available experimental values. Compared with results of the high level EOM-SF-CC method, energy gaps of FN-DMC with RB3LYP orbitals are slightly better than those using ROB3LYP orbitals and results of EOM-SF-CCSD. The present FN-DMC calculations using the simplest ansatz for the trial wavefunction can achieve reasonable results for these diradicals and they can readily be applied to large diradicals.

Excited States with Selected Configuration Interaction-Quantum Monte Carlo: Chemically Accurate Excitation Energies and Geometries

We employ quantum Monte Carlo to obtain chemically accurate vertical and adiabatic excitation energies, and equilibrium excited-state structures for the small, yet challenging, formaldehyde and thioformaldehyde molecules. A key ingredient is a robust protocol to obtain balanced ground- and excited-state Jastrow-Slater wave functions at a given geometry, and to maintain such a balanced description as we relax the structure in the excited state. We use determinantal components generated via a selected configuration interaction scheme which targets the same second-order perturbation energy correction for all states of interest at different geometries, and fully optimize all variational parameters in the resultant Jastrow-Slater wave functions. Importantly, the excitation energies as well as the structural parameters in the ground and excited states are converged with very compact wave functions comprising few thousand determinants in a minimally augmented double-ζ basis set. These results are obtained already at the variational Monte Carlo level, the more accurate diffusion Monte Carlo method yielding only a small improvement in the adiabatic excitation energies. We find that matching Jastrow-Slater wave functions with similar variances can yield excitation energies compatible with our best estimates; however, the variance-matching procedure requires somewhat larger determinantal expansions to achieve the same accuracy, and it is less straightforward to adapt during structural optimization in the excited state.

Accurate Excited-State Geometries: A CASPT2 and Coupled-Cluster Reference Database for Small Molecules

We present an investigation of the excited-state structural parameters determined for a large set of small compounds with the dual goals of defining reference values for further works and assessing the quality of the geometries obtained with relatively cheap computational approaches. In the first stage, we compare the excited-state geometries obtained with ADC(2), CC2, CCSD, CCSDR(3), CC3, and CASPT2 and large atomic basis sets. It is found that CASPT2 and CC3 results are generally in very good agreement with one another (typical differences of ca. 3 × 10^–3 Å) when all electrons are correlated and when the aug-cc-pVTZ atomic basis set is employed with both methods. In a second stage, a statistical analysis reveals that, on the one hand, the excited-state (ES) bond lengths are much more sensitive to the selected level of theory than their ground-state (GS) counterparts and, on the other hand, that CCSDR(3) is probably the most cost-effective method delivering accurate structures. Indeed, CCSD tends to provide too compact multiple bond lengths on an almost systematic basis, whereas both CC2 and ADC(2) tend to exaggerate these bond distances, with more erratic error patterns, especially for the latter method. The deviations are particularly marked for the polarized CO and CN bonds, as well as for the puckering angle in formaldehyde homologues. In the last part of this contribution, we provide a series of CCSDR(3) GS and ES geometries of medium-sized molecules to be used as references in further investigations.

Quantum Monte Carlo Calculations on a Benchmark Molecule-Metal Surface Reaction: H2 + Cu(111)

Accurate modeling of heterogeneous catalysis requires the availability of highly accurate potential energy surfaces. Within density functional theory, these can—unfortunately—depend heavily on the exchange-correlation functional. High-level ab initio calculations, on the other hand, are challenging due to the system size and the metallic character of the metal slab. Here, we present a quantum Monte Carlo (QMC) study for the benchmark system H₂ + Cu(111), focusing on the dissociative chemisorption barrier height. These computationally extremely challenging ab initio calculations agree to within 1.6 ± 1.0 kcal/mol with a chemically accurate semiempirical value. Remaining errors, such as time-step errors and locality errors, are analyzed in detail in order to assess the reliability of the results. The benchmark studies presented here are at the cutting edge of what is computationally feasible at the present time. Illustrating not only the achievable accuracy but also the challenges arising within QMC in such a calculation, our study presents a clear picture of where we stand at the moment and which approaches might allow for even more accurate results in the future.

Molecular structure and spectroscopic signatures of acrolein: theory meets experiment

A comprehensive study of the molecular structure and IR spectrum of cis and trans acrolein has been performed by an integrated computational approach coupling methods rooted in the coupled-cluster ansatz and the density functional theory. From the one side, DFT anharmonic force fields allow us to determine very reliable semiexperimental structures for both isomers, which are in remarkable agreement with the geometries issuing from CCSD(T) computations accounting for the extrapolation to the complete basis set and core correlation. The same kind of coupled-cluster computations provide dipole moment, relative energies, and interconversion barrier in remarkable agreement with experiments. Finally, harmonic CCSD(T) results coupled to DFT evaluation of mechanical and electrical anharmonicity allow us, in the framework of second-order perturbative vibrational theory, to confirm most of the experimental assignments of IR spectra, and to suggest some additional interpretations for congested regions including fundamental bands together with overtones and combination bands.