Showing 1–50 of 78 results for author: Head-Gordon, M

The Role of Spin-Orbit Coupling on the Linear Absorption Spectrum and Intersystem Crossing Rate Coefficients of Ruthenium Polypyridyl Dyes

Authors: Justin J. Talbot, Thomas P. Cheshire, Stephen J. Cotton, Frances A. Houle, Martin Head-Gordon

Abstract: The successful use of molecular dyes for solar energy conversion requires efficient charge injection, which in turn requires the formation of states with sufficiently long lifetimes (e.g. triplets). The molecular structure elements that confer this property can be found empirically, however computational predictions using $\textit{ab initio}$ electronic structure methods are invaluable to identify… ▽ More The successful use of molecular dyes for solar energy conversion requires efficient charge injection, which in turn requires the formation of states with sufficiently long lifetimes (e.g. triplets). The molecular structure elements that confer this property can be found empirically, however computational predictions using $\textit{ab initio}$ electronic structure methods are invaluable to identify structure-property relations for dye sensitizers. The primary challenge for simulations to elucidate the electronic and nuclear origins of these properties is a spin-orbit interaction which drives transitions between electronic states. In this work, we present a computational analysis of the spin-orbit corrected linear absorption cross sections and intersystem crossing rate coefficients for a derivative set of phosphonated tris(2,2'-bipyridine)ruthenium(2+) dye molecules. After sampling the ground state vibrational distributions, the predicted linear absorption cross sections indicate that the mixture between singlet and triplet states plays a crucial role in defining the line shape of the metal-to-ligand charge transfer bands in these derivatives. Additionally, an analysis of the intersystem crossing rate coefficients suggests that transitions from the singlet into the triplet manifolds are ultrafast with rate coefficients on the order of $10^{13}$ s$^{-1}$ for each dye molecule. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 38 pages, 7 figures

Spin parameter optimization for spin-polarized extended tight-binding methods

Authors: Siyavash Moradi, Rebecca Tomann, Josie Hendrix, Martin Head-Gordon, Christopher J. Stein

Abstract: We present an optimization strategy for atom-specific spin-polarization constants within the spin-polarized GFN2-xTB framework, aiming to enhance the accuracy of molecular simulations. We compare a sequential and global optimization of spin parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. Sensitivity analysis using Sobol indices guides the identification of the most influential par… ▽ More We present an optimization strategy for atom-specific spin-polarization constants within the spin-polarized GFN2-xTB framework, aiming to enhance the accuracy of molecular simulations. We compare a sequential and global optimization of spin parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. Sensitivity analysis using Sobol indices guides the identification of the most influential parameters for a given reference dataset, allowing for a nuanced understanding of their impact on diverse molecular properties. In the case of the W4-11 dataset, substantial error reduction was achieved, demonstrating the potential of the optimization. Transferability of the optimized spin-polarization constants over different properties, however, is limited, as we demonstrate by applying the optimized parameters on a set of singlet-triplet gaps in carbenes. Further studies on ionization potentials and electron affinities highlight some inherent limitations of current extended tight-binding methods that can not be resolved by simple parameter optimization. We conclude that the significantly improved accuracy strongly encourages the present re-optimization of the spin-polarization constants, whereas the limited transferability motivates a property-specific optimization strategy. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Measurement of Coherent Vibrational Dynamics with X-ray Transient Absorption Spectroscopy Simultaneously at the Carbon K- and Chlorine L$_{2,3}$- Edges

Authors: Andrew D. Ross, Diptarka Hait, Valeriu Scutelnic, Daniel M. Neumark, Martin Head-Gordon, Stephen R. Leone

Abstract: X-ray Transient Absorption Spectroscopy near the carbon K-edge (1s, $\sim$ 285 eV) and chlorine L$_{2,3}$ edges (2p, $\sim$ 200 eV) is used to study the nuclear dynamics of CCl$_4$ vibrationally activated by impulsive stimulated Raman scattering with a few-cycle 800 nm pump pulse. The totally symmetric stretching mode leads to a strong response in the inner-shell spectra, with the concerted elonga… ▽ More X-ray Transient Absorption Spectroscopy near the carbon K-edge (1s, $\sim$ 285 eV) and chlorine L$_{2,3}$ edges (2p, $\sim$ 200 eV) is used to study the nuclear dynamics of CCl$_4$ vibrationally activated by impulsive stimulated Raman scattering with a few-cycle 800 nm pump pulse. The totally symmetric stretching mode leads to a strong response in the inner-shell spectra, with the concerted elongation (contraction) in bond lengths leading to a red (blue) shift in the X-ray absorption energies associated with core-to-antibonding excitations. The relative slopes of the potential energy surfaces associated with the relevant core-excited states along the symmetric stretching mode are experimentally measured and compared to results from restricted open-shell Kohn-Sham calculations. A combination of experiment and theory indicates that the slope of the core-excited potential energy surface vs totally symmetric bond elongation is $-11.1 \pm 0.8$ eV/Å for the Cl 2p$\to7a_1^*$ excitation, $-9.0\pm0.6$ eV/Å for the Cl 2p$\to8t_2^*$ excitation and $-5.2\pm 0.4$ eV/Å for the C 1s$\to8t_2^*$ excitation, to 95% confidence. The much larger slopes for the Cl 2p excitations compared to the C 1s state are attributed to greater contributions from Cl to the $7a_1^*$ or $8t_2^*$ antibonding orbitals to which the inner-shell electrons are being excited. No net displacement of the center of the vibrational wavefunction along the other vibrational modes is induced by the pump pulse, leading to absence of transient signal. The results highlight the ability of X-ray Transient Absorption Spectroscopy to reveal nuclear dynamics involving tiny ($<0.01$ Å) atomic displacements and also provide direct measurement of forces on core-excited potential energy surfaces. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: ADR and DH contributed equally to this work

A static quantum embedding scheme based on coupled cluster theory

Authors: Avijit Shee, Fabian M. Faulstich, Birgitta Whaley, Lin Lin, Martin Head-Gordon

Abstract: We develop a static quantum embedding scheme, utilizing projection equations to solve coupled cluster (CC) amplitudes. To reduce the computational cost (for example, of a large basis set calculation), we solve the local fragment problem using a high-level coupled cluster method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach i… ▽ More We develop a static quantum embedding scheme, utilizing projection equations to solve coupled cluster (CC) amplitudes. To reduce the computational cost (for example, of a large basis set calculation), we solve the local fragment problem using a high-level coupled cluster method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach is consistently formulated within the coupled cluster framework and will be called MP-CC. We demonstrate the effectiveness of our method through several prototypical molecular examples by analyzing a global quantity, that is, the total correlation energy of the system in the study of potential energy curves (PEC) and thermochemical reaction energies. We have shown that our method can achieve comparable accuracy both with a small and large basis set when the fragment Hilbert space size remains the same. Additionally, our results indicate that increasing the fragment size can systematically enhance the accuracy of observables, approaching the precision of the full coupled cluster solver. △ Less

Submitted 13 April, 2024; originally announced April 2024.

The Staggered Mesh Method: Accurate Exact Exchange towards the Thermodynamic Limit for Solids

Authors: Stephen Jon Quiton, Hamlin Wu, Xin Xing, Lin Lin, Martin Head-Gordon

Abstract: In periodic systems, the Hartree-Fock (HF) exchange energy exhibits the slowest convergence of all HF energy components as the system size approaches the thermodynamic limit. We demonstrate that the recently proposed staggered mesh method for Fock exchange energy [Xing, Li, and Lin, Math. Comp., 2024], which is specifically designed to sidestep certain singularities in exchange energy evaluation,… ▽ More In periodic systems, the Hartree-Fock (HF) exchange energy exhibits the slowest convergence of all HF energy components as the system size approaches the thermodynamic limit. We demonstrate that the recently proposed staggered mesh method for Fock exchange energy [Xing, Li, and Lin, Math. Comp., 2024], which is specifically designed to sidestep certain singularities in exchange energy evaluation, can expedite the finite-size convergence rate for the exact exchange energy across a range of insulators and semiconductors when compared to the regular and truncated Coulomb methods. This remains true even for two computationally cheaper versions of this new method, which we call Non-SCF and Split-SCF staggered mesh. Additionally, a sequence of numerical tests on simple solids showcases the staggered mesh method's ability to improve convergence towards the thermodynamic limit for band gaps, bulk moduli, equilibrium lattice dimensions, energies, and phonon force constants. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 31 pages, 12 figures, submitted to the Journal of Chemical Theory and Computation

Programmable Simulations of Molecules and Materials with Reconfigurable Quantum Processors

Authors: Nishad Maskara, Stefan Ostermann, James Shee, Marcin Kalinowski, Abigail McClain Gomez, Rodrigo Araiza Bravo, Derek S. Wang, Anna I. Krylov, Norman Y. Yao, Martin Head-Gordon, Mikhail D. Lukin, Susanne F. Yelin

Abstract: Simulations of quantum chemistry and quantum materials are believed to be among the most important potential applications of quantum information processors, but realizing practical quantum advantage for such problems is challenging. Here, we introduce a simulation framework for strongly correlated quantum systems that can be represented by model spin Hamiltonians. Our approach leverages reconfigur… ▽ More Simulations of quantum chemistry and quantum materials are believed to be among the most important potential applications of quantum information processors, but realizing practical quantum advantage for such problems is challenging. Here, we introduce a simulation framework for strongly correlated quantum systems that can be represented by model spin Hamiltonians. Our approach leverages reconfigurable qubit architectures to programmably simulate real-time dynamics and introduces an algorithm for extracting chemically relevant spectral properties via classical co-processing of quantum measurement results. We develop a digital-analog simulation toolbox for efficient Hamiltonian time evolution utilizing digital Floquet engineering and hardware-optimized multi-qubit operations to accurately realize complex spin-spin interactions, and as an example present an implementation proposal based on Rydberg atom arrays. Then, we show how detailed spectral information can be extracted from these dynamics through snapshot measurements and single-ancilla control, enabling the evaluation of excitation energies and finite-temperature susceptibilities from a single-dataset. To illustrate the approach, we show how this method can be used to compute key properties of a polynuclear transition-metal catalyst and 2D magnetic materials. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 21 pages and 11 figures, plus supplementary information

Generalization of one-center non orthogonal configuration interaction singles to open shell singlet reference states: Theory and application to valence-core pump-probe states in acetylacetone

Authors: Juan E. Arias-Martinez, Hamlin Wu, Martin Head-Gordon

Abstract: We formulate a one-center non-orthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell mole… ▽ More We formulate a one-center non-orthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With inclusion of core hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMO). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely the core-to-SOMO transition as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Finally, simulations of the TR-NEXAS of acetylacetone after excitation onto its pi to pi-star singlet excited state at the carbon K-edge - an experiment carried out recently - showcases the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on non-adiabatic molecular dynamics simulations. △ Less

