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

Excited States via Coupled Cluster Theory without Equation-of-Motion Methods: Seeking Higher Roots with Application to Doubly Excited States and Double Core Hole States

Authors: Joonho Lee, David W. Small, Martin Head-Gordon

Abstract: In this work, we revisited the idea of using the coupled-cluster ground state formalism to target excited states. Our main focus was targeting doubly excited states and double core hole states. Typical equation-of-motion (EOM) approaches for obtaining these states struggle without higher-order excitations than doubles. We showed that by using a non-aufbau determinant optimized via the maximum over… ▽ More In this work, we revisited the idea of using the coupled-cluster ground state formalism to target excited states. Our main focus was targeting doubly excited states and double core hole states. Typical equation-of-motion (EOM) approaches for obtaining these states struggle without higher-order excitations than doubles. We showed that by using a non-aufbau determinant optimized via the maximum overlap method the CC ground state solver can target higher energy states. Furthermore, just with singles and doubles (i.e., CCSD), we demonstrated that the accuracy of $Δでるた$CCSD and $Δでるた$CCSD(T) far surpasses that of EOM-CCSD for doubly excited states. The accuracy of $Δでるた$CCSD(T) is nearly exact for doubly excited states considered in this work. For double core hole states, we used an improved ansatz for greater numerical stability by freezing core hole orbitals. The improved methods, core valence separation (CVS)-$Δでるた$CCSD and CVS-$Δでるた$CCSD(T), were applied to the calculation of the double ionization potential of small molecules. Even without relativistic corrections, we observed qualitatively accurate results with CVS-$Δでるた$CCSD and CVS-$Δでるた$CCSD(T). Remaining challenges in $Δでるた$CC include the description of open-shell singlet excited states with the single-reference CC ground state formalism as well as excited states with genuine multi-reference character. The tools and intuition developed in this work may serve as a stepping stone towards directly targeting arbitrary excited states using ground state CC methods. △ Less

Submitted 22 September, 2019; originally announced September 2019.

Journal ref: J. Chem. Phys. 151, 214103 (2019)

Making Many-Body Interactions Nearly Pairwise Additive: The Polarized Many-Body Expansion Approach

Authors: Srimukh Prasad Veccham, Joonho Lee, Martin Head-Gordon

Abstract: The Many-Body Expansion (MBE) is a useful tool to simulate condensed phase chemical systems, often avoiding the steep computational cost of usual electronic structure methods. However, it often requires higher than 2-body terms to achieve quantitative accuracy. In this work, we propose the Polarized MBE (PolBE) method where each MBE energy contribution is treated as an embedding problem. In each e… ▽ More The Many-Body Expansion (MBE) is a useful tool to simulate condensed phase chemical systems, often avoiding the steep computational cost of usual electronic structure methods. However, it often requires higher than 2-body terms to achieve quantitative accuracy. In this work, we propose the Polarized MBE (PolBE) method where each MBE energy contribution is treated as an embedding problem. In each energy term, a smaller fragment is embedded into a larger, polarized environment and only a small region is treated at the high-level of theory using embedded mean-field theory. The role of polarized environment was found to be crucial in providing quantitative accuracy at the 2-body level. PolBE accurately predicts non-covalent interaction energies for a number of systems, including CO$_2$, water, and hydrated ion clusters, with a variety of interaction mechanisms, from weak dispersion to strong electrostatics considered in this work. We further demonstrate that the PolBE interaction energy is predominantly pairwise unlike the usual vacuum MBE which requires higher-order terms to achieve similar accuracy. We numerically show that PolBE often performs better than other widely used embedded MBE methods such as the electrostatically embedded MBE. Owing to the lack of expensive diagonalization of Fock matrices and its embarrassingly parallel nature, PolBE is a promising way to access condensed phase systems with hybrid density functionals that are difficult to treat with currently available methods. △ Less

Submitted 18 October, 2019; v1 submitted 4 September, 2019; originally announced September 2019.

Journal ref: J. Chem. Phys. 151, 194101 (2019)

Beyond the Coulson-Fischer point: Characterizing single excitation CI and TDDFT for excited states in single bond dissociations

Authors: Diptarka Hait, Adam Rettig, Martin Head-Gordon

Abstract: Linear response time dependent density functional theory (TDDFT), which builds upon configuration interaction singles (CIS) and TD-Hartree-Fock (TDHF), is the most widely used class of excited state quantum chemistry methods and is often employed to study photochemical processes. This paper studies the behavior of the resulting excited state potential energy surfaces beyond the Coulson-Fisher (CF)… ▽ More Linear response time dependent density functional theory (TDDFT), which builds upon configuration interaction singles (CIS) and TD-Hartree-Fock (TDHF), is the most widely used class of excited state quantum chemistry methods and is often employed to study photochemical processes. This paper studies the behavior of the resulting excited state potential energy surfaces beyond the Coulson-Fisher (CF) point in single bond dissociations, when the optimal reference determinant is spin-polarized. Many excited states exhibit sharp kinks at the CF point, and connect to different dissociation limits via a zone of unphysical concave curvature. In particular, the unrestricted M$_S=0$ lowest triplet T$_1$ state changes character, and does not dissociate into ground state fragments. The unrestricted $M_S=\pm 1$ T$_1$ CIS states better approximate the physical dissociation limit, but their degeneracy is broken beyond the CF point for most single bond dissociations. On the other hand, the $M_S=\pm 1$ T$_1$ TDHF states reach the asymptote too soon, by merging with the ground state from the CF point onwards. Use of local exchange-correlation functionals causes $M_S=\pm 1$ T$_1$ TDDFT states to resemble their unphysical $M_S= 0$ counterpart. The 2 orbital, 2-electron model system of minimal basis H$_2$ is analytically treated to understand the origin of these issues, revealing that the lack of double excitations is at the root of these remarkable observations. The behavior of excited state surfaces is also numerically examined for species like H$_2$, NH$_3$, C$_2$H$_6$ and LiH in extended basis sets. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Journal ref: Phys.Chem. Chem. Phys. 21, 21761-21775. 2019

Third-Order Møller-Plesset Perturbation Theory Made Useful? Choice of Orbitals and Scaling Greatly Improves Accuracy for Thermochemistry, Kinetics and Intermolecular Interactions

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

