1
|
Li HE, Li X, Huang JC, Zhang GZ, Shen ZP, Zhao C, Li J, Hu HS. Variational quantum imaginary time evolution for matrix product state Ansatz with tests on transcorrelated Hamiltonians. J Chem Phys 2024; 161:144104. [PMID: 39377325 DOI: 10.1063/5.0228731] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2024] [Accepted: 09/19/2024] [Indexed: 10/09/2024] Open
Abstract
The matrix product state (MPS) Ansatz offers a promising approach for finding the ground state of molecular Hamiltonians and solving quantum chemistry problems. Building on this concept, the proposed technique of quantum circuit MPS (QCMPS) enables the simulation of chemical systems using a relatively small number of qubits. In this study, we enhance the optimization performance of the QCMPS Ansatz by employing the variational quantum imaginary time evolution (VarQITE) approach. Guided by McLachlan's variational principle, the VarQITE method provides analytical metrics and gradients, resulting in improved convergence efficiency and robustness of the QCMPS. We validate these improvements numerically through simulations of H2, H4, and LiH molecules. In addition, given that VarQITE is applicable to non-Hermitian Hamiltonians, we evaluate its effectiveness in preparing the ground state of transcorrelated Hamiltonians. This approach yields energy estimates comparable to the complete basis set (CBS) limit while using even fewer qubits. In particular, we perform simulations of the beryllium atom and LiH molecule using only three qubits, maintaining high fidelity with the CBS ground state energy of these systems. This qubit reduction is achieved through the combined advantages of both the QCMPS Ansatz and transcorrelation. Our findings demonstrate the potential practicality of this quantum chemistry algorithm on near-term quantum devices.
Collapse
Affiliation(s)
- Hao-En Li
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - 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
| | - Zhu-Ping Shen
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Chen Zhao
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Jun Li
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
- Department of Chemistry and Guangdong Provincial Key Laboratory of Catalytic Chemistry, Southern University of Science and Technology, Shenzhen 518055, China
- Fundamental Science Center of Rare Earths, Ganjiang Innovation Academy, Chinese Academy of Sciences, Ganzhou 341000, 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
|
2
|
Mitchell E, Turney JM, Schaefer HF. Automatic Differentiation for Explicitly Correlated MP2. J Chem Theory Comput 2024; 20:8529-8538. [PMID: 39311755 PMCID: PMC11465469 DOI: 10.1021/acs.jctc.4c00818] [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/23/2024] [Revised: 08/26/2024] [Accepted: 09/04/2024] [Indexed: 10/09/2024]
Abstract
Automatic differentiation (AD) offers a route to achieve arbitrary-order derivatives of challenging wave function methods without the use of analytic gradients or response theory. Currently, AD has been predominantly used in methods where first- and/or second-order derivatives are available, but it has not been applied to methods lacking available derivatives. The most robust approximation of explicitly correlated MP2, MP2-F12/3C(FIX)+CABS, is one such method. By comparing the results of MP2-F12 computed with AD versus finite-differences, it is shown that (a) optimized geometries match to about 10-3 Å for bond lengths and a 10-6 degree for angles, and (b) dipole moments match to about 10-6 D. Hessians were observed to have poorer agreement with numerical results (10-5), which is attributed to deficiencies in AD implementations currently. However, it is notable that vibrational frequencies match within 10-2 cm-1. The use of AD also allowed the prediction of MP2-F12/3C(FIX)+CABS IR intensities for the first time.
Collapse
Affiliation(s)
- Erica
C. Mitchell
- Department
of Chemistry, University of Georgia, 302 East Campus Road, Athens, Georgia 30602, United States
- Center
for Computational Quantum Chemistry, University
of Georgia, Athens, Georgia 30602, United States
| | - Justin M. Turney
- Center
for Computational Quantum Chemistry, University
of Georgia, Athens, Georgia 30602, United States
| | - Henry F. Schaefer
- Department
of Chemistry, University of Georgia, 302 East Campus Road, Athens, Georgia 30602, United States
- Center
for Computational Quantum Chemistry, University
of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
3
|
Di Grande S, Lazzari F, Barone V. Accurate Geometries of Large Molecules at DFT Cost by Semiexperimental and Coupled Cluster Templating Fragments. J Chem Theory Comput 2024. [PMID: 39373710 DOI: 10.1021/acs.jctc.4c00900] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/08/2024]
Abstract
Accurate geometries of small semirigid molecules in the gas phase are available thanks to high-resolution spectroscopy and accurate quantum chemical approaches. These results can be employed for validating cheaper low-level quantum chemical models or correcting the corresponding structures of large molecules. On these grounds, in this work, a large panel of semiexperimental equilibrium structures already available in the literature is used to confirm the average error (1 mÅ for bond lengths and 2 mrad for valence angles) of a version of the Pisa composite schemes (PCS2), which is applicable to molecules containing up to about 20 atoms. Then, the geometries of 30 additional medium-sized systems were optimized at the PCS2 level to cover a more balanced chemical space containing moieties poorly represented in SE compilations. The final database is available on a public domain Web site (https://www.skies-village.it/databases/) and can be employed for correcting structures of larger molecules obtained by hybrid or double-hybrid density functionals in the framework of the templating molecule approach. Several examples show that corrections based on the structures of building blocks taken from this database reduce the error of the B3LYP geometrical parameters of large molecules by nearly an order of magnitude without increasing the computational cost. Furthermore, the results of different density functional theory (DFT) or wave function (e.g., MP2) models can be improved in the same way by simply computing both the whole molecule and suitable building blocks at the chosen level. Then, whenever reference structures of some building blocks containing up to about 20 atoms are not available, they can be purposely optimized at the PCS2 level by employing reasonable computer resources. Therefore, a new DFT-cost tool is now available for the accurate characterization of large molecules by experiment-oriented scientists.
Collapse
Affiliation(s)
- Silvia Di Grande
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Scuola Superiore Meridionale, Largo San Marcellino 10, 80138 Napoli, Italy
| | - Federico Lazzari
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | |
Collapse
|
4
|
Liepuoniute I, Motta M, Pellegrini T, Rice JE, Gujarati TP, Gil S, Jones GO. Simulation of a Diels-Alder reaction on a quantum computer. Phys Chem Chem Phys 2024; 26:25181-25191. [PMID: 39314194 DOI: 10.1039/d4cp01314j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/25/2024]
Abstract
The simulation of chemical reactions is an anticipated application of quantum computers. Using a Diels-Alder reaction as a test case, in this study we explore the potential applications of quantum algorithms and hardware in investigating chemical reactions. Our specific goal is to calculate the activation barrier of a reaction between ethylene and cyclopentadiene forming a transition state. To achieve this goal, we use quantum algorithms for near-term quantum hardware (entanglement forging and quantum subspace expansion) and classical post-processing (many-body perturbation theory) in concert. We conduct simulations on IBM quantum hardware using up to 8 qubits, and compute accurate activation barrier in the reaction between cyclopentadiene and ethylene by accounting for both static and dynamic electronic correlation. This work illustrates a hybrid quantum-classical computational workflow to study chemical reactions on near-term quantum devices, showcasing the potential for performing quantum chemistry simulations on quantum hardware to predict activation barriers in agreement with those predicted by CASCI.
Collapse
Affiliation(s)
- Ieva Liepuoniute
- IBM Quantum, IBM Research - Almaden, 650 Harry Road, San Jose, CA 95120, USA.
| | - Mario Motta
- IBM Quantum, IBM Research - Almaden, 650 Harry Road, San Jose, CA 95120, USA.
- IBM Quantum, T. J. Watson Research Center, Yorktown Heights, NY 10598, USA
| | | | - Julia E Rice
- IBM Quantum, IBM Research - Almaden, 650 Harry Road, San Jose, CA 95120, USA.
| | - Tanvi P Gujarati
- IBM Quantum, IBM Research - Almaden, 650 Harry Road, San Jose, CA 95120, USA.
| | - Sofia Gil
- Cornell University, Ithaca, NY 14850, USA
| | - Gavin O Jones
- IBM Quantum, IBM Research - Almaden, 650 Harry Road, San Jose, CA 95120, USA.
| |
Collapse
|
5
|
Chan GKL. Spiers Memorial Lecture: Quantum chemistry, classical heuristics, and quantum advantage. Faraday Discuss 2024. [PMID: 39258407 DOI: 10.1039/d4fd00141a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/12/2024]
Abstract
We describe the problems of quantum chemistry, the intuition behind classical heuristic methods used to solve them, a conjectured form of the classical complexity of quantum chemistry problems, and the subsequent opportunities for quantum advantage. This article is written for both quantum chemists and quantum information theorists. In particular, we attempt to summarize the domain of quantum chemistry problems as well as the chemical intuition that is applied to solve them within concrete statements (such as a classical heuristic cost conjecture) in the hope that this may stimulate future analysis.
Collapse
Affiliation(s)
- Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California, USA
| |
Collapse
|
6
|
Luu D, Corminboeuf C, Patkowski K. Range Separation of the Interaction Potential in Intermolecular and Intramolecular Symmetry-Adapted Perturbation Theory. J Chem Theory Comput 2024. [PMID: 39255506 DOI: 10.1021/acs.jctc.4c00608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/12/2024]
Abstract
Symmetry-adapted perturbation theory (SAPT) is a popular and versatile tool to compute and decompose noncovalent interaction energies between molecules. The intramolecular SAPT (ISAPT) variant provides a similar energy decomposition between two nonbonded fragments of the same molecule, covalently connected by a third fragment. In this work, we explore an alternative approach where the noncovalent interaction is singled out by a range separation of the Coulomb potential. We investigate two common splittings of the 1/r potential into long-range and short-range parts based on the Gaussian and error functions, and approximate either the entire intermolecular/interfragment interaction or only its attractive terms by the long-range contribution. These range separation schemes are tested for a number of intermolecular and intramolecular complexes. We find that the energy corrections from range-separated SAPT or ISAPT are in reasonable agreement with complete SAPT/ISAPT data. This result should be contrasted with the inability of the long-range multipole expansion to describe crucial short-range charge penetration and exchange effects; it shows that the long-range interaction potential does not just recover the asymptotic interaction energy but also provides a useful account of short-range terms. The best consistency is attained for the error-function separation applied to all interaction terms, both attractive and repulsive. This study is the first step toward a fragmentation-free decomposition of intramolecular nonbonded energy.
Collapse
Affiliation(s)
- Du Luu
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| | - Clemence Corminboeuf
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne, Lausanne CH-1015, Switzerland
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
7
|
Cioslowski J, Strasburger K. Plain Capping for Improved Accuracy of Approximate One- and Two-Electron Densities at Two-Particle Coalescence Points. J Chem Theory Comput 2024. [PMID: 39225495 PMCID: PMC11428139 DOI: 10.1021/acs.jctc.4c00533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
The values of the one-electron and intracule densities at two-particle coalescence points that enter the expressions for relativistic corrections to energies of Coulombic systems cannot be efficiently computed with sufficient accuracy from approximate wave functions expressed in terms of cuspless basis functions such as the explicitly correlated Gaussians. A new approach to alleviation of this problem, called plain capping, is proposed. Unlike those offered by the previously published formalisms, such as the expectation value identities and integral transforms, the accuracy improvements effected by the plain capping are attained with negligible computational effort and minimum programming. In the case of the on-top two-electron densities, whose accurate computation is particularly costly, the plain capping constitutes the only viable means of error reduction available at present.
Collapse
Affiliation(s)
- Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
- Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Straaae 38, 01187 Dresden, Germany
| | - Krzysztof Strasburger
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| |
Collapse
|
8
|
Díaz Soto LJ, Oliveira RR, Baptista L, da Silveira EF, Nascimento MAC. Energy and spectroscopic parameters of neutral and cations isomers of the C nH 2 (n = 2-6) families using high-level ab-initio approaches. J Comput Chem 2024. [PMID: 39177429 DOI: 10.1002/jcc.27485] [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/10/2024] [Revised: 07/27/2024] [Accepted: 08/02/2024] [Indexed: 08/24/2024]
Abstract
Cationic species, previously detected from ion-induced desorption of solid methane by plasma desorption mass spectrometry (PDMS), and neutral species, are investigated using high-level ab-initio approaches. From a set of 25 cationic and 26 neutral structures belonging to CnH2 (n = 2-6) families, it was obtained the energy, rotational constants, harmonic vibrational frequency, charge distribution and excitation energies. The ZPVE-corrected energies, at CCSD(T)-F12; CCSD(T)-F12/RI/(cc-pVTZ-F12, cc-pVTZ-F12-CABS, cc-pVQZ/C) (n = 2-5) and CCSD(T)/cc-pVTZ (n = 6) levels, reveal that the topology of the most stable isomer vary with n and the charge. Out of 674 harmonic frequencies, those with maximum intensity are generally in the 3000-3500 cm-1 range. Analysis of 169 vertical transition energies calculated with the EOM-CCSD approach, suggest three C6H2 species as potential carriers of the diffuse interstellar bands (DIB). Systematic comparison of properties between neutral and cationic species can assist in the structural description of complex matrices.
Collapse
Affiliation(s)
- Lenin J Díaz Soto
- Instituto de Química, Universidade Federal do Rio de Janeiro, Cidade Universitária, CT, Rio de Janeiro, Brazil
| | - Ricardo R Oliveira
- Instituto de Química, Universidade Federal do Rio de Janeiro, Cidade Universitária, CT, Rio de Janeiro, Brazil
| | - Leonardo Baptista
- Departamento de Química e Ambiental, Faculdade de Tecnologia, Universidade do Estado do Rio de Janeiro, Resende, Brazil
| | - Enio F da Silveira
- Departamento de Física, Pontifícia Universidade Católica do Rio de Janeiro, Rio de Janeiro, Brazil
| | | |
Collapse
|
9
|
Szenes K, Mörchen M, Fischill P, Reiher M. Striking the right balance of encoding electron correlation in the Hamiltonian and the wavefunction ansatz. Faraday Discuss 2024. [PMID: 39092888 DOI: 10.1039/d4fd00060a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
Multi-configurational electronic structure theory delivers the most versatile approximations to many-electron wavefunctions, flexible enough to deal with all sorts of transformations, ranging from electronic excitations, to open-shell molecules and chemical reactions. Multi-configurational models are therefore essential to establish universally applicable, predictive ab initio methods for chemistry. Here, we present a discussion of explicit correlation approaches which address the nagging problem of dealing with static and dynamic electron correlation in multi-configurational active-space approaches. We review the latest developments and then point to their key obstacles. Our discussion is supported by new data obtained with tensor network methods. We argue in favor of simple electron-only correlator expressions that may allow one to define transcorrelated models in which the correlator does not bear a dependence on molecular structure.
Collapse
Affiliation(s)
- Kalman Szenes
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland.
| | - Maximilian Mörchen
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland.
| | - Paul Fischill
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland.
| | - Markus Reiher
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland.
| |
Collapse
|
10
|
Kats D, Christlmaier EMC, Schraivogel T, Alavi A. Orbital optimisation in xTC transcorrelated methods. Faraday Discuss 2024. [PMID: 39072553 DOI: 10.1039/d4fd00036f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
We present a combination of the bi-orthogonal orbital optimisation framework with the recently introduced xTC version of transcorrelation. This allows us to implement non-iterative perturbation based methods on top of the transcorrelated Hamiltonian. Additionally, the orbital optimisation influences results of other truncated methods, such as the distinguishable cluster with singles and doubles. The accuracy of these methods in comparison to standard xTC methods is demonstrated, and the advantages and disadvantages of the orbital optimisation are discussed.
Collapse
Affiliation(s)
- Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany.
| | | | - Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany.
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany.
| |
Collapse
|
11
|
Masios N, Hummel F, Grüneis A, Irmler A. Investigating the Basis Set Convergence of Diagrammatically Decomposed Coupled-Cluster Correlation Energy Contributions for the Uniform Electron Gas. J Chem Theory Comput 2024; 20:5937-5950. [PMID: 38976839 PMCID: PMC11270826 DOI: 10.1021/acs.jctc.4c00224] [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/22/2024] [Revised: 05/28/2024] [Accepted: 05/29/2024] [Indexed: 07/10/2024]
Abstract
We investigate the convergence of coupled-cluster (CC) correlation energies and related quantities with respect to the employed basis set size for the uniform electron gas (UEG) to gain a better understanding of the basis set incompleteness error (BSIE). To this end, coupled-cluster doubles (CCD) theory is applied to the three-dimensional UEG for a range of densities, basis set sizes, and electron numbers. We present a detailed analysis of individual diagrammatically decomposed contributions to the amplitudes at the level of CCD theory. In particular, we show that only two terms from the amplitude equations contribute to the asymptotic large-momentum behavior of the transition structure factor, corresponding to the cusp region at short interelectronic distances. However, due to the coupling present in the amplitude equations, all decomposed correlation energy contributions show the same asymptotic convergence behavior to the complete basis set limit. These findings provide an additional rationale for the success of a recently proposed correction to the BSIE of CC theory. Lastly, we examine the BSIE in the CCD plus perturbative triples [CCD(T)] method, as well as in the newly proposed CCD plus complete perturbative triples [CCD(cT)] method.
Collapse
Affiliation(s)
- Nikolaos Masios
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Felix Hummel
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Andreas Grüneis
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Andreas Irmler
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| |
Collapse
|
12
|
Shang Y, Luo SN. Insights into the role of the H-abstraction reaction kinetics of amines in understanding their degeneration fates under atmospheric and combustion conditions. Phys Chem Chem Phys 2024. [PMID: 39028293 DOI: 10.1039/d4cp02187h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
Amines, a class of prototypical volatile organic compounds, have garnered considerable interest within the context of atmospheric and combustion chemistry due to their substantial contributions to the formation of hazardous pollutants in the atmosphere. In the current energy landscape, the implementation of carbon-neutral energy and strategic initiatives leads to generation of new amine sources that cannot be overlooked in terms of the emission scale. To reduce the emission level of amines from their sources and mitigate their impact on the formation of harmful substances, a comprehensive understanding of the fundamental reaction kinetics during the degeneration process of amines is imperative. This perspective article first presents an overview of both traditional amine sources and emerging amine sources within the context of carbon peaking and carbon neutrality and then highlights the importance of H-abstraction reactions in understanding the atmospheric and combustion chemistry of amines from the perspective of reaction kinetics. Subsequently, the current experimental and theoretical techniques for investigating the H-abstraction reactions of amines are introduced, and a concise summary of research endeavors made in this field over the past few decades is provided. In order to provide accurate kinetic parameters of the H-abstraction reactions of amines, advanced kinetic calculations are performed using the multi-path canonical variational theory combined with the small-curvature tunneling and specific-reaction parameter methods. By comparing with the literature data, current kinetic calculations are comprehensively evaluated, and these validated data are valuable for the development of the reaction mechanism of amines.
Collapse
Affiliation(s)
- Yanlei Shang
- Energy Research Institute, Qilu University of Technology (Shandong Academy of Sciences), Jinan, Shandong, 250014, P. R. China.
- School of Materials Science and Engineering, Southwest Jiaotong University, Chengdu, Sichuan, 610031, P. R. China
- Key Laboratory of Extreme Material Dynamics Technology, Chengdu, Sichuan 610031, P. R. China
| | - S N Luo
- School of Materials Science and Engineering, Southwest Jiaotong University, Chengdu, Sichuan, 610031, P. R. China
- Key Laboratory of Extreme Material Dynamics Technology, Chengdu, Sichuan 610031, P. R. China
| |
Collapse
|
13
|
Liao K, Ding L, Schilling C. Quantum Information Orbitals (QIO): Unveiling Intrinsic Many-Body Complexity by Compressing Single-Body Triviality. J Phys Chem Lett 2024; 15:6782-6790. [PMID: 38913404 DOI: 10.1021/acs.jpclett.4c01105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
The simultaneous treatment of static and dynamic correlations in strongly correlated electron systems is a critical challenge. In particular, finding a universal scheme for identifying a single-particle orbital basis that minimizes the representational complexity of the many-body wave function is a formidable and longstanding problem. As a contribution toward its solution, we show that the total orbital correlation actually reveals and quantifies the intrinsic complexity of the wave function, once it is minimized via orbital rotations. To demonstrate the power of this concept in practice, an iterative scheme is proposed to optimize the orbitals by minimizing the total orbital correlation calculated by the tailored coupled cluster singles and doubles (TCCSD) ansatz. The optimized orbitals enable the limited TCCSD ansatz to capture more nontrivial information on the many-body wave function, indicated by the improved wave function and energy. An initial application of this scheme shows great improvement of TCCSD in predicting the singlet ground state potential energy curves of the strongly correlated C2 and Cr2 molecule.
Collapse
Affiliation(s)
- Ke Liao
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Lexin Ding
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Christian Schilling
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| |
Collapse
|
14
|
Carlson CD, Ma J, Al-Jabiri MH, Insausti A, Xu Y. Conformational adaptation and large amplitude motions of 1-phenyl-2,2,2-trifluoroethanol with two water molecules: a rotational spectroscopic and ab initio investigation. Phys Chem Chem Phys 2024; 26:18067-18075. [PMID: 38895791 DOI: 10.1039/d4cp01516a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
The 1 : 2 adduct of 1-phenyl-2,2,2-trifluoroethanol (PhTFE), a chiral fluoroalcohol, with two water molecules (PhTFE⋯2H2O) was investigated via chirped pulse Fourier-transform microwave (CP-FTMW) spectroscopy and theoretical calculations. A systematic search of the PhTFE⋯2H2O conformational landscape identified 38 stable minima at the B3LYP-D3BJ/def2-TZVPPD level of theory, 27 of which are within an energy window of 10 kJ mol-1 after applying zero-point energy corrections. Rotational spectra of a single PhTFE⋯2H2O conformer along with eight deuterated and three oxygen-18 isotopologues were assigned. Interestingly, the observed PhTFE⋯2H2O conformer contains PhTFE II, the second most stable monomer conformer, and the most stable PhTFE I dihydrate is ca. 4 kJ mol-1 higher in energy. In contrast, PhTFE I⋯H2O was identified experimentally and theoretically as the most stable 1 : 1 conformer. Furthermore, the observed dihydrate structure experiences large amplitude motions connecting three theoretical minima which differ only in which water oxygen lone pairs are involved in the hydrogen-bonds, i.e., the free OH pointing directions. Additionally, the ortho and para-H2O tunnelling splittings were detected and attributed to the interchange water hydrogen atoms which interact with the aromatic part of PhTFE but not for the water interacting with PhTFE hydroxy group. Extensive theoretical modelling was carried out to gain insight into the associated large amplitude motions including tunnelling, supported by the experimental isotopic and tunnelling splitting data.
Collapse
Affiliation(s)
- Colton D Carlson
- Department of Chemistry, University of Alberta, Edmonton, AB, T6G 2G2, Canada.
| | - Jiarui Ma
- Department of Chemistry, University of Alberta, Edmonton, AB, T6G 2G2, Canada.
| | - Mohamad H Al-Jabiri
- Department of Chemistry, University of Alberta, Edmonton, AB, T6G 2G2, Canada.
| | - Aran Insausti
- Department of Chemistry, University of Alberta, Edmonton, AB, T6G 2G2, Canada.
- Departamento de Quimica Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco (UPV-EHU), Bilbao 48080, Spain
- Biofisika Institute (CSIC, UPV/EHU), University of the Basque Country (UPV/EHU), Leioa 48940, Spain
| | - Yunjie Xu
- Department of Chemistry, University of Alberta, Edmonton, AB, T6G 2G2, Canada.
| |
Collapse
|
15
|
Plamper D, Vincent A, Fujioka K, Sun R, Weitzel KM. Ion molecule reactions in the HBr + + CH 4 system: a combined experimental and theoretical study. Phys Chem Chem Phys 2024; 26:16732-16746. [PMID: 38814257 DOI: 10.1039/d4cp01121j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2024]
Abstract
Reactions in the system HBr+ + CH4 have been investigated inside a guided ion-beam apparatus under single-collision conditions. The HBr+ is vibrational and rotational state selected in the electronic X2Π1/2 state created by (2+1)-REMPI. Due to the exitation scheme employed different rotational states of the HBr+ are accessible. Four reaction channels have been observed. The cross section, σ, for the exothermic proton transfer channel (PT) decreases with increasing collision energy, steeper than predicted by the Langevin model. The cross section also decreases with increasing rotational energy in the HBr+, with the effect of the rotational energy being stronger than that of translational energy. The cross section for the endothermic charge transfer (CT) increased with increasing collision energy. The energy dependence is well reproduced by a simple line of center (loc) model. Although the bromine transfer (BT) is exothermic the observed cross section increased with increasing collision energy due to an activation barrier on the potential energy surface (PES). Analysis by a modified loc model suggest the relevance of an angle dependence of σ. The cross section for the endothermic hydrogen atom abstraction (HA) exhibits a maximum at 2 eV Ecm. The measured cross sections are rationalized by means of reaction dynamics simulations which show good agreement with the experimental cross sections. The dynamics simulations are carried out with a machine learning potential that is developed and benchmarked with ab initio molecular dynamics simulation. The absolute cross sections predicted by reaction dynamics simulations are well within the same order of magnitude while reproducing the trends over three different collision energies for all four reaction channels. Furthermore, the simulations demonstrate various reaction mechanisms for these reaction channels, including a very interesting HBr+ orientation selectivity for the BT reaction channel.
Collapse
Affiliation(s)
- Dominik Plamper
- Philipps-Universität Marburg, Fachbereich Chemie, 35032 Marburg, Germany.
| | - Allen Vincent
- Department of Chemistry, University of Hawai'i at Manoa, Honolulu, Hawaii 96822, USA.
| | - Kazuumi Fujioka
- Department of Chemistry, University of Hawai'i at Manoa, Honolulu, Hawaii 96822, USA.
| | - Rui Sun
- Department of Chemistry, University of Hawai'i at Manoa, Honolulu, Hawaii 96822, USA.
| | | |
Collapse
|
16
|
Heßelmann A, Giner E, Reinhardt P, Knowles PJ, Werner HJ, Toulouse J. A density-fitting implementation of the density-based basis-set correction method. J Comput Chem 2024; 45:1247-1253. [PMID: 38348951 DOI: 10.1002/jcc.27325] [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: 09/19/2023] [Revised: 12/13/2023] [Accepted: 01/26/2024] [Indexed: 04/19/2024]
Abstract
This work reports an efficient density-fitting implementation of the density-based basis-set correction (DBBSC) method in the MOLPRO software. This method consists in correcting the energy calculated by a wave-function method with a given basis set by an adapted basis-set correction density functional incorporating the short-range electron correlation effects missing in the basis set, resulting in an accelerated convergence to the complete-basis-set limit. Different basis-set correction density-functional approximations are explored and the complementary-auxiliary-basis-set single-excitation correction is added. The method is tested on a benchmark set of reaction energies at the second-order Møller-Plesset (MP2) level and a comparison with the explicitly correlated MP2-F12 method is provided. The results show that the DBBSC method greatly accelerates the basis convergence of MP2 reaction energies, without reaching the accuracy of the MP2-F12 method but with a lower computational cost.
Collapse
Affiliation(s)
- Andreas Heßelmann
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Emmanuel Giner
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, Paris, France
| | - Peter Reinhardt
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, Paris, France
| | | | - Hans-Joachim Werner
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Julien Toulouse
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, Paris, France
- Institut Universitaire de France, Paris, France
| |
Collapse
|
17
|
Dobrautz W, Sokolov IO, Liao K, Ríos PL, Rahm M, Alavi A, Tavernelli I. Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method. J Chem Theory Comput 2024; 20:4146-4160. [PMID: 38723159 PMCID: PMC11137825 DOI: 10.1021/acs.jctc.4c00070] [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/19/2024] [Revised: 04/15/2024] [Accepted: 04/16/2024] [Indexed: 05/29/2024]
Abstract
Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.
Collapse
Affiliation(s)
- Werner Dobrautz
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, 41296 Gothenburg, Sweden
| | - Igor O. Sokolov
- IBM
Quantum, IBM Research Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| | - Ke Liao
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Pablo López Ríos
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Martin Rahm
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, 41296 Gothenburg, Sweden
| | - Ali Alavi
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, U.K.
| | - Ivano Tavernelli
- IBM
Quantum, IBM Research Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| |
Collapse
|
18
|
Urban L, Laqua H, Thompson TH, Ochsenfeld C. Efficient Exploitation of Numerical Quadrature with Distance-Dependent Integral Screening in Explicitly Correlated F12 Theory: Linear Scaling Evaluation of the Most Expensive RI-MP2-F12 Term. J Chem Theory Comput 2024; 20:3706-3718. [PMID: 38626443 PMCID: PMC11099969 DOI: 10.1021/acs.jctc.4c00193] [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/15/2024] [Revised: 03/23/2024] [Accepted: 03/25/2024] [Indexed: 04/18/2024]
Abstract
We present a linear scaling atomic orbital based algorithm for the computation of the most expensive exchange-type RI-MP2-F12 term by employing numerical quadrature in combination with CABS-RI to avoid six-center-three-electron integrals. Furthermore, a robust distance-dependent integral screening scheme, based on integral partition bounds [Thompson, T. H.; Ochsenfeld, C. J. Chem. Phys. 2019, 150, 044101], is used to drastically reduce the number of the required three-center-one-electron integrals substantially. The accuracy of our numerical quadrature/CABS-RI approach and the corresponding integral screening is thoroughly assessed for interaction and isomerization energies across a variety of numerical integration grids. Our method outperforms the standard density fitting/CABS-RI approach with errors below 1 μEh even for small grid sizes and moderate screening thresholds. The choice of the grid size and screening threshold allows us to tailor our ansatz to a desired accuracy and computational efficiency. We showcase the approach's effectiveness for the chemically relevant system valinomycin, employing a triple-ζ F12 basis set combination (C54H90N6O18, 5757 AO basis functions, 10,266 CABS basis functions, 735,783 grid points). In this context, our ansatz achieves higher accuracy combined with a 135× speedup compared to the classical density fitting based variant, requiring notably less computation time than the corresponding RI-MP2 calculation. Additionally, we demonstrate near-linear scaling through calculations on linear alkanes. We achieved an 817-fold acceleration for C80H162 and an extrapolated 28,765-fold acceleration for C200H402, resulting in a substantially reduced computational time for the latter─from 229 days to just 11.5 min. Our ansatz may also be adapted to the remaining MP2-F12 terms, which will be the subject of future work.
Collapse
Affiliation(s)
- Lars Urban
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
- Max
Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| | - Henryk Laqua
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Travis H. Thompson
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
- Max
Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| |
Collapse
|
19
|
Drabik G, Radoń M. Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes. J Chem Theory Comput 2024; 20:3199-3217. [PMID: 38574194 PMCID: PMC11044276 DOI: 10.1021/acs.jctc.4c00092] [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/23/2024] [Revised: 03/19/2024] [Accepted: 03/20/2024] [Indexed: 04/06/2024]
Abstract
Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.
Collapse
Affiliation(s)
- Gabriela Drabik
- Jagiellonian
University, Doctoral School
of Exact and Natural Sciences, Łojasiewicza 11, 30-348 Kraków, Poland
- Jagiellonian
University, Faculty of Chemistry, Gronostajowa 2, 30-387, Kraków Poland
| | - Mariusz Radoń
- Jagiellonian
University, Faculty of Chemistry, Gronostajowa 2, 30-387, Kraków Poland
| |
Collapse
|
20
|
Bierman J, Li Y, Lu J. Qubit Count Reduction by Orthogonally Constrained Orbital Optimization for Variational Quantum Excited-State Solvers. J Chem Theory Comput 2024; 20:3131-3143. [PMID: 38598683 DOI: 10.1021/acs.jctc.3c01297] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2024]
Abstract
We propose a state-averaged orbital optimization scheme for improving the accuracy of excited states of the electronic structure Hamiltonian for use on near-term quantum computers. Instead of parameterizing the orbital rotation operator in the conventional fashion as an exponential of an antihermitian matrix, we parameterize the orbital rotation as a general partial unitary matrix. Whereas conventional orbital optimization methods minimize the state-averaged energy using successive Newton steps of the second-order Taylor expansion of the energy, the method presented here optimizes the state-averaged energy using an orthogonally constrained gradient projection method that does not require any expansion approximations. Through extensive benchmarking of the method on various small molecular systems, we find that the method is capable of producing more accurate results than fixed basis FCI while simultaneously using fewer qubits. In particular, we show that for H2, the method is capable of matching the accuracy of FCI in the cc-pVTZ basis (56 qubits) while only using 14 qubits.
Collapse
Affiliation(s)
- Joel Bierman
- Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - Yingzhou Li
- School of Mathematical Sciences, Fudan University, Shanghai 200433, China
- Shanghai Key Laboratory for Contemporary Applied Mathematics, Shanghai 200433, China
- Key Laboratory of Computational Physical Sciences (MOE), Shanghai 200433, China
| | - Jianfeng Lu
- Department of Physics, Duke University, Durham, North Carolina 27708, United States
- Department of Mathematics, Duke University, Durham, North Carolina 27708, United States
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
21
|
Fortenberry RC. Quantum Chemistry and Astrochemistry: A Match Made in the Heavens. J Phys Chem A 2024; 128:1555-1565. [PMID: 38381079 DOI: 10.1021/acs.jpca.3c07601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2024]
Abstract
Quantum chemistry can uniquely answer astrochemical questions that no other technique can provide. Computations can be parallelized, automated, and left to run continuously providing exceptional molecular throughput that cannot be done through experimentation. Additionally, the granularity of the individual computations that are required of potential energy surfaces, reaction mechanism pathways, or other quantum chemically derived observables produces a unique mosaic that make up the larger whole. These pieces can be dissected for their individual contributions or evaluated in an ad hoc fashion for each of their roles in generating the larger whole. No other scientific approach is capable of reporting such fine-grained insights. Quantum chemistry also works from a bottom-up approach in providing properties directly from the desired molecule instead of a top-down perspective as required of experiment where molecules have to be linked to observed phenomena. Furthermore, modern quantum chemistry is well within the range of "chemical accuracy" and is approaching "spectroscopic accuracy." As such, the seemingly difficult questions asked by astrochemistry that would not be asked initially for any other application require quantum chemical reference data. While the results of quantum chemical computations are needed to interpret astrochemical observation, modeling, or laboratory experimentation, such hard questions, regardless of the original need to answer them, produce unique solutions. While questions in astrochemistry often require novel developments in and implementations of quantum chemistry as outlined herein, the applications of these solutions will stretch beyond astrochemistry and may yet impact fields much closer to Earth.
Collapse
Affiliation(s)
- Ryan C Fortenberry
- Department of Chemistry & Biochemistry, University of Mississippi, Oxford, Mississippi 38677-1848, United States
| |
Collapse
|
22
|
Pathirage PDVS, Phillips JT, Vogiatzis KD. Exploration of the Two-Electron Excitation Space with Data-Driven Coupled Cluster. J Phys Chem A 2024. [PMID: 38422511 DOI: 10.1021/acs.jpca.3c06600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2024]
Abstract
Computational cost limits the applicability of post-Hartree-Fock methods such as coupled-cluster on larger molecular systems. The data-driven coupled-cluster (DDCC) method applies machine learning to predict the coupled-cluster two-electron amplitudes (t2) using data from second-order perturbation theory (MP2). One major limitation of the DDCC models is the size of training sets that increases exponentially with the system size. Effective sampling of the amplitude space can resolve this issue. Five different amplitude selection techniques that reduce the amount of data used for training were evaluated, an approach that also prevents model overfitting and increases the portability of data-driven coupled-cluster singles and doubles to more complex molecules or larger basis sets. In combination with a localized orbital formalism to predict the CCSD t2 amplitudes, we have achieved a 10-fold error reduction for energy calculations.
Collapse
Affiliation(s)
- P D Varuna S Pathirage
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996-1600, United States
| | - Justin T Phillips
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996-1600, United States
| | - Konstantinos D Vogiatzis
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996-1600, United States
| |
Collapse
|
23
|
Vaddypally S, Byrne AN, Goldsmith CF, Zdilla MJ, Kiselev VG. Metal-Free Reversible Double Cyclization of Cyanuric Diazide to an Asymmetric Bitetrazolate via Cleavage of the Six-Membered Aromatic Ring. Inorg Chem 2024; 63:2322-2326. [PMID: 38262914 DOI: 10.1021/acs.inorgchem.3c04338] [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/2024]
Abstract
Crystallization of the reaction mixture of 2-amino-4,6-diazido-1,3,5-triazine and excess tert-butylamine results in the isolation of tert-butylammonium N,N-[1'H-(1,5'-bitetrazol)-5-yl]cyanamidate, suggesting a complex decyclization/cyclization rearrangement involving breakage of the six-membered aromatic ring and the formation of two new five-membered azole rings mediated by deprotonation of the precursor by the amine. The addition of tert-butylamine to 2-amino-4,6-diazido-1,3,5-triazine gives spectroscopic indication of thermodynamically unfavorable reactivity in low-dielectric solvents, and high-level quantum chemical computations also suggest its formation to be unfavorable. A computed interconversion pathway describes the likely reaction mechanism and supports the general thermodynamic unfavorability of the reaction and the requirement for a high-dielectric environment to template formation of the ionic product and its trapping by crystallization.
Collapse
Affiliation(s)
- Shivaiah Vaddypally
- Department of Chemistry, Temple University, 1901 North 13th Street, Philadelphia, Pennsylvania 19122, United States
| | - Alex N Byrne
- Department of Chemistry, Temple University, 1901 North 13th Street, Philadelphia, Pennsylvania 19122, United States
| | - C Franklin Goldsmith
- School of Engineering, Brown University, 184 Hope Street, Providence, Rhode Island 02912, United States
| | - Michael J Zdilla
- Department of Chemistry, Temple University, 1901 North 13th Street, Philadelphia, Pennsylvania 19122, United States
| | - Vitaly G Kiselev
- School of Engineering, Brown University, 184 Hope Street, Providence, Rhode Island 02912, United States
- Novosibirsk State University, 1 Pirogova Street, 630090 Novosibirsk, Russia
- Institute of Chemical Kinetics and Combustion, Siberian Branch of the Russian Academy of Sciences, 3 Institutskaya Street, 630090 Novosibirsk, Russia
- Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 4 Kosygina Street, 119991 Moscow, Russia
| |
Collapse
|
24
|
Yang S, Zhang IY, Ren X. Developing correlation-consistent numeric atom-centered orbital basis sets for krypton: Applications in RPA-based correlated calculations. J Chem Phys 2024; 160:024112. [PMID: 38193553 DOI: 10.1063/5.0174952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Accepted: 12/04/2023] [Indexed: 01/10/2024] Open
Abstract
Localized atomic orbitals are the preferred basis set choice for large-scale explicit correlated calculations, and high-quality hierarchical correlation-consistent basis sets are a prerequisite for correlated methods to deliver numerically reliable results. At present, numeric atom-centered orbital (NAO) basis sets with valence correlation consistency (VCC), designated as NAO-VCC-nZ, are only available for light elements from hydrogen (H) to argon (Ar) [Zhang et al., New J. Phys. 15, 123033 (2013)]. In this work, we extend this series by developing NAO-VCC-nZ basis sets for krypton (Kr), a prototypical element in the fourth row of the periodic table. We demonstrate that NAO-VCC-nZ basis sets facilitate the convergence of electronic total-energy calculations using the Random Phase Approximation (RPA), which can be used together with a two-point extrapolation scheme to approach the complete basis set limit. Notably, the Basis Set Superposition Error (BSSE) associated with the newly generated NAO basis sets is minimal, making them suitable for applications where BSSE correction is either cumbersome or impractical to do. After confirming the reliability of NAO basis sets for Kr, we proceed to calculate the Helmholtz free energy for Kr crystal at the theoretical level of RPA plus renormalized single excitation correction. From this, we derive the pressure-volume (P-V) diagram, which shows excellent agreement with the latest experimental data. Our work demonstrates the capability of correlation-consistent NAO basis sets for heavy elements, paving the way toward numerically reliable correlated calculations for bulk materials.
Collapse
Affiliation(s)
- Sixian Yang
- Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, China
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Igor Ying Zhang
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai, Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Xinguo Ren
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
25
|
Volkmann H, Sathyanarayanan R, Saenz A, Jansen K, Kühn S. Chemically Accurate Potential Curves for H 2 Molecules Using Explicitly Correlated Qubit-ADAPT. J Chem Theory Comput 2024. [PMID: 38215397 DOI: 10.1021/acs.jctc.3c01281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2024]
Abstract
With the recent advances in the development of devices capable of performing quantum computations, a growing interest in finding near-term applications has emerged in many areas of science. In the era of nonfault tolerant quantum devices, algorithms that only require comparably short circuits accompanied by high repetition rates are considered to be a promising approach for assisting classical machines with finding a solution on computationally hard problems. The ADAPT approach previously introduced in Nat. Commun. 10, 3007 (2019) extends the class of variational quantum eigensolver algorithms with dynamically growing ansätze in order to find approximations to the ground and excited state energies of molecules. In this work, the ADAPT algorithm has been combined with a first-quantized formulation for the hydrogen molecule in the Born-Oppenheimer approximation, employing the explicitly correlated basis functions introduced in J. Chem. Phys. 43, 2429 (1965). By the virtue of their explicit electronic correlation properties, it is shown in classically performed simulations that chemical accuracy (<1.6 mHa) can be reached for ground and excited state potential curves using reasonably short circuits.
Collapse
Affiliation(s)
- Hakon Volkmann
- AG Moderne Optik, Institut für Physik, Humboldt-Universität zu Berlin, Newtonstraße 15, 12489 Berlin, Germany
| | - Raamamurthy Sathyanarayanan
- AG Moderne Optik, Institut für Physik, Humboldt-Universität zu Berlin, Newtonstraße 15, 12489 Berlin, Germany
| | - Alejandro Saenz
- AG Moderne Optik, Institut für Physik, Humboldt-Universität zu Berlin, Newtonstraße 15, 12489 Berlin, Germany
| | - Karl Jansen
- CQTA, DESY Zeuthen, Platanenallee 6, 15738 Zeuthen, Germany
- Computation-Based Science and Technology Research Center, The Cyprus Institute, 20 Kavafi Street, 2121 Nicosia, Cyprus
| | - Stefan Kühn
- CQTA, DESY Zeuthen, Platanenallee 6, 15738 Zeuthen, Germany
- Computation-Based Science and Technology Research Center, The Cyprus Institute, 20 Kavafi Street, 2121 Nicosia, Cyprus
| |
Collapse
|
26
|
Pitman SJ, Evans AK, Ireland RT, Lempriere F, McKemmish LK. Benchmarking Basis Sets for Density Functional Theory Thermochemistry Calculations: Why Unpolarized Basis Sets and the Polarized 6-311G Family Should Be Avoided. J Phys Chem A 2023; 127:10295-10306. [PMID: 37982604 DOI: 10.1021/acs.jpca.3c05573] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2023]
Abstract
Basis sets are a crucial but often largely overlooked choice in setting up quantum chemistry calculations. The choice of the basis set can be critical in determining the accuracy and calculation time of your quantum chemistry calculations. Clear recommendations based on thorough benchmarking are essential but not readily available currently. This study investigates the relative quality of basis sets for general properties by benchmarking basis set performance for a diverse set of 139 reactions (from the diet-150-GMTKN55 data set). In our analysis, we find the distributions of errors are often significantly non-Gaussian, meaning that the joint consideration of median errors, mean absolute errors, and outlier statistics is helpful to provide a holistic understanding of basis set performance. Our direct comparison of performance between most modern basis sets provides quantitative evidence for basis set recommendations that broadly align with the established understanding of basis set experts and is evident in the design of modern basis sets. For example, while zeta is a good measure of quality, it is not the only determining factor for an accurate calculation with unpolarized double- and triple-ζ basis sets (like 6-31G and 6-311G) having very poor performance. Appropriate use of polarization functions (e.g., 6-31G*) is essential to obtain the accuracy offered by double- or triple-ζ basis sets. In our study, the best performances for double- and triple-ζ basis sets are 6-31++G** and pcseg-2, respectively. However, the performances of singly polarized double-ζ and doubly polarized triple-ζ basis sets are quite similar with one key exception: the polarized 6-311G basis set family has poor parametrization, which means its performance is more like a double-ζ than a triple-ζ basis set. All versions of the 6-311G basis set family should be avoided entirely for valence chemistry calculations moving forward.
Collapse
Affiliation(s)
- Samuel J Pitman
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Alicia K Evans
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Robbie T Ireland
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Felix Lempriere
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Laura K McKemmish
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
27
|
Le Sech C, Sarsa A. Partial Separability of the Schrödinger Equation Combined with a Jastrow Factor. J Chem Theory Comput 2023; 19:8090-8096. [PMID: 37967284 DOI: 10.1021/acs.jctc.3c00635] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2023]
Abstract
Describing the Coulomb interactions between electrons in atomic or molecular systems is an important step to help us obtain accurate results for the different observables in the system. One convenient approach is to separate the dynamic electronic correlation, i.e., Coulomb electron-electron repulsion, from the motion of the electrons in the nucleus electric field. The wave function is written as the product of two terms, one accounting for the electron-electron interactions, which is symmetric under identical particle exchange, and the other is antisymmetric and represents the dynamics and exchange of electrons within the nuclear electric field. In this work, we present a novel computational scheme based on this idea that leads to an expression of the energy as the sum of two terms. To illustrate the method, we look into few-body Coulombic systems, H2, H3+, and Li(1s2,2s), and discuss the possible extension to larger systems. A simple correlation factor, based on the Jastrow exponential term, is employed to represent the dynamics of the electron pairs, leading to simple analytical forms and accurate results. We also present and illustrate a different approach with the Li atom based on the partial separability applied to a portion of the atom.
Collapse
Affiliation(s)
- Claude Le Sech
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay-ISMO (UMR 8214), Orsay Cedex 91405, France
| | - Antonio Sarsa
- Departamento de Física, Campus de Rabanales Edif. C2, Universidad de Córdoba, Córdoba E-14071, Spain
| |
Collapse
|
28
|
Wang Y, Guo Y, Neese F, Valeev EF, Li W, Li S. Cluster-in-Molecule Approach with Explicitly Correlated Methods for Large Molecules. J Chem Theory Comput 2023; 19:8076-8089. [PMID: 37920973 DOI: 10.1021/acs.jctc.3c00627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
In this article, we present a series of explicitly correlated local correlation methods developed under the cluster-in-molecule (CIM) framework, including explicitly correlated second-order Møller-Plesset perturbation (MP2), coupled-cluster singles and doubles (CCSD), domain-based local pair natural orbital CCSD (DLPNO-CCSD), and DLPNO-CCSD with perturbative triples (DLPNO-CCSD(T)). In these methods, F12 correction is decomposed into contributions from each occupied local molecular orbital and then evaluated independently in a given cluster, which consists of a subset of localized orbitals. These newly developed methods allow F12 calculations of large molecules (up to 145 atoms for quasi-one-dimensional systems) on a single node. We use these methods to investigate the relative stability between extended and folded alkane C30H62, the relative stability of four secondary structures of a polyglycine Ace(Gly)10NH2, and the binding energies of two host-guest complexes. The results demonstrate that the combination of CIM with F12 methods is a promising way to investigate large molecules with small basis set errors.
Collapse
Affiliation(s)
- Yuqi Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
29
|
Briccolani-Bandini L, Brémond E, Pagliai M, Cardini G, Ciofini I, Adamo C. Concerted versus stepwise proton transfer reactions in the [2, 2'-bipyridyl]-3-3'-diol molecule: A static and dynamic ab-initio investigation. J Comput Chem 2023; 44:2308-2318. [PMID: 37584183 DOI: 10.1002/jcc.27198] [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: 04/11/2023] [Revised: 07/09/2023] [Accepted: 07/11/2023] [Indexed: 08/17/2023]
Abstract
The double proton transfer (PT) reaction has been investigated in the [2,2'-bipyridyl]-3-3'-diol, a complex molecule where the proton movements is coupled to significant rearrangement of the electronic structure. Moreover, the reaction could be concerted, that is the two protons are exchanged simultaneously, or stepwise, where the two protons are transferred sequentially. To this end, a static exploration of the potential energy surface (PES) was carried together with the analysis of the free-energy surface (FES), both surfaces being evaluated at density functional theory level and different exchange-correlation functionals. While the concerted mechanism has been clearly discharged, the characteristics of the stepwise PT significantly depends on the chosen functionals, some suggesting a clear stepwise mechanism characterized by a stable reaction intermediates and two transitions states, whereas other approaches propend for a asynchronous PT, with a single TS. These features appear on both PES and FES, albeit some differences appears due to their different nature.
Collapse
Affiliation(s)
| | - Eric Brémond
- ITODYS, CNRS, Université Paris Cité, Paris, France
| | - Marco Pagliai
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Sesto Fiorentino, Italy
| | - Gianni Cardini
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Sesto Fiorentino, Italy
| | - Ilaria Ciofini
- CNRS, i-CLeHS, UMR 8060, Chimie ParisTech, PSL University, Paris, France
| | - Carlo Adamo
- CNRS, i-CLeHS, UMR 8060, Chimie ParisTech, PSL University, Paris, France
- Institut Universitaire de France, Paris, France
| |
Collapse
|
30
|
Derbali E, Ajili Y, Mehnen B, Żuchowski PS, Kędziera D, Al-Mogren MM, Jaidane NE, Hochlaf M. Towards the generation of potential energy surfaces of weakly bound medium-sized molecular systems: the case of benzonitrile-He complex. Phys Chem Chem Phys 2023; 25:30198-30210. [PMID: 37807943 DOI: 10.1039/d3cp02720a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/10/2023]
Abstract
Currently, the explicitly correlated coupled cluster method is used routinely to generate the multi-dimensional potential energy surfaces (mD-PESs) of van der Waals complexes of small molecular systems relevant for atmospheric, astrophysical and industrial applications. Although very accurate, this method is computationally prohibitive for medium and large molecules containing clusters. For instance, the recent detections of complex organic molecules (COMs) in the interstellar medium, such as benzonitrile, revealed the need to establish an accurate enough electronic structure approach to map the mD-PESs of these species interacting with the surrounding gases. As a benchmark, we have treated the case of the polar molecule benzonitrile interacting with helium, where we use post-Hartree-Fock and symmetry-adapted perturbation theory (SAPT) techniques. Accordingly, we show that MP2 and distinguishable-cluster approximation (DCSD) cannot be used for this purpose, whereas accurate enough PESs may be obtained using the corresponding explicitly correlated versions (MP2-F12 or DCSD-F12) with a reduction in computational costs. Alternatively, computations revealed that SAPT(DFT) is as performant as CCSD(T)-F12/aug-cc-pVTZ, making it the method of choice for mapping the mD-PESs of COMs containing clusters. Therefore, we have used this approach to generate the 3D-PES of the benzonitrile-He complex along the intermonomer Jacobi coordinates. As an application, we have incorporated the analytic form of this PES into quantum dynamical computations to determine the cross sections of the rotational (de-)excitation of benzonitrile colliding with helium at a collision energy of 10 cm-1.
Collapse
Affiliation(s)
- Eya Derbali
- Laboratoire de Spectroscopie Atomique, Moléculaire et Applications LSAMA, Université de Tunis Al Manar, Tunis, Tunisia.
| | - Yosra Ajili
- Laboratoire de Spectroscopie Atomique, Moléculaire et Applications LSAMA, Université de Tunis Al Manar, Tunis, Tunisia.
| | - Bilel Mehnen
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadz Street 5, 87-100 Toruń, Poland
| | - Piotr S Żuchowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadz Street 5, 87-100 Toruń, Poland
| | - Dariusz Kędziera
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, ul. Gagarina 7, PL 87-100 Toruń, Poland.
| | - Muneerah Mogren Al-Mogren
- Department of Chemistry, College of Sciences, King Saud University, PO Box 2455, Riyadh 11451, Saudi Arabia
| | - Nejm-Edine Jaidane
- Laboratoire de Spectroscopie Atomique, Moléculaire et Applications LSAMA, Université de Tunis Al Manar, Tunis, Tunisia.
| | - Majdi Hochlaf
- Université Gustave Eiffel, COSYS/IMSE, 5 Bd Descartes 77454, Champs sur Marne, France.
| |
Collapse
|
31
|
Crisci L, Di Grande S, Cavallotti C, Barone V. Toward an Accurate Black-Box Tool for the Kinetics of Gas-Phase Reactions Involving Barrier-less Elementary Steps. J Chem Theory Comput 2023; 19:7626-7639. [PMID: 37880932 PMCID: PMC10653117 DOI: 10.1021/acs.jctc.3c00857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Revised: 10/06/2023] [Accepted: 10/13/2023] [Indexed: 10/27/2023]
Abstract
An enhanced computational protocol has been devised for the accurate characterization of gas-phase barrier-less reactions in the framework of the reaction-path (RP) and variable reaction coordinate variational transition-state theory. In particular, the synergistic combination of density functional theory and Monte Carlo sampling to optimize reactive fluxes led to a reliable yet effective computational workflow. A black-box strategy has been developed for selecting the most suited density functional with reference to a high-level one-dimensional reference potential. At the same time, different descriptions of hindered rotations are automatically selected, depending on the corresponding harmonic frequencies along the RP. The performance of the new tool is investigated by means of two prototypical reactions involving different degrees of static and dynamic correlation, namely, H2S + Cl and CH3 + CH3. The remarkable agreement of the computed kinetic parameters with the available experimental data confirms the accuracy and robustness of the proposed approach. Together with their intrinsic interest, these results also pave the way toward systematic investigations of gas-phase reactions involving barrier-less elementary steps by a reliable, user-friendly tool, which can be confidently used also by nonspecialists.
Collapse
Affiliation(s)
- Luigi Crisci
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di M.S. Angelo, via Cintia 21, I-80126 Napoli, Italy
| | - Silvia Di Grande
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Carlo Cavallotti
- Department
of Chemistry, Materials and Chemical Engineering “G. Natta”, Politecnico di Milano, I-20131 Milano, Italy
| | - Vincenzo Barone
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
32
|
Ten-No SL. Nonunitary projective transcorrelation theory inspired by the F12 ansatz. J Chem Phys 2023; 159:171103. [PMID: 37921247 DOI: 10.1063/5.0175337] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023] Open
Abstract
An alternative nonunitary transcorrelation, inspired by the F12 ansatz, is investigated. In contrast to the Jastrow transcorrelation of Boys-Handy, the effective Hamiltonian of this projective transcorrelation features: 1. a series terminating formally at four-body interactions. 2. no spin-contamination within the non-relativistic framework. 3. simultaneous satisfaction of the singlet and triplet first-order cusp conditions. 4. arbitrary choices of pairs for correlation including frozen-core approximations. We discuss the connection between the projective transcorrelation and F12 theory with applications to small molecules, to show that the cusp conditions play an important role to reduce the uncertainty arising from the nonunitary transformation.
Collapse
Affiliation(s)
- Seiichiro L Ten-No
- Graduate School of System Informatics, Kobe University, Nada-ku, Kobe 657-8501, Japan
| |
Collapse
|
33
|
Valeev EF, Harrison RJ, Holmes AA, Peterson CC, Penchoff DA. Direct Determination of Optimal Real-Space Orbitals for Correlated Electronic Structure of Molecules. J Chem Theory Comput 2023; 19:7230-7241. [PMID: 37791808 DOI: 10.1021/acs.jctc.3c00732] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
We demonstrate how to determine numerically nearly exact orthonormal orbitals that are optimal for the evaluation of the energy of arbitrary (correlated) states of atoms and molecules by minimization of the energy Lagrangian. Orbitals are expressed in real space using a multiresolution spectral element basis that is refined adaptively to achieve the user-specified target precision while avoiding the ill-conditioning issues that plague AO basis set expansions traditionally used for correlated models of molecular electronic structure. For light atoms, the orbital solver, in conjunction with a variational electronic structure model [selected Configuration Interaction (CI)] provides energies of comparable precision to a state-of-the-art atomic CI solver. The computed electronic energies of atoms and molecules are significantly more accurate than the counterparts obtained with the orbital sets of the same rank expanded in Gaussian AO bases, and can be determined even when linear dependence issues preclude the use of the AO bases. It is feasible to optimize more than 100 fully correlated numerical orbitals on a single computer node, and significant room exists for additional improvement. These findings suggest that real-space orbital representations might be the preferred alternative to AO representations for high-end models of correlated electronic states of molecules and materials.
Collapse
Affiliation(s)
- Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Robert J Harrison
- Department of Applied Mathematics & Statistics, Stony Brook University, Stony Brook, New York 11794, United States
| | - Adam A Holmes
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Charles C Peterson
- Office of Advanced Research Computing, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Deborah A Penchoff
- UT Innovative Computing Laboratory, University of Tennessee, Knoxville, Tennessee 37996, United States
| |
Collapse
|
34
|
Di Felice R, Mayes ML, Richard RM, Williams-Young DB, Chan GKL, de Jong WA, Govind N, Head-Gordon M, Hermes MR, Kowalski K, Li X, Lischka H, Mueller KT, Mutlu E, Niklasson AMN, Pederson MR, Peng B, Shepard R, Valeev EF, van Schilfgaarde M, Vlaisavljevich B, Windus TL, Xantheas SS, Zhang X, Zimmerman PM. A Perspective on Sustainable Computational Chemistry Software Development and Integration. J Chem Theory Comput 2023; 19:7056-7076. [PMID: 37769271 PMCID: PMC10601486 DOI: 10.1021/acs.jctc.3c00419] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Indexed: 09/30/2023]
Abstract
The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.
Collapse
Affiliation(s)
- Rosa Di Felice
- Departments
of Physics and Astronomy and Quantitative and Computational Biology, University of Southern California, Los Angeles, California 90089, United States
- CNR-NANO
Modena, Modena 41125, Italy
| | - Maricris L. Mayes
- Department
of Chemistry and Biochemistry, University
of Massachusetts Dartmouth, North Dartmouth, Massachusetts 02747, United States
| | | | | | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Wibe A. de Jong
- Lawrence
Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Niranjan Govind
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Martin Head-Gordon
- Pitzer Center
for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Karol Kowalski
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409, United States
| | - Karl T. Mueller
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced
Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Anders M. N. Niklasson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Mark R. Pederson
- Department
of Physics, The University of Texas at El
Paso, El Paso, Texas 79968, United States
| | - Bo Peng
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, United States
| | - Edward F. Valeev
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | | | - Bess Vlaisavljevich
- Department
of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Theresa L. Windus
- Department
of Chemistry, Iowa State University and
Ames Laboratory, Ames, Iowa 50011, United States
| | - Sotiris S. Xantheas
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Advanced
Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Xing Zhang
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Paul M. Zimmerman
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
35
|
Peng H, Yang S, Jiang H, Weng H, Ren X. Basis-Set-Error-Free Random-Phase Approximation Correlation Energies for Atoms Based on the Sternheimer Equation. J Chem Theory Comput 2023; 19:7199-7214. [PMID: 37811855 DOI: 10.1021/acs.jctc.3c00668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/10/2023]
Abstract
The finite basis set errors for all-electron random-phase approximation (RPA) correlation energy calculations are analyzed for isolated atomic systems. We show that, within the resolution-of-identity (RI) RPA framework, the major source of the basis set errors is the incompleteness of the single-particle atomic orbitals used to expand the Kohn-Sham eigenstates, instead of the auxiliary basis set (ABS) to represent the density response function χ0 and the bare Coulomb operator v. By solving the Sternheimer equation for the first-order wave function on a dense radial grid, we are able to eliminate the major error─the incompleteness error of the single-particle atomic basis set─for atomic RPA calculations. The error stemming from a finite ABS can be readily rendered vanishingly small by increasing the size of the ABS, or by iteratively determining the eigenmodes of the χ0v operator. The variational property of the RI-RPA correlation energy can be further exploited to optimize the ABS in order to achieve fast convergence of the RI-RPA correlation energy. These numerical techniques enable us to obtain basis-set-error-free RPA correlation energies for atoms, and in this work, such energies for atoms from H to Kr are presented. The implications of the numerical techniques developed in the present work for addressing the basis set issue for molecules and solids are discussed.
Collapse
Affiliation(s)
- Hao Peng
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Sixian Yang
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, China
| | - Hong Jiang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Hongming Weng
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Songshan Lake Materials Laboratory, Dongguan 523808, Guangdong, China
| | - Xinguo Ren
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Songshan Lake Materials Laboratory, Dongguan 523808, Guangdong, China
| |
Collapse
|
36
|
Di Grande S, Kállay M, Barone V. Accurate thermochemistry at affordable cost by means of an improved version of the JunChS-F12 model chemistry. J Comput Chem 2023; 44:2149-2157. [PMID: 37432050 DOI: 10.1002/jcc.27187] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2023] [Accepted: 06/23/2023] [Indexed: 07/12/2023]
Abstract
The junChS-F12 composite method has been improved by means of the latest implementation of the CCSD(F12*)(T+) ansatz and validated for the thermochemistry of molecules containing atoms of the first three rows of the periodic table. A thorough benchmark showed that this model, in conjunction with cost-effective revDSD-PBEP86-D3(BJ) reference geometries, offers an optimal compromise between accuracy and computational cost. If improved geometries are sought, the most effective option is to add MP2-F12 core-valence correlation corrections to CCSD(T)-F12b/jun-cc-pVTZ geometries without the need of performing any extrapolation to the complete basis set limit. In the same vein, CCSD(T)-F12b/jun-cc-pVTZ harmonic frequencies are remarkably accurate without any additional contribution. Pilot applications to noncovalent intermolecular interactions, conformational landscapes, and tautomeric equilibria confirm the effectiveness and reliability of the model.
Collapse
Affiliation(s)
- S Di Grande
- Classe di Scienze, Scuola Normale Superiore di Pisa, Pisa, Italy
- Scuola Superiore Meridionale, Napoli, Italy
| | - M Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Budapest, Hungary
| | - V Barone
- Classe di Scienze, Scuola Normale Superiore di Pisa, Pisa, Italy
| |
Collapse
|
37
|
Semidalas E, Martin JML. Correlation Consistent Basis Sets for Explicitly Correlated Theory: The Transition Metals. J Chem Theory Comput 2023; 19:5806-5820. [PMID: 37540641 PMCID: PMC10500978 DOI: 10.1021/acs.jctc.3c00506] [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/16/2023] [Indexed: 08/06/2023]
Abstract
We present correlation consistent basis sets for explicitly correlated (F12) calculations, denoted VnZ(-PP)-F12-wis (n = D,T), for the d-block elements. The cc-pVDZ-F12-wis basis set is contracted to [8s7p5d2f] for the 3d-block, while its ECP counterpart for the 4d and 5d-blocks, cc-pVDZ-PP-F12-wis, is contracted to [6s6p5d2f]. The corresponding contracted sizes for cc-pVTZ(-PP)-F12-wis are [9s8p6d3f2g] for the 3d-block elements and [7s7p6d3f2g] for the 4d and 5d-block elements. Our VnZ(-PP)-F12-wis basis sets are evaluated on challenging test sets for metal-organic barrier heights (MOBH35) and group-11 metal clusters (CUAGAU-2). In F12 calculations, they are found to be about as close to the complete basis set limit as the combination of standard cc-pVnZ-F12 on main-group elements with the standard aug-cc-pV(n+1)Z(-PP) basis sets on the transition metal(s). While our basis sets are somewhat more compact than aug-cc-pV(n+1)Z(-PP), the CPU time benefit is negligible for catalytic complexes that contain only one or two transition metals among dozens of main-group elements; however, it is somewhat more significant for metal clusters.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Molecular Chemistry
and Materials Science, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Molecular Chemistry
and Materials Science, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
38
|
Lee N, Thom AJW. Studies on the Transcorrelated Method. J Chem Theory Comput 2023; 19:5743-5759. [PMID: 37640393 PMCID: PMC10500994 DOI: 10.1021/acs.jctc.3c00046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2023] [Indexed: 08/31/2023]
Abstract
We investigate the possibility of using a transcorrelated (TC) Hamiltonian to describe electron correlation. A method to obtain TC wavefunctions was developed based on the mathematical framework of the bi-variational principle. This involves the construction of an effective TC Hamiltonian matrix, which can be solved in a self-consistent manner. This was optimized using a method we call second-order-moment minimization and demonstrate that it is possible to obtain highly accurate energies for some closed-shell atoms and helium-like ions. The effects of certain correlator terms on the description of electron-electron and electron-nuclear cusps were also examined graphically, and some TC wavefunctions were compared against near-exact Hylleraas wavefunctions.
Collapse
Affiliation(s)
- Nicholas Lee
- Department
of Chemistry, Physical and Theoretical Chemistry
Laboratory, South Parks
Road, Oxford OX1 3QZ, U.K.
| | - Alex J. W. Thom
- Yusuf
Hamied Department of Chemistry, Lensfield Road, Cambridge CB2 1EW, U.K.
| |
Collapse
|
39
|
Barone V, Di Grande S, Lazzari F, Mendolicchio M. Accurate Structures and Spectroscopic Parameters of Guanine Tautomers in the Gas Phase by the Pisa Conventional and Explicitly Correlated Composite Schemes (PCS and PCS-F12). J Phys Chem A 2023; 127:6771-6778. [PMID: 37535450 PMCID: PMC10440789 DOI: 10.1021/acs.jpca.3c03999] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Revised: 07/21/2023] [Indexed: 08/05/2023]
Abstract
A general strategy for the accurate computation of structural and spectroscopic properties of biomolecule building blocks in the gas phase is proposed and validated for tautomeric equilibria. The main features of the new model are the inclusion of core-valence correlation in geometry optimizations by a double hybrid functional and the systematic use of wave-function composite methods in conjunction with cc-pVnZ-F12 basis sets with separate extrapolation of MP2 and post-MP2 contributions. The resulting Pisa composite scheme employing conventional (PCS) or explicitly correlated (PCS-F12) approaches is applied to the challenging problem of guanine tautomers in the gas phase. The results are in remarkable agreement with the experimental structures, relative stabilities, and spectroscopic signatures of different tautomers. The accuracy of the results obtained at reasonable cost by means of black-box parameter-free approaches paves the way toward systematic investigations of other molecular bricks of life also by non-specialists.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | - Silvia Di Grande
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, Napoli 80138, Italy
| | - Federico Lazzari
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | | |
Collapse
|
40
|
Haupt JP, Hosseini SM, López Ríos P, Dobrautz W, Cohen A, Alavi A. Optimizing Jastrow factors for the transcorrelated method. J Chem Phys 2023; 158:2895246. [PMID: 37290083 DOI: 10.1063/5.0147877] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Accepted: 05/25/2023] [Indexed: 06/10/2023] Open
Abstract
We investigate the optimization of flexible tailored real-space Jastrow factors for use in the transcorrelated (TC) method in combination with highly accurate quantum chemistry methods, such as initiator full configuration interaction quantum Monte Carlo (FCIQMC). Jastrow factors obtained by minimizing the variance of the TC reference energy are found to yield better, more consistent results than those obtained by minimizing the variational energy. We compute all-electron atomization energies for the challenging first-row molecules C2, CN, N2, and O2 and find that the TC method yields chemically accurate results using only the cc-pVTZ basis set, roughly matching the accuracy of non-TC calculations with the much larger cc-pV5Z basis set. We also investigate an approximation in which pure three-body excitations are neglected from the TC-FCIQMC dynamics, saving storage and computational costs, and show that it affects relative energies negligibly. Our results demonstrate that the combination of tailored real-space Jastrow factors with the multi-configurational TC-FCIQMC method provides a route to obtaining chemical accuracy using modest basis sets, obviating the need for basis-set extrapolation and composite techniques.
Collapse
Affiliation(s)
- J Philip Haupt
- Max-Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | | | - Pablo López Ríos
- Max-Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Werner Dobrautz
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology, 41296 Gothenburg, Sweden
| | - Aron Cohen
- DeepMind, 6 Pancras Square, London N1C 4AG, United Kingdom
| | - Ali Alavi
- Max-Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
41
|
Rzhevskiy SA, Minaeva LI, Topchiy MA, Melnikov IN, Kiselev VG, Pivkina AN, Fomenkov IV, Asachenko AF. Synthesis, Characterization, and Properties of High-Energy Fillers Derived from Nitroisobutylglycerol. Int J Mol Sci 2023; 24:8541. [PMID: 37239887 PMCID: PMC10218491 DOI: 10.3390/ijms24108541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Revised: 05/03/2023] [Accepted: 05/08/2023] [Indexed: 05/28/2023] Open
Abstract
Herein we report a comprehensive laboratory synthesis of a series of energetic azidonitrate derivatives (ANDP, SMX, AMDNNM, NIBTN, NPN, 2-nitro-1,3-dinitro-oxypropane) starting from the readily available nitroisobutylglycerol. This simple protocol allows obtaining the high-energy additives from the available precursor in yields higher than those reported using safe and simple operations not presented in previous works. A detailed characterization of the physical, chemical, and energetic properties including impact sensitivity and thermal behavior of these species was performed for the systematic evaluation and comparison of the corresponding class of energetic compounds.
Collapse
Affiliation(s)
- Sergey A. Rzhevskiy
- A.V. Topchiev Institute of Petrochemical Synthesis, Russian Academy of Sciences, Leninsky Prospect 29, 119991 Moscow, Russia; (S.A.R.); (L.I.M.); (M.A.T.)
| | - Lidiya I. Minaeva
- A.V. Topchiev Institute of Petrochemical Synthesis, Russian Academy of Sciences, Leninsky Prospect 29, 119991 Moscow, Russia; (S.A.R.); (L.I.M.); (M.A.T.)
| | - Maxim A. Topchiy
- A.V. Topchiev Institute of Petrochemical Synthesis, Russian Academy of Sciences, Leninsky Prospect 29, 119991 Moscow, Russia; (S.A.R.); (L.I.M.); (M.A.T.)
| | - Igor N. Melnikov
- Semenov Federal Research Center for Chemical Physics RAS, 4 Kosygina Str., 119991 Moscow, Russia; (I.N.M.); (A.N.P.)
| | - Vitaly G. Kiselev
- Institute of Chemical Kinetics and Combustion SB RAS, 3 Institutskaya Str., 630090 Novosibirsk, Russia;
| | - Alla N. Pivkina
- Semenov Federal Research Center for Chemical Physics RAS, 4 Kosygina Str., 119991 Moscow, Russia; (I.N.M.); (A.N.P.)
| | - Igor V. Fomenkov
- Zelinsky Institute of Organic Chemistry RAS, 47 Leninsky Ave., 119991 Moscow, Russia;
| | - Andrey F. Asachenko
- A.V. Topchiev Institute of Petrochemical Synthesis, Russian Academy of Sciences, Leninsky Prospect 29, 119991 Moscow, Russia; (S.A.R.); (L.I.M.); (M.A.T.)
| |
Collapse
|
42
|
Nash HW, Shaw RA, Hill JG. Correlation consistent auxiliary basis sets in density fitting Hartree-Fock: The atoms sodium through argon revisited. J Comput Chem 2023; 44:1119-1128. [PMID: 36636897 PMCID: PMC10947126 DOI: 10.1002/jcc.27069] [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: 08/30/2022] [Revised: 11/29/2022] [Accepted: 12/14/2022] [Indexed: 01/14/2023]
Abstract
We present a series of auxiliary basis sets, for the elements Na to Ar, for use in density-fitted Hartree-Fock calculations with the correlation consistent cc-pV(n + d)Z orbital basis sets. Benchmarking on total molecular energies, reaction energies and the spectroscopic constants of the SO molecule demonstrate that the new sets address the deficiencies of using existing auxiliary sets in combination with these orbital basis sets. We also report auxiliary basis sets for Na and Mg matched to cc-pVnZ, along with recommendations for pairing auxiliary basis sets to the cc-pVnZ-F12 basis sets for Hartree-Fock calculations.
Collapse
Affiliation(s)
- Harry W. Nash
- Department of ChemistryUniversity of SheffieldSheffieldUK
| | - Robert A. Shaw
- Department of ChemistryUniversity of SheffieldSheffieldUK
| | | |
Collapse
|
43
|
Ma H, Liu J, Shang H, Fan Y, Li Z, Yang J. Multiscale quantum algorithms for quantum chemistry. Chem Sci 2023; 14:3190-3205. [PMID: 36970085 PMCID: PMC10034224 DOI: 10.1039/d2sc06875c] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Accepted: 02/15/2023] [Indexed: 02/18/2023] Open
Abstract
Exploring the potential applications of quantum computers in material design and drug discovery is attracting more and more attention after quantum advantage has been demonstrated using Gaussian boson sampling. However, quantum resource requirements in material and (bio)molecular simulations are far beyond the capacity of near-term quantum devices. In this work, multiscale quantum computing is proposed for quantum simulations of complex systems by integrating multiple computational methods at different scales of resolution. In this framework, most computational methods can be implemented in an efficient way on classical computers, leaving the critical portion of the computation to quantum computers. The simulation scale of quantum computing strongly depends on available quantum resources. As a near-term scheme, we integrate adaptive variational quantum eigensolver algorithms, second-order Møller-Plesset perturbation theory and Hartree-Fock theory within the framework of the many-body expansion fragmentation approach. This new algorithm is applied to model systems consisting of hundreds of orbitals with decent accuracy on the classical simulator. This work should encourage further studies on quantum computing for solving practical material and biochemistry problems.
Collapse
Affiliation(s)
- Huan Ma
- Hefei National Laboratory, University of Science and Technology of China Hefei 230088 China
| | - Jie Liu
- Hefei National Laboratory, University of Science and Technology of China Hefei 230088 China
| | - Honghui Shang
- State Key Laboratory of Computer Architecture, Institute of Computing Technology, Chinese Academy of Sciences Beijing 100190 China
| | - Yi Fan
- Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China Hefei 230026 China
| | - Zhenyu Li
- Hefei National Laboratory, University of Science and Technology of China Hefei 230088 China
- Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China Hefei 230026 China
| | - Jinlong Yang
- Hefei National Laboratory, University of Science and Technology of China Hefei 230088 China
- Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China Hefei 230026 China
| |
Collapse
|
44
|
Melnikov IN, Kiselev VG, Dalinger IL, Starosotnikov AM, Muravyev NV, Pivkina AN. Thermochemistry, Tautomerism, and Thermal Stability of 5,7-Dinitrobenzotriazoles. Int J Mol Sci 2023; 24:ijms24065330. [PMID: 36982405 PMCID: PMC10049112 DOI: 10.3390/ijms24065330] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2022] [Revised: 02/09/2023] [Accepted: 03/08/2023] [Indexed: 03/14/2023] Open
Abstract
Nitro derivatives of benzotriazoles are safe energetic materials with remarkable thermal stability. In the present study, we report on the kinetics and mechanism of thermal decomposition for 5,7-dinitrobenzotriazole (DBT) and 4-amino-5,7-dinitrobenzotriazole (ADBT). The pressure differential scanning calorimetry was employed to study the decomposition kinetics of DBT experimentally because the measurements under atmospheric pressure are disturbed by competing evaporation. The thermolysis of DBT in the melt is described by a kinetic scheme with two global reactions. The first stage is a strong autocatalytic process that includes the first-order reaction (Ea1I = 173.9 ± 0.9 kJ mol−1, log(A1I/s−1) = 12.82 ± 0.09) and the catalytic reaction of the second order with Ea2I = 136.5 ± 0.8 kJ mol−1, log(A2I/s−1) = 11.04 ± 0.07. The experimental study was complemented by predictive quantum chemical calculations (DLPNO-CCSD(T)). The calculations reveal that the 1H tautomer is the most energetically preferable form for both DBT and ADBT. Theory suggests the same decomposition mechanisms for DBT and ADBT, with the most favorable channels being nitro-nitrite isomerization and C–NO2 bond cleavage. The former channel has lower activation barriers (267 and 276 kJ mol−1 for DBT and ADBT, respectively) and dominates at lower temperatures. At the same time, due to the higher preexponential factor, the radical bond cleavage, with reaction enthalpies of 298 and 320 kJ mol−1, dominates in the experimental temperature range for both DBT and ADBT. In line with the theoretical predictions of C–NO2 bond energies, ADBT is more thermally stable than DBT. We also determined a reliable and mutually consistent set of thermochemical values for DBT and ADBT by combining the theoretically calculated (W1-F12 multilevel procedure) gas-phase enthalpies of formation and experimentally measured sublimation enthalpies.
Collapse
Affiliation(s)
- Igor N. Melnikov
- Semenov Federal Research Center for Chemical Physics RAS, 4 Kosygina Str., 119991 Moscow, Russia
| | - Vitaly G. Kiselev
- Semenov Federal Research Center for Chemical Physics RAS, 4 Kosygina Str., 119991 Moscow, Russia
- Physics Department, Novosibirsk State University, 1 Pirogova Str., 630090 Novosibirsk, Russia
- Institute of Chemical Kinetics and Combustion SB RAS, 3 Institutskaya Str., 630090 Novosibirsk, Russia
| | - Igor L. Dalinger
- Zelinsky Institute of Organic Chemistry RAS, 47 Leninsky Ave., 119991 Moscow, Russia
| | | | - Nikita V. Muravyev
- Semenov Federal Research Center for Chemical Physics RAS, 4 Kosygina Str., 119991 Moscow, Russia
- Correspondence: ; Tel.: +7-499-137-8203
| | - Alla N. Pivkina
- Semenov Federal Research Center for Chemical Physics RAS, 4 Kosygina Str., 119991 Moscow, Russia
| |
Collapse
|
45
|
Steenbakkers K, Marimuthu AN, Redlich B, Groenenboom GC, Brünken S. A vibrational action spectroscopic study of the Renner-Teller- and spin-orbit-affected cyanoacetylene radical cation HC 3N . J Chem Phys 2023; 158:084305. [PMID: 36859081 DOI: 10.1063/5.0135000] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
The linear radical cation of cyanoacetylene, HC3N+ (2Π), is not only of astrophysical interest, being the, so far undetected, cationic counterpart of the abundant cyanoaceteylene, but also of fundamental spectroscopic interest due to its strong spin-orbit and Renner-Teller interactions. Here, we present the first broadband vibrational action spectroscopic investigation of this ion through the infrared pre-dissociation (IRPD) method using a Ne tag. Experiments have been performed using the FELion cryogenic ion-trap instrument in combination with the FELIX free-electron lasers and a Laservision optical parametric oscillator/optical parametric amplifier system. The vibronic splitting patterns of the three interacting bending modes (ν5, ν6, ν7), ranging from 180 to 1600 cm-1, could be fully resolved revealing several bands that were previously unobserved. The associated Renner-Teller and intermode coupling constants have been determined by fitting an effective Hamiltonian to the experimental data, and the obtained spectroscopic constants are in reasonable agreement with previous photoelectron spectroscopy (PES) studies and ab initio calculations on the HC3N+ ion. The influence of the attached Ne atom on the infrared spectrum has been investigated by ab initio calculations at the RCCSD(T)-F12a level of theory, which strongly indicates that the discrepancies between the IRPD and PES data are a result of the effects of the Ne attachment.
Collapse
Affiliation(s)
- Kim Steenbakkers
- FELIX Laboratory, Institute for Molecules and Materials, Radboud University, Nijmegen, The Netherlands
| | - Aravindh N Marimuthu
- FELIX Laboratory, Institute for Molecules and Materials, Radboud University, Nijmegen, The Netherlands
| | - Britta Redlich
- FELIX Laboratory, Institute for Molecules and Materials, Radboud University, Nijmegen, The Netherlands
| | - Gerrit C Groenenboom
- Radboud University, Institute for Molecules and Materials, Nijmegen, The Netherlands
| | - Sandra Brünken
- FELIX Laboratory, Institute for Molecules and Materials, Radboud University, Nijmegen, The Netherlands
| |
Collapse
|
46
|
Minenkov Y, Cavallo L, Peterson KA. Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds. J Comput Chem 2023; 44:687-696. [PMID: 36399072 DOI: 10.1002/jcc.27033] [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: 07/13/2022] [Revised: 09/23/2022] [Accepted: 09/28/2022] [Indexed: 11/19/2022]
Abstract
The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted-core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller-Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol-1 relative to the experiment. The two-point QZ/5Z and TZ/QZ CBS 1 / l max 3 extrapolation schemes are inferior to their inverse power counterpart ( 1 / l max + 1 / 2 4 ) by 1.3 and 4.3 kJ mol-1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol-1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double-ζ/triple-ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.
Collapse
Affiliation(s)
- Yury Minenkov
- N.N. Semenov Federal Research Center for Chemical Physics RAS, Moscow, Russian Federation.,Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow, Russian Federation
| | - Luigi Cavallo
- KAUST Catalysis Center (KCC), King Abdullah University of Science and Technology, Thuwal, Saudi Arabia
| | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington, USA
| |
Collapse
|
47
|
Samsonova I, Tucker GB, Alaal N, Brorsen KR. Hydrogen-Atom Electronic Basis Sets for Multicomponent Quantum Chemistry. ACS OMEGA 2023; 8:5033-5041. [PMID: 36777583 PMCID: PMC9910068 DOI: 10.1021/acsomega.2c07782] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Accepted: 01/13/2023] [Indexed: 06/18/2023]
Abstract
Multicomponent methods are a conceptually simple way to include nuclear quantum effects into quantum chemistry calculations. In multicomponent methods, the electronic molecular orbitals are described using the linear combination of atomic orbitals approximation. This requires the selection of a one-particle electronic basis set which, in practice, is commonly a correlation-consistent basis set. In multicomponent method studies, it has been demonstrated that large electronic basis sets are required for quantum hydrogen nuclei to accurately describe electron-nuclear correlation. However, as we show in this study, much of the need for large electronic basis sets is due to the correlation-consistent electronic basis sets not being optimized to describe nuclear properties and electron-nuclear correlation. Herein, we introduce a series of correlation-consistent electronic basis sets for hydrogen atoms called cc-pVnZ-mc with additional basis functions optimized to reproduce multicomponent density functional theory protonic densities. These new electronic basis sets are shown to yield better protonic densities with fewer electronic basis functions than the standard correlation-consistent basis sets and reproduce other protonic properties such as proton affinities and protonic excitation energies, even though they were not optimized for these purposes. The cc-pVnZ-mc basis sets should enable multicomponent many-body calculations on larger systems due to the improved computational efficiency they provide for a given level of accuracy.
Collapse
|
48
|
Nasiri S, Shomenov T, Bubin S, Adamowicz L. Dissociation energy and the lowest vibrational transition in LiH without assuming the non-Born–Oppenheimer approximation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2147105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Affiliation(s)
- Saeed Nasiri
- Department of Physics, Nazarbayev University, Nur-Sultan, Kazakhstan
| | | | - Sergiy Bubin
- Department of Physics, Nazarbayev University, Nur-Sultan, Kazakhstan
| | - Ludwik Adamowicz
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, AZ, USA
- Department of Physics, University of Arizona, Tucson, AZ, USA
- Centre for Advanced Study (CAS), the Norwegian Academy of Science and Letters, Oslo, Norway
| |
Collapse
|
49
|
Mehta N, Martin JML. Reduced-Scaling Double Hybrid Density Functional Theory with Rapid Basis Set Convergence through Localized Pair Natural Orbital F12. J Phys Chem Lett 2022; 13:9332-9338. [PMID: 36178852 PMCID: PMC9575149 DOI: 10.1021/acs.jpclett.2c02620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Accepted: 09/26/2022] [Indexed: 06/16/2023]
Abstract
Following earlier work [Mehta, N.; Martin, J. M. L. J. Chem. Theory Comput.2022, 10.1021/acs.jctc.2c00426] that showed how the slow basis set convergence of the double hybrid density functional theory can be obviated by the use of F12 explicit correlation in the GLPT2 step (second order Görling-Levy perturbation theory), we demonstrate here for the very large and chemically diverse GMTKN55 benchmark suite that the CPU time scaling of this step can be reduced (asymptotically linearized) using the localized pair natural orbital (PNO-L) approximation at negligible cost in accuracy.
Collapse
Affiliation(s)
- Nisha Mehta
- Department of Molecular Chemistry and
Materials Science, Weizmann Institute of
Science, Reḥovot7610001, Israel
| | - Jan M. L. Martin
- Department of Molecular Chemistry and
Materials Science, Weizmann Institute of
Science, Reḥovot7610001, Israel
| |
Collapse
|
50
|
Mehta N, Martin JML. Explicitly Correlated Double-Hybrid DFT: A Comprehensive Analysis of the Basis Set Convergence on the GMTKN55 Database. J Chem Theory Comput 2022; 18:5978-5991. [PMID: 36099641 PMCID: PMC9558368 DOI: 10.1021/acs.jctc.2c00426] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2022] [Indexed: 11/28/2022]
Abstract
Double-hybrid density functional theory (DHDFT) offers a pathway to accuracy approaching composite wavefunction approaches such as G4 theory. However, the Görling-Levy second-order perturbation theory (GLPT2) term causes them to partially inherit the slow ∝L-3 (with L the maximum angular momentum) basis set convergence of correlated wavefunction methods. This could potentially be remedied by introducing F12 explicit correlation: we investigate the basis set convergence of both DHDFT and DHDFT-F12 (where GLPT2 is replaced by GLPT2-F12) for the large and chemically diverse general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) benchmark suite. The B2GP-PLYP-D3(BJ) and revDSD-PBEP86-D4 DHDFs are investigated as test cases, together with orbital basis sets as large as aug-cc-pV5Z and F12 basis sets as large as cc-pVQZ-F12. We show that F12 greatly accelerates basis set convergence of DHDFs, to the point that even the modest cc-pVDZ-F12 basis set is closer to the basis set limit than cc-pV(Q+d)Z or def2-QZVPPD in orbital-based approaches, and in fact comparable in quality to cc-pV(5+d)Z. Somewhat surprisingly, aug-cc-pVDZ-F12 is not required even for the anionic subsets. In conclusion, DHDF-F12/VDZ-F12 eliminates concerns about basis set convergence in both the development and applications of double-hybrid functionals. Mass storage and I/O bottlenecks for larger systems can be circumvented by localized pair natural orbital approximations, which also exhibit much gentler system size scaling.
Collapse
Affiliation(s)
- Nisha Mehta
- Department of Molecular Chemistry and
Materials Science, Weizmann Institute of
Science, 7610001 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Molecular Chemistry and
Materials Science, Weizmann Institute of
Science, 7610001 Reḥovot, Israel
| |
Collapse
|