Submitted 14 October, 2023; originally announced October 2023.

Comments: Manuscript: 42 pages, 6 figures. SI: 48 pages, 10 figures

Optimizing the Regularization in Size-Consistent Second-Order Brillouin-Wigner Perturbation Theory

Authors: Kevin Carter-Fenk, James Shee, Martin Head-Gordon

Abstract: Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improv… ▽ More Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter ($αあるふぁ$) was set to recover the exact dissociation limit of H2 in a minimal basis set. Alternatively $αあるふぁ$ can be viewed as a regularization parameter, where each value of $αあるふぁ$ represents a valid variant of BW-s2, which we denote as BW-s2($αあるふぁ$). In this work, we semi-empirically optimize $αあるふぁ$ for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the $\textit{ab initio}$ parameterization of BW-s2 and MP2. We demonstrate that the optimal $αあるふぁ$ parameter ($αあるふぁ = 4$) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2($αあるふぁ$) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2($αあるふぁ$) is iterative $O(N^5)$, this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost. △ Less

Submitted 4 October, 2023; v1 submitted 4 September, 2023; originally announced September 2023.

Comments: 7 pages main text, 7 pages supporting information, 11 figures

The Doubles Connected Moments Expansion: A Tractable Approximate Horn-Weinstein Approach for Quantum Chemistry

Authors: Brad Ganoe, Martin Head-Gordon

Abstract: Ab initio methods based on the second-order and higher connected moments, or cumulants, of a reference function have seen limited use in the determination of correlation energies of chemical systems throughout the years. Moment-based methods have remained unattractive relative to more ubiquitous methods, such as perturbation theory and coupled cluster theory, due in part to the intractable cost of… ▽ More Ab initio methods based on the second-order and higher connected moments, or cumulants, of a reference function have seen limited use in the determination of correlation energies of chemical systems throughout the years. Moment-based methods have remained unattractive relative to more ubiquitous methods, such as perturbation theory and coupled cluster theory, due in part to the intractable cost of assembling moments of high-order and poor performance of low-order expansions. Many of the traditional quantum chemical methodologies can be recast as a selective summation of perturbative contributions to their energy; using this familiar structure as a guide in selecting terms, we develop a scheme to approximate connected moments limited to double excitations. The tractable Double Connected Moments (DCM(N)) approximation is developed and tested against a multitude of common single-reference methods to determine its efficacy in the determination of the correlation energy of model systems and small molecules. The DCM(N) sequence of energies exhibits smooth convergence, with compute costs that scale as a non-iterative O(N^6) with molecule size, M. Numerical tests on correlation energy recovery for 55 small molecules comprising the G1 test set in the cc-pVDZ basis show that DCM(N) strongly outperforms MP2 and even CCD with a Hartree-Fock reference. When using an approximate Brueckner reference from orbital-optimized (oo) MP2, the resulting oo:DCM(N) energies converge to values more accurate than CCSD for 49 of 55 molecules. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Earth Mover's Distance as a metric to evaluate the extent of charge transfer in excitations using discretized real-space densities

Authors: Zhe Wang, Jiashu Liang, Martin Head-Gordon

Abstract: This paper presents a novel theoretical measure, $μみゅー^{\text{EMD}}$, based on the Earth Mover's Distance, for quantifying the density shift caused by electronic excitations in molecules. As input, the EMD metric uses only the discretized ground and excited state electron densities in real space, rendering it compatible with almost all electronic structure methods used to calculate excited states. Th… ▽ More This paper presents a novel theoretical measure, $μみゅー^{\text{EMD}}$, based on the Earth Mover's Distance, for quantifying the density shift caused by electronic excitations in molecules. As input, the EMD metric uses only the discretized ground and excited state electron densities in real space, rendering it compatible with almost all electronic structure methods used to calculate excited states. The EMD metric is compared against other popular theoretical metrics for describing the extent of electron-hole separation in a wide range of excited states (valence, Rydberg, charge-transfer, etc). The results showcase the EMD metric's effectiveness across all excitation types and suggest that it is useful as an additional tool to characterize electronic excitations. The study also reveals that $μみゅー^{\text{EMD}}$ can function as a promising diagnostic tool for predicting the failure of pure exchange-correlation functionals. Specifically, we show statistical relationships between the functional-driven errors, the exact exchange content within the functional, and the magnitude of $μみゅー^{\text{EMD}}$ values. △ Less

Submitted 6 October, 2023; v1 submitted 14 August, 2023; originally announced August 2023.

Comments: 31 pages, 7 figures

Fantastical Excited State Optimized Structures and Where to Find Them

Authors: Justin J. Talbot, Juan E. Arias-Martinez, Stephen J. Cotton, Martin Head-Gordon

Abstract: The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximatio… ▽ More The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximation, where the nuclear and electronic degrees of freedom separate. This work demonstrates through ab initio calculations and simple nonadiabatic dynamics that some excited state minimum structures are fantastical: they appear to exist as stable configurations only as a consequence of the PES construct, rather than being physically observable. One such case is the S2 excited state of phosphine and a second case are local minima of a number of states of tris(bipyridine)ruthenium(II). Each fantastical structure exhibits an unphysically high predicted harmonic frequency and associated force constant. This fact can serve as a valuable diagnostic of when an optimized excited state structure is non-observable. Their origin lies in the coupling between different electronic states, and the resulting avoided crossings. The upper state may exhibit a minimum very close to the crossing, where the force constant relates to the strength of the electronic coupling rather than to any characteristic excited state vibration. Nonadiabatic dynamics results using a Landau-Zener model illustrate that fantastical excited state structures have extremely short lifetimes on the order of a few femtoseconds. Their appearance in a calculation signals the presence of a nearby avoided crossing or conical intersection through which the system will rapidly cross to a lower surface. △ Less

Submitted 12 August, 2023; originally announced August 2023.

Comments: 7 pages, 4 figures

Greater Transferability and Accuracy of Norm-conserving Pseudopotentials using Nonlinear Core Corrections

Authors: Wan-Lu Li, Kaixuan Chen, Elliot Rossomme, Martin Head-Gordon, Teresa Head-Gordon

Abstract: We present an investigation into the transferability of pseudopotentials (PPs) with a nonlinear core correction (NLCC) using the Goedecker, Teter, and Hutter (GTH) protocol across a range of pure GGA, meta-GGA and hybrid functionals, and their impact on thermochemical and non-thermochemical properties. The GTH-NLCC PP for the PBE density functional demonstrates remarkable transferability to the PB… ▽ More We present an investigation into the transferability of pseudopotentials (PPs) with a nonlinear core correction (NLCC) using the Goedecker, Teter, and Hutter (GTH) protocol across a range of pure GGA, meta-GGA and hybrid functionals, and their impact on thermochemical and non-thermochemical properties. The GTH-NLCC PP for the PBE density functional demonstrates remarkable transferability to the PBE0 and $ωおめが$B97X-V exchange-correlation functionals, and relative to no NLCC, improves agreement significantly for thermochemical benchmarks compared to all-electron calculations. On the other hand, the B97M-rV meta-GGA functional performs poorly with the PBE-derived GTH-NLCC PP, which is mitigated by reoptimizing the NLCC parameters for this specific functional. The findings reveal that atomization energies exhibit the greatest improvements from use of the NLCC, which thus provides an important correction needed for covalent interactions relevant to applications involving chemical reactivity. Finally we test the NLCC-GTH PPs when combined with medium-size TZV2P molecularly optimized (MOLOPT) basis sets which are typically utilized in condensed phase simulations, and show that they lead to consistently good results when compared to all-electron calculations for atomization energies, ionization potentials, barrier heights, and non-covalent interactions, but lead to somewhat larger errors for electron affinities. △ Less

Submitted 18 July, 2023; originally announced July 2023.

A New View on Density Corrected DFT: Can One Get a Better Answer for a Good Reason?

Authors: Devin J. Hernandez, Adam Rettig, Martin Head-Gordon