Abstract: We develop and test methods that include second and third-order perturbation theory (MP3) using orbitals obtained from regularized orbital-optimized second-order perturbation theory, $κかっぱ$-OOMP2, denoted as MP3:$κかっぱ$-OOMP2. Testing MP3:$κかっぱ$-OOMP2 shows RMS errors that are 1.7 to 5 times smaller than MP3 across 7 data sets. To do still better, empirical training of the scaling factors for the second- an… ▽ More We develop and test methods that include second and third-order perturbation theory (MP3) using orbitals obtained from regularized orbital-optimized second-order perturbation theory, $κかっぱ$-OOMP2, denoted as MP3:$κかっぱ$-OOMP2. Testing MP3:$κかっぱ$-OOMP2 shows RMS errors that are 1.7 to 5 times smaller than MP3 across 7 data sets. To do still better, empirical training of the scaling factors for the second- and third-order correlation energies and the regularization parameter on one of those data sets led to an unregularized scaled ($c_2=1.0$; $c_3=0.8$) denoted as MP2.8:$κかっぱ$-OOMP2. MP2.8:$κかっぱ$-OOMP2 yields significant additional improvement over MP3:$κかっぱ$-OOMP2 in 4 of 6 test data sets on thermochemistry, kinetics, and noncovalent interactions. Remarkably, these two methods outperform coupled cluster with singles and doubles in 5 of the 7 data sets considered, at greatly reduced cost (no $\mathcal{O}(N^6)$ iterations). △ Less

Submitted 12 October, 2019; v1 submitted 6 June, 2019; originally announced June 2019.

Journal ref: J. Phys. Chem. Lett. 2019, 10, 15, 4170-4176

A Multi-Resolution 3D-DenseNet for Chemical Shift Prediction in NMR Crystallography

Authors: Shuai Liu, Jie Li, Kochise C. Bennett, Brad Ganoe, Tim Stauch, Martin Head-Gordon, Alexander Hexemer, Daniela Ushizima, Teresa Head-Gordon

Abstract: We have developed a deep learning algorithm for chemical shift prediction for atoms in molecular crystals that utilizes an atom-centered Gaussian density model for the 3D data representation of a molecule. We define multiple channels that describe different spatial resolutions for each atom type that utilizes cropping, pooling, and concatenation to create a multi-resolution 3D-DenseNet architectur… ▽ More We have developed a deep learning algorithm for chemical shift prediction for atoms in molecular crystals that utilizes an atom-centered Gaussian density model for the 3D data representation of a molecule. We define multiple channels that describe different spatial resolutions for each atom type that utilizes cropping, pooling, and concatenation to create a multi-resolution 3D-DenseNet architecture (MR-3D-DenseNet). Because the training and testing time scale linearly with the number of samples, the MR-3D-DenseNet can exploit data augmentation that takes into account the property of rotational invariance of the chemical shifts, thereby also increasing the size of the training dataset by an order of magnitude without additional cost. We obtain very good agreement for 13C, 15N, and 17O chemical shifts, with the highest accuracy found for 1H chemical shifts that is equivalent to the best predictions using ab initio quantum chemistry methods. △ Less

Submitted 31 May, 2019; originally announced June 2019.

Development of an Advanced Force Field for Water using Variational Energy Decomposition Analysis

Authors: A. K. Das, L. Urban, I. Leven, M. Loipersberger, A. Aldossary, M. Head-Gordon, T. Head-Gordon

Abstract: Given the piecewise approach to modeling intermolecular interactions for force fields, they can be difficult to parameterize since they are fit to data like total energies that only indirectly connect to their separable functional forms. Furthermore, by neglecting certain types of molecular interactions such as charge penetration and charge transfer, most classical force fields must rely on, but d… ▽ More Given the piecewise approach to modeling intermolecular interactions for force fields, they can be difficult to parameterize since they are fit to data like total energies that only indirectly connect to their separable functional forms. Furthermore, by neglecting certain types of molecular interactions such as charge penetration and charge transfer, most classical force fields must rely on, but do not always demonstrate, how cancellation of errors occurs among the remaining molecular interactions accounted for such as exchange repulsion, electrostatics, and polarization. In this work we present the first generation of the (many-body) MB-UCB force field that explicitly accounts for the decomposed molecular interactions commensurate with a variational energy decomposition analysis, including charge transfer, with force field design choices that reduce the computational expense of the MB-UCB potential while remaining accurate. We optimize parameters using only single water molecule and water cluster data up through pentamers, with no fitting to condensed phase data, and we demonstrate that high accuracy is maintained when the force field is subsequently validated against conformational energies of larger water cluster data sets, radial distribution functions of the liquid phase, and the temperature dependence of thermodynamic and transport water properties. We conclude that MB-UCB is comparable in performance to MB-Pol, but is less expensive and more transferable by eliminating the need to represent short-ranged interactions through large parameter fits to high order polynomials. △ Less

Submitted 27 July, 2019; v1 submitted 19 May, 2019; originally announced May 2019.

Kohn-Sham Density Functional Theory with Complex, Spin-Restricted Orbitals: Accessing a New Class of Densities without the Symmetry Dilemma

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

Abstract: We show that using complex, spin-restricted orbitals (cR) in Kohn-Sham density functional theory (KS-DFT) allows one to access a new class of densities that is not accessible by either spin-restricted (RKS) or spin-unrestricted (UKS) orbitals. We further show that the real part of a cRKS density matrix can be non-idempotent when the imaginary part of the density matrix is not zero. Using cRKS orbi… ▽ More We show that using complex, spin-restricted orbitals (cR) in Kohn-Sham density functional theory (KS-DFT) allows one to access a new class of densities that is not accessible by either spin-restricted (RKS) or spin-unrestricted (UKS) orbitals. We further show that the real part of a cRKS density matrix can be non-idempotent when the imaginary part of the density matrix is not zero. Using cRKS orbitals shows significant improvements in the triplet-singlet gaps of the TS12 benchmark set for the SPW92, PBE, BLYP, and TPSS functionals. Moreover, it was shown that RKS and UKS yield qualitatively wrong charge densities and spin densities, respectively, leading to worse energetics. We demonstrate that the modern functionals SCAN, MN15-L, and B97M-V show surprisingly no improvement even with a qualitatively more accurate density from cRKS orbitals. To this end, our work not only paves the way to escape the symmetry dilemma whenever there exists a cRKS solution, but also suggests a new route to design better approximate XC functionals. △ Less

Submitted 17 April, 2019; originally announced April 2019.

Journal ref: Phys. Rev. Lett. 123, 113001 (2019)

Two Single-Reference Approaches to Singlet Biradicaloid Problems: Complex, Restricted Orbitals and Approximate Spin-Projection Combined With Regularized Orbital-Optimized Møller-Plesset Perturbation Theory

