1
|
Shin H, Gasperich K, Rojas T, Ngo AT, Krogel JT, Benali A. Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction. J Chem Theory Comput 2024. [PMID: 39213586 DOI: 10.1021/acs.jctc.4c00335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. We show that the sCI framework enables accurate prediction of band gaps in LiCoO2, closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t2g and eg orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.
Collapse
Affiliation(s)
- Hyeondeok Shin
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Kevin Gasperich
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Tomas Rojas
- Department of Chemical Engineering, University of Illinois at Chicago, Chicago, Illinois 60608, United States
- Material Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Anh T Ngo
- Department of Chemical Engineering, University of Illinois at Chicago, Chicago, Illinois 60608, United States
- Material Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Anouar Benali
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| |
Collapse
|
2
|
Kossoski F, Boggio-Pasqua M, Loos PF, Jacquemin D. Reference Energies for Double Excitations: Improvement and Extension. J Chem Theory Comput 2024; 20:5655-5678. [PMID: 38885174 DOI: 10.1021/acs.jctc.4c00410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
In the realm of photochemistry, the significance of double excitations (also known as doubly excited states), where two electrons are concurrently elevated to higher energy levels, lies in their involvement in key electronic transitions essential in light-induced chemical reactions as well as their challenging nature from the computational theoretical chemistry point of view. Based on state-of-the-art electronic structure methods (such as high-order coupled-cluster, selected configuration interaction, and multiconfigurational methods), we improve and expand our prior set of accurate reference excitation energies for electronic states exhibiting a substantial amount of double excitations [Loos et al. J. Chem. Theory Comput. 2019, 15, 1939]. This extended collection encompasses 47 electronic transitions across 26 molecular systems that we separate into two distinct subsets: (i) 28 "genuine" doubly excited states where the transitions almost exclusively involve doubly excited configurations and (ii) 19 "partial" doubly excited states which exhibit a more balanced character between singly and doubly excited configurations. For each subset, we assess the performance of high-order coupled-cluster (CC3, CCSDT, CC4, and CCSDTQ) and multiconfigurational methods (CASPT2, CASPT3, PC-NEVPT2, and SC-NEVPT2). Using as a probe the percentage of single excitations involved in a given transition (%T1) computed at the CC3 level, we also propose a simple correction that reduces the errors of CC3 by a factor of 3, for both sets of excitations. We hope that this more complete and diverse compilation of double excitations will help future developments of electronic excited-state methodologies.
Collapse
Affiliation(s)
- Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| |
Collapse
|
3
|
Zhou H, Kincaid B, Wang G, Annaberdiyev A, Ganesh P, Mitas L. A new generation of effective core potentials: Selected lanthanides and heavy elements. J Chem Phys 2024; 160:084302. [PMID: 38391016 DOI: 10.1063/5.0180057] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Accepted: 11/05/2023] [Indexed: 02/24/2024] Open
Abstract
We construct correlation-consistent effective core potentials (ccECPs) for a selected set of heavy atoms and f elements that are currently of significant interest in materials and chemical applications, including Y, Zr, Nb, Rh, Ta, Re, Pt, Gd, and Tb. As is customary, ccECPs consist of spin-orbit (SO) averaged relativistic effective potential (AREP) and effective SO terms. For the AREP part, our constructions are carried out within a relativistic coupled-cluster framework while also taking into account objective function one-particle characteristics for improved convergence in optimizations. The transferability is adjusted using binding curves of hydride and oxide molecules. We address the difficulties encountered with f elements, such as the presence of large cores and multiple near-degeneracies of excited levels. For these elements, we construct ccECPs with core-valence partitioning that includes 4f subshell in the valence space. The developed ccECPs achieve an excellent balance between accuracy, size of the valence space, and transferability and are also suitable to be used in plane wave codes with reasonable energy cutoffs.
Collapse
Affiliation(s)
- Haihan Zhou
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Benjamin Kincaid
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Guangming Wang
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Abdulgani Annaberdiyev
- Center for Nanophase Materials Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Panchapakesan Ganesh
- Center for Nanophase Materials Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Lubos Mitas
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| |
Collapse
|
4
|
Jacquemin D, Kossoski F, Gam F, Boggio-Pasqua M, Loos PF. Reference Vertical Excitation Energies for Transition Metal Compounds. J Chem Theory Comput 2023. [PMID: 37965941 DOI: 10.1021/acs.jctc.3c01080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
To enrich and enhance the diversity of the quest database of highly accurate excitation energies [Véril, M.; et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with a singlet or doublet ground state containing a fourth-row transition metal (CuCl, CuF, CuH, ScF, ScH, ScO, ScS, TiN, ZnH, ZnO, and ZnS) are considered, and the corresponding excitation energies are computed using high-level coupled-cluster (CC) methods, namely, CC3, CCSDT, CC4, and CCSDTQ, as well as multiconfigurational methods such as CASPT2 and NEVPT2. In many cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. Based on these calculations, theoretical best estimates of the transition energies are established in both the aug-cc-pVDZ and aug-cc-pVTZ basis sets. This allows us to accurately assess the performance of the CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported.
Collapse
Affiliation(s)
- Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Franck Gam
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| |
Collapse
|
5
|
Kossoski F, Loos PF. Seniority and Hierarchy Configuration Interaction for Radicals and Excited States. J Chem Theory Comput 2023. [PMID: 37965728 DOI: 10.1021/acs.jctc.3c00946] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
Hierarchy configuration interaction (hCI) has recently been introduced as an alternative configuration interaction (CI) route combining excitation degree and seniority number and has been shown to efficiently recover both dynamic and static correlations for closed-shell molecular systems [ J. Phys. Chem. Lett. 2022, 13, 4342]. Here we generalize hCI for an arbitrary reference determinant, allowing calculations for radicals and excited states in a state-specific way. We gauge this route against excitation-based CI (eCI) and seniority-based CI (sCI) by evaluating how different ground-state properties of radicals converge to the full CI limit. We find that hCI outperforms or matches eCI, whereas sCI is far less accurate, in line with previous observations for closed-shell molecules. Employing second-order Epstein-Nesbet (EN2) perturbation theory as a correction significantly accelerates the convergence of hCI and eCI. We further explore various hCI and sCI models to calculate the excitation energies of closed- and open-shell systems. Our results underline that the choice of both the reference determinant and the set of orbitals drives the fine balance between correlation of ground and excited states. State-specific hCI2 and higher-order models perform similarly to their eCI counterparts, whereas lower orders of hCI deliver poor results unless supplemented by the EN2 correction, which substantially improves their accuracy. In turn, sCI1 produces decent excitation energies for radicals, encouraging the development of related seniority-based coupled-cluster methods.
Collapse
Affiliation(s)
- Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| |
Collapse
|
6
|
Kossoski F, Loos PF. State-Specific Configuration Interaction for Excited States. J Chem Theory Comput 2023; 19:2258-2269. [PMID: 37024102 PMCID: PMC10134430 DOI: 10.1021/acs.jctc.3c00057] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
Abstract
We introduce and benchmark a systematically improvable route for excited-state calculations, labeled state-specific configuration interaction (ΔCI), which is a particular realization of multiconfigurational self-consistent field and multireference configuration interaction. Starting with a reference built from optimized configuration state functions, separate CI calculations are performed for each targeted state (hence, state-specific orbitals and determinants). Accounting for single and double excitations produces the ΔCISD model, which can be improved with second-order Epstein-Nesbet perturbation theory (ΔCISD+EN2) or a posteriori Davidson corrections (ΔCISD+Q). These models were gauged against a vast and diverse set of 294 reference excitation energies. We have found that ΔCI is significantly more accurate than standard ground-state-based CI, whereas close performances were found between ΔCISD and EOM-CC2 and between ΔCISD+EN2 and EOM-CCSD. For larger systems, ΔCISD+Q delivers more accurate results than EOM-CC2 and EOM-CCSD. The ΔCI route can handle challenging multireference problems, singly and doubly excited states, from closed- and open-shell species, with overall comparable accuracy and thus represents a promising alternative to more established methodologies. In its current form, however, it is reliable only for relatively low-lying excited states.
Collapse
Affiliation(s)
- Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| |
Collapse
|
7
|
Wheeler WA, Pathak S, Kleiner KG, Yuan S, Rodrigues JNB, Lorsung C, Krongchon K, Chang Y, Zhou Y, Busemeyer B, Williams KT, Muñoz A, Chow CY, Wagner LK. PyQMC: An all-Python real-space quantum Monte Carlo module in PySCF. J Chem Phys 2023; 158:114801. [PMID: 36948839 DOI: 10.1063/5.0139024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023] Open
Abstract
We describe a new open-source Python-based package for high accuracy correlated electron calculations using quantum Monte Carlo (QMC) in real space: PyQMC. PyQMC implements modern versions of QMC algorithms in an accessible format, enabling algorithmic development and easy implementation of complex workflows. Tight integration with the PySCF environment allows for a simple comparison between QMC calculations and other many-body wave function techniques, as well as access to high accuracy trial wave functions.
Collapse
Affiliation(s)
- William A Wheeler
- Department of Materials Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Shivesh Pathak
- Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Kevin G Kleiner
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Shunyue Yuan
- Department of Applied Physics and Materials Science, California Institute of Technology, Pasadena, California 91125, USA
| | - João N B Rodrigues
- Centro de Ciências Naturais e Humanas, Universidade Federal do ABC-UFABC, Santo André, São Paulo 09210-580, Brazil
| | - Cooper Lorsung
- Department of Mechanical Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Kittithat Krongchon
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Yueqing Chang
- Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Yiqing Zhou
- Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York 14853, USA
| | | | | | - Alexander Muñoz
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Chun Yu Chow
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Lucas K Wagner
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
8
|
Otis L, Neuscamman E. Optimization Stability in Excited-State-Specific Variational Monte Carlo. J Chem Theory Comput 2023; 19:767-782. [PMID: 36662538 DOI: 10.1021/acs.jctc.2c00642] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
We investigate the issue of optimization stability in variance-based state-specific variational Monte Carlo, discussing the roles of the objective function, the complexity of wave function ansatz, the amount of sampling effort, and the choice of minimization algorithm. Using a small cyanine dye molecule as a test case, we systematically perform minimizations using variants of the linear method as both a standalone algorithm and in a hybrid combination with accelerated descent. We demonstrate that adaptive step control is crucial for maintaining the linear method's stability when optimizing complicated wave functions and that the hybrid method enjoys both greater stability and minimization performance.
Collapse
Affiliation(s)
- Leon Otis
- Department of Physics, University of California Berkeley, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
9
|
Otis L, Neuscamman E. A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1659] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Affiliation(s)
- Leon Otis
- Department of Physics University of California Berkeley Berkeley California USA
| | - Eric Neuscamman
- Department of Chemistry University of California Berkeley Berkeley California USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley California USA
| |
Collapse
|
10
|
Entwistle MT, Schätzle Z, Erdman PA, Hermann J, Noé F. Electronic excited states in deep variational Monte Carlo. Nat Commun 2023; 14:274. [PMID: 36650151 PMCID: PMC9845370 DOI: 10.1038/s41467-022-35534-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Accepted: 12/08/2022] [Indexed: 01/19/2023] Open
Abstract
Obtaining accurate ground and low-lying excited states of electronic systems is crucial in a multitude of important applications. One ab initio method for solving the Schrödinger equation that scales favorably for large systems is variational quantum Monte Carlo (QMC). The recently introduced deep QMC approach uses ansatzes represented by deep neural networks and generates nearly exact ground-state solutions for molecules containing up to a few dozen electrons, with the potential to scale to much larger systems where other highly accurate methods are not feasible. In this paper, we extend one such ansatz (PauliNet) to compute electronic excited states. We demonstrate our method on various small atoms and molecules and consistently achieve high accuracy for low-lying states. To highlight the method's potential, we compute the first excited state of the much larger benzene molecule, as well as the conical intersection of ethylene, with PauliNet matching results of more expensive high-level methods.
Collapse
Affiliation(s)
- M T Entwistle
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 12, 14195, Berlin, Germany
| | - Z Schätzle
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 12, 14195, Berlin, Germany
| | - P A Erdman
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 12, 14195, Berlin, Germany
| | - J Hermann
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 12, 14195, Berlin, Germany.
| | - F Noé
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 12, 14195, Berlin, Germany.
- Microsoft Research AI4Science, Berlin, Germany.
- Department of Physics, FU Berlin, Arnimallee 14, 14195, Berlin, Germany.
- Department of Chemistry, Rice University, Houston, TX, 77005, USA.
| |
Collapse
|
11
|
Hanscam R, Neuscamman E. Applying Generalized Variational Principles to Excited-State-Specific Complete Active Space Self-consistent Field Theory. J Chem Theory Comput 2022; 18:6608-6621. [PMID: 36215108 DOI: 10.1021/acs.jctc.2c00639] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We employ a generalized variational principle to improve the stability, reliability, and precision of fully excited-state-specific complete active space self-consistent field theory. Compared to previous approaches that similarly seek to tailor this ansatz's orbitals and configuration interaction expansion for an individual excited state, we find the present approach to be more resistant to root flipping and better at achieving tight convergence to an energy stationary point. Unlike state-averaging, this approach allows orbital shapes to be optimal for individual excited states, which is especially important for charge-transfer states and some doubly excited states. We demonstrate the convergence and state-targeting abilities of this method in LiH, ozone, and MgO, showing in the latter that it is capable of finding three excited-state energy stationary points that no previous method has been able to locate.
Collapse
Affiliation(s)
- Rebecca Hanscam
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
12
|
Loos PF, Lipparini F, Matthews DA, Blondel A, Jacquemin D. A Mountaineering Strategy to Excited States: Revising Reference Values with EOM-CC4. J Chem Theory Comput 2022; 18:4418-4427. [PMID: 35737466 DOI: 10.1021/acs.jctc.2c00416] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In the framework of the computational determination of highly accurate vertical excitation energies in small organic compounds, we explore the possibilities offered by the equation-of-motion formalism relying on the approximate fourth-order coupled-cluster (CC) method, CC4. We demonstrate, using an extended set of more than 200 reference values based on CC including up to quadruples excitations (CCSDTQ), that CC4 is an excellent approximation to CCSDTQ for excited states with a dominant contribution from single excitations with an average deviation as small as 0.003 eV. We next assess the accuracy of several additive basis set correction schemes, in which vertical excitation energies obtained with a compact basis set and a high-order CC method are corrected with lower-order CC calculations performed in a larger basis set. Such strategies are found to be overall very beneficial, though their accuracy depends significantly on the actual scheme. Finally, CC4 is employed to improve several theoretical best estimates of the QUEST database for molecules containing between four and six (nonhydrogen) atoms, for which previous estimates were computed at the CCSDT level.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Devin A Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Aymeric Blondel
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
13
|
Burton HGA. Energy Landscape of State-Specific Electronic Structure Theory. J Chem Theory Comput 2022; 18:1512-1526. [PMID: 35179023 PMCID: PMC9082508 DOI: 10.1021/acs.jctc.1c01089] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Indexed: 11/29/2022]
Abstract
State-specific approximations can provide a more accurate representation of challenging electronic excitations by enabling relaxation of the electron density. While state-specific wave functions are known to be local minima or saddle points of the approximate energy, the global structure of the exact electronic energy remains largely unexplored. In this contribution, a geometric perspective on the exact electronic energy landscape is introduced. On the exact energy landscape, ground and excited states form stationary points constrained to the surface of a hypersphere, and the corresponding Hessian index increases at each excitation level. The connectivity between exact stationary points is investigated, and the square-magnitude of the exact energy gradient is shown to be directly proportional to the Hamiltonian variance. The minimal basis Hartree-Fock and excited-state mean-field representations of singlet H2 (STO-3G) are then used to explore how the exact energy landscape controls the existence and properties of state-specific approximations. In particular, approximate excited states correspond to constrained stationary points on the exact energy landscape, and their Hessian index also increases for higher energies. Finally, the properties of the exact energy are used to derive the structure of the variance optimization landscape and elucidate the challenges faced by variance optimization algorithms, including the presence of unphysical saddle points or maxima of the variance.
Collapse
Affiliation(s)
- Hugh G. A. Burton
- Physical and Theoretical
Chemistry Laboratory, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
14
|
Corzo HH, Abou Taka A, Pribram-Jones A, Hratchian HP. Using projection operators with maximum overlap methods to simplify challenging self-consistent field optimization. J Comput Chem 2021; 43:382-390. [PMID: 34936117 DOI: 10.1002/jcc.26797] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2021] [Revised: 10/30/2021] [Accepted: 11/28/2021] [Indexed: 02/04/2023]
Abstract
Maximum overlap methods are effective tools for optimizing challenging ground- and excited-state wave functions using self-consistent field models such as Hartree-Fock and Kohn-Sham density functional theory. Nevertheless, such models have shown significant sensitivity to the user-defined initial guess of the target wave function. In this work, a projection operator framework is defined and used to provide a metric for non-aufbau orbital selection in maximum-overlap-methods. The resulting algorithms, termed the Projection-based Maximum Overlap Method (PMOM) and Projection-based Initial Maximum Overlap Method (PIMOM), are shown to perform exceptionally well when using simple user-defined target solutions based on occupied/virtual molecular orbital permutations. This work also presents a new metric that provides a simple and conceptually convenient measure of agreement between the desired target and the current or final SCF results during a calculation employing a maximum-overlap method.
Collapse
Affiliation(s)
- Hector H Corzo
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, California, USA
| | - Ali Abou Taka
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, California, USA
| | - Aurora Pribram-Jones
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, California, USA
| | - Hrant P Hratchian
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, California, USA
| |
Collapse
|
15
|
|