Abstract: Despite its widespread use, density functional theory (DFT) has several notable areas of failure; perhaps the most well-studied of these failures is self-interaction error (SIE). Density corrected DFT (DC-DFT) was proposed as a potential solution to systems where SIE causes traditional DFT to fail. The Hartree-Fock (HF) density is then used for cases where the DFT energy is suitable but the self-c… ▽ More Despite its widespread use, density functional theory (DFT) has several notable areas of failure; perhaps the most well-studied of these failures is self-interaction error (SIE). Density corrected DFT (DC-DFT) was proposed as a potential solution to systems where SIE causes traditional DFT to fail. The Hartree-Fock (HF) density is then used for cases where the DFT energy is suitable but the self-consistent density is erroneous. In this study, we investigate the utility of the higher quality orbital optimized MP2 densities in DC-DFT for barrier heights and halogen bonded complexes. For functionals such as PBE and r$^2$SCAN, find that these densities yield worse results than the HF density due to favorable cancellation between the density-driven and functional-driven errors, confirming a recent study. Error decomposition reveals functional driven error, not density driven error, to be the primary cause of inaccuracy in DFT calculations where SIE is prominent. We therefore advise caution when using HF-DFT, because the only rigorous way to remove large functional-driven errors in lower rungs of Jacob's ladder is by climbing to higher rungs that include exact exchange. We recommend that better functionals be improved by using a better density in SIE-sensitive cases. Examples support the value of this variant of DC-DFT. We also emphasize that DC-DFT potential energy surfaces have first derivative discontinuities at Coulson-Fischer points, in contrast to the second derivative discontinuities in SCF solutions. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Highly Accurate Prediction of NMR Chemical Shifts from Low-Level Quantum Mechanics Calculations Using Machine Learning

Authors: Jie Li, Jiashu Liang, Zhe Wang, Aleksandra L. Ptaszek, Xiao Liu, Brad Ganoe, Martin Head-Gordon, Teresa Head-Gordon

Abstract: Theoretical predictions of NMR chemical shifts from first-principles can greatly facilitate experimental interpretation and structure identification. However, accurate prediction of chemical shifts using the best coupled cluster methods can be prohibitively expensive for systems larger than ten to twenty non-hydrogen atoms on today's computers. By contrast machine learning methods offer inexpensiv… ▽ More Theoretical predictions of NMR chemical shifts from first-principles can greatly facilitate experimental interpretation and structure identification. However, accurate prediction of chemical shifts using the best coupled cluster methods can be prohibitively expensive for systems larger than ten to twenty non-hydrogen atoms on today's computers. By contrast machine learning methods offer inexpensive alternatives but are hampered by generalization to molecules outside the original training set. Here we propose a novel machine learning feature representation informed by intermediate calculations of atomic chemical shielding tensors within a molecular environment using an inexpensive quantum mechanics method, and training it to predict NMR chemical shieldings of a high-level composite theory that is comparable to CCSD(T) in the complete basis set limit. The inexpensive shift machine learning (iShiftML) algorithm is trained through a new progressive active learning workflow that reduces the total number of expensive calculations required when constructing the dataset, while allowing the model to continuously improve on data it has never seen. Furthermore, we show that the error estimations from our model correlate quite well with actual errors to provide confidence values on new predictions. We illustrate the predictive capacity of iShiftML across gas phase experimental chemical shifts for small organic molecules and much larger and more complex natural products in which we can accurately differentiate between subtle diastereomers based on chemical shift assignments. △ Less

Submitted 14 June, 2023; originally announced June 2023.

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Authors: Adam Rettig, Joonho Lee, Martin Head-Gordon

Abstract: Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals. This is done to lower the computational scaling with respect to the number of basis func… ▽ More Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals. This is done to lower the computational scaling with respect to the number of basis functions ($N$), and $\mathbf k$-points ($N_k$). Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e. a set of points, $N_\text{ISDF}$) to further reduce computation time and memory usage. This algorithm has linear scaling cost with $\mathbf k$-points, no explicit dependence of $N_\text{ISDF}$ on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized systems, with additional gains for large systems. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. We perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Journal ref: J. Chem. Theory Comput. 2023, 19, 17, 5773-5784

Probing C-I Bond Fission in the UV Photochemistry of 2-Iodothiophene with Core-to-Valence Transient Absorption Spectroscopy

Authors: Benjamin W. Toulson, Diptarka Hait, Davide Faccialà, Daniel M. Neumark, Stephen R. Leone, Martin Head-Gordon, Oliver Gessner

Abstract: The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photoinduced transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268 nm UV light-induced dynamics in 2-iodothiophene (C$_4$H$_3$IS) is performed. The dynamics are experimentally m… ▽ More The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photoinduced transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268 nm UV light-induced dynamics in 2-iodothiophene (C$_4$H$_3$IS) is performed. The dynamics are experimentally monitored with a femtosecond XUV probe pulse that measures iodine N-edge 4d core-to-valence transitions. Experiments are complemented by density functional theory calculations of both the pump-pulse induced valence excitations as well as the XUV probe-induced core-to-valence transitions. Possible intramolecular relaxation dynamics are investigated by ab initio molecular dynamics simulations. Gradual absorption changes up to ~0.5-1 ps after excitation are observed for both the parent molecular species and emerging iodine fragments, with the latter appearing with a characteristic rise time of 160$\pm$30 fs. Comparison of spectral intensities and energies with the calculations identify an iodine dissociation pathway initiated by a predominant $πぱい\toπぱい^*$ excitation. In contrast, initial excitation to a nearby n$_\perp\toσしぐま^*$ excited state appears unlikely based on a significantly smaller oscillator strength and the absence of any corresponding XUV absorption signatures. Excitation to the $πぱい\toπぱい^*$ state is followed by contraction of the C-I bond, enabling a nonadiabatic transition to a dissociative $πぱい\toσしぐま_\textrm{C-I}^*$ state. For the subsequent fragmentation, a narrow bond-length region along the C-I stretch coordinate between 230 and 280 pm is identified, where the transition between the parent molecule and the thienyl radical + iodine atom products becomes prominent in the XUV spectrum due to rapid localization of two singly-occupied molecular orbitals on the two fragments. △ Less

Submitted 30 June, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Journal ref: J. Chem. Phys. 159, 034304 (2023)

Repartitioned Brillouin-Wigner Perturbation Theory with a Size-Consistent Second-Order Correlation Energy

Authors: Kevin Carter-Fenk, Martin Head-Gordon

Abstract: Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regul… ▽ More Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size-consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through second order, is size-extensive, size-consistent (provided its Hartree-Fock reference is also), and orbital invariant. Our second-order size-consistent Brillouin-Wigner (BW-s2) approach is capable of describing the exact dissociation limit of H$_2$ in a minimal basis set regardless of the spin-polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, while rivaling coupled-cluster with single and double substitutions (CCSD) for thermochemical properties. △ Less

Submitted 24 May, 2023; v1 submitted 10 March, 2023; originally announced March 2023.

Comments: 15 pages, 7 figures

Femtosecond Symmetry Breaking and Coherent Relaxation of Methane Cations at the Carbon K-Edge

Authors: Enrico Ridente, Diptarka Hait, Eric A. Haugen, Andrew D. Ross, Daniel M. Neumark, Martin Head-Gordon, Stephen R. Leone

Abstract: Understanding the relaxation pathways of photoexcited molecules is essential to gain atomistic level insight into photochemistry. Herein, we perform a time-resolved study of ultrafast molecular symmetry breaking via geometric relaxation (Jahn-Teller distortion) on the methane cation. Attosecond transient absorption spectroscopy with soft X-rays at the carbon K-edge reveals that the distortion occu… ▽ More Understanding the relaxation pathways of photoexcited molecules is essential to gain atomistic level insight into photochemistry. Herein, we perform a time-resolved study of ultrafast molecular symmetry breaking via geometric relaxation (Jahn-Teller distortion) on the methane cation. Attosecond transient absorption spectroscopy with soft X-rays at the carbon K-edge reveals that the distortion occurs within $10\pm 2$ femtoseconds after few-femtosecond strong-field ionization of methane. The distortion activates coherent oscillations in the scissoring vibrational mode of the symmetry broken cation, which are detected in the X-ray signal. These oscillations are damped within $58\pm13$ femtoseconds, as vibrational coherence is lost with the energy redistributing into lower-frequency vibrational modes. This study completely reconstructs the molecular relaxation dynamics of this prototypical example and opens new avenues for exploring complex systems. △ Less

Submitted 25 December, 2022; originally announced December 2022.

Journal ref: Science 380,713-717(2023)

Ultrafast X-ray Spectroscopy of Intersystem Crossing in Hexafluoroacetylacetone: Chromophore Photophysics and Spectral Changes in the Face of Electron Withdrawing Groups

Authors: Eric A. Haugen, Diptarka Hait, Valeriu Scutelnic, Tian Xue, Martin Head-Gordon, Stephen R. Leone

Abstract: Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. We observe the excited state dynamics following excitation with 266 nm UV light to the $^1πぱいπぱい^{*}$ (S$_2$) state with element and… ▽ More Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. We observe the excited state dynamics following excitation with 266 nm UV light to the $^1πぱいπぱい^{*}$ (S$_2$) state with element and site-specificity using a broadband soft X-ray pulse produced by high harmonic generation. These results are compared to X-ray spectra computed from orbital optimized density functional theory methods. It is found that the electron withdrawing fluorine atoms decongest the X-ray absorption spectrum by enhancing separation between features originating from different carbon atoms. This facilitates the elucidation of structural and electronic dynamics at the chromophore. The evolution of the core-to-valence resonances at the carbon K-edge reveals an ultrafast population transfer between the $^1nπぱい^{*}$ (S$_1$) and $^3πぱいπぱい^{*}$ (T$_1$) states on a $1.6\pm0.4$ ps timescale, which is similar to the 1.5 ps timescale earlier observed for acetylacetone [J. Am. Chem. Soc. 139, 16576 (2017)]. It therefore appears that terminal fluorination has little influence on the intersystem crossing rate of the acetylacetone chromophore. In addition, the significant role of hydrogen-bond opened and twisted rotational isomers is elucidated in the excited state dynamics by comparison of the experimental transient X-ray spectra with theory. △ Less

