1
|
Damour Y, Scemama A, Kossoski F, Loos PF. Selected Configuration Interaction for Resonances. J Phys Chem Lett 2024; 15:8296-8305. [PMID: 39107252 DOI: 10.1021/acs.jpclett.4c02060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/09/2024]
Abstract
Electronic resonances are metastable states that can decay by electron loss. They are ubiquitous across various fields of science, such as chemistry, physics, and biology. However, current theoretical and computational models for resonances cannot yet rival the level of accuracy achieved by bound-state methodologies. Here, we generalize selected configuration interaction (SCI) to treat resonances by using the complex absorbing potential (CAP) technique. By modifying the selection procedure and the extrapolation protocol of standard SCI, the resulting CAP-SCI method yields resonance positions and widths of full configuration interaction quality. Initial results for the shape resonances of N2- and CO- reveal the important effect of high-order correlation, which shifts the values obtained with CAP-augmented equation-of-motion coupled-cluster with singles and doubles by more than 0.1 eV. The present CAP-SCI approach represents a cornerstone in the development of highly accurate methodologies for resonances.
Collapse
Affiliation(s)
- Yann Damour
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - 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
|
2
|
Evangelista FA, Li C, Verma P, Hannon KP, Schriber JB, Zhang T, Cai C, Wang S, He N, Stair NH, Huang M, Huang R, Misiewicz JP, Li S, Marin K, Zhao Z, Burns LA. Forte: A suite of advanced multireference quantum chemistry methods. J Chem Phys 2024; 161:062502. [PMID: 39132791 DOI: 10.1063/5.0216512] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2024] [Accepted: 06/24/2024] [Indexed: 08/13/2024] Open
Abstract
Forte is an open-source library specialized in multireference electronic structure theories for molecular systems and the rapid prototyping of new methods. This paper gives an overview of the capabilities of Forte, its software architecture, and examples of applications enabled by the methods it implements.
Collapse
Affiliation(s)
- Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Prakash Verma
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Kevin P Hannon
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Jeffrey B Schriber
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
- Department of Chemistry and Biochemistry, Iona University, New Rochelle, New York 10801, USA
| | - Tianyuan Zhang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Chenxi Cai
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Shuhe Wang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Nan He
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Nicholas H Stair
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Meng Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Renke Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Jonathon P Misiewicz
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Shuhang Li
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Kevin Marin
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Zijun Zhao
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
3
|
Loos PF, Jacquemin D. A mountaineering strategy to excited states: Accurate vertical transition energies and benchmarks for substituted benzenes. J Comput Chem 2024; 45:1791-1805. [PMID: 38661240 DOI: 10.1002/jcc.27358] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2024] [Revised: 02/09/2024] [Accepted: 02/15/2024] [Indexed: 04/26/2024]
Abstract
In an effort to expand the existing QUEST database of accurate vertical transition energies [Véril et al. WIREs Comput. Mol. Sci. 2021, 11, e1517], we have modeled more than 100 electronic excited states of different natures (local, charge-transfer, Rydberg, singlet, and triplet) in a dozen of mono- and di-substituted benzenes, including aniline, benzonitrile, chlorobenzene, fluorobenzene, nitrobenzene, among others. To establish theoretical best estimates for these vertical excitation energies, we have employed advanced coupled-cluster methods including iterative triples (CC3 and CCSDT) and, when technically possible, iterative quadruples (CC4). These high-level computational approaches provide a robust foundation for benchmarking a series of popular wave function methods. The evaluated methods all include contributions from double excitations (ADC(2), CC2, CCSD, CIS(D), EOM-MP2, STEOM-CCSD), along with schemes that also incorporate perturbative or iterative triples (ADC(3), CCSDR(3), CCSD(T)(a) ⋆ , and CCSDT-3). This systematic exploration not only broadens the scope of the QUEST database but also facilitates a rigorous assessment of different theoretical approaches in the framework of a homologous chemical series, offering valuable insights into the accuracy and reliability of these methods in such cases. We found that both ADC(2.5) and CCSDT-3 can provide very consistent estimates, whereas among less expensive methods SCS-CC2 is likely the most effective approach. Importantly, we show that some lower order methods may offer reasonable trends in the homologous series while providing quite large average errors, and vice versa. Consequently, benchmarking the accuracy of a model based solely on absolute transition energies may not be meaningful for applications involving a series of similar compounds.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes, France
- Institut Universitaire de France (IUF), Paris, France
| |
Collapse
|
4
|
Johnson PA. Beyond a Richardson-Gaudin Mean-Field: Slater-Condon Rules and Perturbation Theory. J Phys Chem A 2024; 128:6033-6045. [PMID: 39007410 DOI: 10.1021/acs.jpca.4c02857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/16/2024]
Abstract
Richardson-Gaudin states provide a basis of the Hilbert space for strongly correlated electrons. In this study, optimal expressions for the transition density matrix elements between Richardson-Gaudin states are obtained with a cost comparable with the corresponding reduced density matrix elements. Analogues of the Slater-Condon rules are identified based on the number of near-zero singular values of the RG state overlap matrix. Finally, a perturbative approach is shown to be close in quality to a configuration interaction of Richardson-Gaudin states while being feasible to compute.
Collapse
Affiliation(s)
- Paul A Johnson
- Département de Chimie, Université Laval, Québec, Québec G1V 0A6, Canada
| |
Collapse
|
5
|
Zhao Z, Evangelista FA. Toward Accurate Spin-Orbit Splittings from Relativistic Multireference Electronic Structure Theory. J Phys Chem Lett 2024; 15:7103-7110. [PMID: 38954768 PMCID: PMC11261625 DOI: 10.1021/acs.jpclett.4c01372] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2024] [Revised: 06/25/2024] [Accepted: 06/27/2024] [Indexed: 07/04/2024]
Abstract
Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. However, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.
Collapse
Affiliation(s)
- Zijun Zhao
- Department of Chemistry and
Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Francesco A. Evangelista
- Department of Chemistry and
Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
6
|
Scheurer M, Anselmetti GLR, Oumarou O, Gogolin C, Rubin NC. Tailored and Externally Corrected Coupled Cluster with Quantum Inputs. J Chem Theory Comput 2024; 20:5068-5093. [PMID: 38829984 DOI: 10.1021/acs.jctc.4c00037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
We propose to use wave function overlaps obtained from a quantum computer as inputs for the classical split-amplitude techniques, tailored and externally corrected coupled cluster, to achieve balanced treatment of static and dynamic correlation effects in molecular electronic structure simulations. By combining insights from statistical properties of matchgate shadows, which are used to measure quantum trial state overlaps, with classical correlation diagnostics, we can provide quantum resource estimates well into the classically no longer exactly solvable regime. We find that rather imperfect wave functions and remarkably low shot counts are sufficient to cure qualitative failures of plain coupled cluster singles doubles and to obtain chemically precise dynamic correlation energy corrections. We provide insights into which wave function preparation schemes have a chance of yielding quantum advantage, and we test our proposed method using overlaps measured on Google's Sycamore device.
Collapse
|
7
|
Marie A, Loos PF. Reference Energies for Valence Ionizations and Satellite Transitions. J Chem Theory Comput 2024; 20:4751-4777. [PMID: 38776293 PMCID: PMC11171335 DOI: 10.1021/acs.jctc.4c00216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 04/10/2024] [Accepted: 04/11/2024] [Indexed: 05/24/2024]
Abstract
Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.
Collapse
Affiliation(s)
- Antoine Marie
- Laboratoire de Chimie et Physique
Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse 31062, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique
Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse 31062, France
| |
Collapse
|
8
|
Hutton L, Moreno Carrascosa A, Prentice AW, Simmermacher M, Runeson JE, Paterson MJ, Kirrander A. Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone. J Chem Phys 2024; 160:204307. [PMID: 38814011 DOI: 10.1063/5.0203667] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 05/05/2024] [Indexed: 05/31/2024] Open
Abstract
Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α-CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.
Collapse
Affiliation(s)
- Lewis Hutton
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Andrés Moreno Carrascosa
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Andrew W Prentice
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, United Kingdom
| | - Mats Simmermacher
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Johan E Runeson
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Martin J Paterson
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, United Kingdom
| | - Adam Kirrander
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
9
|
Damour Y, Scemama A, Jacquemin D, Kossoski F, Loos PF. State-Specific Coupled-Cluster Methods for Excited States. J Chem Theory Comput 2024; 20:4129-4145. [PMID: 38749498 PMCID: PMC11137840 DOI: 10.1021/acs.jctc.4c00034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Revised: 03/06/2024] [Accepted: 03/06/2024] [Indexed: 05/29/2024]
Abstract
We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.
Collapse
Affiliation(s)
- Yann Damour
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| | - Anthony Scemama
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| | - 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, 31000 Toulouse, France
| | - Pierre-François Loos
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| |
Collapse
|
10
|
Li C, Mao S, Huang R, Evangelista FA. Frozen Natural Orbitals for the State-Averaged Driven Similarity Renormalization Group. J Chem Theory Comput 2024; 20:4170-4181. [PMID: 38747709 DOI: 10.1021/acs.jctc.4c00152] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
We present a reduced-cost implementation of the state-averaged driven similarity renormalization group (SA-DSRG) based on the frozen natural orbital (FNO) approach. The natural orbitals (NOs) are obtained by diagonalizing the one-body reduced density matrix from SA-DSRG second-order perturbation theory (SA-DSRG-PT2). We consider three criteria to truncate the virtual NOs for the subsequent electron correlation treatment beyond SA-DSRG-PT2. An additive second-order correction is applied to the SA-DSRG Hamiltonian to reintroduce correlation effects from the discarded orbitals. The FNO SA-DSRG method is benchmarked on 35 small organic molecules in the QUEST database. When keeping 98-99% of the cumulative occupation numbers, the mean absolute error in the vertical transition energies due to FNO is less than 0.01 eV. Using the same FNO threshold, we observe a speedup of 9 times compared to the conventional SA-DSRG implementation for nickel carbonyl with a quadruple-ζ basis set. The FNO approach enables nonperturbative SA-DSRG computations on chloroiron corrole [FeCl(C19H11N4)] with more than 1000 basis functions, surpassing the current limit of a conventional implementation.
Collapse
Affiliation(s)
- Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Shuxian Mao
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Renke Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
11
|
Kulahlioglu AH, Dreuw A. Dense-sparse quantum Monte Carlo algebraic diagrammatic construction and importance ranking. J Chem Phys 2024; 160:204111. [PMID: 38785284 DOI: 10.1063/5.0209137] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2024] [Accepted: 05/07/2024] [Indexed: 05/25/2024] Open
Abstract
Quantum Monte Carlo Algebraic Diagrammatic Construction (QMCADC) has been proposed as a reformulation of the second-order ADC scheme for the polarization propagator within the projection quantum Monte Carlo formalism. Dense-sparse partitioning and importance ranking filtering strategies are now exploited to accelerate its convergence and to alleviate the sign problem inherent in such calculations. By splitting the configuration space into dense and sparse subsets, the corresponding projection operator is decomposed into four distinct blocks. Deterministic calculations handle the dense-to-dense and sparse-to-dense blocks, while the remaining blocks, dense-to-sparse and sparse-to-sparse, are stochastically evaluated. The dense set is efficiently stored in a fixed-size array, and the sparse set is represented through conventional floating random Monte Carlo walks. The stochastic projection is further refined through importance ranking criteria, enabling a reduction in the required number of walkers with a controllable bias. Our results demonstrate the integration of dense-sparse partitioning with importance ranking filtering to significantly enhance the efficiency of QMCADC, enabling large-scale molecular excited-state calculations. Furthermore, this novel approach maximizes the utilization of the sparsity of ADC(2), transforming QMCADC into a tailored framework for ADC calculations.
Collapse
Affiliation(s)
- Adem Halil Kulahlioglu
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
12
|
Tang CL, Heide AG, Heide AD, Douberly GE, Turney JM, Schaefer HF. Exploring the Tl 2 H 2 potential energy surface: A comparative analysis with group 13 systems and experiment. J Comput Chem 2024; 45:985-994. [PMID: 38197269 DOI: 10.1002/jcc.27293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2023] [Revised: 12/04/2023] [Accepted: 12/10/2023] [Indexed: 01/11/2024]
Abstract
Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet Tl 2 H 2 potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug-cc-pwCVQZ-PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero-point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 M 2 H 2 (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.
Collapse
Affiliation(s)
- Carson L Tang
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia, USA
| | - Alexander G Heide
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia, USA
| | - Alexandra D Heide
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia, USA
| | - Gary E Douberly
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia, USA
| | - Justin M Turney
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia, USA
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia, USA
| |
Collapse
|
13
|
Miranda-Quintana RA, Kim TD, Lokhande RA, Richer M, Sánchez-Díaz G, Gaikwad PB, Ayers PW. Flexible Ansatz for N-Body Perturbation Theory. J Phys Chem A 2024; 128:3458-3467. [PMID: 38651558 DOI: 10.1021/acs.jpca.4c00855] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
We propose a new perturbation theory framework that can be used to help with the projective solution of the Schrödinger equation for arbitrary wave functions. This Flexible Ansatz for N-body Perturbation Theory (FANPT) is based on our previously proposed Flexible Ansatz for the N-body Configuration Interaction (FANCI). We derive recursive FANPT expressions, including arbitrary orders in the perturbation hierarchy. We show that the FANPT equations are well-behaved across a wide range of conditions, including static correlation-dominated configurations and highly nonlinear wave functions.
Collapse
Affiliation(s)
- Ramón Alain Miranda-Quintana
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32603, United States
| | - Taewon D Kim
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32603, United States
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, Ontario L8S 4M1, Canada
| | - Rugwed A Lokhande
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32603, United States
| | - M Richer
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, Ontario L8S 4M1, Canada
| | - Gabriela Sánchez-Díaz
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, Ontario L8S 4M1, Canada
| | - Pratiksha B Gaikwad
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32603, United States
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, Ontario L8S 4M1, Canada
| |
Collapse
|
14
|
Burton HGA, Loos PF. Rationale for the extrapolation procedure in selected configuration interaction. J Chem Phys 2024; 160:104102. [PMID: 38456526 DOI: 10.1063/5.0192458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2023] [Accepted: 02/18/2024] [Indexed: 03/09/2024] Open
Abstract
Selected configuration interaction (SCI) methods have emerged as state-of-the-art methodologies for achieving high accuracy and generating benchmark reference data for ground and excited states in small molecular systems. However, their precision relies heavily on extrapolation procedures to produce a final estimate of the exact result. Using the structure of the exact electronic energy landscape, we provide a rationale for the common linear extrapolation of the variational energy as a function of the second-order perturbative correction. In particular, we demonstrate that the energy gap and the coupling between the so-called internal and external spaces are the key factors determining the rate at which the linear regime is reached. Starting from the first principles, we also derive a new non-linear extrapolation formula that improves the post-processing of data generated from SCI methods and can be applied to both ground- and excited-state energies.
Collapse
Affiliation(s)
- Hugh G A Burton
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
15
|
Gao H, Imamura S, Kasagi A, Yoshida E. Distributed Implementation of Full Configuration Interaction for One Trillion Determinants. J Chem Theory Comput 2024; 20:1185-1192. [PMID: 38314701 PMCID: PMC10867839 DOI: 10.1021/acs.jctc.3c01190] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Revised: 12/24/2023] [Accepted: 01/12/2024] [Indexed: 02/07/2024]
Abstract
Full configuration interaction (FCI) can provide an exact molecular ground-state energy within a given basis set and serve as a benchmark for approximate methods in quantum chemical calculations, including the emerging variational quantum eigensolver. However, its exponential computational and memory requirements easily exceed the capability of a single server and limit its applicability to large molecules. In this paper, we present a distributed FCI implementation employing a hybrid parallelization scheme with multithreading and multiprocessing to expand FCI's applicability. We optimize this scheme to minimize the bottlenecks arising from interprocess communications and interthread data management. Our implementation achieves higher scalability than the naive combination of prior works and successfully calculates the exact energy of C3H8/STO-3G with 1.31 trillion determinants, which is the largest FCI calculation to the best of our knowledge. Furthermore, we provide a comprehensive list of FCI results with 136 combinations of molecules and basis sets for future evaluation and development of approximate methods.
Collapse
Affiliation(s)
- Hong Gao
- Computing Laboratory, Fujitsu Laboratories, Fujitsu Limited, Kawasaki City, Kanagawa, 211-0053, Japan
| | - Satoshi Imamura
- Computing Laboratory, Fujitsu Laboratories, Fujitsu Limited, Kawasaki City, Kanagawa, 211-0053, Japan
| | - Akihiko Kasagi
- Computing Laboratory, Fujitsu Laboratories, Fujitsu Limited, Kawasaki City, Kanagawa, 211-0053, Japan
| | - Eiji Yoshida
- Computing Laboratory, Fujitsu Laboratories, Fujitsu Limited, Kawasaki City, Kanagawa, 211-0053, Japan
| |
Collapse
|
16
|
Ren X, Zou J, Zhang H, Li W, Li S. Block-Correlated Coupled Cluster Theory with up to Four-Pair Correlation for Accurate Static Correlation of Strongly Correlated Systems. J Phys Chem Lett 2024; 15:693-700. [PMID: 38207241 DOI: 10.1021/acs.jpclett.3c03373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
A block-correlated coupled cluster method with up to four-pair correlation based on the generalized valence bond wave function (GVB-BCCC4) is first implemented, which offers an alternative method for electronic structure calculations of strongly correlated systems. We developed some techniques to derive a set of compact and cost-effective equations for GVB-BCCC4, which include the definition of n-block (n = 1-4) Hamiltonian matrices, the combination of excitation operators, and the definition of independent amplitudes. We then applied the GVB-BCCC4 method to investigate several potential energy surfaces of strongly correlated systems with singlet ground states. Our calculations demonstrate that the GVB-BCCC4 method can provide nearly exact static correlation energies as the density matrix renormalization group method (on the basis of the same GVB orbitals). This work highlights the significance of four-pair correlation in quantitative descriptions of static correlation energy for strongly correlated systems.
Collapse
Affiliation(s)
- Xiaochuan Ren
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, People's Republic of China
| | - Jingxiang Zou
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, People's Republic of China
| | - Haodong Zhang
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, People's Republic of China
| | - Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, People's Republic of China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, People's Republic of China
| |
Collapse
|
17
|
Prentice AW, Coe JP, Paterson MJ. Modular Approach to Selected Configuration Interaction in an Arbitrary Spin Basis: Implementation and Comparison of Approaches. J Chem Theory Comput 2023; 19:9161-9176. [PMID: 38061390 PMCID: PMC10753805 DOI: 10.1021/acs.jctc.3c00897] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Revised: 11/24/2023] [Accepted: 11/27/2023] [Indexed: 12/27/2023]
Abstract
A modular selected configuration interaction (SCI) code has been developed that is based on the existing Monte-Carlo configuration interaction code (MCCI). The modularity allows various selection protocols to be implemented with ease and allows for fair comparison between wave functions built via different criteria. We have initially implemented adaptations of existing SCI theories, which are based on either energy- or coefficient-driven selection schemes. These codes have been implemented not only in the basis of Slater determinants (SDs) but also in the basis of configuration state functions (CSFs) and extended to state-averaged regimes. This allows one to take advantage of the reduced dimensionality of the wave function in the CSF basis and also the guarantee of pure spin states. All SCI methods were found to be able to predict potential energy surfaces to high accuracy, producing compact wave functions, when compared to full configuration interaction (FCI) for a variety of bond-breaking potential energy surfaces. The compactness of the error-controlled adaptive configuration interaction approach, particularly in the CSF basis, was apparent with nonparallelity errors within chemical accuracy while containing as little as 0.02% of the FCI CSF space. The size-to-accuracy was also extended to FCI spaces approaching one billion configurations.
Collapse
Affiliation(s)
- Andrew W. Prentice
- Institute of Chemical Sciences, School
of Engineering and Physical Sciences, Heriot-Watt
University, Edinburgh EH14 4AS, U.K.
| | - Jeremy P. Coe
- Institute of Chemical Sciences, School
of Engineering and Physical Sciences, Heriot-Watt
University, Edinburgh EH14 4AS, U.K.
| | - Martin J. Paterson
- Institute of Chemical Sciences, School
of Engineering and Physical Sciences, Heriot-Watt
University, Edinburgh EH14 4AS, U.K.
| |
Collapse
|
18
|
Ugandi M, Roemelt M. A configuration-based heatbath-CI for spin-adapted multireference electronic structure calculations with large active spaces. J Comput Chem 2023; 44:2374-2390. [PMID: 37589287 DOI: 10.1002/jcc.27203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Revised: 07/24/2023] [Accepted: 07/27/2023] [Indexed: 08/18/2023]
Abstract
This work reports on a spin-pure configuration-based implementation of the heatbath configuration interaction (HCI) algorithm for selective configuration interaction. Besides the obvious advantage of being spin-pure, the presented method combines the compactness of the configurational ansatz with the known efficiency of the HCI algorithm and a variety of algorithmic and conceptual ideas to achieve a high level of performance. In particular, through pruning of the selected configurational space after HCI selection by means of a more strict criterion, a more compact wavefunction representation is obtained. Moreover, the underlying logic of the method allows us to minimize the number of redundant matrix-matrix multiplications while making use of just-in-time compilation to achieve fast diagonalization of the Hamiltonian. The critical search for 2-electron connections within the configurational space is facilitated by a tree-based representation thereof as suggested previously by Gopal et al. Usage of a prefix-based parallelization and batching during the calculation of the PT2-correction leads to a good load balancing and significantly reduced memory requirements for these critical steps of the calculation. In this way, the need for a semistochastic approach to the PT2 correction is avoided even for large configurational spaces. Finally, several test-cases are discussed to demonstrate the strengths and weaknesses of the presented method.
Collapse
Affiliation(s)
- Mihkel Ugandi
- Institut für Chemie, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Michael Roemelt
- Institut für Chemie, Humboldt-Universität zu Berlin, Berlin, Germany
| |
Collapse
|
19
|
Johnson PA, DePrince AE. Single Reference Treatment of Strongly Correlated H 4 and H 10 Isomers with Richardson-Gaudin States. J Chem Theory Comput 2023; 19:8129-8146. [PMID: 37955440 DOI: 10.1021/acs.jctc.3c00807] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2023]
Abstract
Richardson-Gaudin (RG) states are employed as a variational wave function ansatz for strongly correlated isomers of H4 and H10. In each case, a single RG state describes the seniority-zero sector quite well. Simple natural orbital functionals offer a cheap and reasonable approximation of the outstanding weak correlation in the seniority-zero sector, while systematic improvement is achieved by performing a configuration interaction in terms of RG states.
Collapse
Affiliation(s)
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| |
Collapse
|
20
|
Li X, Huang JC, Zhang GZ, Li HE, Cao CS, Lv D, Hu HS. A Nonstochastic Optimization Algorithm for Neural-Network Quantum States. J Chem Theory Comput 2023; 19:8156-8165. [PMID: 37962975 DOI: 10.1021/acs.jctc.3c00831] [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
Neural-network quantum states (NQS) employ artificial neural networks to encode many-body wave functions in a second quantization through variational Monte Carlo (VMC). They have recently been applied to accurately describe electronic wave functions of molecules and have shown the challenges in efficiency compared with traditional quantum chemistry methods. Here, we introduce a general nonstochastic optimization algorithm for NQS in chemical systems, which deterministically generates a selected set of important configurations simultaneously with energy evaluation of NQS. This method bypasses the need for Markov-chain Monte Carlo within the VMC framework, thereby accelerating the entire optimization process. Furthermore, this newly developed nonstochastic optimization algorithm for NQS offers comparable or superior accuracy compared to its stochastic counterpart and ensures more stable convergence. The application of this model to test molecules exhibiting strong electron correlations provides further insight into the performance of NQS in chemical systems and opens avenues for future enhancements.
Collapse
Affiliation(s)
- Xiang Li
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Jia-Cheng Huang
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Guang-Ze Zhang
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Hao-En Li
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Chang-Su Cao
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
- ByteDance Research, Zhonghang Plaza, No. 43, North Third Ring West Road, Haidian District, Beijing 100089, China
| | - Dingshun Lv
- ByteDance Research, Zhonghang Plaza, No. 43, North Third Ring West Road, Haidian District, Beijing 100089, China
| | - Han-Shi Hu
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| |
Collapse
|
21
|
Coe JP. Analytic Non-adiabatic Couplings for Selected Configuration Interaction via Approximate Degenerate Coupled Perturbed Hartree-Fock. J Chem Theory Comput 2023; 19:8053-8065. [PMID: 37939698 PMCID: PMC10687870 DOI: 10.1021/acs.jctc.3c00601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2023] [Revised: 10/24/2023] [Accepted: 10/25/2023] [Indexed: 11/10/2023]
Abstract
We use degenerate perturbation theory and assume that for degenerate pairs of orbitals, the coupled perturbed Hartree-Fock coefficients are symmetric in the degenerate basis to show [Formula: see text] is the only modification needed in the original molecular orbital basis. This enables us to develop efficient and accurate analytic nonadiabatic couplings between electronic states for selected configuration interactions (CIs). Even when the states belong to different irreducible representations, degenerate orbital pairs cannot be excluded by symmetry. For various excited states of carbon monoxide and trigonal planar ammonia, we benchmark the method against the full CI and find it to be accurate. We create a semi-numerical approach and use it to show that the analytic approach is correct even when a high-symmetry structure is distorted to break symmetry so that near degeneracies in orbitals occur. For a range of geometries of trigonal planar ammonia, we find that the analytic non-adiabatic couplings for selected CI can achieve sufficient accuracy using a small fraction of the full CI space.
Collapse
Affiliation(s)
- Jeremy P. Coe
- Institute of Chemical Sciences, School
of Engineering and Physical Sciences, Heriot-Watt
University, Edinburgh EH14 4AS, U.K.
| |
Collapse
|
22
|
Wang Z, Zhang Z, Lu J, Li Y. Coordinate Descent Full Configuration Interaction for Excited States. J Chem Theory Comput 2023; 19:7731-7739. [PMID: 37870778 DOI: 10.1021/acs.jctc.3c00452] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2023]
Abstract
An efficient excited state method, named xCDFCI, in the configuration interaction framework is proposed. xCDFCI extends the unconstrained nonconvex optimization problem in coordinate descent full configuration interaction (CDFCI) to a multicolumn version for low-lying excited states computation. The optimization problem is addressed via a tailored coordinate descent method. In each iteration, a determinant is selected based on an approximated gradient, and coefficients of all states associated with the selected determinant are updated. A deterministic compression is applied to limit memory usage. We test xCDFCI applied to H2O and N2 molecules under the cc-pVDZ basis set. For both systems, five low-lying excited states in the same symmetry sector are calculated, together with the ground state. xCDFCI also produces accurate binding curves of the carbon dimer in the cc-pVDZ basis with chemical accuracy, where the ground state and four excited states in the same symmetry sector are benchmarked.
Collapse
Affiliation(s)
- Zhe Wang
- Department of Mathematics, Duke University, Durham, North Carolina 27708-0187, United States
| | - Zhiyuan Zhang
- School of Future Technology, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jianfeng Lu
- Department of Mathematics, Duke University, Durham, North Carolina 27708-0187, United States
- Department of Chemistry and Department of Physics, Duke University, Durham, North Carolina 27708-0187, United States
| | - Yingzhou Li
- School of Mathematical Sciences, Fudan University, Shanghai 200433, China
| |
Collapse
|
23
|
Li Manni G, Fdez. Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry. J Chem Theory Comput 2023; 19:6933-6991. [PMID: 37216210 PMCID: PMC10601490 DOI: 10.1021/acs.jctc.3c00182] [Citation(s) in RCA: 63] [Impact Index Per Article: 63.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Indexed: 05/24/2023]
Abstract
The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.
Collapse
Affiliation(s)
- Giovanni Li Manni
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ignacio Fdez. Galván
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Ali Alavi
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Yusuf Hamied
Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Flavia Aleotti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Francesco Aquilante
- Theory and
Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jochen Autschbach
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
| | - Davide Avagliano
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Alberto Baiardi
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jie J. Bao
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Stefano Battaglia
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Letitia Birnoschi
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Alejandro Blanco-González
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Sergey I. Bokarev
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
- Chemistry
Department, School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| | - Ria Broer
- Theoretical
Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Roberto Cacciari
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Paul B. Calio
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Rebecca K. Carlson
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Rafael Carvalho Couto
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luis Cerdán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
- Instituto
de Óptica (IO−CSIC), Consejo
Superior de Investigaciones Científicas, 28006, Madrid, Spain
| | - Liviu F. Chibotaru
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Nicholas F. Chilton
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | | | - Irene Conti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Sonia Coriani
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Juliana Cuéllar-Zuquin
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Razan E. Daoud
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Nike Dattani
- HPQC Labs, Waterloo, N2T 2K9 Ontario Canada
- HPQC College, Waterloo, N2T 2K9 Ontario Canada
| | - Piero Decleva
- Istituto
Officina dei Materiali IOM-CNR and Dipartimento di Scienze Chimiche
e Farmaceutiche, Università degli
Studi di Trieste, I-34121 Trieste, Italy
| | - Coen de Graaf
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
- ICREA, Pg. Lluís
Companys 23, 08010 Barcelona, Spain
| | - Mickaël
G. Delcey
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luca De Vico
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Werner Dobrautz
- Chalmers
University of Technology, Department of Chemistry
and Chemical Engineering, 41296 Gothenburg, Sweden
| | - Sijia S. Dong
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Chemical Biology, Department of Physics, and Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Rulin Feng
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Department
of Chemistry, Fudan University, Shanghai 200433, China
| | - Nicolas Ferré
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | | | - Laura Gagliardi
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Garavelli
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Yafu Guan
- State Key
Laboratory of Molecular Reaction Dynamics and Center for Theoretical
Computational Chemistry, Dalian Institute
of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Meiyuan Guo
- SSRL, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Matthew R. Hennefarth
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Chad E. Hoyer
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Miquel Huix-Rotllant
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | - Vishal Kumar Jaiswal
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Andy Kaiser
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Danil S. Kaliakin
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Marjan Khamesian
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Daniel S. King
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Vladislav Kochetov
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Marek Krośnicki
- Institute
of Theoretical Physics and Astrophysics, Faculty of Mathematics, Physics
and Informatics, University of Gdańsk, ul Wita Stwosza 57, 80-952, Gdańsk, Poland
| | | | - Ernst D. Larsson
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Susi Lehtola
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 University of Helsinki, Finland
| | - Marie-Bernadette Lepetit
- Condensed
Matter Theory Group, Institut Néel, CNRS UPR 2940, 38042 Grenoble, France
- Theory
Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409-1061, United States
| | - Pablo López Ríos
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Marcus Lundberg
- Department
of Chemistry − Ångström Laboratory, Uppsala University, SE-75120 Uppsala, Sweden
| | - Dongxia Ma
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | | | - Francesco Montorsi
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Maximilian Mörchen
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Artur Nenov
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Vu Ha Anh Nguyen
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Yoshio Nishimoto
- Graduate
School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Meagan S. Oakley
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Markus Oppel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Daniele Padula
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Riddhish Pandharkar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Quan Manh Phung
- Department
of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
- Institute
of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
| | - Felix Plasser
- Department
of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K.
| | - Gerardo Raggi
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Quantum
Materials and Software LTD, 128 City Road, London, EC1V 2NX, United Kingdom
| | - Elisa Rebolini
- Scientific
Computing Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Markus Reiher
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Ivan Rivalta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Daniel Roca-Sanjuán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Thies Romig
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Arta Anushirwan Safari
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Aitor Sánchez-Mansilla
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
| | - Andrew M. Sand
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Biochemistry, Butler University, Indianapolis, Indiana 46208, United States
| | - Igor Schapiro
- Institute
of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Thais R. Scott
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Javier Segarra-Martí
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Francesco Segatta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Dumitru-Claudiu Sergentu
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Laboratory
RA-03, RECENT AIR, A. I. Cuza University of Iaşi, RA-03 Laboratory (RECENT AIR), Iaşi 700506, Romania
| | - Prachi Sharma
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, USA
| | - Yinan Shu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Jakob K. Staab
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Tjerk P. Straatsma
- National
Center for Computational Sciences, Oak Ridge
National Laboratory, Oak Ridge, Tennessee 37831-6373, United States
- Department
of Chemistry and Biochemistry, University
of Alabama, Tuscaloosa, Alabama 35487-0336, United States
| | | | - Bruno Nunes Cabral Tenorio
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Donald G. Truhlar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Liviu Ungur
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Morgane Vacher
- Nantes
Université, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| | - Valera Veryazov
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Torben Arne Voß
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Oskar Weser
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Dihua Wu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Xuchun Yang
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - David Yarkony
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Chen Zhou
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - J. Patrick Zobel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Roland Lindh
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Uppsala
Center for Computational Chemistry (UC3), Uppsala University, PO Box 576, SE-751 23 Uppsala. Sweden
| |
Collapse
|
24
|
Park JW. Dynamic Correlation on the Adaptive Sampling Configuration Interaction (ASCI) Reference Function: ASCI-DSRG-MRPT2. J Chem Theory Comput 2023; 19:6263-6272. [PMID: 37611192 DOI: 10.1021/acs.jctc.3c00688] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/25/2023]
Abstract
A balanced description of static and dynamic electron correlations is at the heart of quantum chemical methods. To obtain accurate results in strongly correlated systems using wave-function-based methods, a large active space is necessary to ensure correct descriptions of static correlations. Correcting the results for dynamic correlations is also necessary. In this work, we present implementations of second-order perturbation theory for dynamic correlations based on the adaptive sampling configuration interaction self-consistent field (ASCI-SCF) method. In particular, we implemented spin-free driven similarity renormalization group second-order multireference perturbation theory (DSRG-MRPT2). The extrapolation of the ASCI + PT2 energy based on the relaxed Hamiltonian in DSRG-MRPT2 gives a reasonable approximation of DSRG-MRPT2 based on CASSCF. We demonstrate the application of the ASCI-DSRG-MRPT2 method in evaluations of the spin-state energy gaps in iron porphyrins, polyacenes, and periacenes along with the reaction energies of methane oxidation by FeO+ and electrocyclic ring formation in cethrene.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
25
|
Li S, Misiewicz JP, Evangelista FA. Intruder-free cumulant-truncated driven similarity renormalization group second-order multireference perturbation theory. J Chem Phys 2023; 159:114106. [PMID: 37712785 DOI: 10.1063/5.0159403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 08/14/2023] [Indexed: 09/16/2023] Open
Abstract
Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet-triplet gap of CH2, the potential energy curves of the XΣg+1 and AΣu+3 states of N2, and the singlet-triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5 kcal mol-1 or less. At the same time, the amount of memory required is reduced from O(NA6) to O(NA5), where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods.
Collapse
Affiliation(s)
- Shuhang Li
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Jonathon P Misiewicz
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
26
|
Gururangan K, Piecuch P. Converging high-level coupled-cluster energetics via adaptive selection of excitation manifolds driven by moment expansions. J Chem Phys 2023; 159:084108. [PMID: 37610021 DOI: 10.1063/5.0162873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Accepted: 08/03/2023] [Indexed: 08/24/2023] Open
Abstract
A novel approach to rapidly converging high-level coupled-cluster (CC) energetics in an automated fashion is proposed. The key idea is an adaptive selection of excitation manifolds defining higher--than--two-body components of the cluster operator inspired by CC(P;Q) moment expansions. The usefulness of the resulting methodology is illustrated by molecular examples where the goal is to recover the electronic energies obtained using the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT) when the noniterative triples corrections to CCSD fail.
Collapse
Affiliation(s)
- Karthik Gururangan
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
27
|
Magoulas I, Evangelista FA. Unitary Coupled Cluster: Seizing the Quantum Moment. J Phys Chem A 2023; 127:6567-6576. [PMID: 37523485 PMCID: PMC10424243 DOI: 10.1021/acs.jpca.3c02781] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2023] [Revised: 07/08/2023] [Indexed: 08/02/2023]
Abstract
Shallow, CNOT-efficient quantum circuits are crucial for performing accurate computational chemistry simulations on current noisy quantum hardware. Here, we explore the usefulness of noniterative energy corrections, based on the method of moments of coupled-cluster theory, for accelerating convergence toward full configuration interaction. Our preliminary numerical results relying on iteratively constructed ansätze suggest that chemically accurate energies can be obtained with substantially more compact circuits, implying enhanced resilience to gate and decoherence noise.
Collapse
Affiliation(s)
- Ilias Magoulas
- Department of Chemistry and
Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Francesco A. Evangelista
- Department of Chemistry and
Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
28
|
Malrieu JP, Heully JL. From complete to selected model spaces in determinant-based multi-reference second-order perturbation treatments. J Chem Phys 2023; 159:034110. [PMID: 37470425 DOI: 10.1063/5.0153416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2023] [Accepted: 06/20/2023] [Indexed: 07/21/2023] Open
Abstract
The present paper reformulates and improves a previously proposed determinant-based second-order multi-reference perturbative formalism. Through a rather simple modification of the energy denominators, this formalism takes into account the interactions between the model space determinants, which are repeated in outer space. The method has been shown to be size-consistent when the model space is a complete active space, which is a severe limit. It is shown here that the completeness of the model space is not necessary to keep this property, provided that the zero-order function satisfies some conditions. For instance, size consistency may be obtained from truncated complete active spaces. It may even be satisfied from Singles and Doubles Configuration Interactions, provided that a coupled electron pair approximation is used in the definition of the model space wave function. The physical content of the method is illustrated by a series of model problems, showing its robustness. A major benefit of the fact that the perturbers are single determinants is the possibility to revise with full flexibility the model-space component of the wave function, i.e., to treat the feedback effect of the dynamic correlation on the valence component of the wave function.
Collapse
|
29
|
Jones GM, Li RR, DePrince AE, Vogiatzis KD. Data-Driven Refinement of Electronic Energies from Two-Electron Reduced-Density-Matrix Theory. J Phys Chem Lett 2023:6377-6385. [PMID: 37418691 DOI: 10.1021/acs.jpclett.3c01382] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2023]
Abstract
The exponential computational cost of describing strongly correlated electrons can be mitigated by adopting a reduced-density matrix (RDM)-based description of the electronic structure. While variational two-electron RDM (v2RDM) methods can enable large-scale calculations on such systems, the quality of the solution is limited by the fact that only a subset of known necessary N-representability constraints can be applied to the 2RDM in practical calculations. Here, we demonstrate that violations of partial three-particle (T1 and T2) N-representability conditions, which can be evaluated with knowledge of only the 2RDM, can serve as physics-based features in a machine-learning (ML) protocol for improving energies from v2RDM calculations that consider only two-particle (PQG) conditions. Proof-of-principle calculations demonstrate that the model yields substantially improved energies relative to reference values from configuration-interaction-based calculations.
Collapse
Affiliation(s)
- Grier M Jones
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Run R Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | | |
Collapse
|
30
|
Filip MA, Thom AJW. A hybrid stochastic configuration interaction-coupled cluster approach for multireference systems. J Chem Phys 2023; 158:2889005. [PMID: 37154279 DOI: 10.1063/5.0145767] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Accepted: 04/19/2023] [Indexed: 05/10/2023] Open
Abstract
The development of multireference coupled cluster (MRCC) techniques has remained an open area of study in electronic structure theory for decades due to the inherent complexity of expressing a multiconfigurational wavefunction in the fundamentally single-reference coupled cluster framework. The recently developed multireference-coupled cluster Monte Carlo (mrCCMC) technique uses the formal simplicity of the Monte Carlo approach to Hilbert space quantum chemistry to avoid some of the complexities of conventional MRCC, but there is room for improvement in terms of accuracy and, particularly, computational cost. In this paper, we explore the potential of incorporating ideas from conventional MRCC-namely, the treatment of the strongly correlated space in a configuration interaction formalism-to the mrCCMC framework, leading to a series of methods with increasing relaxation of the reference space in the presence of external amplitudes. These techniques offer new balances of stability and cost against accuracy, as well as a means to better explore and better understand the structure of solutions to the mrCCMC equations.
Collapse
Affiliation(s)
- Maria-Andreea Filip
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
- Peterhouse, University of Cambridge, Cambridge, United Kingdom
| | - Alex J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
31
|
Haidar M, Rančić MJ, Ayral T, Maday Y, Piquemal J. Open source variational quantum eigensolver extension of the quantum learning machine for quantum chemistry. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/17/2023]
Affiliation(s)
- Mohammad Haidar
- Sorbonne Université Laboratoire de Chimie Théorique (UMR‐7616‐CNRS) Paris France
- Sorbonne Université, CNRS Université Paris Cité, Laboratoire Jacques Louis Lions (LJLL) Paris France
- TotalEnergies, Tour Coupole La Défense Paris France
| | | | - Thomas Ayral
- Atos Quantum Laboratory Les Clayes‐sous‐Bois France
| | - Yvon Maday
- Sorbonne Université, CNRS Université Paris Cité, Laboratoire Jacques Louis Lions (LJLL) Paris France
- Institut Universitaire de France Paris France
| | - Jean‐Philip Piquemal
- Sorbonne Université Laboratoire de Chimie Théorique (UMR‐7616‐CNRS) Paris France
| |
Collapse
|
32
|
Wang X, Sharma S. Relativistic Semistochastic Heat-Bath Configuration Interaction. J Chem Theory Comput 2023; 19:848-855. [PMID: 36700783 DOI: 10.1021/acs.jctc.2c01025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
In this work we present the extension of semistochastic heat-bath configuration interaction (SHCI) to work with any two-component and four-component Hamiltonian. The vertical detachment energy (VDE) of AuH2- and zero-field splitting (ZFS) of NpO22+ is calculated by correlating more than 100 spinors in both cases. This work demonstrates the capability of SHCI to treat problems where both relativistic effect and electron correlation are important.
Collapse
Affiliation(s)
- Xubo Wang
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado80309, United States
| | - Sandeep Sharma
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado80309, United States
| |
Collapse
|
33
|
Ghosh SK, Ghosh D. Machine learning matrix product state ansatz for strongly correlated systems. J Chem Phys 2023; 158:064108. [PMID: 36792489 DOI: 10.1063/5.0133399] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
Machine learning (ML) has been used to optimize the matrix product state (MPS) ansatz for the wavefunction of strongly correlated systems. The ML optimization of MPS has been tested for the Heisenberg Hamiltonian on one-dimensional and ladder lattices, which correspond to conjugated molecular systems. The input descriptors and output for the supervised ML are lattice configurations and configuration interaction coefficients, respectively. Efficient learning can be achieved from data over the full Hilbert space via exact diagonalization or full configuration interaction, as well as over a much smaller sub-space via Monte Carlo Configuration Interaction. We show that this circumvents the need to calculate energy and operator expectation values and is therefore a computationally efficient alternative to variational optimization.
Collapse
Affiliation(s)
- Sumanta K Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A and 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, India
| | - Debashree Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A and 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, India
| |
Collapse
|
34
|
Iino T, Shiozaki T, Yanai T. Algorithm for analytic nuclear energy gradients of state averaged DMRG-CASSCF theory with newly derived coupled-perturbed equations. J Chem Phys 2023; 158:054107. [PMID: 36754810 DOI: 10.1063/5.0130636] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
We present an algorithm for evaluating analytic nuclear energy gradients of the state-averaged density matrix renormalization group complete-active-space self-consistent field (SA-DMRG-CASSCF) theory based on the newly derived coupled-perturbed (CP) DMRG-CASSCF equations. The Lagrangian for the conventional SA-CASSCF analytic gradient theory is extended to the SA-DMRG-CASSCF variant that can fully consider a whole set of constraints on the parameters of multi-root canonical matrix product states formed at all the DMRG block configurations. An efficient algorithm to solve the CP-DMRG-CASSCF equations for determining the multipliers was developed. The complexity of the resultant analytic gradient algorithm is overall the same as that of the unperturbed SA-DMRG-CASSCF algorithm. In addition, a reduced-scaling approach was developed to directly compute the SA reduced density matrices (SA-RDMs) and their perturbed ones without calculating separate state-specific RDMs. As part of our implementation scheme, we neglect the term associated with the constraint on the active orbitals in terms of the active-active rotation in the Lagrangian. Thus, errors from the true analytic gradients may be caused in this scheme. The proposed gradient algorithm was tested with the spin-adapted implementation by checking how accurately the computed analytic energy gradients reproduce numerical gradients of the SA-DMRG-CASSCF energies using a common number of renormalized bases. The illustrative applications show that the errors are sufficiently small when using a typical number of the renormalized bases, which is required to attain adequate accuracy in DMRG's total energies.
Collapse
Affiliation(s)
- Tsubasa Iino
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| | - Toru Shiozaki
- Quantum Simulation Technologies, Inc., Boston, Massachusetts 02135, USA
| | - Takeshi Yanai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| |
Collapse
|
35
|
Abstract
We develop analytic gradients for selected configuration interaction wave functions. Despite all pairs of molecular orbitals now potentially having to be considered for the coupled perturbed Hartree-Fock equations, we show that degenerate orbital pairs belonging to different irreducible representations in the largest abelian subgroup do not need to be included and instabilities due to degeneracies are avoided. We introduce seminumerical gradients and use them to validate the analytic approach even when near degeneracies are present due to high-symmetry geometries being slightly distorted to break symmetry. The method is applied to carbon monoxide, ammonia, square planar H4, hexagonal planar H6, and methane for a range of bond lengths where we demonstrate that analytic gradients for selected configuration interaction can approach the quality of full configuration interaction yet only use a very small fraction of its determinants.
Collapse
|
36
|
Dang DK, Kammeraad JA, Zimmerman PM. Advances in Parallel Heat Bath Configuration Interaction. J Phys Chem A 2023; 127:400-411. [PMID: 36580361 DOI: 10.1021/acs.jpca.2c07949] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Heat-bath configuration interaction (HCI) is a deterministic method that approaches the full CI limit at greatly reduced computational cost. In this work, computational improvements to the HCI algorithm are introduced targeting speed, parallel efficiency, and memory requirements. The new implementation introduces a hash function to distribute determinants and takes advantage of MPI and OpenMP for parallelism allowing for a (22e,168o) active space to be studied, which explicitly includes 2.39 × 107 variational determinants and 8.95 × 1010 perturbative determinants. Benchmarks show up to 86% parallel efficiency of the perturbative step on 32 nodes (4096 cores) and a total efficiency of 74%. The new HCI implementation is benchmarked for accuracy against prior results and applied to study the triplet-quintet gap in the challenging [FeO(NH3)5]2+ complex.
Collapse
Affiliation(s)
- Duy-Khoi Dang
- Department of Chemistry, University of Michigan 930 North University Avenue Ann Arbor, Michigan 48109, United States
| | - Joshua A Kammeraad
- Department of Chemistry, University of Michigan 930 North University Avenue Ann Arbor, Michigan 48109, United States
| | - Paul M Zimmerman
- Department of Chemistry, University of Michigan 930 North University Avenue Ann Arbor, Michigan 48109, United States
| |
Collapse
|
37
|
Damour Y, Quintero-Monsebaiz R, Caffarel M, Jacquemin D, Kossoski F, Scemama A, Loos PF. Ground- and Excited-State Dipole Moments and Oscillator Strengths of Full Configuration Interaction Quality. J Chem Theory Comput 2023; 19:221-234. [PMID: 36548519 DOI: 10.1021/acs.jctc.2c01111] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We report ground- and excited-state dipole moments and oscillator strengths (computed in different "gauges" or representations) of full configuration interaction (FCI) quality using the selected configuration interaction method known as Configuration Interaction using a Perturbative Selection made Iteratively (CIPSI). Thanks to a set encompassing 35 ground- and excited-state properties computed in 11 small molecules, the present near-FCI estimates allow us to assess the accuracy of high-order coupled-cluster (CC) calculations including up to quadruple excitations. In particular, we show that incrementing the excitation degree of the CC expansion (from CC with singles and doubles (CCSD) to CC with singles, doubles, and triples (CCSDT) or from CCSDT to CC with singles, doubles, triples, and quadruples (CCSDTQ)) reduces the average error with respect to the near-FCI reference values by approximately 1 order of magnitude.
Collapse
Affiliation(s)
- Yann Damour
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Raúl Quintero-Monsebaiz
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Michel Caffarel
- 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
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- 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
|
38
|
Cotton SJ. A truncated Davidson method for the efficient “chemically accurate” calculation of full configuration interaction wavefunctions without any large matrix diagonalization. J Chem Phys 2022; 157:224105. [DOI: 10.1063/5.0115796] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
This work develops and illustrates a new method of calculating “chemically accurate” electronic wavefunctions (and energies) via a truncated full configuration interaction (CI) procedure, which arguably circumvents the large matrix diagonalization that is the core problem of full CI and is also central to modern selective CI approaches. This is accomplished simply by following the standard/ubiquitous Davidson method in its “direct” form—wherein, in each iteration, the electronic Hamiltonian operator is applied directly in second quantization to the Ritz vector/wavefunction from the prior iteration—except that (in this work) only a small portion of the resultant expansion vector is actually even computed (through the application of only a similarly small portion of the Hamiltonian). Specifically, at each iteration of this truncated Davidson approach, the new expansion vector is taken to be twice as large as that from the prior iteration. In this manner, a small set of highly truncated expansion vectors (say 10–30) of increasing precision is incrementally constructed, forming a small subspace within which diagonalization of the Hamiltonian yields clear, consistent, and monotonically variational convergence to the approximate full CI limit. The good efficiency in which convergence to the level of chemical accuracy (1.6 mhartree) is achieved suggests, at least for the demonstrated problem sizes—Hilbert spaces of 1018 and wavefunctions of 108 determinants—that this truncated Davidson methodology can serve as a replacement of standard CI and complete-active space approaches in circumstances where only a few chemically significant digits of accuracy are required and/or meaningful in view of ever-present basis set limitations.
Collapse
Affiliation(s)
- Stephen J. Cotton
- Quantum Artificial Intelligence Laboratory (QuAIL), NASA Ames Research Center, Moffett Field, California 94035, USA and Employed By KBR, 601 Jefferson St., Houston, Texas 77002, USA
| |
Collapse
|
39
|
Pandharkar R, Hermes MR, Cramer CJ, Gagliardi L. Localized Active Space-State Interaction: a Multireference Method for Chemical Insight. J Chem Theory Comput 2022; 18:6557-6566. [PMID: 36257065 DOI: 10.1021/acs.jctc.2c00536] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Multireference electronic structure methods, like the complete active space (CAS) self-consistent field model, have long been used to characterize chemically interesting processes. Important work has been done in recent years to develop modifications having a lower computational cost than CAS, but typically these methods offer no more chemical insight than that from the CAS solution being approximated. In this paper, we present the localized active space-state interaction (LASSI) method that can be used not only to lower the intrinsic cost of the multireference calculation but also to improve interpretability. The localized active space (LAS) approach utilizes the local nature of the electron-electron correlation to express a composite wave function as an antisymmetrized product of unentangled wave functions in local active subspaces. LASSI then uses these LAS states as a basis from which to express complete molecular wave functions. This not only makes the molecular wave function more compact but also permits flexibility in choosing those states to be included in the basis. Such selective inclusion of states translates to the selective inclusion of specific types of interactions, thereby allowing a quantitative analysis of these interactions. We demonstrate the use of LASSI to study charge migration and spin-flip excitations in multireference organic molecules. We also compute the J coupling parameter for a bimetallic compound using various LAS bases to construct the Hamiltonian to provide insights into the coupling mechanism.
Collapse
Affiliation(s)
- Riddhish Pandharkar
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois60637, United States.,Argonne National Laboratory, Lemont, Illinois60439, USA
| | - Matthew R Hermes
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois60637, United States
| | - Christopher J Cramer
- Underwriters Laboratories Inc., 333 Pfingsten Road., Northbrook, Illinois60062, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois60637, United States.,Argonne National Laboratory, Lemont, Illinois60439, USA
| |
Collapse
|
40
|
Lunghi A, Sanvito S. Computational design of magnetic molecules and their environment using quantum chemistry, machine learning and multiscale simulations. Nat Rev Chem 2022; 6:761-781. [PMID: 37118096 DOI: 10.1038/s41570-022-00424-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/15/2022] [Indexed: 11/09/2022]
Abstract
Having served as a playground for fundamental studies on the physics of d and f electrons for almost a century, magnetic molecules are now becoming increasingly important for technological applications, such as magnetic resonance, data storage, spintronics and quantum information. All of these applications require the preservation and control of spins in time, an ability hampered by the interaction with the environment, namely with other spins, conduction electrons, molecular vibrations and electromagnetic fields. Thus, the design of a novel magnetic molecule with tailored properties is a formidable task, which does not only concern its electronic structures but also calls for a deep understanding of the interaction among all the degrees of freedom at play. This Review describes how state-of-the-art ab initio computational methods, combined with data-driven approaches to materials modelling, can be integrated into a fully multiscale strategy capable of defining design rules for magnetic molecules.
Collapse
|
41
|
Coe JP, Moreno Carrascosa A, Simmermacher M, Kirrander A, Paterson MJ. Efficient Computation of Two-Electron Reduced Density Matrices via Selected Configuration Interaction. J Chem Theory Comput 2022; 18:6690-6699. [PMID: 36198067 PMCID: PMC9648180 DOI: 10.1021/acs.jctc.2c00738] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We create an approach to efficiently calculate two-electron
reduced
density matrices (2-RDMs) using selected configuration interaction
wavefunctions. This is demonstrated using the specific example of
Monte Carlo configuration interaction (MCCI). The computation of the
2-RDMs is accelerated by using ideas from fast implementations of
full configuration interaction (FCI) and recent advances in implementing
the Slater–Condon rules using hardware bitwise operations.
This method enables a comparison of MCCI and truncated CI 2-RDMs with
FCI values for a range of molecules, which includes stretched bonds
and excited states. The accuracy in energies, wavefunctions, and 2-RDMs
is seen to exhibit a similar behavior. We find that MCCI can reach
sufficient accuracy of the 2-RDM using significantly fewer configurations
than truncated CI, particularly for systems with strong multireference
character.
Collapse
Affiliation(s)
- Jeremy P Coe
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, EdinburghEH14 4AS, U.K
| | | | - Mats Simmermacher
- EaStCHEM, School of Chemistry, University of Edinburgh, EdinburghEH9 3FJ, U.K
| | - Adam Kirrander
- EaStCHEM, School of Chemistry, University of Edinburgh, EdinburghEH9 3FJ, U.K
| | - Martin J Paterson
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, EdinburghEH14 4AS, U.K
| |
Collapse
|
42
|
Reliably assessing the electronic structure of cytochrome P450 on today's classical computers and tomorrow's quantum computers. Proc Natl Acad Sci U S A 2022; 119:e2203533119. [PMID: 36095200 PMCID: PMC9499570 DOI: 10.1073/pnas.2203533119] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
Chemical simulation is one of the most promising applications for future quantum computers. It is thought that quantum computers may enable accurate simulation for complex molecules that are otherwise impossible to simulate classically; that is, it displays quantum advantage. To better understand quantum advantage in chemical simulation, we explore what quantum and classical resources are required to simulate a series of pharmaceutically relevant molecules. Using classical methods, we show that reliable classical simulation of these molecules requires significant resources and therefore is a promising candidate for quantum simulation. We estimate the quantum resources, both in overall simulation time and the size. The insights from this study pave the way for future quantum simulation of complex molecules. An accurate assessment of how quantum computers can be used for chemical simulation, especially their potential computational advantages, provides important context on how to deploy these future devices. To perform this assessment reliably, quantum resource estimates must be coupled with classical computations attempting to answer relevant chemical questions and to define the classical algorithms simulation frontier. Herein, we explore the quantum computation and classical computation resources required to assess the electronic structure of cytochrome P450 enzymes (CYPs) and thus define a classical–quantum advantage boundary. This is accomplished by analyzing the convergence of density matrix renormalization group plus n-electron valence state perturbation theory (DMRG+NEVPT2) and coupled-cluster singles doubles with noniterative triples [CCSD(T)] calculations for spin gaps in models of the CYP catalytic cycle that indicate multireference character. The quantum resources required to perform phase estimation using qubitized quantum walks are calculated for the same systems. Compilation into the surface code provides runtime estimates to compare directly to DMRG runtimes and to evaluate potential quantum advantage. Both classical and quantum resource estimates suggest that simulation of CYP models at scales large enough to balance dynamic and multiconfigurational electron correlation has the potential to be a quantum advantage problem and emphasizes the important interplay between classical computations and quantum algorithms development for chemical simulation.
Collapse
|
43
|
Zou J, Wang Q, Ren X, Wang Y, Zhang H, Li S. Efficient Implementation of Block-Correlated Coupled Cluster Theory Based on the Generalized Valence Bond Reference for Strongly Correlated Systems. J Chem Theory Comput 2022; 18:5276-5285. [PMID: 35922401 DOI: 10.1021/acs.jctc.2c00445] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
An optimized implementation of block-correlated coupled cluster theory based on the generalized valence bond wave function (GVB-BCCC) for the singlet ground state of strongly correlated systems is presented. The GVB-BCCC method with two-pair correlation (GVB-BCCC2b) or up to three-pair correlation (GVB-BCCC3b) will be the focus of this work. Three major techniques have been adopted to dramatically accelerate GVB-BCCC2b and GVB-BCCC3b calculations. First, the GVB-BCCC2b and GVB-BCCC3b codes are noticeably optimized by removing redundant calculations. Second, independent amplitudes are identified by constraining excited configurations to be pure singlet states and only independent amplitudes need to be solved. Third, an incremental updating scheme for the amplitudes in solving the GVB-BCCC equations is adopted. With these techniques, accurate GVB-BCCC3b calculations are now accessible for systems with relatively large active spaces (50 electrons in 50 orbitals) and GVB-BCCC2b calculations are affordable for systems with much larger active spaces. We have applied GVB-BCCC methods to investigate three typical kinds of systems: polyacenes, pentaprismane, and [Cu2O2]2+ isomers. For polyacenes, we demonstrate that GVB-BCCC3b can capture more than 94% of the total correlation energy even for 12-acene with 50 π electrons. For the potential energy curve of simultaneously stretching 15 C-C bonds in pentaprismane, our calculations show that the GVB-BCCC3b results are quite close to the results from the density matrix renormalization group (DMRG) over the whole range. For two dinuclear copper oxide isomers, their relative energy predicted by GVB-BCCC3b is also in good accord with the DMRG result. All calculations show that the inclusion of three-pair correlation in GVB-BCCC is critical for accurate descriptions of strongly correlated systems.
Collapse
Affiliation(s)
- Jingxiang Zou
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| | - Qingchun Wang
- Key Laboratory of Quantum Information, Chinese Academy of Sciences, School of Physics, University of Science and Technology of China, Hefei, Anhui 230026, People's Republic of China
| | - Xiaochuan Ren
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| | - Yuqi Wang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| | - Haodong Zhang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| | - Shuhua Li
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| |
Collapse
|
44
|
Duan C, Ladera AJ, Liu JCL, Taylor MG, Ariyarathna IR, Kulik HJ. Exploiting Ligand Additivity for Transferable Machine Learning of Multireference Character across Known Transition Metal Complex Ligands. J Chem Theory Comput 2022; 18:4836-4845. [PMID: 35834742 DOI: 10.1021/acs.jctc.2c00468] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Accurate virtual high-throughput screening (VHTS) of transition metal complexes (TMCs) remains challenging due to the possibility of high multireference (MR) character that complicates property evaluation. We compute MR diagnostics for over 5,000 ligands present in previously synthesized octahedral mononuclear transition metal complexes in the Cambridge Structural Database (CSD). To accomplish this task, we introduce an iterative approach for consistent ligand charge assignment for ligands in the CSD. Across this set, we observe that the MR character correlates linearly with the inverse value of the averaged bond order over all bonds in the molecule. We then demonstrate that ligand additivity of the MR character holds in TMCs, which suggests that the TMC MR character can be inferred from the sum of the MR character of the ligands. Encouraged by this observation, we leverage ligand additivity and develop a ligand-derived machine learning representation to train neural networks to predict the MR character of TMCs from properties of the constituent ligands. This approach yields models with excellent performance and superior transferability to unseen ligand chemistry and compositions.
Collapse
Affiliation(s)
- Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Adriana J Ladera
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Julian C-L Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Michael G Taylor
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Isuru R Ariyarathna
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
45
|
Smith JET, Lee J, Sharma S. Near-Exact Nuclear Gradients of Complete Active Space Self-Consistent Field Wave Functions. J Chem Phys 2022; 157:094104. [DOI: 10.1063/5.0085515] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this paper, we study the nuclear gradients of heat bath configuration interaction self-consistent field (HCISCF) wave functions and use them to optimize molecular geometries for various molecules.We show that the HCISCF nuclear gradients are fairly insensitive to the size of the "selected" variational space, which allows us to reduce the computational cost without introducing significant error.The ability of HCISCF to treat larger active spaces combined with the flexibility for users to control the computational cost makes the method very attractive for studying strongly correlated systems which require a larger active space than possible with complete active space self-consistent field (CASSCF).Finally, we study the realistic catalyst, Fe(PDI), and highlight some of the challenges this system poses for density functional theory (DFT).We demonstrate how HCISCF can clarify the energetic stability of geometries obtained from DFT when the results are strongly dependent on the functional.We also use the HCISCF gradients to optimize geometries for this species and study the adiabatic singlet-triplet gap. During geometry optimization, we find that multiple near-degenerate local minima exist on the triplet potential energy surface.
Collapse
Affiliation(s)
- James E. T. Smith
- Center for Computational Quantum Phyics, Flatiron Institute, United States of America
| | - Joonho Lee
- Chemistry, Columbia University Department of Chemistry, United States of America
| | - Sandeep Sharma
- University of Colorado at Boulder, United States of America
| |
Collapse
|
46
|
Baiardi A, Lesiuk M, Reiher M. Explicitly Correlated Electronic Structure Calculations with Transcorrelated Matrix Product Operators. J Chem Theory Comput 2022; 18:4203-4217. [PMID: 35666238 DOI: 10.1021/acs.jctc.2c00167] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, we present the first implementation of the transcorrelated electronic Hamiltonian in an optimization procedure for matrix product states by the density matrix renormalization group (DMRG) algorithm. In the transcorrelation ansatz, the electronic Hamiltonian is similarity-transformed with a Jastrow factor to describe the cusp in the wave function at electron-electron coalescence. As a result, the wave function is easier to approximate accurately with the conventional expansion in terms of one-particle basis functions and Slater determinants. The transcorrelated Hamiltonian in first quantization comprises up to three-body interactions, which we deal with in the standard way by applying robust density fitting to two- and three-body integrals entering the second-quantized representation of this Hamiltonian. The lack of hermiticity of the transcorrelated Hamiltonian is taken care of along the lines of the first work on transcorrelated DMRG [ J. Chem. Phys. 2020, 153, 164115] by encoding it as a matrix product operator and optimizing the corresponding ground state wave function with imaginary-time time-dependent DMRG. We demonstrate our quantum chemical transcorrelated DMRG approach at the example of several atoms and first-row diatomic molecules. We show that transcorrelation improves the convergence rate to the complete basis set limit in comparison to conventional DMRG. Moreover, we study extensions of our approach that aim at reducing the cost of handling the matrix product operator representation of the transcorrelated Hamiltonian.
Collapse
Affiliation(s)
- Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Michał Lesiuk
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland.,Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
47
|
Chen BJ, Chai JD. TAO-DFT fictitious temperature made simple. RSC Adv 2022; 12:12193-12210. [PMID: 35481082 PMCID: PMC9026342 DOI: 10.1039/d2ra01632j] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 04/04/2022] [Indexed: 11/21/2022] Open
Abstract
Over the past few years, thermally-assisted-occupation density functional theory (TAO-DFT) [J.-D. Chai, J. Chem. Phys., 2012, 136, 154104] has been proved to be an efficient electronic structure method for investigating the ground-state properties of large electronic systems with strong static correlation effects. In TAO-DFT, the strength of static correlation in an electronic system at zero temperature is closely related to the so-called fictitious temperature (i.e., the temperature of the corresponding noninteracting reference system). In this work, we propose a simple model to define the optimal system-independent fictitious temperature of a given energy functional in TAO-DFT. Besides, we employ this model to determine the optimal system-independent fictitious temperature of a global hybrid functional in TAO-DFT as a function of the fraction of exact exchange. In addition, we adopt TAO-DFT with various global hybrid functionals and system-independent fictitious temperatures to explore the ground-state properties of several electronic systems with strong static correlation effects, such as the linear acenes and cyclic carbon chains. Furthermore, we discuss the role of exact exchange and an optimal system-independent fictitious temperature in TAO-DFT. Owing to the much reduced self-interaction error, TAO-DFT with exact exchange and an optimal system-independent fictitious temperature can accurately predict the radical character and bond length alternation of cyclic carbon chains (with even number of carbon atoms), which are challenging problems for traditional electronic structure methods. Optimal system-independent fictitious temperature θ of TAO-GH as a function of the fraction of exact exchange ax.![]()
Collapse
Affiliation(s)
- Bo-Jyun Chen
- Department of Physics, National Taiwan University Taipei 10617 Taiwan
| | - Jeng-Da Chai
- Department of Physics, National Taiwan University Taipei 10617 Taiwan .,Center for Theoretical Physics, Center for Quantum Science and Engineering, National Taiwan University Taipei 10617 Taiwan.,Physics Division, National Center for Theoretical Sciences Taipei 10617 Taiwan
| |
Collapse
|
48
|
Magoulas I, Shen J, Piecuch P. Addressing strong correlation by approximate coupled-pair methods with active-space and full treatments of three-body clusters. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2057365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Ilias Magoulas
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
- Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA
| |
Collapse
|
49
|
Zhang N, Xiao Y, Liu W. SOiCI and iCISO: combining iterative configuration interaction with spin-orbit coupling in two ways. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:224007. [PMID: 35287124 DOI: 10.1088/1361-648x/ac5db4] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Accepted: 03/14/2022] [Indexed: 06/14/2023]
Abstract
The near-exact iCIPT2 approach for strongly correlated systems of electrons, which stems from the combination of iterative configuration interaction (iCI, an exact solver of full CI) with configuration selection for static correlation and second-order perturbation theory (PT2) for dynamic correlation, is extended to the relativistic domain. In the spirit of spin separation, relativistic effects are treated in two steps: scalar relativity is treated by the infinite-order, spin-free part of the exact two-component (X2C) relativistic Hamiltonian, whereas spin-orbit coupling (SOC) is treated by the first-order, Douglas-Kroll-Hess-like SOC operator derived from the same X2C Hamiltonian. Two possible combinations of iCIPT2 with SOC are considered, i.e., SOiCI and iCISO. The former treats SOC and electron correlation on an equal footing, whereas the latter treats SOC in the spirit of state interaction, by constructing and diagonalizing an effective spin-orbit Hamiltonian matrix in a small number of correlated scalar states. Both double group and time reversal symmetries are incorporated to simplify the computation. Pilot applications reveal that SOiCI is very accurate for the spin-orbit splitting (SOS) of heavy atoms, whereas the computationally very cheap iCISO can safely be applied to the SOS of light atoms and even of systems containing heavy atoms when SOC is largely quenched by ligand fields.
Collapse
Affiliation(s)
- Ning Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, People's Republic of China
| | - Yunlong Xiao
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, People's Republic of China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, People's Republic of China
| |
Collapse
|
50
|
Rask AE, Zimmerman PM. The many-body electronic interactions of Fe(II)–porphyrin. J Chem Phys 2022; 156:094110. [DOI: 10.1063/5.0079310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Fe(II)–porphyrin complexes exhibit a diverse range of electronic interactions between the metal and macrocycle. Herein, the incremental full configuration interaction method is applied to the entire space of valence orbitals of a Fe(II)–porphyrin model using a modest basis set. A novel visualization framework is proposed to analyze individual many-body contributions to the correlation energy, providing detailed maps of this complex’s highly correlated electronic structure. This technique is used to parse the numerous interactions of two low-lying triplet states (3A2g and 3Eg) and to show that strong metal d–d and macrocycle π–π orbital interactions preferentially stabilize the 3A2g state. d–π interactions, on the other hand, preferentially stabilize the 3Eg state and primarily appear when correlating six electrons at a time. Ultimately, the Fe(II)–porphyrin model’s full set of 88 valence electrons are correlated in 275 orbitals, showing the interactions up to the 4-body level, which covers the great majority of correlations in this system.
Collapse
Affiliation(s)
- A. E. Rask
- Department of Chemistry, University of Michigan, 930 N. University Ave., Ann Arbor, Michigan 48109, USA
| | - P. M. Zimmerman
- Department of Chemistry, University of Michigan, 930 N. University Ave., Ann Arbor, Michigan 48109, USA
| |
Collapse
|