Authors: Joonho Lee, Martin Head-Gordon

Abstract: We present a comprehensive study of two single-reference approaches to singlet biradicaloids. These two approaches are based on the recently developed regularized orbital-optimized Møller-Plesset method ($κかっぱ$-OOMP2). The first approach is to combine the Yamaguchi's approximate projection (AP) scheme and $κかっぱ$-OOMP2 with unrestricted (U) orbitals ($κかっぱ$-UOOMP2). By capturing only essential symmetry brea… ▽ More We present a comprehensive study of two single-reference approaches to singlet biradicaloids. These two approaches are based on the recently developed regularized orbital-optimized Møller-Plesset method ($κかっぱ$-OOMP2). The first approach is to combine the Yamaguchi's approximate projection (AP) scheme and $κかっぱ$-OOMP2 with unrestricted (U) orbitals ($κかっぱ$-UOOMP2). By capturing only essential symmetry breaking, $κかっぱ$-UOOMP2 can serve as a suitable basis for AP. The second approach is $κかっぱ$-OOMP2 with complex, restricted (cR) orbitals ($κかっぱ$-cROOMP2). Though its applicability is more limited due to the comparative rarity of cR solutions, $κかっぱ$-cROOMP2 offers a simple framework for describing singlet biradicaloids with complex polarization while removing artificial spatial symmetry breaking. We compare the scope of these two methods with numerical studies. We show that AP+$κかっぱ$-UOOMP2 and $κかっぱ$-cROOMP2 can perform similarly well in the TS12 set, a data set that includes 12 data points for triplet-singlet gaps of several atoms and diatomic molecules with a triplet ground state. This was also found to be true for the barrier height of a reaction involving attack on a cysteine ion by a singlet oxygen molecule. However, we also demonstrate that in highly symmetric systems like $\text{C}_{30}$ ($\text{D}_{5h}$) $κかっぱ$-cROOMP2 is more suitable as it conserves spatial symmetry. Lastly, we present an organic biradicaloid that does not have a $κかっぱ$-cROOMP2 solution in which case only AP+$κかっぱ$-UOOMP2 is applicable. We recommend $κかっぱ$-cROOMP2 whenever complex polarization is essential and AP+$κかっぱ$-UOOMP2 for biradicaloids without essential complex polarization but with essential spin-polarization. △ Less

Submitted 27 June, 2019; v1 submitted 26 March, 2019; originally announced March 2019.

Journal ref: J. Chem. Phys. 150, 244106 (2019)

Distinguishing Artificial and Essential Symmetry Breaking in a Single Determinant: Approach and Application to the C$_{60}$, C$_{36}$, and C$_{20}$ Fullerenes

Authors: Joonho Lee, Martin Head-Gordon

Abstract: We present a thorough analysis of symmetry breaking observed in Hartree-Fock (HF) solutions of fullerenes C$_{60}$, C$_{36}$, and C$_{20}$ in order to characterize the nature of electron correlation in them. Our analysis is based on (1) the critical regularization strength to restore symmetry breaking in the recently developed regularized orbital optimized second-order Møller-Plesset perturbation… ▽ More We present a thorough analysis of symmetry breaking observed in Hartree-Fock (HF) solutions of fullerenes C$_{60}$, C$_{36}$, and C$_{20}$ in order to characterize the nature of electron correlation in them. Our analysis is based on (1) the critical regularization strength to restore symmetry breaking in the recently developed regularized orbital optimized second-order Møller-Plesset perturbation theory ($κかっぱ$-OOMP2), (2) singlet-triplet gaps from various MP2 methods, and (3) natural orbital occupation numbers from restricted coupled-cluster with singles and doubles (RCCSD) and coupled-cluster valence bond with singles and doubles (CCVB-SD). Based on these three independent probes, we conclude that C$_{36}$ (D$_\text{6h}$) exhibits genuine strong correlation and symmetry breaking whereas C$_{60}$ exhibits {\it artificial} HF symmetry breaking and is not strongly correlated. Investigating the critical regularization strength, we discuss strong correlation in C$_{20}$ at the Jahn-Teller distorted geometries (C$_\text{2h}$, D$_\text{2h}$, C$_\text{i}$, and D$_\text{3h}$) and the I$_\text{h}$ geometry. Only C$_{20}$ (I$_\text{h}$) was found to be strongly correlated while others exhibit {\it artificial} HF symmetry breaking. This analysis highlights a useful feature of the recommended $κかっぱ$-OOMP2 method. It is an electronic structure method that describes dynamic correlation, and attenuates strong correlation in MP2 towards zero by regularization. Therefore, $κかっぱ$-OOMP2 will exhibit symmetry breaking in its reference determinant only when correlation is strong (i.e., essential symmetry breaking). Artificial symmetry breaking (arising in HF due to neglect of dynamic correlation) is thus removed in $κかっぱ$-OOMP2. △ Less

Submitted 13 December, 2018; originally announced December 2018.

Well-behaved versus ill-behaved density functionals for single bond dissociation: Separating success from disaster functional by functional for stretched H$_2$

Authors: Diptarka Hait, Adam Rettig, Martin Head-Gordon