Submitted 12 December, 2022; v1 submitted 27 September, 2022; originally announced September 2022.

Journal ref: J. Phys. Chem. A 2023

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm-Dancoff approximation

Authors: Justin J. Talbot, Martin Head-Gordon, Stephen J. Cotton

Abstract: The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavorable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm-Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realisti… ▽ More The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavorable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm-Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately-sized basis sets. In this work, some of the benefits and limitations of TDDFT/TDA are discussed and analyzed with regard to its applicability as a "back-end" electronic structure method for the symmetric quasi-classical Meyer-Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyze a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 35 pages, 8 figures, 2 tables

Efficient Calculation of NMR Shielding Constants Using Composite Method Approximations and Locally Dense Basis Sets

Authors: Jiashu Liang, Zhe Wang, Jie Li, Jonathan Wong, Xiao Liu, Brad Ganoe, Teresa Head-Gordon, Martin Head-Gordon

Abstract: This paper presents a systematic study of applying composite method approximations with locally dense basis sets (LDBS) to efficiently calculate NMR shielding constants in small and medium-sized molecules. The pcSseg-n series of basis sets are shown to have similar accuracy to the pcS-n series when n $\geq1$ and can slightly reduce compute costs. We identify two different LDBS partition schemes th… ▽ More This paper presents a systematic study of applying composite method approximations with locally dense basis sets (LDBS) to efficiently calculate NMR shielding constants in small and medium-sized molecules. The pcSseg-n series of basis sets are shown to have similar accuracy to the pcS-n series when n $\geq1$ and can slightly reduce compute costs. We identify two different LDBS partition schemes that perform very effectively for density functional calculations. We select a large subset of the recent NS372 database containing 290 H, C, N, and O shielding values evaluated by reference methods on 106 molecules to carefully assess methods of the high, medium, and low compute costs to make practical recommendations. Our assessment covers conventional electronic structure methods (DFT and wavefunction) with global basis calculations, as well as their use in one of the satisfactory LDBS approaches, and a range of composite approaches, also with and without LDBS. Altogether 99 methods are evaluated. On this basis, we recommend different methods to reach three different levels of accuracy and time requirements across the four nuclei considered. △ Less

Submitted 11 November, 2022; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: 29 pages, 7 figures, 2 tables

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

Authors: Jiashu Liang, Xintian Feng, Martin Head-Gordon

Abstract: VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωおめがB97X-V and the hybrid mGGA, ωおめがB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and eff… ▽ More VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωおめがB97X-V and the hybrid mGGA, ωおめがB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. This study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωおめがB97M-V, and ωおめがB97X-V perform very well. The convergence of frequencies with respect to grid size and atomic orbital basis set size is studied and recommendations reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωおめがB97X-V, M06-SX, and ωおめがB97M-V). △ Less

Submitted 2 April, 2023; v1 submitted 31 August, 2022; originally announced August 2022.

Comments: 30 pages, 5 figures, 4 tables

Faster Exact Exchange for Solids via occ-RI-K: Application to Combinatorially Optimized Range-Separated Hybrid Functionals for Simple Solids with Pseudopotentials Near the Basis Set Limit

Authors: Joonho Lee, Adam Rettig, Xintian Feng, Evgeny Epifanovsky, Martin Head-Gordon

Abstract: In this work, we developed and showcased the occ-RI-K algorithm to compute the exact exchange contribution in density functional calculations of solids near the basis set limit. Within the gaussian planewave (GPW) density fitting, our algorithm achieves a 1-2 orders of magnitude speedup compared to conventional GPW algorithms. Since our algorithm is well-suited for simulations with large basis set… ▽ More In this work, we developed and showcased the occ-RI-K algorithm to compute the exact exchange contribution in density functional calculations of solids near the basis set limit. Within the gaussian planewave (GPW) density fitting, our algorithm achieves a 1-2 orders of magnitude speedup compared to conventional GPW algorithms. Since our algorithm is well-suited for simulations with large basis sets, we applied it to 12 hybrid density functionals with pseudopotentials and a large uncontracted basis set to assess their performance on band gaps of 25 simple solids near the basis set limit. The largest calculation performed in this work involves 16 electrons and 350 basis functions in the unit cell utilizing a 6x6x6 k-mesh. With 20-27\% exact exchange, global hybrid functionals (B3LYP, PBE0, revPBE0, B97-3, SCAN0) perform similarly with a root-mean-square-deviation (RMSD) of 0.61-0.77 eV while other global hybrid functionals such as M06-2X (2.02 eV) and MN15 (1.05 eV) show higher RMSD due to their increased fraction of exact exchange. A short-range hybrid functional, HSE achieves a similar RMSD (0.76 eV) but shows a noticeable underestimation of band gaps due to the complete lack of long-range exchange. We found that two combinatorially optimized range-separated hybrid functionals, $ωおめが$B97X-rV (3.94 eV) and $ωおめが$B97M-rV (3.40 eV), and the two other range separated hybrid functionals, CAM-B3LYP (2.41 eV) and CAM-QTP01 (4.16 eV), significantly overestimate the band gap because of their high fraction of long-range exact exchange. Given the failure of $ωおめが$B97X-rV and $ωおめが$B97M-rV, we have yet to find a density functional that offers consistent performance for both molecules and solids. Our algorithm development and density functional assessment will serve as a stepping stone towards developing more accurate hybrid functionals and applying them to practical applications. △ Less

Submitted 8 November, 2022; v1 submitted 18 July, 2022; originally announced July 2022.

Say NO to Optimization: A Non-Orthogonal Quantum Eigensolver

Authors: Unpil Baek, Diptarka Hait, James Shee, Oskar Leimkuhler, William J. Huggins, Torin F. Stetina, Martin Head-Gordon, K. Birgitta Whaley

Abstract: A balanced description of both static and dynamic correlations in electronic systems with nearly degenerate low-lying states presents a challenge for multi-configurational methods on classical computers. We present here a quantum algorithm utilizing the action of correlating cluster operators to provide high-quality wavefunction ansätze employing a non-orthogonal multireference basis that captures… ▽ More A balanced description of both static and dynamic correlations in electronic systems with nearly degenerate low-lying states presents a challenge for multi-configurational methods on classical computers. We present here a quantum algorithm utilizing the action of correlating cluster operators to provide high-quality wavefunction ansätze employing a non-orthogonal multireference basis that captures a significant portion of the exact wavefunction in a highly compact manner, and that allows computation of the resulting energies and wavefunctions at polynomial cost with a quantum computer. This enables a significant improvement over the corresponding classical non-orthogonal solver, which incurs an exponential cost when evaluating off-diagonal matrix elements between the ansatz states, and is therefore intractable. We implement the non-orthogonal quantum eigensolver (NOQE) here with an efficient ansatz parameterization inspired by classical quantum chemistry methods that succeed in capturing significant amounts of electronic correlation accurately. By taking advantage of classical methods for chemistry, NOQE provides a flexible, compact, and rigorous description of both static and dynamic electronic correlation, making it an attractive method for the calculation of electronic states of a wide range of molecular systems. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 27+11 pages, 14+4 figures

Electron-Affinity Time-Dependent Density Functional Theory: Formalism and Applications to Core-Excited States

Authors: Kevin Carter-Fenk, Leonardo A. Cunha, Juan E. Arias-Martinez, Martin Head-Gordon

Abstract: The particle-hole interaction problem is longstanding within time-dependent density functional theory (TDDFT) and leads to extreme errors in the prediction of K-edge X-ray absorption spectra (XAS). We derive a linear-response formalism that uses optimized orbitals of the n-1-electron system as reference, building orbital relaxation and a proper hole into the initial density. Our approach is an exa… ▽ More The particle-hole interaction problem is longstanding within time-dependent density functional theory (TDDFT) and leads to extreme errors in the prediction of K-edge X-ray absorption spectra (XAS). We derive a linear-response formalism that uses optimized orbitals of the n-1-electron system as reference, building orbital relaxation and a proper hole into the initial density. Our approach is an exact generalization of the static-exchange approximation that ameliorates particle-hole interaction error associated with the adiabatic approximation and reduces errors in TDDFT XAS by orders of magnitude. With a statistical performance of just 0.5 eV root-mean-square error and the same computational scaling as TDDFT under the core-valence separation approximation, we anticipate that this approach will be of great utility in XAS calculations of large systems. △ Less

Submitted 19 September, 2022; v1 submitted 17 May, 2022; originally announced May 2022.

Comments: 11 pages (main text), 12 pages (supporting information), 9 figures

Jahn-Teller Distortion and Dissociation of CCl$_4^+$ by Transient X-ray Spectroscopy Simultaneously at the Carbon K- and Chlorine L-Edge

Authors: Andrew D. Ross, Diptarka Hait, Valeriu Scutelnic, Eric A. Haugen, Enrico Ridente, Mikias B. Balkew, Daniel M. Neumark, Martin Head-Gordon, Stephen R. Leone

Abstract: X-ray Transient Absorption Spectroscopy (XTAS) and theoretical calculations are used to study CCl$_4^+$ prepared by 800 nm strong-field ionization. XTAS simultaneously probes atoms at the carbon K-edge (280-300 eV) and chlorine L-edge (195-220 eV). Comparison of experiment to X-ray spectra computed by orbital-optimized density functional theory (OO-DFT) indicates that after ionization, CCl$_4^+$ u… ▽ More X-ray Transient Absorption Spectroscopy (XTAS) and theoretical calculations are used to study CCl$_4^+$ prepared by 800 nm strong-field ionization. XTAS simultaneously probes atoms at the carbon K-edge (280-300 eV) and chlorine L-edge (195-220 eV). Comparison of experiment to X-ray spectra computed by orbital-optimized density functional theory (OO-DFT) indicates that after ionization, CCl$_4^+$ undergoes symmetry breaking driven by Jahn-Teller distortion away from the initial tetrahedral structure (T$_d$) in 6$\pm$2 fs. The resultant symmetry-broken covalently bonded form subsequently separates to a noncovalently bound complex between CCl$_3^+$ and Cl over 90$\pm$10 fs, which is again predicted by theory. Finally, after more than 800 fs, L-edge signals for atomic Cl are observed, indicating dissociation to free CCl$_3^+$ and Cl. The results for Jahn-Teller distortion to the symmetry-broken form of CCl$_4^+$ and formation of the Cl -- CCl$_3^+$ complex characterize previously unobserved new species along the route to dissociation. △ Less

Submitted 6 June, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Journal ref: Chem. Sci., 2022,13, 9310-9320

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

Authors: Juan E. Arias-Martinez, Leonardo A. Cunha, Katherine J. Oosterbaan, Joonho Lee, Martin Head-Gordon

Abstract: We investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without any use of response theory or equation-of-motion formalisms. Three schemes are employed to successfully address the convergence difficulties associated… ▽ More We investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without any use of response theory or equation-of-motion formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for Delta-CCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed Delta-CCSD schemes seem to fare better than the widely used equation-of-motion CCSD (EOM-CCSD) with core-valence separation protocol, with statistical errors being reduced by more than a factor of two when compared to FC-CVS-EOM-CCSD. △ Less

Submitted 11 July, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

Computing X-ray absorption spectra from linear-response particles atop optimized holes

Authors: Diptarka Hait, Katherine J. Oosterbaan, Kevin Carter-Fenk, Martin Head-Gordon

Abstract: State specific orbital optimized density functional theory (OO-DFT) methods like restricted open-shell Kohn-Sham (ROKS) can attain semiquantitative accuracy for predicting X-ray absorption spectra of closed-shell molecules. OO-DFT methods however require that each state be individually optimized. In this work, we present an approach to generate an approximate core-excited state density for use wit… ▽ More State specific orbital optimized density functional theory (OO-DFT) methods like restricted open-shell Kohn-Sham (ROKS) can attain semiquantitative accuracy for predicting X-ray absorption spectra of closed-shell molecules. OO-DFT methods however require that each state be individually optimized. In this work, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, that is capable of giving reasonable accuracy without requiring state-specific optimization. This is achieved by fully optimizing the core-hole through the core-ionized state, followed by use of electron-addition configuration interaction singles (EA-CIS) to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method, and can attain $\sim 0.6$ eV RMS error for the K-edges of C-F through the use of local functionals like PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment, and can thus help reduce the computational cost for obtaining OO-DFT quality spectra. ROKS(STEX) therefore appears to be a useful technique for efficient prediction of X-ray absorption spectra. △ Less

Submitted 25 March, 2022; originally announced March 2022.

Revisiting the performance of time-dependent density functional theory for electronic excitations: Assessment of 43 popular and recently developed functionals from rungs one to four

Authors: Jiashu Liang, Xintian Feng, Diptarka Hait, Martin Head-Gordon

Abstract: In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between performance and accuracy. The functionals $ωおめが$B97X-D and BMK are found to offer the best performance ov… ▽ More In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between performance and accuracy. The functionals $ωおめが$B97X-D and BMK are found to offer the best performance overall with a Root-Mean Square Error (RMSE) of 0.28 eV, better than the computationally more demanding CIS(D) wavefunction method with a RMSE of 0.36 eV. The results also suggest that Jacob's ladder still holds for TDDFT excitation energies, though hybrid meta-GGAs are not generally better than hybrid GGAs. Effects of basis set convergence, gauge invariance correction to meta-GGAs, and nonlocal correlation (VV10) are also studied, and practical basis set recommendations are provided. △ Less

Submitted 18 May, 2022; v1 submitted 26 February, 2022; originally announced February 2022.

Comments: 12 pages, 8 figures; 8 pages, 7 figures for supporting info

Changes in polarization dictate necessary approximations for modeling electronic de-excitation intensity: an application to X-ray emission

Authors: Subhayan Roychoudhury, Leonardo A. Cunha, Martin Head-Gordon, David Prendergast

Abstract: We systematically investigate the underlying relations among different levels of approximation for simulating electronic de-excitations, with a focus on modeling X-ray emission spectroscopy (XES). Using Fermi's golden rule and explicit modeling of the initial, core-excited state and the final, valence-hole state, we show that XES can be accurately modeled by using orbital optimization for the vari… ▽ More We systematically investigate the underlying relations among different levels of approximation for simulating electronic de-excitations, with a focus on modeling X-ray emission spectroscopy (XES). Using Fermi's golden rule and explicit modeling of the initial, core-excited state and the final, valence-hole state, we show that XES can be accurately modeled by using orbital optimization for the various final states within a Slater-determinant framework. However, in this paper, we introduce a much cheaper approach reliant only on a single self-consistent field for all the final states, and show that it is typically sufficient. Further approximations reveal that these fundamentally many-body transitions can be reasonably approximated by projections of ground state orbitals, but that the ground state alone is insufficient. Furthermore, except in cases where the core-ionization induces negligible changes in polarization, linear-response approaches within the adiabatic approximation will have difficulty in accurately modeling de-excitation to the core level. Therefore, change in the net dipole moment of the valence electrons can serve as a metric for the validity of the linear-response approximation. △ Less

Submitted 18 December, 2021; originally announced December 2021.

Relativistic Orbital Optimized Density Functional Theory for Accurate Core-Level Spectroscopy

Authors: Leonardo A. Cunha, Diptarka Hait, Richard Kang, Yuezhi Mao, Martin Head-Gordon

Abstract: Core-level spectra of 1s electrons of elements heavier than Ne show significant relativistic effects. We combine advances in orbital optimized DFT (OO-DFT) with the spin-free exact two-component (X2C) model for scalar relativistic effects, to study K-edge spectra of third period elements. OO-DFT/X2C is found to be quite accurate at predicting energies, yielding $\sim 0.5$ eV RMS error vs experimen… ▽ More Core-level spectra of 1s electrons of elements heavier than Ne show significant relativistic effects. We combine advances in orbital optimized DFT (OO-DFT) with the spin-free exact two-component (X2C) model for scalar relativistic effects, to study K-edge spectra of third period elements. OO-DFT/X2C is found to be quite accurate at predicting energies, yielding $\sim 0.5$ eV RMS error vs experiment with the modern SCAN (and related) functionals. This marks a significant improvement over the $>50$ eV deviations that are typical for the popular time-dependent DFT (TDDFT) approach. Consequently, experimental spectra are quite well reproduced by OO-DFT/X2C, sans empirical shifts for alignment. OO-DFT/X2C combines high accuracy with ground state DFT cost and is thus a promising route for computing core-level spectra of third period elements. We also explored K and L edges of 3d transition metals to identify limitations of the OO-DFT/X2C approach in modeling the spectra of heavier atoms. △ Less

Submitted 30 March, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Journal ref: J. Phys. Chem. Lett. 2022, 13, 3438-3449

Approaching the Basis Set Limit in Gaussian-Orbital-Based Periodic Calculations with Transferability: Performance of Pure Density Functionals for Simple Semiconductors

Authors: Joonho Lee, Xintian Feng, Leonardo A. Cunha, Jerome F. Gonthier, Evgeny Epifanovsky, Martin Head-Gordon

Abstract: Simulating solids with quantum chemistry methods and Gaussian-type orbitals (GTOs) has been gaining popularity. Nonetheless, there are few systematic studies that assess the basis set incompleteness error (BSIE) in these GTO-based simulations over a variety of solids. In this work, we report a GTO-based implementation for solids, and apply it to address the basis set convergence issue. We employ a… ▽ More Simulating solids with quantum chemistry methods and Gaussian-type orbitals (GTOs) has been gaining popularity. Nonetheless, there are few systematic studies that assess the basis set incompleteness error (BSIE) in these GTO-based simulations over a variety of solids. In this work, we report a GTO-based implementation for solids, and apply it to address the basis set convergence issue. We employ a simple strategy to generate large uncontracted (unc) GTO basis sets, that we call the unc-def2-GTH sets. These basis sets exhibit systematic improvement towards the basis set limit as well as good transferability based on application to a total of 43 simple semiconductors. Most notably, we found the BSIE of unc-def2-QZVP-GTH to be smaller than 0.7 m$E_h$ per atom in total energies and 20 meV in band gaps for all systems considered here. Using unc-def2-QZVP-GTH, we report band gap benchmarks of a combinatorially designed meta generalized gradient functional (mGGA), B97M-rV, and show that B97M-rV performs similarly (a root-mean-square-deviation (RMSD) of 1.18 eV) to other modern mGGA functionals, M06-L (1.26 eV), MN15-L (1.29 eV), and SCAN (1.20 eV). This represents a clear improvement over older pure functionals such as LDA (1.71 eV) and PBE (1.49 eV) though all these mGGAs are still far from being quantitatively accurate. We also provide several cautionary notes on the use of our uncontracted bases and on future research on GTO basis set development for solids. △ Less