Abstract: Unrestricted DFT methods are typically expected to describe the homolytic dissociation of nonpolar single bonds in neutral species with qualitative accuracy, due to the lack of significant delocalization error. We however find that many widely used density functional approximations fail to describe features along the dissociation curve of the simple H$_2$ molecule. This is not an universal failure… ▽ More Unrestricted DFT methods are typically expected to describe the homolytic dissociation of nonpolar single bonds in neutral species with qualitative accuracy, due to the lack of significant delocalization error. We however find that many widely used density functional approximations fail to describe features along the dissociation curve of the simple H$_2$ molecule. This is not an universal failure of DFT in the sense that many classic functionals like PBE and B3LYP give very reasonable results, as do some more modern methods like MS2. However, some other widely used functionals like B97-D (empirically fitted) and TPSS (non-empirically constrained) predict qualitatively wrong static polarizabilities, force constants and some even introduce an artificial barrier against association of independent H atoms to form H$_2$. The polarizability and force constant prediction failures appear to stem from incomplete spin localization into individual H atoms beyond the Coulson-Fisher point, resulting in `fractionally bonded' species where the ionic contributions to the Slater determinant are not completely eliminated, unlike the case of unrestricted Hartree-Fock. These errors therefore appear to be a consequence of poor self-consistent density prediction by the problematic functional. The same reasons could potentially lead to spurious barriers towards H atom association, indirectly also leading to incorrect forces. These unphysicalities suggest that the use of problematic functionals is probably unwise in \textit{ab initio} dynamics calculations, especially if strong electrostatic interactions are possible. △ Less

Submitted 22 October, 2018; originally announced October 2018.

Journal ref: J. Chem. Phys. 150, 094115 (2019)

Generalized Unitary Coupled Cluster Wavefunctions for Quantum Computation

Authors: Joonho Lee, William J. Huggins, Martin Head-Gordon, K. Birgitta Whaley

Abstract: We introduce a unitary coupled-cluster (UCC) ansatz termed $k$-UpCCGSD that is based on a family of sparse generalized doubles (D) operators which provides an affordable and systematically improvable unitary coupled-cluster wavefunction suitable for implementation on a near-term quantum computer. $k$-UpCCGSD employs $k$ products of the exponential of pair coupled-cluster double excitation operator… ▽ More We introduce a unitary coupled-cluster (UCC) ansatz termed $k$-UpCCGSD that is based on a family of sparse generalized doubles (D) operators which provides an affordable and systematically improvable unitary coupled-cluster wavefunction suitable for implementation on a near-term quantum computer. $k$-UpCCGSD employs $k$ products of the exponential of pair coupled-cluster double excitation operators (pCCD), together with generalized single (S) excitation operators. We compare its performance in both efficiency of implementation and accuracy with that of the generalized UCC ansatz employing the full generalized SD excitation operators (UCCGSD), as well as with the standard ansatz employing only SD excitations (UCCSD). $k$-UpCCGSD is found to show the best scaling for quantum computing applications, requiring a circuit depth of $\mathcal O(kN)$, compared with $\mathcal O(N^3)$ for UCCGSD and $\mathcal O((N-ηいーた)^2 ηいーた)$ for UCCSD where $N$ is the number of spin orbitals and $ηいーた$ is the number of electrons. We analyzed the accuracy of these three ansätze by making classical benchmark calculations on the ground state and the first excited state of H$_4$ (STO-3G, 6-31G), H$_2$O (STO-3G), and N$_2$ (STO-3G), making additional comparisons to conventional coupled cluster methods. The results for ground states show that $k$-UpCCGSD offers a good tradeoff between accuracy and cost, achieving chemical accuracy for lower cost of implementation on quantum computers than both UCCGSD and UCCSD. Excited states are calculated with an orthogonally constrained variational quantum eigensolver approach. This is seen to generally yield less accurate energies than for the corresponding ground states. We demonstrate that using a specialized multi-determinantal reference state constructed from classical linear response calculations allows these excited state energetics to be improved. △ Less

Submitted 3 January, 2019; v1 submitted 4 October, 2018; originally announced October 2018.

Postponing the orthogonality catastrophe: efficient state preparation for electronic structure simulations on quantum devices

Authors: Norm M. Tubman, Carlos Mejuto-Zaera, Jeffrey M. Epstein, Diptarka Hait, Daniel S. Levine, William Huggins, Zhang Jiang, Jarrod R. McClean, Ryan Babbush, Martin Head-Gordon, K. Birgitta Whaley

Abstract: Despite significant work on resource estimation for quantum simulation of electronic systems, the challenge of preparing states with sufficient ground state support has so far been largely neglected. In this work we investigate this issue in several systems of interest, including organic molecules, transition metal complexes, the uniform electron gas, Hubbard models, and quantum impurity models ar… ▽ More Despite significant work on resource estimation for quantum simulation of electronic systems, the challenge of preparing states with sufficient ground state support has so far been largely neglected. In this work we investigate this issue in several systems of interest, including organic molecules, transition metal complexes, the uniform electron gas, Hubbard models, and quantum impurity models arising from embedding formalisms such as dynamical mean-field theory. Our approach uses a state-of-the-art classical technique for high-fidelity ground state approximation. We find that easy-to-prepare single Slater determinants such as the Hartree-Fock state often have surprisingly robust support on the ground state for many applications of interest. For the most difficult systems, single-determinant reference states may be insufficient, but low-complexity reference states may suffice. For this we introduce a method for preparation of multi-determinant states on quantum computers. △ Less

Submitted 14 September, 2018; originally announced September 2018.

Comments: 8 pages + SI, 5 figures

Open-Shell Coupled-Cluster Valence-Bond Theory Augmented with an Independent Amplitude Approximation for Three-Pair Correlations: Application to a Model Oxygen-Evolving Complex and Single Molecular Magnet

Authors: Joonho Lee, David W. Small, Martin Head-Gordon

Abstract: We report the failure of coupled-cluster valence-bond (CCVB) theory with two-pair configurations [J. Chem. Phys. 2009, 130, 084103 (2009)] for open-shell (OS) spin-frustrated systems where including three-pair configurations is necessary to properly describe strong spin-correlations. We extend OS-CCVB by augmenting the model with three-pair configurations within the independent amplitude approxima… ▽ More We report the failure of coupled-cluster valence-bond (CCVB) theory with two-pair configurations [J. Chem. Phys. 2009, 130, 084103 (2009)] for open-shell (OS) spin-frustrated systems where including three-pair configurations is necessary to properly describe strong spin-correlations. We extend OS-CCVB by augmenting the model with three-pair configurations within the independent amplitude approximation (IAA). The resulting new electronic structure model, OS-CCVB+i3, involves only a quadratic number of independent wavefunction parameters. It includes the recently reported closed-shell CCVB+i3 as a special case. Its cost is dominated by integral transformations and it is capable of breaking multiple bonds exactly for all systems examined so far. The strength of OS-CCVB+i3 is highlighted in realistic systems including the [CaMn$_3$O$_4$] cubane subunit of the oxygen-evolving complex and a molecular magnet with the [Cr$_9$] core unit as well as model systems such as N$_3$, V$_3$O$_3$, and P$_5$. We show that OS-CCVB+i3 is only slightly dependent on the underlying perfect-pairing reference while OS-CCVB shows a stronger dependence. We also emphasize the compactness of the OS-CCVB+i3 wavefunction compared to the heat-bath configuration interaction wavefunction, a recently introduced soft exponential-scaling approach. △ Less

Submitted 31 December, 2018; v1 submitted 20 August, 2018; originally announced August 2018.

Delocalization errors in density functional theory are essentially quadratic in fractional occupation number

Authors: Diptarka Hait, Martin Head-Gordon

Abstract: Approximate functionals used in practical density functional theory (DFT) deviate from the piecewise linear behavior of the exact functional for fractional charges. This deviation causes excess charge delocalization, which leads to incorrect densities, molecular properties, barrier heights, band gaps and excitation energies. We present a simple delocalization function for characterizing this error… ▽ More Approximate functionals used in practical density functional theory (DFT) deviate from the piecewise linear behavior of the exact functional for fractional charges. This deviation causes excess charge delocalization, which leads to incorrect densities, molecular properties, barrier heights, band gaps and excitation energies. We present a simple delocalization function for characterizing this error and find it to be almost perfectly linear vs the fractional electron number for systems spanning in size from the H atom to the C$_{12}$H$_{14}$ polyene. This causes the delocalization energy error to be a quadratic polynomial in the fractional electron number, which permits us to assess the comparative performance of 47 popular and recent functionals through the curvature. The quadratic form further suggests that information about a single fractional charge is sufficient to eliminate the principal source of delocalization error. Generalizing traditional two-point information like ionization potentials or electron affinities to account for a third, fractional charge based data point could therefore permit fitting/tuning of functionals with lower delocalization error. △ Less

Submitted 13 September, 2018; v1 submitted 7 August, 2018; originally announced August 2018.

Comments: Discussion about fractional binding issues in anions have been added, with other minor fixes/elaborations

Journal ref: J. Phys. Chem. Lett. 2018, 9, 6280-6288

An efficient deterministic perturbation theory for selected configuration interaction methods

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

Abstract: The interplay between advances in stochastic and deterministic algorithms has recently led to development of interesting new selected configuration interaction (SCI) methods for solving the many body Schrödinger equation. The performance of these SCI methods can be greatly improved with a second order perturbation theory (PT2) correction, which is often evaluated in a stochastic or hybrid-stochast… ▽ More The interplay between advances in stochastic and deterministic algorithms has recently led to development of interesting new selected configuration interaction (SCI) methods for solving the many body Schrödinger equation. The performance of these SCI methods can be greatly improved with a second order perturbation theory (PT2) correction, which is often evaluated in a stochastic or hybrid-stochastic manner. In this work, we present a highly efficient, fully deterministic PT2 algorithm for SCI methods and demonstrate that our approach is orders of magnitude faster than recent proposals for stochastic SCI+PT2. We also show that it is important to have a compact reference SCI wave function, in order to obtain optimal SCI+PT2 energies. This indicates that it advantageous to use accurate search algorithms such as 'ASCI search' rather than more approximate approaches. Our deterministic PT2 algorithm is based on sorting techniques that have been developed for modern computing architectures and is inherently straightforward to use on parallel computing architectures. Related architectures such as GPU implementations can be also used to further increase the efficiency. Overall, we demonstrate that the algorithms presented in this work allow for efficient evaluation of trillions of PT2 contributions with modest computing resources. △ Less

Submitted 6 August, 2018; originally announced August 2018.

Comments: 13 pages, 2 figures

Regularized Orbital-Optimized Second-Order Møller-Plesset Perturbation Theory: A Reliable Fifth-Order Scaling Electron Correlation Model with Orbital Energy Dependent Regularizers

Authors: Joonho Lee, Martin Head-Gordon

Abstract: We derive and assess two new classes of regularizers that cope with offending denominators in the single-reference second-order Møller-Plesset perturbation theory (MP2). In particular, we discuss the use of two types of orbital energy dependent regularizers, $κかっぱ$ and $σしぐま$, in conjunction with orbital-optimized MP2 (OOMP2). The resulting fifth-order scaling methods, $κかっぱ$-OOMP2 and $σしぐま$-OOMP2, have been… ▽ More We derive and assess two new classes of regularizers that cope with offending denominators in the single-reference second-order Møller-Plesset perturbation theory (MP2). In particular, we discuss the use of two types of orbital energy dependent regularizers, $κかっぱ$ and $σしぐま$, in conjunction with orbital-optimized MP2 (OOMP2). The resulting fifth-order scaling methods, $κかっぱ$-OOMP2 and $σしぐま$-OOMP2, have been examined for bond-breaking, thermochemistry, and biradical problems. Both methods with strong enough regularization restore restricted to unrestricted instability (i.e. Coulson-Fischer points) that unregularized OOMP2 lacks when breaking bonds in $\text{H}_{2}$, $\text{C}_2\text{H}_6$, $\text{C}_2\text{H}_4$, and $\text{C}_2\text{H}_2$. The training of the $κかっぱ$ and $σしぐま$ regularization parameters was performed with the W4-11 set. We further developed scaled correlation energy variants, $κかっぱ$-S-OOMP2 and $σしぐま$-S-OOMP2, by training on the TAE140 subset of the W4-11 set. Those new OOMP2 methods were tested on the RSE43 set and the TA14 set where unmodified OOMP2 itself performs very well. The modifications we made were found insignificant in these data sets. Furthermore, we tested the new OOMP2 methods on singlet biradicaloids using Yamaguchi's approximate spin-projection. Unlike the unregularized OOMP2, which fails to converge these systems due to the singularity, we show that regularized OOMP2 methods successfully capture strong biradicaloid characters. While further assessment on larger datasets is desirable, $κかっぱ$-OOMP2 with $κかっぱ$ = 1.45 $E_{h}^{-1}$ appears to combine favorable recovery of Coulson-Fischer points with good numerical performance. △ Less

Submitted 25 August, 2018; v1 submitted 16 July, 2018; originally announced July 2018.

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method

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

Abstract: Recent advances in selected CI, including the adaptive sampling configuration interaction (ASCI) algorithm and its heat bath extension, have made the ASCI approach competitive with the most accurate techniques available, and hence an increasingly powerful tool in solving quantum Hamiltonians. In this work, we show that a useful paradigm for generating efficient selected CI/exact diagonalization al… ▽ More Recent advances in selected CI, including the adaptive sampling configuration interaction (ASCI) algorithm and its heat bath extension, have made the ASCI approach competitive with the most accurate techniques available, and hence an increasingly powerful tool in solving quantum Hamiltonians. In this work, we show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques we present ASCI results applied to a large range of systems and basis sets in order to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-zeta benchmark data for the G1 dataset. The largest of these calculations is Si$_{2}$H$_{6}$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets. △ Less

Submitted 28 December, 2019; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: 22 pages,8 figures, 15 tables (added supplemental information on Cr2 in the svp basis)

Journal ref: J. Chem. Theory Comput. 2020, 16, 4, 2139-2159

xDH double hybrid functionals can be qualitatively incorrect for non-equilibrium geometries: Dipole moment inversion and barriers to radical-radical association using XYG3 and XYGJ-OS

Authors: Diptarka Hait, Martin Head-Gordon

Abstract: Double hybrid (DH) density functionals are amongst the most accurate density functional approximations developed so far, largely due to incorporation of correlation effects from unoccupied orbitals via second order perturbation theory (PT2). The xDH family of DH functionals calculate energy directly from orbitals optimized by a lower level approach like B3LYP, without self-consistent optimization.… ▽ More Double hybrid (DH) density functionals are amongst the most accurate density functional approximations developed so far, largely due to incorporation of correlation effects from unoccupied orbitals via second order perturbation theory (PT2). The xDH family of DH functionals calculate energy directly from orbitals optimized by a lower level approach like B3LYP, without self-consistent optimization. XYG3 and XYGJ-OS are two widely used xDH functionals that are known to be quite accurate at equilibrium geometries. Here, we show that the XYG3 and XYGJ-OS functionals can be ill behaved for stretched bonds well beyond the Coulson-Fischer point, predicting unphysical dipole moments and humps in potential energy curves for some simple systems like the HF molecule. Numerical experiments and analysis show these failures are not due to PT2. Instead, a large mismatch at stretched bond-lengths between the reference B3LYP orbitals and the optimized orbitals associated with the non-PT2 part of XYG3 lead to an unphysically large non-Hellman-Feynman contribution to first order properties like forces and electron densities. △ Less

Submitted 27 March, 2018; originally announced March 2018.

Journal ref: Journal of Chemical Physics 2018 148 (17), 171102

Reaction Mechanism of the Selective Reduction of CO$_2$ to CO by a Tetraaza [Co$^\text{II}$N$_4$H]$^{2+}$ Complex in the Presence of Protons

Authors: Alejandro J. Garza, Srimanta Pakhira, Alexis T. Bell, Jose L. Mendoza-Cortes, Martin Head-Gordon

Abstract: The tetraaza [Co$^\text{II}$N$_4$H]$^{2+}$ complex (\textbf{1}) is remarkable for its ability to selectively reduce CO$_2$ to CO with 45\% Faradaic efficiency and a CO to H$_2$ ratio of 3:2. We employ density functional theory (DFT) to determine the reasons behind the unusual catalytic properties of \textbf{1} and the most likely mechanism for CO$_2$ reduction. The selectivity for CO$_2$ over prot… ▽ More The tetraaza [Co$^\text{II}$N$_4$H]$^{2+}$ complex (\textbf{1}) is remarkable for its ability to selectively reduce CO$_2$ to CO with 45\% Faradaic efficiency and a CO to H$_2$ ratio of 3:2. We employ density functional theory (DFT) to determine the reasons behind the unusual catalytic properties of \textbf{1} and the most likely mechanism for CO$_2$ reduction. The selectivity for CO$_2$ over proton reduction is explained by analyzing the catalyst's affinity for the possible ligands present under typical reaction conditions: acetonitrile, water, CO$_2$, and bicarbonate. After reduction of the catalyst by two electrons, formation of [Co$^\text{I}$N$_4$H]$^{+}$-CO$_{2}^{-}$ is strongly favored. Based on thermodynamic and kinetic data, we establish that the only likely route for producing CO from here consists of a protonation step to yield [Co$^\text{I}$N$_4$H]$^{+}$-CO$_{2}$H, followed by reaction with CO$_2$ to form [Co$^\text{II}$N$_4$H]$^{2+}$-CO and bicarbonate. This conclusion corroborates the idea of a direct role of CO$_2$ as a Lewis acid to assist in {C-O} bond dissociation, a conjecture put forward by other authors to explain recent experimental observations. The pathway to formic acid is predicted to be forbidden by high activation barriers, in accordance with the products that are known to be generated by \textbf{1}. Calculated physical observables such as standard reduction potentials and the turnover frequency for our proposed catalytic cycle are in agreement with available experimental data reported in the literature. The mechanism also makes a prediction that may be experimentally verified: that the rate of CO formation should increase linearly with the partial pressure of CO$_2$. △ Less

Submitted 27 March, 2018; originally announced March 2018.

Comments: 8 pages, 2 figures

How accurate is density functional theory at predicting dipole moments? An assessment using a new database of 200 benchmark values

Authors: Diptarka Hait, Martin Head-Gordon

Abstract: Dipole moments are a simple, global measure of the accuracy of the electron density of a polar molecule. Dipole moments also affect the interactions of a molecule with other molecules as well as electric fields. To directly assess the accuracy of modern density functionals for calculating dipole moments, we have developed a database of 200 benchmark dipole moments, using coupled cluster theory thr… ▽ More Dipole moments are a simple, global measure of the accuracy of the electron density of a polar molecule. Dipole moments also affect the interactions of a molecule with other molecules as well as electric fields. To directly assess the accuracy of modern density functionals for calculating dipole moments, we have developed a database of 200 benchmark dipole moments, using coupled cluster theory through triple excitations, extrapolated to the complete basis set limit. This new database is used to assess the performance of 88 popular or recently developed density functionals. The results suggest that double hybrid functionals perform the best, yielding dipole moments within about 3.6-4.5% regularized RMS error versus the reference values---which is not very different from the 4% regularized RMS error produced by coupled cluster singles and doubles. Many hybrid functionals also perform quite well, generating regularized RMS errors in the 5-6% range. Some functionals however exhibit large outliers and local functionals in general perform less well than hybrids or double hybrids. △ Less

Submitted 15 February, 2018; v1 submitted 15 September, 2017; originally announced September 2017.

Comments: Added several double hybrid functionals, most of which turned out to be better than any functional from Rungs 1-4 of Jacob's ladder and are actually competitive with CCSD

Journal ref: J. Chem. Theory Comput., 2018, 14 (4), 1969-1981

Cluster decomposition of full configuration interaction wave functions: a tool for chemical interpretation of systems with strong correlation

Authors: Susi Lehtola, Norm M. Tubman, K. Birgitta Whaley, Martin Head-Gordon

Abstract: Approximate full configuration interaction (FCI) calculations have recently become tractable for systems of unforeseen size thanks to stochastic and adaptive approximations to the exponentially scaling FCI problem. The result of an FCI calculation is a weighted set of electronic configurations, which can also be expressed in terms of excitations from a reference configuration. The excitation ampli… ▽ More Approximate full configuration interaction (FCI) calculations have recently become tractable for systems of unforeseen size thanks to stochastic and adaptive approximations to the exponentially scaling FCI problem. The result of an FCI calculation is a weighted set of electronic configurations, which can also be expressed in terms of excitations from a reference configuration. The excitation amplitudes contain information on the complexity of the electronic wave function, but this information is contaminated by contributions from disconnected excitations, i.e. those excitations that are just products of independent lower-level excitations. The unwanted contributions can be removed via a cluster decomposition procedure, making it possible to examine the importance of connected excitations in complicated multireference molecules which are outside the reach of conventional algorithms. We present an implementation of the cluster decomposition analysis and apply it to both true FCI wave functions, as well as wave functions generated from the adaptive sampling CI (ASCI) algorithm. The cluster decomposition is useful for interpreting calculations in chemical studies, as a diagnostic for the convergence of various excitation manifolds, as well as as a guidepost for polynomially scaling electronic structure models. Applications are presented for (i) the double dissociation of water, (ii) the carbon dimer, (iii) the πぱい space of polyacenes, as well as (iv) the chromium dimer. While the cluster amplitudes exhibit rapid decay with increasing rank for the first three systems, even connected octuple excitations still appear important in Cr$_2$, suggesting that spin-restricted single-reference coupled-cluster approaches may not be tractable for some problems in transition metal chemistry. △ Less

Submitted 24 October, 2017; v1 submitted 13 July, 2017; originally announced July 2017.

Comments: 15 pages, 5 figures

Journal ref: J. Chem. Phys. 147, 154105 (2017)

Orbital optimization in the perfect pairing hierarchy. Applications to full-valence calculations on linear polyacenes

Authors: Susi Lehtola, John Parkhill, Martin Head-Gordon

Abstract: We describe the implementation of orbital optimization for the models in the perfect pairing hierarchy [Lehtola et al, J. Chem. Phys. 145, 134110 (2016)]. Orbital optimization, which is generally necessary to obtain reliable results, is pursued at perfect pairing (PP) and perfect quadruples (PQ) levels of theory for applications on linear polyacenes, which are believed to exhibit strong correlatio… ▽ More We describe the implementation of orbital optimization for the models in the perfect pairing hierarchy [Lehtola et al, J. Chem. Phys. 145, 134110 (2016)]. Orbital optimization, which is generally necessary to obtain reliable results, is pursued at perfect pairing (PP) and perfect quadruples (PQ) levels of theory for applications on linear polyacenes, which are believed to exhibit strong correlation in the πぱい space. While local minima and σしぐま-πぱい symmetry breaking solutions were found for PP orbitals, no such problems were encountered for PQ orbitals. The PQ orbitals are used for single-point calculations at PP, PQ and perfect hextuples (PH) levels of theory, both only in the πぱい subspace, as well as in the full σしぐまπぱい valence space. It is numerically demonstrated that the inclusion of single excitations is necessary also when optimized orbitals are used. PH is found to yield good agreement with previously published density matrix renormalization group (DMRG) data in the πぱい space, capturing over 95% of the correlation energy. Full-valence calculations made possible by our novel, efficient code reveal that strong correlations are weaker when larger bases or active spaces are employed than in previous calculations. The largest full-valence PH calculations presented correspond to a (192e,192o) problem. △ Less

Submitted 7 March, 2018; v1 submitted 3 May, 2017; originally announced May 2017.

Comments: 19 pages, 4 figures

Journal ref: Mol. Phys. 116, 547 (2018)

Ultrafast Triplet Pair Formation and Subsequent Thermally Activated Dissociation Control Efficient Endothermic Singlet Exciton Fission

Authors: Hannah L. Stern, Alexandre Cheminal, Shane R. Yost, Katharina Broch, Sam L. Bayliss, Kai Chen, Maxim Tabachnyk, Karl Thorley, Neil Greenham, Justin Hodgkiss, John Anthony, Martin Head-Gordon, Andrew J. Musser, Akshay Rao, Richard H. Friend

Abstract: Singlet exciton fission (SF), the conversion of one spin-singlet exciton (S1) into two spin-triplet excitons (T1), could provide a means to overcome the Shockley-Queisser limit in photovoltaics. SF as measured by the decay of S1 has been shown to occur efficiently and independently of temperature even when the energy of S1 is as much as 200 meV less than 2T1. Here, we study films of TIPS-tetracene… ▽ More Singlet exciton fission (SF), the conversion of one spin-singlet exciton (S1) into two spin-triplet excitons (T1), could provide a means to overcome the Shockley-Queisser limit in photovoltaics. SF as measured by the decay of S1 has been shown to occur efficiently and independently of temperature even when the energy of S1 is as much as 200 meV less than 2T1. Here, we study films of TIPS-tetracene using transient optical spectroscopy and show that the initial rise of the triplet pair state (TT) occurs in 300 fs, matched by rapid loss of S1 stimulated emission, and that this process is mediated by the strong coupling of electronic and vibrational degrees of freedom. This is followed by a slower 10 ps morphology-dependent phase of S1 decay and TT growth. We observe the TT to be thermally dissociated on 10-100 ns timescales to form free triplets. This provides a model for "temperature independent", efficient TT formation and thermally activated TT separation. △ Less

Submitted 5 April, 2017; originally announced April 2017.

Cost-effective description of strong correlation: efficient implementations of the perfect quadruples and perfect hextuples models

Authors: Susi Lehtola, John Parkhill, Martin Head-Gordon

Abstract: Novel implementations based on dense tensor storage are presented for the singlet-reference perfect quadruples (PQ) [Parkhill, Lawler, and Head-Gordon, J. Chem. Phys. 130, 084101 (2009)] and perfect hextuples (PH) [Parkhill and Head-Gordon, J. Chem. Phys. 133, 024103 (2010)] models. The methods are obtained as block decompositions of conventional coupled-cluster theory that are exact for four elec… ▽ More Novel implementations based on dense tensor storage are presented for the singlet-reference perfect quadruples (PQ) [Parkhill, Lawler, and Head-Gordon, J. Chem. Phys. 130, 084101 (2009)] and perfect hextuples (PH) [Parkhill and Head-Gordon, J. Chem. Phys. 133, 024103 (2010)] models. The methods are obtained as block decompositions of conventional coupled-cluster theory that are exact for four electrons in four orbitals (PQ) and six electrons in six orbitals (PH), but that can also be applied to much larger systems. PQ and PH have storage requirements that scale as the square, and as the cube of the number of active electrons, respectively, and exhibit quartic scaling of the computational effort for large systems. Applications of the new implementations are presented for full-valence calculations on linear polyenes (C n H n+2 ), which highlight the excellent computational scaling of the present implementations that can routinely handle active spaces of hundreds of electrons. The accuracy of the models is studied in the πぱい space of the polyenes, in hydrogen chains (H 50 ), and in the πぱい space of polyacene molecules. In all cases, the results compare favorably to density matrix renormalization group values. With the novel implementation of PQ, active spaces of 140 electrons in 140 orbitals can be solved in a matter of minutes on a single core workstation, and the relatively low polynomial scaling means that very large systems are also accessible using parallel computing. △ Less

Submitted 31 August, 2016; originally announced September 2016.

Comments: 13 pages, 7 figures

Journal ref: J. Chem. Phys. 145, 134110 (2016)

A deterministic alternative to the full configuration interaction quantum Monte Carlo method

Authors: Norm M. Tubman, Joonho Lee, Tyler Y. Takeshita, Martin Head-Gordon, K. Birgitta Whaley

Abstract: Development of exponentially scaling methods has seen great progress in tackling larger systems than previously thought possible. One such technique, full configuration interaction quantum Monte Carlo, is a useful algorithm that allows exact diagonalization through stochastically sampling determinants. The method derives its utility from the information in the matrix elements of the Hamiltonian, a… ▽ More Development of exponentially scaling methods has seen great progress in tackling larger systems than previously thought possible. One such technique, full configuration interaction quantum Monte Carlo, is a useful algorithm that allows exact diagonalization through stochastically sampling determinants. The method derives its utility from the information in the matrix elements of the Hamiltonian, along with a stochastic projected wave function, to find the important parts of Hilbert space. However, the stochastic representation of the wave function is not required to search Hilbert space efficiently, and here we describe a highly efficient deterministic method to achieve chemical accuracy for a wide range of systems, including the difficult Cr$_{2}$ dimer. In addition our method also allows efficient calculation of excited state energies, for which we illustrate with benchmark results for the excited states of C$_{2}$. △ Less

Submitted 8 March, 2016; originally announced March 2016.

Comments: 4 pages, 2 figures

Journal ref: J. Chem. Phys. 145, 044112 (2016)

Improving approximate-optimized effective potentials by imposing exact conditions: Theory and applications to electronic statics and dynamics

Authors: Yair Kurzweil Martin Head-Gordon

Abstract: We develop a method that can constrain any local exchange-correlation potential to preserve ba-sic exact conditions. Using the method of Lagrange multipliers, we calculate for each set of given Kohn-Sham orbitals, a constraint-preserving potential which is closest to the given exchange-correlation potential. The method is applicable to both the time-dependent (TD) and independent cases. The exac… ▽ More We develop a method that can constrain any local exchange-correlation potential to preserve ba-sic exact conditions. Using the method of Lagrange multipliers, we calculate for each set of given Kohn-Sham orbitals, a constraint-preserving potential which is closest to the given exchange-correlation potential. The method is applicable to both the time-dependent (TD) and independent cases. The exact conditions that are enforced for the time-independent case are: Galilean covariance, zero net force and torque, and Levy-Perdew virial theorem. For the TD case we enforce translational covariance, zero net force, the Levy-Perdew virial theorem and energy balance. We test our method on the exchange Krieger-Li-Iafrate (xKLI) approximate-optimized effective potential (AOEP), for both cases. For the time-independent case, we calculated the ground state properties of some hydrogen chains and small sodium clusters, for some constrained xKLI potentials and Hartree-Fock (HF) exchange. The results indicate that enforcing the exact conditions is not important for these cases. On the other hand, in the case, constraining both energy balance and zero net force yields improved results relative to TDHF calculations. We explored the electron dynamics in small sodium clusters driven by CW laser pulses. For each laser pulse we compared calculations from TD constrained xKLI, TD partially constrained xKLI and TDHF. We found that electron dynamics like elec-tron ionization and moment of inertia dynamics for the constrained xKLI are most similar to the TDHF results. Also, energy conservation is better by at least an order of magnitude with respect to the unconstrained xKLI. We also discuss the problems that arise in satisfying constraints in the TD case with a non CW driving force. △ Less

Submitted 8 May, 2009; originally announced May 2009.

Analysis of multi-configuration Kohn-Sham methods

Authors: Yair Kurzweil, M. Head-Gordon

Abstract: We consider the extension of the standard single-determinant Kohn-Sham method to the case of a multiconfiguration trial wavefunction. By applying the rigorous Kohn-Sham method to this case, we construct the proper interacting and non-interacting energy functionals. Following the Hohenberg-Kohn theorem for both energy functionals, we derive the corre-sponding multiconfiguration Kohn-Sham equation… ▽ More We consider the extension of the standard single-determinant Kohn-Sham method to the case of a multiconfiguration trial wavefunction. By applying the rigorous Kohn-Sham method to this case, we construct the proper interacting and non-interacting energy functionals. Following the Hohenberg-Kohn theorem for both energy functionals, we derive the corre-sponding multiconfiguration Kohn-Sham equations. At the end of the analysis we show that, at the ground state, the multiconfiguration wavefunction must collapse into a single-determinant wavefunction, equal to the regular KS wavefunction. We also discuss the non-collapse of the wavefunction in other multiconfiguration density functional theory methods where the auxiliary system is partially interacting. △ Less

Submitted 5 December, 2008; originally announced December 2008.

Simulated Quantum Computation of Molecular Energies

Authors: Alán Aspuru-Guzik, Anthony D. Dutoi, Peter J. Love, Martin Head-Gordon

Abstract: The calculation time for the energy of atoms and molecules scales exponentially with system size on a classical computer but polynomially using quantum algorithms. We demonstrate that such algorithms can be applied to problems of chemical interest using modest numbers of quantum bits. Calculations of the water and lithium hydride molecular ground-state energies have been carried out on a quantum… ▽ More The calculation time for the energy of atoms and molecules scales exponentially with system size on a classical computer but polynomially using quantum algorithms. We demonstrate that such algorithms can be applied to problems of chemical interest using modest numbers of quantum bits. Calculations of the water and lithium hydride molecular ground-state energies have been carried out on a quantum computer simulator using a recursive phase-estimation algorithm. The recursive algorithm reduces the number of quantum bits required for the readout register from about 20 to 4. Mappings of the molecular wave function to the quantum bits are described. An adiabatic method for the preparation of a good approximate ground-state wave function is described and demonstrated for a stretched hydrogen molecule. The number of quantum bits required scales linearly with the number of basis functions, and the number of gates required grows polynomially with the number of quantum bits. △ Less

Submitted 26 April, 2006; originally announced April 2006.

Comments: 20 pages, 3 figures

Journal ref: Science 309 p. 1704, 2005