Submitted 7 October, 2021; v1 submitted 29 August, 2021; originally announced August 2021.

NewtonNet: A Newtonian message passing network for deep learning of interatomic potentials and forces

Authors: Mojtaba Haghighatlari, Jie Li, Xingyi Guan, Oufan Zhang, Akshaya Das, Christopher J. Stein, Farnaz Heidar-Zadeh, Meili Liu, Martin Head-Gordon, Luke Bertels, Hongxia Hao, Itai Leven, Teresa Head-Gordon

Abstract: We report a new deep learning message passing network that takes inspiration from Newton's equations of motion to learn interatomic potentials and forces. With the advantage of directional information from trainable latent force vectors, and physics-infused operators that are inspired by the Newtonian physics, the entire model remains rotationally equivariant, and many-body interactions are inferr… ▽ More We report a new deep learning message passing network that takes inspiration from Newton's equations of motion to learn interatomic potentials and forces. With the advantage of directional information from trainable latent force vectors, and physics-infused operators that are inspired by the Newtonian physics, the entire model remains rotationally equivariant, and many-body interactions are inferred by more interpretable physical features. We test NewtonNet on the prediction of several reactive and non-reactive high quality ab initio data sets including single small molecule dynamics, a large set of chemically diverse molecules, and methane and hydrogen combustion reactions, achieving state-of-the-art test performance on energies and forces with far greater data and computational efficiency than other deep learning models. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Exploring Spin Symmetry-Breaking Effects for Static Field Ionization of Atoms: Is There an Analog to the Coulson-Fischer Point in Bond Dissociation?

Authors: Leonardo A. Cunha, Joonho Lee, Diptarka Hait, C. William McCurdy, Martin Head-Gordon

Abstract: Löwdin's symmetry dilemma is an ubiquitous issue in approximate quantum chemistry. In the context of Hartree-Fock (HF) theory, the use of Slater determinants with some imposed constraints to preserve symmetries of the exact problem may lead to physically unreasonable potential energy surfaces. On the other hand, lifting these constraints leads to the so-called broken symmetry solutions that usuall… ▽ More Löwdin's symmetry dilemma is an ubiquitous issue in approximate quantum chemistry. In the context of Hartree-Fock (HF) theory, the use of Slater determinants with some imposed constraints to preserve symmetries of the exact problem may lead to physically unreasonable potential energy surfaces. On the other hand, lifting these constraints leads to the so-called broken symmetry solutions that usually provide better energetics, at the cost of losing information about good quantum numbers that describe the state of the system. This behavior has been previously extensively studied in the context of bond dissociation. This paper studies the behavior of different classes of Hartree-Fock spin polarized solutions (restricted, unrestricted, generalized) in the context of ionization by strong static electric fields. We find that, for simple two-electron systems, UHF is able to provide a qualitatively good description of states involved during the ionization process (neutral, singly-ionized and doubly ionized states), whereas RHF fails to describe the singly ionized state. For more complex systems, even though UHF is able to capture some of the expected characteristics of the ionized states, it is constrained to a single $M_s$ (diabatic) manifold in the energy surface as a function of field intensity. In this case a better qualitative picture can be painted by GHF as it is able to explore different spin manifolds and follow the lowest solution due to lack of collinearity constraints on the spin quantization axis. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Journal ref: J. Chem. Phys. 155, 014309 (2021)

Real time evolution for ultracompact Hamiltonian eigenstates on quantum hardware

Authors: Katherine Klymko, Carlos Mejuto-Zaera, Stephen J. Cotton, Filip Wudarski, Miroslav Urbanek, Diptarka Hait, Martin Head-Gordon, K. Birgitta Whaley, Jonathan Moussa, Nathan Wiebe, Wibe A. de Jong, Norm M. Tubman

Abstract: In this work we present a detailed analysis of variational quantum phase estimation (VQPE), a method based on real-time evolution for ground and excited state estimation on near-term hardware. We derive the theoretical ground on which the approach stands, and demonstrate that it provides one of the most compact variational expansions to date for solving strongly correlated Hamiltonians. At the cen… ▽ More In this work we present a detailed analysis of variational quantum phase estimation (VQPE), a method based on real-time evolution for ground and excited state estimation on near-term hardware. We derive the theoretical ground on which the approach stands, and demonstrate that it provides one of the most compact variational expansions to date for solving strongly correlated Hamiltonians. At the center of VQPE lies a set of equations, with a simple geometrical interpretation, which provides conditions for the time evolution grid in order to decouple eigenstates out of the set of time evolved expansion states, and connects the method to the classical filter diagonalization algorithm. Further, we introduce what we call the unitary formulation of VQPE, in which the number of matrix elements that need to be measured scales linearly with the number of expansion states, and we provide an analysis of the effects of noise which substantially improves previous considerations. The unitary formulation allows for a direct comparison to iterative phase estimation. Our results mark VQPE as both a natural and highly efficient quantum algorithm for ground and excited state calculations of general many-body systems. We demonstrate a hardware implementation of VQPE for the transverse field Ising model. Further, we illustrate its power on a paradigmatic example of strong correlation (Cr2 in the SVP basis set), and show that it is possible to reach chemical accuracy with as few as ~50 timesteps. △ Less

Submitted 7 April, 2021; v1 submitted 15 March, 2021; originally announced March 2021.

Orbital Optimized Density Functional Theory for Electronic Excited States

Authors: Diptarka Hait, Martin Head-Gordon

Abstract: Density functional theory (DFT) based modeling of electronic excited states is of importance for investigation of the photophysical/photochemical properties and spectroscopic characterization of large systems. The widely used linear response time-dependent DFT (TDDFT) approach is however not effective at modeling many types of excited states, including (but not limited to) charge-transfer states,… ▽ More Density functional theory (DFT) based modeling of electronic excited states is of importance for investigation of the photophysical/photochemical properties and spectroscopic characterization of large systems. The widely used linear response time-dependent DFT (TDDFT) approach is however not effective at modeling many types of excited states, including (but not limited to) charge-transfer states, doubly excited states and core-level excitations. In this perspective, we discuss state-specific orbital optimized (OO) DFT approaches as an alterative to TDDFT for electronic excited states. We motivate the use of OO-DFT methods and discuss reasons behind their relatively restricted historical usage (vs TDDFT). We subsequently highlight modern developments that address these factors and allow efficient and reliable OO-DFT computations. Several successful applications of OO-DFT for challenging electronic excitations are also presented, indicating their practical efficacy. OO-DFT approaches are thus increasingly becoming a useful route for computing excited states of large chemical systems. We conclude by discussing the limitations and challenges still facing OO-DFT methods, as well as some potential avenues for addressing them. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Journal ref: J. Phys. Chem. Lett. 2021, 12, 4517-4529

Assessment of the Performance of Density Functionals for Predicting Potential Energy Curves in Hydrogen Storage Applications

Authors: Srimukh Prasad Veccham, Martin Head-Gordon

Abstract: The availability of accurate computational tools for modeling and simulation is vital to accelerate the discovery of materials capable of storing hydrogen (H2) under given parameters of pressure swing and temperature. Previously, we compiled the H2Bind275 dataset consisting of equilibrium geometries and assessed the performance of 55 density functionals over this dataset (Veccham, S. P.; Head-Gord… ▽ More The availability of accurate computational tools for modeling and simulation is vital to accelerate the discovery of materials capable of storing hydrogen (H2) under given parameters of pressure swing and temperature. Previously, we compiled the H2Bind275 dataset consisting of equilibrium geometries and assessed the performance of 55 density functionals over this dataset (Veccham, S. P.; Head-Gordon, M. J. Chem. Theory Comput., 2020, 16, 4963--4982). As it is crucial for computational tools to accurately model the entire potential energy curve (PEC), in addition to the equilibrium geometry, we have extended this dataset with 389 new data points to include two compressed and three elongated geometries along 78 PECs for H2 binding, forming the H2Bind78x7 dataset. Assessing the performance of 55 density functionals on this significantly larger and more comprehensive H2Bind78x7 dataset, we have identified the best performing density functionals for H2 binding applications: PBE0-DH, $ωおめが$B97X-V, $ωおめが$B97M-V, and DSD-PBEPBE-D3(BJ). Addition of Hartree Fock exchange improves the performance of density functionals, albeit not uniformly throughout the PEC. We recommend the usage of wB97X-V and wB97M-V density functionals as they give good performance for both geometries and energies. In addition, we have also identified B97M-V and B97M-rV as the best semi-local density functionals for predicting H2 binding energy at its equilibrium geometry. △ Less

Submitted 4 February, 2021; originally announced February 2021.

Too big, too small or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Authors: Diptarka Hait, Yu Hsuan Liang, Martin Head-Gordon

Abstract: Multipole moments are the first order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory (DFT) prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero… ▽ More Multipole moments are the first order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory (DFT) prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher order moments origin-dependent. Therefore, instead of using the $3 \times 3$ quadrupole moment matrix, we utilize the translationally invariant $3 \times 3$ matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by $\mathcal{K}$). The principal components of ${\mathcal{K}}$ are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the prinicpal components of ${\mathcal{K}}$ for 100 small molecules at the coupled cluster singles and doubles with perturbative triples (CCSD(T)) at the complete basis set (CBS) limit is developed, resulting in 213 independent ${\mathcal{K}}$ components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 yield reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H and Be atoms in particular are challenging for nearly all methods, indicating that future functional development could benefit from inclusion of their density information in training or testing protocols. △ Less

Submitted 20 January, 2021; v1 submitted 25 November, 2020; originally announced November 2020.

Journal ref: J. Chem. Phys. 154 (7), 074109 (2021)

A Non-Perturbative Pairwise-Additive Analysis of Charge Transfer Contributions to Intermolecular Interaction Energies

Authors: Srimukh Prasad Veccham, Joonho Lee, Yuezhi Mao, Paul R. Horn, Martin Head-Gordon

Abstract: Energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs) decomposes the interaction energy between molecules into physically interpretable components like geometry distortion, frozen interactions, polarization, and charge transfer (CT, also sometimes called charge delocalization) interactions. In this work, a numerically exact scheme to decompose the CT interact… ▽ More Energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs) decomposes the interaction energy between molecules into physically interpretable components like geometry distortion, frozen interactions, polarization, and charge transfer (CT, also sometimes called charge delocalization) interactions. In this work, a numerically exact scheme to decompose the CT interaction energy into pairwise additive terms is introduced for the ALMO-EDA using density functional theory. Unlike perturbative pairwise charge-decomposition analysis, the new approach does not break down for strongly interacting systems, or show significant exchange-correlation functional dependence in the decomposed energy components. Both the energy lowering and the charge flow associated with CT can be decomposed. Complementary occupied-virtual orbital pairs (COVPs) that capture the dominant donor and acceptor CT orbitals are obtained for the new decomposition. It is applied to systems with different types of interactions including DNA base-pairs, borane-ammonia adducts, and transition metal hexacarbonyls. While consistent with most existing understanding of the nature of CT in these systems, the results also reveal some new insights into the origin of trends in donor-acceptor interactions. △ Less

Submitted 10 November, 2020; originally announced November 2020.

Journal ref: Phys. Chem. Chem. Phys., 2021, 23, 928-943

Third order Møller-Plesset theory made more useful? The role of density functional theory orbitals

Authors: Adam Rettig, Diptarka Hait, Luke W. Bertels, Martin Head-Gordon

Abstract: The practical utility of Møller-Plesset (MP) perturbation theory is severely constrained by the use of Hartree-Fock (HF) orbitals. It has recently been shown that use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error (J. Phys. Chem. Lett. 10, 4170, 2019). In this work we examine whether density functional theory (DFT) optimiz… ▽ More The practical utility of Møller-Plesset (MP) perturbation theory is severely constrained by the use of Hartree-Fock (HF) orbitals. It has recently been shown that use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error (J. Phys. Chem. Lett. 10, 4170, 2019). In this work we examine whether density functional theory (DFT) optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, non-covalent interactions, and dipole moments relative to standard HF based MP theory. Indeed MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled cluster singles and doubles (CCSD) for several datasets. We also found that the results are not particularly functional sensitive in most cases, (although range-separated hybrid functionals with low delocalization error perform the best). MP3 based on DFT orbitals thus appears to be an efficient, non-iterative $O(N^6)$ scaling wave function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate. △ Less

Submitted 4 January, 2021; v1 submitted 22 September, 2020; originally announced September 2020.

Journal ref: J. Chem. Theory Comput. 2020, 16, 12, 7473-7489

The Ground State Electronic Energy of Benzene

Authors: Janus J. Eriksen, Tyler A. Anderson, J. Emiliano Deustua, Khaldoon Ghanem, Diptarka Hait, Mark R. Hoffmann, Seunghoon Lee, Daniel S. Levine, Ilias Magoulas, Jun Shen, Norman M. Tubman, K. Birgitta Whaley, Enhua Xu, Yuan Yao, Ning Zhang, Ali Alavi, Garnet Kin-Lic Chan, Martin Head-Gordon, Wenjian Liu, Piotr Piecuch, Sandeep Sharma, Seiichiro L. Ten-no, C. J. Umrigar, Jürgen Gauss

Abstract: We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-$ζぜーた$ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire of… ▽ More We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-$ζぜーた$ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around $-863$ m$E_{\text{H}}$. However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 m$E_{\text{H}}$), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward. △ Less

Submitted 7 October, 2020; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: 29 pages, 1 figure, 2 tables. SI as an ancillary file

Journal ref: J. Phys. Chem. Lett., 11, 8922 (2020)

Polishing the Gold Standard: The Role of Orbital Choice in CCSD(T) Vibrational Frequency Prediction

Authors: Luke W. Bertels, Joonho Lee, Martin Head-Gordon

Abstract: While CCSD(T) with spin-restricted Hartree-Fock (RHF) orbitals has long been lauded for its ability to accurately describe closed-shell interactions, the performance of CCSD(T) on open-shell species is much more erratic, especially when using a spin-unrestricted HF (UHF) reference. Previous studies have shown improved treatment of open-shell systems when a non-HF set of molecular orbitals, like Br… ▽ More While CCSD(T) with spin-restricted Hartree-Fock (RHF) orbitals has long been lauded for its ability to accurately describe closed-shell interactions, the performance of CCSD(T) on open-shell species is much more erratic, especially when using a spin-unrestricted HF (UHF) reference. Previous studies have shown improved treatment of open-shell systems when a non-HF set of molecular orbitals, like Brueckner or Kohn-Sham density functional theory (DFT) orbitals, is used as a reference. Inspired by the success of regularized orbital-optimized second-order Møller-Plesset perturbation theory ($κかっぱ$-OOMP2) orbitals as reference orbitals for MP3, we investigate the use of $κかっぱ$-OOMP2 orbitals and various DFT orbitals as reference orbitals for CCSD(T) calculations of the corrected ground-state harmonic vibrational frequencies of a set of 36 closed-shell (29 neutrals, 6 cations, 1 anion) and 59 open-shell diatomic species (38 neutrals, 15 cations, 6 anions). The use of $κかっぱ$-OOMP2 orbitals in this context alleviates difficult cases observed for both UHF orbitals and OOMP2 orbitals. Removing 2 multireference systems and 12 systems with ambiguous experimental data leaves a pruned data set. Overall performance on the pruned data set highlights CCSD(T) with a B97 orbital reference (CCSD(T):B97), CCSD(T) with a $κかっぱ$-OOMP2 orbital reference (CCSD(T):$κかっぱ$-OOMP2), and CCSD(T) with a B97M-rV orbital reference (CCSD(T):B97M-rV) with RMSDs of 8.48 cm$^{-1}$ and 8.50 cm$^{-1}$, and 8.75 cm$^{-1}$ respectively, outperforming CCSD(T):UHF by nearly a factor of 5. Moreover, the performance on the closed- and open-shell subsets show these methods are able to treat open-shell and closed-shell systems with comparable accuracy and robustness. The use of $κかっぱ$-OOMP2 orbitals has also proven useful in diagnosing multireference character that can hinder the reliability of CCSD(T). △ Less

Submitted 16 July, 2020; v1 submitted 16 July, 2020; originally announced July 2020.

Journal ref: J. Chem. Theory Comput. 2021, 17, 2, 742-755

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations

Authors: Diptarka Hait, Eric A. Haugen, Zheyue Yang, Katherine J. Oosterbaan, Stephen R. Leone, Martin Head-Gordon

Abstract: State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted on account of the risk of `variational collapse' down to the ground state. We employ the recently developed square gradient minimization (SGM, J. Chem. Theory Comput. 16, 1699-1710, 2020) algorithm to reliably avoid variati… ▽ More State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted on account of the risk of `variational collapse' down to the ground state. We employ the recently developed square gradient minimization (SGM, J. Chem. Theory Comput. 16, 1699-1710, 2020) algorithm to reliably avoid variational collapse and study the effectiveness of orbital optimized density functional theory (DFT) at predicting second period element 1s core-level spectra of open-shell systems. Several density functionals (including SCAN, B3LYP and $ωおめが$B97X-D3) are found to predict excitation energies from the core to singly occupied levels to high accuracy ($\le 0.3$ eV RMS error), against available experimental data. Higher excited states are however more challenging by virtue of being intrinsically multiconfigurational. We thus present a CI inspired route to self-consistently recouple single determinant mixed configurations obtained from DFT, in order to obtain approximate doublet states. This recoupling scheme is used to predict the C K-edge spectra of the allyl radical, the O K-edge spectra of CO$^+$ and the N K-edge of NO$_2$ to high accuracy relative to experiment, indicating substantial promise in using this approach for computation of core-level spectra for doublet species (vs more traditional time dependent DFT, EOM-CCSD or using unrecoupled mixed configurations). We also present general guidelines for computing core-excited states from orbital optimized DFT. △ Less

Submitted 21 September, 2020; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: Added more data

Journal ref: J. Chem. Phys. 153 (13), 134108 (2020)

The Critical Role of Thermal Fluctuations for Electrocatalytic Metal Surface Properties and CO Binding Trends

Authors: Wan-Lu Li, Christianna N. Lininger, Valerie Vaissier Welborn, Elliot Rossomme, Alexis T. Bell, Martin Head-Gordon, Teresa Head-Gordon

Abstract: This work addresses a longstanding theoretical discrepancy using Density Functional Theory (DFT) with experimental observations of CO binding trends on electrocatalytically relevant metals for the CO2 reduction reaction (CO2RR). By introducing thermal fluctuations using appropriate statistical mechanical NVT and NPT ensembles, we show that DFT with universal dispersion interactions yields qualitat… ▽ More This work addresses a longstanding theoretical discrepancy using Density Functional Theory (DFT) with experimental observations of CO binding trends on electrocatalytically relevant metals for the CO2 reduction reaction (CO2RR). By introducing thermal fluctuations using appropriate statistical mechanical NVT and NPT ensembles, we show that DFT with universal dispersion interactions yields qualitatively better metal surface strain trends and CO binding energetics, consistently predicts the correct site preference for all metals due to thermally induced surface distortions that preferentially exposes the undercoordinated atop site for Cu(111) and Pt(111), and for the weak binding Ag(111) and Au(111) surfaces at finite temperatures shows CO-metal interactions that are a mixture of chemisorbed and physisorbed species. This study better places theory as an equal partner to experimental heterogeneous catalysis by demonstrating the need to fully account for finite temperature fluctuations to make contact with surface science experiments. △ Less

Submitted 11 November, 2020; v1 submitted 6 June, 2020; originally announced June 2020.

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Authors: Srimukh Prasad Veccham, Martin Head-Gordon

Abstract: Efficient and high capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density func… ▽ More Efficient and high capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a dataset (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this dataset using coupled cluster theory. Secondly, we have assessed the performance of 55 density functional approximations for predicting H$_2$ interaction energies and have identified two hybrid density functionals ($ωおめが$B97X-V and $ωおめが$B97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semi-local density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals like revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1kJ/mol. △ Less

Submitted 24 June, 2020; v1 submitted 25 March, 2020; originally announced March 2020.

Comments: 63 pages, 11 figures

CASSCF with Extremely Large Active Spaces using the Adaptive Sampling Configuration Interaction Method

Authors: Daniel S. Levine, Diptarka Hait, Norm M. Tubman, Susi Lehtola, K. Birgitta Whaley, Martin Head-Gordon

Abstract: The complete active space self-consistent field (CASSCF) method is the principal approach employed for studying strongly correlated systems. However, exact CASSCF can only be performed on small active spaces of ~20 electrons in ~20 orbitals due to exponential growth in the computational cost. We show that employing the Adaptive Sampling Configuration Interaction (ASCI) method as an approximate Ful… ▽ More The complete active space self-consistent field (CASSCF) method is the principal approach employed for studying strongly correlated systems. However, exact CASSCF can only be performed on small active spaces of ~20 electrons in ~20 orbitals due to exponential growth in the computational cost. We show that employing the Adaptive Sampling Configuration Interaction (ASCI) method as an approximate Full CI solver in the active space allows CASSCF-like calculations within chemical accuracy (<1 kcal/mol for relative energies) in active spaces with more than ~50 active electrons in ~50 active orbitals, significantly increasing the sizes of systems amenable to accurate multiconfigurational treatment. The main challenge with using any selected CI-based approximate CASSCF is the orbital optimization problem; they tend to exhibit large numbers of local minima in orbital space due to their lack of invariance to active-active rotations (in addition to the local minima that exist in exact CASSCF). We highlight methods that can avoid spurious local extrema as a practical solution to the orbital optimization problem. We employ ASCI-SCF to demonstrate lack of polyradical character in moderately sized periacenes with up to 52 correlated electrons and compare against heat-bath CI on an iron porphyrin system with more than 40 correlated electrons. △ Less

Submitted 4 February, 2020; v1 submitted 18 December, 2019; originally announced December 2019.

Comments: 44 pages, 5 figures

Journal ref: J. Chem. Theory Comput. 16, 2340 (2020)

Highly Accurate Prediction of Core Spectra of Molecules at Density Functional Theory Cost: Attaining sub eV Error from a Restricted Open-Shell Kohn-Sham Approach

Authors: Diptarka Hait, Martin Head-Gordon

Abstract: We present the use of the recently developed Square Gradient Minimization (SGM) algorithm for excited state orbital optimization, to obtain spin-pure Restricted Open-Shell Kohn-Sham (ROKS) energies for core excited states of molecules. The SGM algorithm is robust against variational collapse, and offers a reliable route to converging orbitals for target excited states at only 2-3 times the cost of… ▽ More We present the use of the recently developed Square Gradient Minimization (SGM) algorithm for excited state orbital optimization, to obtain spin-pure Restricted Open-Shell Kohn-Sham (ROKS) energies for core excited states of molecules. The SGM algorithm is robust against variational collapse, and offers a reliable route to converging orbitals for target excited states at only 2-3 times the cost of ground state orbital optimization (per iteration). ROKS/SGM with the modern SCAN/$ωおめが$B97X-V functionals is found to predict the K edge of C,N,O and F to a root mean squared error of $\sim$0.3 eV. ROKS/SGM is equally effective at predicting L edge spectra of third period elements, provided a perturbative spin-orbit correction is employed. This high accuracy can be contrasted with traditional TDDFT, which typically has greater than 10 eV error and requires translation of computed spectra to align with experiment. ROKS is computationally affordable (having the same scaling as ground state DFT, and a slightly larger prefactor) and can be applied to geometry optimizations/ab-initio molecular dynamics of core excited states, as well as condensed phase simulations. ROKS can also model doubly excited/ionized states with one broken electron pair, which are beyond the ability of linear response based methods. △ Less

Submitted 9 January, 2020; v1 submitted 11 December, 2019; originally announced December 2019.

Comments: Added spectra and L edge calculations

Journal ref: J. Phys. Chem. Lett. 2020, 11, 3, 775-786

Excited state orbital optimization via minimizing the square of the gradient: General approach and application to singly and doubly excited states via density functional theory

Authors: Diptarka Hait, Martin Head-Gordon

Abstract: We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting Square Gradient Minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite dif… ▽ More We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting Square Gradient Minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant $Δでるた$SCF and spin-purified Restricted Open-Shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the Maximum Overlap Method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that $Δでるた$SCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems. △ Less

Submitted 12 November, 2019; originally announced November 2019.

Journal ref: J. Chem. Theory Comput. 2020, 16, 3, 1699-1710

Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller-Plesset Perturbation Theory

Authors: Joonho Lee, Lin Lin, Martin Head-Gordon

Abstract: We present a systematically improvable tensor hypercontraction (THC) factorization based on interpolative separable density fitting (ISDF). We illustrate algorithmic details to achieve this within the framework of Becke's atom-centered quadrature grid. A single ISDF parameter $c_\text{ISDF}$ controls the tradeoff between accuracy and cost. In particular, $c_\text{ISDF}$ sets the number of interpol… ▽ More We present a systematically improvable tensor hypercontraction (THC) factorization based on interpolative separable density fitting (ISDF). We illustrate algorithmic details to achieve this within the framework of Becke's atom-centered quadrature grid. A single ISDF parameter $c_\text{ISDF}$ controls the tradeoff between accuracy and cost. In particular, $c_\text{ISDF}$ sets the number of interpolation points used in THC, $N_\text{IP} = c_\text{ISDF}\times N_\text{X}$ with $N_\text{X}$ being the number of auxiliary basis functions. In conjunction with the resolution-of-the-identity (RI) technique, we develop and investigate the THC-RI algorithms for cubic-scaling exact exchange for Hartree-Fock and range-separated hybrids (e.g., $ωおめが$B97X-V) and quartic-scaling second- and third-order Møller-Plesset theory (MP2 and MP3). These algorithms were evaluated over the W4-11 thermochemistry (atomization energy) set and A24 non-covalent interaction benchmark set with standard Dunning basis sets (cc-pVDZ, cc-pVTZ, aug-cc-pVDZ, and aug-cc-pVTZ). We demonstrate the convergence of THC-RI algorithms to numerically exact RI results using ISDF points. Based on these, we make recommendations on $c_\text{ISDF}$ for each basis set and method. We also demonstrate the utility of THC-RI exact exchange and MP2 for larger systems such as water clusters and $\text{C}_{20}$. We stress that more challenges await in obtaining accurate and numerically stable THC factorization for wavefunction amplitudes as well as the space spanned by virtual orbitals in large basis sets and implementing sparsity-aware THC-RI algorithms. △ Less

Submitted 1 November, 2019; originally announced November 2019.

Journal ref: J. Chem. Theory Comput. 2020, 16, 1, 243-263

The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients

Authors: Christopher J. Stein, John M. Herbert, Martin Head-Gordon

Abstract: We present the theory and implementation of a Poisson-Boltzmann implicit solvation model for electrolyte solutions. This model can be combined with arbitrary electronic structure methods that provide an accurate charge density of the solute. A hierarchy of approximations for this model includes a linear approximation for weak electrostatic potentials, finite size of the mobile electrolyte ions and… ▽ More We present the theory and implementation of a Poisson-Boltzmann implicit solvation model for electrolyte solutions. This model can be combined with arbitrary electronic structure methods that provide an accurate charge density of the solute. A hierarchy of approximations for this model includes a linear approximation for weak electrostatic potentials, finite size of the mobile electrolyte ions and a Stern-layer correction. Recasting the Poisson-Boltzmann equations into Euler-Lagrange equations then significantly simplifies the derivation of the free energy of solvation for these approximate models. The parameters of the model are then either fit directly to experimental observables, e.g. the finite ion size, or optimized for agreement with experimental results. Experimental data for this optimization is available in the form of Sechenov coefficients that describe the linear dependence of the salting-out effect of solutes with respect to the electrolyte concentration. In the final part we rationalize the qualitative disagreement of the finite ion size modification to the Poisson-Boltzmann model with experimental observations by taking into account the electrolyte concentration dependence of the Stern layer. A route towards a revised model that captures the experimental observations while including the finite ion size effects is then outlined. This implementation paves the way for the study of electrochemical and electrocatalytic processes of molecules and cluster models with accurate electronic structure methods. △ Less

Submitted 8 October, 2019; originally announced October 2019.

Comments: 21 pages, 5 figures, 3 tables