1
|
Larsson HR, Viel A. 2500 vibronic eigenstates of the NO 3 radical. Phys Chem Chem Phys 2024; 26:24506-24523. [PMID: 39283267 DOI: 10.1039/d4cp02653e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/26/2024]
Abstract
The nitrate radical NO3 plays an important role in atmospheric chemistry, yet many aspects of its coupled and anharmonic vibronic structure remain elusive. Here, using an accurate, coupled full-dimensional diabatic potential that includes five electronic states, we revisit the vibronic spectrum associated with the electronic state. Using recently developed tensor network state methods, we are able to compute more than 2500 vibronic states, thereby increasing the number of computed full-dimensional states by a factor of 50, compared to previous work. While we obtain good agreement with experiment for most of the assigned vibronic levels, for several others, we observe striking disagreement. Further, for the antisymmetric bending motion we find remarkably large symmetry-induced level splittings that are larger than the zero-order reference. We discuss non-negligible nonadiabatic effects and show that the Born-Oppenheimer approximation leads to significant errors in the spectrum.
Collapse
Affiliation(s)
- Henrik R Larsson
- Department of Chemistry and Biochemistry, University of California, Merced, CA 95343, USA. NO3a[at]larsson-research.δe
| | - Alexandra Viel
- Univ. Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France. alexandra.viel[at]univ-rennes1.fr
| |
Collapse
|
2
|
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
|
3
|
Mendolicchio M, Barone V. Accurate Vibrational and Ro-Vibrational Contributions to the Properties of Large Molecules by a New Engine Employing Curvilinear Internal Coordinates and Vibrational Perturbation Theory to Second Order. J Chem Theory Comput 2024. [PMID: 39215708 DOI: 10.1021/acs.jctc.4c00857] [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 unbiased comparison between theory and experiment requires approaches more sophisticated than the basic harmonic-oscillator rigid-rotor model, for taking into account vibrational averaging effects and ro-vibrational couplings in molecules of increasing size. Second-order vibrational perturbation theory based on curvilinear internal coordinates (ICs) offers a remarkable compromise between accuracy and computational cost, thanks to the reduction of mode-mode couplings with respect to their counterparts based on Cartesian coordinates. Therefore, we have developed, implemented, and validated a general engine employing ICs, which allows the accurate evaluation of vibrational averages and ro-vibrational couplings for molecules containing up to about 50 atoms beyond the harmonic approximation. After validation of the new tool for relatively small molecules, the effectiveness of ICs has been demonstrated for some flexible and/or quite large molecular bricks of life.
Collapse
|
4
|
Ammar A, Scemama A, Loos PF, Giner E. Compactification of determinant expansions via transcorrelation. J Chem Phys 2024; 161:084104. [PMID: 39171701 DOI: 10.1063/5.0217650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2024] [Accepted: 07/22/2024] [Indexed: 08/23/2024] Open
Abstract
Although selected configuration interaction (SCI) algorithms can tackle much larger Hilbert spaces than the conventional full CI method, the scaling of their computational cost with respect to the system size remains inherently exponential. In addition, inaccuracies in describing the correlation hole at small interelectronic distances lead to the slow convergence of the electronic energy relative to the size of the one-electron basis set. To alleviate these effects, we show that the non-Hermitian, transcorrelated (TC) version of SCI significantly compactifies the determinant space, allowing us to reach a given accuracy with a much smaller number of determinants. Furthermore, we note a significant acceleration in the convergence of the TC-SCI energy as the basis set size increases. The extent of this compression and the energy convergence rate are closely linked to the accuracy of the correlation factor used for the similarity transformation of the Coulombic Hamiltonian. Our systematic investigation of small molecular systems in increasingly large basis sets illustrates the magnitude of these effects.
Collapse
Affiliation(s)
- Abdallah Ammar
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Emmanuel Giner
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
| |
Collapse
|
5
|
Feldmann R, Reiher M. Renormalized Internally Contracted Multireference Coupled Cluster with Perturbative Triples. J Chem Theory Comput 2024; 20. [PMID: 39158160 PMCID: PMC11360144 DOI: 10.1021/acs.jctc.4c00679] [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/25/2024] [Revised: 07/17/2024] [Accepted: 07/22/2024] [Indexed: 08/20/2024]
Abstract
In this work, we combine the many-body formulation of the internally contracted multireference coupled cluster (ic-MRCC) method with Evangelista's multireference formulation of the driven similarity renormalization group (DSRG). The DSRG method can be viewed as a unitary multireference coupled cluster theory, which renormalizes the amplitudes based on a flow equation approach to eliminate numerical instabilities. We extend this approach by demonstrating that the unitary flow equation approach can be adapted for nonunitary transformations, rationalizing the renormalization of ic-MRCC amplitudes. We denote the new approach, the renormalized ic-MRCC (ric-MRCC) method. To achieve high accuracy with a reasonable computational cost, we introduce a new approximation to the Baker-Campbell-Hausdorff expansion. We fully consider the linear commutator while approximating the quadratic commutator, for which we neglect specific contractions involving amplitudes with active indices. Moreover, we introduce approximate perturbative triples to obtain the ric-MRCCSD[T] method. We demonstrate the accuracy of our approaches in comparison to advanced multireference methods for the potential energy curves of H8, F2, H2O, N2, and Cr2. Additionally, we show that ric-MRCCSD and ric-MRCSSD[T] match the accuracy of CCSD(T) for evaluating spectroscopic constants and of full configuration interaction energies for a set of small molecules.
Collapse
Affiliation(s)
- Robin Feldmann
- 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
|
6
|
Dwivedi A, Lopez-Ruiz MA, Iyengar SS. Resource Optimization for Quantum Dynamics with Tensor Networks: Quantum and Classical Algorithms. J Phys Chem A 2024; 128:6774-6797. [PMID: 39101545 DOI: 10.1021/acs.jpca.4c03407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/06/2024]
Abstract
The exponential scaling of the quantum degrees of freedom with the size of the system is one of the biggest challenges in computational chemistry and particularly in quantum dynamics. We present a tensor network approach for the time-evolution of the nuclear degrees of freedom of multiconfigurational chemical systems at a reduced storage and computational complexity. We also present quantum algorithms for the resultant dynamics. To preserve the compression advantage achieved via tensor network decompositions, we present an adaptive algorithm for the regularization of nonphysical bond dimensions, preventing the potentially exponential growth of these with time. While applicable to any quantum dynamical problem, our method is particularly valuable for dynamical simulations of nuclear chemical systems. Our algorithm is demonstrated using ab initio potentials obtained for a symmetric hydrogen-bonded system, namely, the protonated 2,2'-bipyridine, and compared to exact diagonalization numerical results.
Collapse
Affiliation(s)
- Anurag Dwivedi
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
- Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Miguel Angel Lopez-Ruiz
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
- Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
- Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
7
|
Braunscheidel NM, Bachhar A, Mayhall NJ. Accurate and interpretable representation of correlated electronic structure via Tensor Product Selected CI. Faraday Discuss 2024. [PMID: 39119803 DOI: 10.1039/d4fd00049h] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2024]
Abstract
The task of computing wavefunctions that are accurate, yet simple enough mathematical objects to use for reasoning, has long been a challenge in quantum chemistry. The difficulty in drawing physical conclusions from a wavefunction is often related to the generally large number of configurations with similar weights. In Tensor Product Selected Configuration Interaction (TPSCI), we use a locally correlated tensor product state basis, which has the effect of concentrating the weight of a state onto a smaller number of physically interpretable degrees of freedom. In this paper, we apply TPSCI to a series of three molecular systems ranging in separability, one of which is the first application of TPSCI to an open-shell bimetallic system. For each of these systems, we obtain accurate solutions to large active spaces, and analyze the resulting wavefunctions through a series of different approaches including (i) direct inspection of the TPS basis coefficients, (ii) construction of Bloch effective Hamiltonians, and (iii) computation of cluster correlation functions.
Collapse
Affiliation(s)
| | - Arnab Bachhar
- Department of Chemistry, Virginia Tech, Blacksburg, VA 24060, USA.
| | | |
Collapse
|
8
|
Xu Y, Liu C, Ma H. Kylin-V: An open-source package calculating the dynamic and spectroscopic properties of large systems. J Chem Phys 2024; 161:052501. [PMID: 39087896 DOI: 10.1063/5.0220712] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2024] [Accepted: 07/14/2024] [Indexed: 08/02/2024] Open
Abstract
Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples.
Collapse
Affiliation(s)
- Yihe Xu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Chungen Liu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, School of Chemistry and Chemical Engineering, Shandong University, Qingdao 266237, China
| |
Collapse
|
9
|
Abstract
The cumulant expansion of the Green's function is a computationally efficient beyond-GW approach renowned for its significant enhancement of satellite features in materials. In contrast to the ubiquitous GW approximation of many-body perturbation theory, ab initio cumulant expansions performed on top of GW (GW + C) have demonstrated the capability to handle multi-particle processes by incorporating higher-order correlation effects or vertex corrections, yielding better agreements between experiment and theory for satellite structures. While widely employed in condensed matter physics, very few applications of GW + C have been published on molecular systems. Here, we assess the performance of this scheme on a series of 10-electron molecular systems (Ne, HF, H2O, NH3, and CH4) where full configuration interaction estimates of the outer-valence quasiparticle and satellite energies are available.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France.
| | - Antoine Marie
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France.
| | - Abdallah Ammar
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France.
| |
Collapse
|
10
|
Choi S, Loaiza I, Lang RA, Martínez-Martínez LA, Izmaylov AF. Probing Quantum Efficiency: Exploring System Hardness in Electronic Ground State Energy Estimation. J Chem Theory Comput 2024; 20:5982-5993. [PMID: 38950444 DOI: 10.1021/acs.jctc.4c00298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/03/2024]
Abstract
We consider the question of how correlated the system hardness is between classical algorithms of electronic structure theory in ground state estimation and quantum algorithms. To define the system hardness for classical algorithms, we employ empirical criterion based on the deviation of electronic energies produced by coupled cluster and configuration interaction methods from the exact ones along multiple bonds dissociation in a set of molecular systems. For quantum algorithms, we have selected the Variational Quantum Eigensolver (VQE) and Quantum Phase Estimation (QPE) methods. As characteristics of the system hardness for quantum methods, we analyzed circuit depths for the state preparation, the number of quantum measurements needed for the energy expectation value, and various cost characteristics for the Hamiltonian encodings via Trotter approximation and linear combination of unitaries (LCU). Our results show that the quantum resource requirements are mostly unaffected by classical hardness, with the only exception being the state preparation part, which contributes to both VQE and QPE algorithm costs. However, there are clear indications that constructing the initial state with a significant overlap with the true ground state is easier than obtaining the state with an energy expectation value within chemical precision. These results support optimism regarding the identification of a molecular system where a quantum algorithm excels over its classical counterpart, as quantum methods can maintain efficiency in classically challenging systems.
Collapse
Affiliation(s)
- Seonghoon Choi
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Ignacio Loaiza
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
- Zapata Computing Canada Inc., Toronto, Ontario M5C 3A1, Canada
| | - Robert A Lang
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Luis A Martínez-Martínez
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Artur F Izmaylov
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
11
|
Pauletti T, Sanino M, Gimenes L, Carvalho IM, França VV. Quantum phase transitions in one-dimensional nanostructures: a comparison between DFT and DMRG methodologies. J Mol Model 2024; 30:268. [PMID: 39012396 DOI: 10.1007/s00894-024-06059-z] [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: 02/15/2024] [Accepted: 07/01/2024] [Indexed: 07/17/2024]
Abstract
CONTEXT In the realm of quantum chemistry, the accurate prediction of electronic structure and properties of nanostructures remains a formidable challenge. Density functional theory (DFT) and density matrix renormalization group (DMRG) have emerged as two powerful computational methods for addressing electronic correlation effects in diverse molecular systems. We compare ground-state energies (e 0 ), density profiles ( n ), and average entanglement entropies (S ¯ ) in metals, insulators and at the transition from metal to insulator, in homogeneous, superlattices, and harmonically confined chains described by the fermionic one-dimensional Hubbard model. While for the homogeneous systems, there is a clear hierarchy between the deviations, D % ( S ¯ ) < D % ( e 0 ) < D ¯ % ( n ) , and all the deviations decrease with the chain size; for superlattices and harmonic confinement, the relation among the deviations is less trivial and strongly dependent on the superlattice structure and the confinement strength considered. For the superlattices, in general, increasing the number of impurities in the unit cell represents lower precision in the DFT calculations. For the confined chains, DFT performs better for metallic phases, while the highest deviations appear for the Mott and band-insulator phases. This work provides a comprehensive comparative analysis of these methodologies, shedding light on their respective strengths, limitations, and applications. METHODS The DFT calculations were performed using the standard Kohn-Sham scheme within the BALDA approach. It integrated the numerical Bethe-Ansatz (BA) solution of the Hubbard model as the homogeneous density functional within a local-density approximation (LDA) for the exchange-correlation energy. The DMRG algorithms were implemented using the ITensor library, which is based on the matrix product states (MPS) ansatz. The calculations were performed until the energy reaches convergence of at least 10 - 8 .
Collapse
Affiliation(s)
- T Pauletti
- Institute of Chemistry, São Paulo State University, Francisco Degni 55, Araraquara, 14800-090, São Paulo, Brazil
| | - M Sanino
- Institute of Chemistry, São Paulo State University, Francisco Degni 55, Araraquara, 14800-090, São Paulo, Brazil
| | - L Gimenes
- Institute of Chemistry, São Paulo State University, Francisco Degni 55, Araraquara, 14800-090, São Paulo, Brazil
| | - I M Carvalho
- Institute of Chemistry, São Paulo State University, Francisco Degni 55, Araraquara, 14800-090, São Paulo, Brazil.
| | - V V França
- Institute of Chemistry, São Paulo State University, Francisco Degni 55, Araraquara, 14800-090, São Paulo, Brazil.
| |
Collapse
|
12
|
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
|
13
|
Aldossary A, Campos-Gonzalez-Angulo JA, Pablo-García S, Leong SX, Rajaonson EM, Thiede L, Tom G, Wang A, Avagliano D, Aspuru-Guzik A. In Silico Chemical Experiments in the Age of AI: From Quantum Chemistry to Machine Learning and Back. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2024; 36:e2402369. [PMID: 38794859 DOI: 10.1002/adma.202402369] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Revised: 04/28/2024] [Indexed: 05/26/2024]
Abstract
Computational chemistry is an indispensable tool for understanding molecules and predicting chemical properties. However, traditional computational methods face significant challenges due to the difficulty of solving the Schrödinger equations and the increasing computational cost with the size of the molecular system. In response, there has been a surge of interest in leveraging artificial intelligence (AI) and machine learning (ML) techniques to in silico experiments. Integrating AI and ML into computational chemistry increases the scalability and speed of the exploration of chemical space. However, challenges remain, particularly regarding the reproducibility and transferability of ML models. This review highlights the evolution of ML in learning from, complementing, or replacing traditional computational chemistry for energy and property predictions. Starting from models trained entirely on numerical data, a journey set forth toward the ideal model incorporating or learning the physical laws of quantum mechanics. This paper also reviews existing computational methods and ML models and their intertwining, outlines a roadmap for future research, and identifies areas for improvement and innovation. Ultimately, the goal is to develop AI architectures capable of predicting accurate and transferable solutions to the Schrödinger equation, thereby revolutionizing in silico experiments within chemistry and materials science.
Collapse
Affiliation(s)
- Abdulrahman Aldossary
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
| | | | - Sergio Pablo-García
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, ON, M5S 2E4, Canada
| | - Shi Xuan Leong
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
| | - Ella Miray Rajaonson
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- Vector Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto, ON, M5G 1M1, Canada
| | - Luca Thiede
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, ON, M5S 2E4, Canada
- Vector Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto, ON, M5G 1M1, Canada
| | - Gary Tom
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- Vector Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto, ON, M5G 1M1, Canada
| | - Andrew Wang
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
| | - Davide Avagliano
- Chimie ParisTech, PSL University, CNRS, Institute of Chemistry for Life and Health Sciences (iCLeHS UMR 8060), Paris, F-75005, France
| | - Alán Aspuru-Guzik
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, ON, M5S 2E4, Canada
- Vector Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto, ON, M5G 1M1, Canada
- Department of Materials Science & Engineering, University of Toronto, 184 College St., Toronto, ON, M5S 3E4, Canada
- Department of Chemical Engineering & Applied Chemistry, University of Toronto, 200 College St., Toronto, ON, M5S 3E5, Canada
- Lebovic Fellow, Canadian Institute for Advanced Research (CIFAR), 66118 University Ave., Toronto, M5G 1M1, Canada
- Acceleration Consortium, 80 St George St, Toronto, M5S 3H6, Canada
| |
Collapse
|
14
|
Choudhury A, Santra S, Ghosh D. Understanding the Photoprocesses in Biological Systems: Need for Accurate Multireference Treatment. J Chem Theory Comput 2024; 20:4951-4964. [PMID: 38864715 DOI: 10.1021/acs.jctc.4c00027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2024]
Abstract
Light-matter interaction is crucial to life itself and revolves around many of the central processes in biology. The need for understanding these photochemical and photophysical processes cannot be overemphasized. Interaction of light with biological systems starts with the absorption of light and subsequent phenomena that occur in the excited states of the system. However, excited states are typically difficult to understand within the mean field approximation of quantum chemical methods. Therefore, suitable multireference methods and methodologies have been developed to understand these phenomena. In this Perspective, we will describe a few methods and methodologies suitable for these descriptions and discuss some persisting difficulties.
Collapse
Affiliation(s)
- Arpan Choudhury
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Supriyo Santra
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Debashree Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
15
|
Li X, Huang JC, Zhang GZ, Li HE, Shen ZP, Zhao C, Li J, Hu HS. Improved optimization for the neural-network quantum states and tests on the chromium dimer. J Chem Phys 2024; 160:234102. [PMID: 38884396 DOI: 10.1063/5.0214150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2024] [Accepted: 05/30/2024] [Indexed: 06/18/2024] Open
Abstract
The advent of Neural-network Quantum States (NQS) has significantly advanced wave function ansatz research, sparking a resurgence in orbital space variational Monte Carlo (VMC) exploration. This work introduces three algorithmic enhancements to reduce computational demands of VMC optimization using NQS: an adaptive learning rate algorithm, constrained optimization, and block optimization. We evaluate the refined algorithm on complex multireference bond stretches of H2O and N2 within the cc-pVDZ basis set and calculate the ground-state energy of the strongly correlated chromium dimer (Cr2) in the Ahlrichs SV basis set. Our results achieve superior accuracy compared to coupled cluster theory at a relatively modest CPU cost. This work demonstrates how to enhance optimization efficiency and robustness using these strategies, opening a new path to optimize large-scale restricted Boltzmann machine-based NQS more effectively and marking a substantial advancement in NQS's practical quantum chemistry applications.
Collapse
Affiliation(s)
- Xiang Li
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Jia-Cheng Huang
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Guang-Ze Zhang
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Hao-En Li
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - 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
|
16
|
Weymuth T, Unsleber JP, Türtscher PL, Steiner M, Sobez JG, Müller CH, Mörchen M, Klasovita V, Grimmel SA, Eckhoff M, Csizi KS, Bosia F, Bensberg M, Reiher M. SCINE-Software for chemical interaction networks. J Chem Phys 2024; 160:222501. [PMID: 38857173 DOI: 10.1063/5.0206974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 05/09/2024] [Indexed: 06/12/2024] Open
Abstract
The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge.
Collapse
Affiliation(s)
- Thomas Weymuth
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan P Unsleber
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Paul L Türtscher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Miguel Steiner
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan-Grimo Sobez
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Charlotte H Müller
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Maximilian Mörchen
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Veronika Klasovita
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stephanie A Grimmel
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Marco Eckhoff
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Katja-Sophia Csizi
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Francesco Bosia
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Moritz Bensberg
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
17
|
Jensen PWK, Kjellgren ER, Reinholdt P, Ziems KM, Coriani S, Kongsted J, Sauer SPA. Quantum Equation of Motion with Orbital Optimization for Computing Molecular Properties in Near-Term Quantum Computing. J Chem Theory Comput 2024; 20:3613-3625. [PMID: 38701352 DOI: 10.1021/acs.jctc.4c00069] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2024]
Abstract
Determining the properties of molecules and materials is one of the premier applications of quantum computing. A major question in the field is how to use imperfect near-term quantum computers to solve problems of practical value. Inspired by the recently developed variants of the quantum counterpart of the equation-of-motion (qEOM) approach and the orbital-optimized variational quantum eigensolver (oo-VQE), we present a quantum algorithm (oo-VQE-qEOM) for the calculation of molecular properties by computing expectation values on a quantum computer. We perform noise-free quantum simulations of BeH2 in the series of STO-3G/6-31G/6-31G* basis sets and of H4 and H2O in 6-31G using an active space of four electrons and four spatial orbitals (8 qubits) to evaluate excitation energies, electronic absorption, and, for twisted H4, circular dichroism spectra. We demonstrate that the proposed algorithm can reproduce the results of conventional classical CASSCF calculations for these molecular systems.
Collapse
Affiliation(s)
- Phillip W K Jensen
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen, Denmark
| | - Erik Rosendahl Kjellgren
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense, Denmark
| | - Karl Michael Ziems
- Department of Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
| | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense, Denmark
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen, Denmark
| |
Collapse
|
18
|
Fitzpatrick A, Nykänen A, Talarico NW, Lunghi A, Maniscalco S, García-Pérez G, Knecht S. Self-Consistent Field Approach for the Variational Quantum Eigensolver: Orbital Optimization Goes Adaptive. J Phys Chem A 2024; 128:2843-2856. [PMID: 38547028 DOI: 10.1021/acs.jpca.3c05882] [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 present a self-consistent field (SCF) approach within the adaptive derivative-assembled problem-tailored ansatz variational quantum eigensolver (ADAPT-VQE) framework for efficient quantum simulations of chemical systems on near-term quantum computers. To this end, our ADAPT-VQE-SCF approach combines the idea of generating an ansatz with a small number of parameters, resulting in shallow-depth quantum circuits with a direct minimization of an energy expression that is correct to second order with respect to changes in the molecular orbital basis. Our numerical analysis, including calculations for the transition-metal complex ferrocene [Fe (C5H5)2], indicates that convergence in the self-consistent orbital optimization loop can be reached without a considerable increase in the number of two-qubit gates in the quantum circuit by comparison to a VQE optimization in the initial molecular orbital basis. Moreover, the orbital optimization can be carried out simultaneously within each iteration of the ADAPT-VQE cycle. ADAPT-VQE-SCF thus allows us to implement a routine analogous to the complete active space SCF, a cornerstone of state-of-the-art computational chemistry, in a hardware-efficient manner on near-term quantum computers. Hence, ADAPT-VQE-SCF paves the way toward a paradigm shift for quantitative quantum-chemistry simulations on quantum computers by requiring fewer qubits and opening up for the use of large and flexible atomic orbital basis sets in contrast to earlier methods that are predominantly based on the idea of full active spaces with minimal basis sets.
Collapse
Affiliation(s)
- Aaron Fitzpatrick
- Algorithmiq Ltd, Kanavakatu 3C, Helsinki FI-00160, Finland
- Trinity Quantum Alliance, Unit 16, Trinity Technology and Enterprise Centre, Pearse Street, Dublin 2 D02 YN67, Ireland
| | - Anton Nykänen
- Algorithmiq Ltd, Kanavakatu 3C, Helsinki FI-00160, Finland
| | | | - Alessandro Lunghi
- School of Physics, AMBER and CRANN Institute, Trinity College, Dublin 2, Ireland
| | | | | | - Stefan Knecht
- Algorithmiq Ltd, Kanavakatu 3C, Helsinki FI-00160, Finland
- ETH Zürich, Department of Chemistry and Applied Life Sciences Vladimir-Prelog-Weg 1-5/10, Zürich 8093, Switzerland
| |
Collapse
|
19
|
Kjellgren ER, Reinholdt P, Fitzpatrick A, Talarico WN, Jensen PWK, Sauer SPA, Coriani S, Knecht S, Kongsted J. The variational quantum eigensolver self-consistent field method within a polarizable embedded framework. J Chem Phys 2024; 160:124114. [PMID: 38533884 DOI: 10.1063/5.0190594] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Accepted: 03/10/2024] [Indexed: 03/28/2024] Open
Abstract
We formulate and implement the Variational Quantum Eigensolver Self Consistent Field (VQE-SCF) algorithm in combination with polarizable embedding (PE), thereby extending PE to the regime of quantum computing. We test the resulting algorithm, PE-VQE-SCF, on quantum simulators and demonstrate that the computational stress on the quantum device is only slightly increased in terms of gate counts compared to regular VQE-SCF. On the other hand, no increase in shot noise was observed. We illustrate how PE-VQE-SCF may lead to the modeling of real chemical systems using a simulation of the reaction barrier of the Diels-Alder reaction between furan and ethene as an example.
Collapse
Affiliation(s)
- Erik Rosendahl Kjellgren
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | | | - Walter N Talarico
- Algorithmiq Ltd., Kanavakatu 3C, FI-00160 Helsinki, Finland
- Department of Applied Physics, QTF Centre of Excellence, Center for Quantum Engineering, Aalto University School of Science, FIN-00076 AALTO Espoo, Finland
| | - Phillip W K Jensen
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Sonia Coriani
- DTU Chemistry - Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Stefan Knecht
- Algorithmiq Ltd., Kanavakatu 3C, FI-00160 Helsinki, Finland
- ETH Zürich, Department of Chemistry and Applied Life Sciences, Vladimir-Prelog-Weg 1-5/10, CH-8093 Zürich, Switzerland
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| |
Collapse
|
20
|
Cheng Y, Ma H. Renormalized-Residue-Based Multireference Configuration Interaction Method for Strongly Correlated Systems. J Chem Theory Comput 2024; 20:1988-2009. [PMID: 38380619 DOI: 10.1021/acs.jctc.3c01247] [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
The implementation of multireference configuration interaction (MRCI) methods in quantum systems with large active spaces is hindered by the expansion of configuration bases or the intricate handling of reduced density matrices (RDMs). In this work, we present a spin-adapted renormalized-residue-based MRCI (RR-MRCI) approach that leverages renormalized residues to effectively capture the entanglement between active and inactive orbitals. This approach is reinforced by a novel efficient algorithm, which also facilitates an efficient deployment of spin-adapted matrix product state MRCI (MPS-MRCI). The RR-MRCI framework possesses several advantages: (1) It considers the orbital entanglement and utilizes highly compressed MPS structure, improving computational accuracy and efficiency compared with internally contracted (ic) MRCI. (2) Utilizing small-sized buffer environments of a few external orbitals as probes based on quantum information theory, it enhances computational efficiency over MPS-MRCI and offers potential application to large molecular systems. (3) The RR framework can be implemented in conjunction with ic-MRCI, eliminating the need for high-rank RDMs, by using distinct renormalized residues. We evaluated this method across nine diverse molecular systems, including Cu2O22+ with an active space of (24e,24o) and two complexes of lanthanide and actinide with active space (38e,36o), demonstrating the method's versatility and efficacy.
Collapse
Affiliation(s)
- Yifan Cheng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, School of Chemistry and Chemical Engineering, Shandong University, Qingdao, Shandong 266237, China
| |
Collapse
|
21
|
Christopoulou G, Di Paola C, Elzinga FE, Jallat A, Muñoz Ramo D, Krompiec M. Quantum hardware calculations of the activation and dissociation of nitrogen on iron clusters and surfaces. Phys Chem Chem Phys 2024; 26:5895-5906. [PMID: 38315092 DOI: 10.1039/d3cp05167f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2024]
Abstract
Catalytic processes are the cornerstone of chemical industry, and catalytic conversion of nitrogen to ammonia remains one of the largest industrial processes implemented. Rational design of catalysts and catalytic reactions largely depends on approximate computational chemistry methods, such as density functional theory, which, however, suffer from limited accuracy, especially for strongly-correlated materials. Rigorous ab initio methods which account for static and dynamic electron correlation, while arbitrarily accurate for small systems, are generally too expensive to be applied to modelling of catalytic cycles, due to prohibitive time and space computational complexity with respect to the size of the active space. Recent advances in quantum computing give hope for enabling access to accurate ab initio methods at scale. Herein, we present a prototype hybrid quantum-classical workflow for modeling chemical reactions on surfaces, applied to proof-of-concept models of activation and dissociation of nitrogen on small Fe clusters and a single-layer (221) iron surface. First, we determined the structures of species present in the catalytic cycle at DFT level and studied their electronic structure using CASSCF. We show that it is possible to decouple the half-filled Fe-3d band from the Fe-N and N-N bond orbitals, thereby reducing the active space significantly. Subsequently, we translated the CASSCF wavefunctions into corresponding qubit quantum states, using the Adaptive Variational Quantum Eigensolver, and estimated their energies using a state vector simulator, H1-1E quantum emulator and (for selected systems) H1-1 quantum computer. We demonstrated that if a sufficiently small active orbital space is chosen, ground state energies obtained with classical methods and with the quantum computer are in reasonable agreement. We argue that once quantum computing methods are scaled up so that larger active spaces are accessible, they can offer a tremendous practical advantage to the computational catalysis community.
Collapse
Affiliation(s)
| | - Cono Di Paola
- Quantinuum, Terrington House, 13-15 Hills Road, CB2 1NL, Cambridge, UK.
| | | | | | - David Muñoz Ramo
- Quantinuum, Terrington House, 13-15 Hills Road, CB2 1NL, Cambridge, UK.
| | - Michal Krompiec
- Quantinuum, Terrington House, 13-15 Hills Road, CB2 1NL, Cambridge, UK.
| |
Collapse
|
22
|
Attal L, Falvo C, Calvo F, Parneix P. Modeling the dynamics of quantum systems coupled to large-dimensional baths using effective energy states. J Chem Phys 2024; 160:044107. [PMID: 38270236 DOI: 10.1063/5.0184299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Accepted: 12/14/2023] [Indexed: 01/26/2024] Open
Abstract
The quantum dynamics of a low-dimensional system in contact with a large but finite harmonic bath is theoretically investigated by coarse-graining the bath into a reduced set of effective energy states. In this model, the couplings between the system and the bath are obtained from statistically averaging over the discrete, degenerate effective states. Our model is aimed at intermediate bath sizes in which non-Markovian processes and energy transfer between the bath and the main system are important. The method is applied to a model system of a Morse oscillator coupled to 40 harmonic modes. The results are found to be in excellent agreement with the direct quantum dynamics simulations presented in the work of Bouakline et al. [J. Phys. Chem. A 116, 11118-11127 (2012)], but at a much lower computational cost. Extension to larger baths is discussed in comparison to the time-convolutionless method. We also extend this study to the case of a microcanonical bath with finite initial internal energies. The computational efficiency and convergence properties of the effective bath states model with respect to relevant parameters are also discussed.
Collapse
Affiliation(s)
- Loïse Attal
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Cyril Falvo
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
- Université Grenoble-Alpes, CNRS, LiPhy, 38000 Grenoble, France
| | - Florent Calvo
- Université Grenoble-Alpes, CNRS, LiPhy, 38000 Grenoble, France
| | - Pascal Parneix
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
23
|
Xiang C, Jia W, Fang WH, Li Z. Distributed Multi-GPU Ab Initio Density Matrix Renormalization Group Algorithm with Applications to the P-Cluster of Nitrogenase. J Chem Theory Comput 2024; 20:775-786. [PMID: 38198503 DOI: 10.1021/acs.jctc.3c01228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
The presence of many degenerate d/f orbitals makes polynuclear transition-metal compounds, such as iron-sulfur clusters in nitrogenase, challenging for state-of-the-art quantum chemistry methods. To address this challenge, we present the first distributed multi-graphics processing unit (GPU) ab initio density matrix renormalization group (DMRG) algorithm suitable for modern high-performance computing (HPC) infrastructures. The central idea is to parallelize the most computationally intensive part─the multiplication of O(K2) operators with a trial wave function, where K is the number of spatial orbitals, by combining operator parallelism for distributing the workload with a batched algorithm for performing contractions on GPU. With this new implementation, we are able to reach an unprecedentedly large bond dimension D = 14,000 on 48 GPUs (NVIDIA A100 80 GB SXM) for an active space model (114 electrons in 73 active orbitals) of the P-cluster, which is nearly 3 times larger than the bond dimensions reported in previous DMRG calculations for the same system using only central processing units (CPUs).
Collapse
Affiliation(s)
- Chunyang Xiang
- State Key Lab of Processors, Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100190, China
| | - Weile Jia
- State Key Lab of Processors, Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100190, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
24
|
Friesecke G, Barcza G, Legeza Ö. Predicting the FCI Energy of Large Systems to Chemical Accuracy from Restricted Active Space Density Matrix Renormalization Group Calculations. J Chem Theory Comput 2024; 20:87-102. [PMID: 38109339 DOI: 10.1021/acs.jctc.3c01001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
We theoretically derive and validate with large scale simulations a remarkably accurate power law scaling of errors for the restricted active space density matrix renormalization group (DMRG-RAS) method [J. Phys. Chem. A 126, 9709] in electronic structure calculations. This yields a new extrapolation method, DMRG-RAS-X, which reaches chemical accuracy for strongly correlated systems such as the chromium dimer, dicarbon up to a large cc-pVQZ basis and even a large chemical complex such as the FeMoco with significantly lower computational demands than those of previous methods. The method is free of empirical parameters, performed robustly and reliably in all examples we tested, and has the potential to become a vital alternative method for electronic structure calculations in quantum chemistry and more generally for the computation of strong correlations in nuclear and condensed matter physics.
Collapse
Affiliation(s)
- Gero Friesecke
- Department of Mathematics, Technical University of Munich, München 85748, Germany
| | - Gergely Barcza
- Strongly Correlated Systems Lendület Research Group, Wigner Research Centre for Physics, Budapest H-1525, Hungary
| | - Örs Legeza
- Strongly Correlated Systems Lendület Research Group, Wigner Research Centre for Physics, Budapest H-1525, Hungary
- Institute for Advanced Study, Technical University of Munich, Germany, Lichtenbergstrasse 2a, Garching 85748, Germany
| |
Collapse
|
25
|
Glaser N, Baiardi A, Reiher M. Flexible DMRG-Based Framework for Anharmonic Vibrational Calculations. J Chem Theory Comput 2023; 19:9329-9343. [PMID: 38060309 PMCID: PMC10753801 DOI: 10.1021/acs.jctc.3c00902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 10/16/2023] [Accepted: 10/23/2023] [Indexed: 12/08/2023]
Abstract
We present a novel formulation of the vibrational density matrix renormalization group (vDMRG) algorithm tailored to strongly anharmonic molecules described by general, high-dimensional model representations of potential energy surfaces. For this purpose, we extend the vDMRG framework to support vibrational Hamiltonians expressed in the so-called n-mode second-quantization formalism. The resulting n-mode vDMRG method offers full flexibility with respect to both the functional form of the PES and the choice of the single-particle basis set. We leverage this framework to apply, for the first time, vDMRG based on an anharmonic modal basis set optimized with the vibrational self-consistent field algorithm on an on-the-fly constructed PES. We also extend the n-mode vDMRG framework to include excited-state-targeting algorithms in order to efficiently calculate anharmonic transition frequencies. We demonstrate the capabilities of our novel n-mode vDMRG framework for methyloxirane, a challenging molecule with 24 coupled vibrational modes.
Collapse
Affiliation(s)
- Nina Glaser
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Alberto Baiardi
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
26
|
Zhai H, Larsson HR, Lee S, Cui ZH, Zhu T, Sun C, Peng L, Peng R, Liao K, Tölle J, Yang J, Li S, Chan GKL. Block2: A comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond. J Chem Phys 2023; 159:234801. [PMID: 38108484 DOI: 10.1063/5.0180424] [Citation(s) in RCA: 13] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 11/16/2023] [Indexed: 12/19/2023] Open
Abstract
block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.
Collapse
Affiliation(s)
- Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Henrik R Larsson
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Zhi-Hao Cui
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Tianyu Zhu
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Chong Sun
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Linqing Peng
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Ruojing Peng
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Ke Liao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Johannes Tölle
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Junjie Yang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Shuoxue Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
27
|
Ugandi M, Roemelt M. A configuration-based heatbath-CI for spin-adapted multireference electronic structure calculations with large active spaces. J Comput Chem 2023; 44:2374-2390. [PMID: 37589287 DOI: 10.1002/jcc.27203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Revised: 07/24/2023] [Accepted: 07/27/2023] [Indexed: 08/18/2023]
Abstract
This work reports on a spin-pure configuration-based implementation of the heatbath configuration interaction (HCI) algorithm for selective configuration interaction. Besides the obvious advantage of being spin-pure, the presented method combines the compactness of the configurational ansatz with the known efficiency of the HCI algorithm and a variety of algorithmic and conceptual ideas to achieve a high level of performance. In particular, through pruning of the selected configurational space after HCI selection by means of a more strict criterion, a more compact wavefunction representation is obtained. Moreover, the underlying logic of the method allows us to minimize the number of redundant matrix-matrix multiplications while making use of just-in-time compilation to achieve fast diagonalization of the Hamiltonian. The critical search for 2-electron connections within the configurational space is facilitated by a tree-based representation thereof as suggested previously by Gopal et al. Usage of a prefix-based parallelization and batching during the calculation of the PT2-correction leads to a good load balancing and significantly reduced memory requirements for these critical steps of the calculation. In this way, the need for a semistochastic approach to the PT2 correction is avoided even for large configurational spaces. Finally, several test-cases are discussed to demonstrate the strengths and weaknesses of the presented method.
Collapse
Affiliation(s)
- Mihkel Ugandi
- Institut für Chemie, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Michael Roemelt
- Institut für Chemie, Humboldt-Universität zu Berlin, Berlin, Germany
| |
Collapse
|
28
|
Ding L, Knecht S, Schilling C. Quantum Information-Assisted Complete Active Space Optimization (QICAS). J Phys Chem Lett 2023:11022-11029. [PMID: 38047727 DOI: 10.1021/acs.jpclett.3c02536] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2023]
Abstract
We propose an effective quantum information-assisted complete active space optimization scheme (QICAS). What sets QICAS apart from other correlation-based selection schemes is (i) the use of unique measures from quantum information that assess the correlation in electronic structures in an unambiguous and predictive manner and (ii) an orbital optimization step that minimizes the correlation discarded by the active space approximation. Equipped with these features, QICAS yields, for smaller correlated molecule, sets of optimized orbitals with respect to which the complete active space configuration interaction energy reaches the corresponding complete active space self-consistent field (CASSCF) energy within chemical accuracy. For more challenging systems such as the chromium dimer, QICAS offers an excellent starting point for CASSCF by greatly reducing the number of iterations required for numerical convergence. Accordingly, our study validates a profound empirical conjecture: the energetically optimal nonactive spaces are predominantly those that contain the least entanglement.
Collapse
Affiliation(s)
- 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
| | - Stefan Knecht
- Algorithmiq Ltd., Kanavakatu 3C, FI-00160 Helsinki, Finland
- ETH Zürich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - 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
|
29
|
Zhai H, Lee S, Cui ZH, Cao L, Ryde U, Chan GKL. Multireference Protonation Energetics of a Dimeric Model of Nitrogenase Iron-Sulfur Clusters. J Phys Chem A 2023; 127:9974-9984. [PMID: 37967028 PMCID: PMC10694817 DOI: 10.1021/acs.jpca.3c06142] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Revised: 10/31/2023] [Accepted: 11/01/2023] [Indexed: 11/17/2023]
Abstract
Characterizing the electronic structure of the iron-sulfur clusters in nitrogenase is necessary to understand their role in the nitrogen fixation process. One challenging task is to determine the protonation state of the intermediates in the nitrogen fixing cycle. Here, we use a dimeric iron-sulfur model to study relative energies of protonation at C, S, or Fe. Using a composite method based on coupled cluster and density matrix renormalization group energetics, we converge the relative energies of four protonated configurations with respect to basis set and correlation level. We find that accurate relative energies require large basis sets as well as a proper treatment of multireference and relativistic effects. We have also tested ten density functional approximations for these systems. Most of them give large errors in their relative energies. The best performing functional in this system is B3LYP, which gives mean absolute and maximum deviations of only 10 and 13 kJ/mol with respect to our correlated wave function estimates, respectively, comparable to the uncertainty in our correlated estimates. Our work provides benchmark results for the calibration of new approximate electronic structure methods and density functionals for these problems.
Collapse
Affiliation(s)
- Huanchen Zhai
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Seunghoon Lee
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Zhi-Hao Cui
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Lili Cao
- Department
of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Ulf Ryde
- Department
of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
30
|
King DS, Truhlar DG, Gagliardi L. Variational Active Space Selection with Multiconfiguration Pair-Density Functional Theory. J Chem Theory Comput 2023; 19:8118-8128. [PMID: 37905518 DOI: 10.1021/acs.jctc.3c00792] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
The selection of an adequate set of active orbitals for modeling strongly correlated electronic states is difficult to automate because it is highly dependent on the states and molecule of interest. Although many approaches have shown some success, no single approach has worked well in all cases. In light of this, we present the "discrete variational selection" (DVS) approach to active space selection, in which one generates multiple trial wave functions from a diverse set of systematically constructed active spaces and then selects between these wave functions variationally. We apply this DVS approach to 207 vertical excitations of small-to-medium-sized organic and inorganic molecules (with 3 to 18 atoms) in the QUESTDB database by (i) constructing various sets of active space orbitals through diagonalization of parametrized operators and (ii) choosing the result with the lowest average energy among the states of interest. This approach proves ineffective when variationally selecting between wave functions using the density matrix renormalization group (DMRG) or complete active space self-consistent field (CASSCF) energy but is able to provide good results when variationally selecting between wave functions using the energy of the translated PBE (tPBE) functional from multiconfiguration pair-density functional theory (MC-PDFT). Applying this DVS-tPBE approach to selection among state-averaged DMRG wave functions, we obtain a mean unsigned error of only 0.17 eV using hybrid MC-PDFT. This result matches that of our previous benchmark without the need to filter out poor active spaces and with no further orbital optimization following active space selection of the SA-DMRG wave functions. Furthermore, we find that DVS-tPBE is able to robustly and effectively select between the new SA-DMRG wave functions and our previous SA-CASSCF results.
Collapse
Affiliation(s)
- Daniel S King
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Group, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
31
|
Wang Z, Zhang Z, Lu J, Li Y. Coordinate Descent Full Configuration Interaction for Excited States. J Chem Theory Comput 2023; 19:7731-7739. [PMID: 37870778 DOI: 10.1021/acs.jctc.3c00452] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2023]
Abstract
An efficient excited state method, named xCDFCI, in the configuration interaction framework is proposed. xCDFCI extends the unconstrained nonconvex optimization problem in coordinate descent full configuration interaction (CDFCI) to a multicolumn version for low-lying excited states computation. The optimization problem is addressed via a tailored coordinate descent method. In each iteration, a determinant is selected based on an approximated gradient, and coefficients of all states associated with the selected determinant are updated. A deterministic compression is applied to limit memory usage. We test xCDFCI applied to H2O and N2 molecules under the cc-pVDZ basis set. For both systems, five low-lying excited states in the same symmetry sector are calculated, together with the ground state. xCDFCI also produces accurate binding curves of the carbon dimer in the cc-pVDZ basis with chemical accuracy, where the ground state and four excited states in the same symmetry sector are benchmarked.
Collapse
Affiliation(s)
- Zhe Wang
- Department of Mathematics, Duke University, Durham, North Carolina 27708-0187, United States
| | - Zhiyuan Zhang
- School of Future Technology, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jianfeng Lu
- Department of Mathematics, Duke University, Durham, North Carolina 27708-0187, United States
- Department of Chemistry and Department of Physics, Duke University, Durham, North Carolina 27708-0187, United States
| | - Yingzhou Li
- School of Mathematical Sciences, Fudan University, Shanghai 200433, China
| |
Collapse
|
32
|
Mi W, Luo K, Trickey SB, Pavanello M. Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations. Chem Rev 2023; 123:12039-12104. [PMID: 37870767 DOI: 10.1021/acs.chemrev.2c00758] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2023]
Abstract
Kohn-Sham Density Functional Theory (KSDFT) is the most widely used electronic structure method in chemistry, physics, and materials science, with thousands of calculations cited annually. This ubiquity is rooted in the favorable accuracy vs cost balance of KSDFT. Nonetheless, the ambitions and expectations of researchers for use of KSDFT in predictive simulations of large, complicated molecular systems are confronted with an intrinsic computational cost-scaling challenge. Particularly evident in the context of first-principles molecular dynamics, the challenge is the high cost-scaling associated with the computation of the Kohn-Sham orbitals. Orbital-free DFT (OFDFT), as the name suggests, circumvents entirely the explicit use of those orbitals. Without them, the structural and algorithmic complexity of KSDFT simplifies dramatically and near-linear scaling with system size irrespective of system state is achievable. Thus, much larger system sizes and longer simulation time scales (compared to conventional KSDFT) become accessible; hence, new chemical phenomena and new materials can be explored. In this review, we introduce the historical contexts of OFDFT, its theoretical basis, and the challenge of realizing its promise via approximate kinetic energy density functionals (KEDFs). We review recent progress on that challenge for an array of KEDFs, such as one-point, two-point, and machine-learnt, as well as some less explored forms. We emphasize use of exact constraints and the inevitability of design choices. Then, we survey the associated numerical techniques and implemented algorithms specific to OFDFT. We conclude with an illustrative sample of applications to showcase the power of OFDFT in materials science, chemistry, and physics.
Collapse
Affiliation(s)
- Wenhui Mi
- Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, PR China
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, PR China
- International Center of Future Science, Jilin University, Changchun 130012, PR China
| | - Kai Luo
- Department of Applied Physics, Nanjing University of Science and Technology, Nanjing 210094, PR China
| | - S B Trickey
- Quantum Theory Project, Department of Physics and Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Michele Pavanello
- Department of Physics and Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
33
|
Li Manni G, Fdez. Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry. J Chem Theory Comput 2023; 19:6933-6991. [PMID: 37216210 PMCID: PMC10601490 DOI: 10.1021/acs.jctc.3c00182] [Citation(s) in RCA: 64] [Impact Index Per Article: 64.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Indexed: 05/24/2023]
Abstract
The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.
Collapse
Affiliation(s)
- Giovanni Li Manni
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ignacio Fdez. Galván
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Ali Alavi
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Yusuf Hamied
Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Flavia Aleotti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Francesco Aquilante
- Theory and
Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jochen Autschbach
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
| | - Davide Avagliano
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Alberto Baiardi
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jie J. Bao
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Stefano Battaglia
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Letitia Birnoschi
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Alejandro Blanco-González
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Sergey I. Bokarev
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
- Chemistry
Department, School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| | - Ria Broer
- Theoretical
Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Roberto Cacciari
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Paul B. Calio
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Rebecca K. Carlson
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Rafael Carvalho Couto
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luis Cerdán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
- Instituto
de Óptica (IO−CSIC), Consejo
Superior de Investigaciones Científicas, 28006, Madrid, Spain
| | - Liviu F. Chibotaru
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Nicholas F. Chilton
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | | | - Irene Conti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Sonia Coriani
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Juliana Cuéllar-Zuquin
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Razan E. Daoud
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Nike Dattani
- HPQC Labs, Waterloo, N2T 2K9 Ontario Canada
- HPQC College, Waterloo, N2T 2K9 Ontario Canada
| | - Piero Decleva
- Istituto
Officina dei Materiali IOM-CNR and Dipartimento di Scienze Chimiche
e Farmaceutiche, Università degli
Studi di Trieste, I-34121 Trieste, Italy
| | - Coen de Graaf
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
- ICREA, Pg. Lluís
Companys 23, 08010 Barcelona, Spain
| | - Mickaël
G. Delcey
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luca De Vico
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Werner Dobrautz
- Chalmers
University of Technology, Department of Chemistry
and Chemical Engineering, 41296 Gothenburg, Sweden
| | - Sijia S. Dong
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Chemical Biology, Department of Physics, and Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Rulin Feng
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Department
of Chemistry, Fudan University, Shanghai 200433, China
| | - Nicolas Ferré
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | | | - Laura Gagliardi
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Garavelli
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Yafu Guan
- State Key
Laboratory of Molecular Reaction Dynamics and Center for Theoretical
Computational Chemistry, Dalian Institute
of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Meiyuan Guo
- SSRL, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Matthew R. Hennefarth
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Chad E. Hoyer
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Miquel Huix-Rotllant
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | - Vishal Kumar Jaiswal
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Andy Kaiser
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Danil S. Kaliakin
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Marjan Khamesian
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Daniel S. King
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Vladislav Kochetov
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Marek Krośnicki
- Institute
of Theoretical Physics and Astrophysics, Faculty of Mathematics, Physics
and Informatics, University of Gdańsk, ul Wita Stwosza 57, 80-952, Gdańsk, Poland
| | | | - Ernst D. Larsson
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Susi Lehtola
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 University of Helsinki, Finland
| | - Marie-Bernadette Lepetit
- Condensed
Matter Theory Group, Institut Néel, CNRS UPR 2940, 38042 Grenoble, France
- Theory
Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409-1061, United States
| | - Pablo López Ríos
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Marcus Lundberg
- Department
of Chemistry − Ångström Laboratory, Uppsala University, SE-75120 Uppsala, Sweden
| | - Dongxia Ma
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | | | - Francesco Montorsi
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Maximilian Mörchen
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Artur Nenov
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Vu Ha Anh Nguyen
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Yoshio Nishimoto
- Graduate
School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Meagan S. Oakley
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Markus Oppel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Daniele Padula
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Riddhish Pandharkar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Quan Manh Phung
- Department
of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
- Institute
of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
| | - Felix Plasser
- Department
of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K.
| | - Gerardo Raggi
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Quantum
Materials and Software LTD, 128 City Road, London, EC1V 2NX, United Kingdom
| | - Elisa Rebolini
- Scientific
Computing Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Markus Reiher
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Ivan Rivalta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Daniel Roca-Sanjuán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Thies Romig
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Arta Anushirwan Safari
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Aitor Sánchez-Mansilla
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
| | - Andrew M. Sand
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Biochemistry, Butler University, Indianapolis, Indiana 46208, United States
| | - Igor Schapiro
- Institute
of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Thais R. Scott
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Javier Segarra-Martí
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Francesco Segatta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Dumitru-Claudiu Sergentu
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Laboratory
RA-03, RECENT AIR, A. I. Cuza University of Iaşi, RA-03 Laboratory (RECENT AIR), Iaşi 700506, Romania
| | - Prachi Sharma
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, USA
| | - Yinan Shu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Jakob K. Staab
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Tjerk P. Straatsma
- National
Center for Computational Sciences, Oak Ridge
National Laboratory, Oak Ridge, Tennessee 37831-6373, United States
- Department
of Chemistry and Biochemistry, University
of Alabama, Tuscaloosa, Alabama 35487-0336, United States
| | | | - Bruno Nunes Cabral Tenorio
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Donald G. Truhlar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Liviu Ungur
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Morgane Vacher
- Nantes
Université, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| | - Valera Veryazov
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Torben Arne Voß
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Oskar Weser
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Dihua Wu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Xuchun Yang
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - David Yarkony
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Chen Zhou
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - J. Patrick Zobel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Roland Lindh
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Uppsala
Center for Computational Chemistry (UC3), Uppsala University, PO Box 576, SE-751 23 Uppsala. Sweden
| |
Collapse
|
34
|
Kim SY, Park JW. Approximate Excited-State Geometry Optimization with the State-Averaged Adaptive Sampling Configuration Interaction Algorithm with Large Active Spaces. J Chem Theory Comput 2023; 19:7260-7272. [PMID: 37800852 DOI: 10.1021/acs.jctc.3c00808] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/07/2023]
Abstract
The selected configuration interaction (SCI) wave function is a useful approximation to the full configuration interaction (FCI) one. The adaptive sampling CI (ASCI) method is a deterministic SCI method. By combining ASCI and orbital optimization, the ASCI self-consistent field (ASCI-SCF) method, which is an approximation of the complete active space self-consistent field (CASSCF) method, can be formulated as well. However, their applicability has been tested mainly on the systems in their electronically ground states. In this work, we implement the state-average (SA) ansatz in ASCI-SCF calculations to calculate excited states. We also derive expressions for the approximate analytical gradient and implement them as a computer program. We demonstrate the applicability of the current method for calculating vertical and adiabatic excitation energies and optimizing the molecular geometries of thermally activated delayed fluorescence (TADF) molecules.
Collapse
Affiliation(s)
- So Yeon Kim
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| | - Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
35
|
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: 0] [Impact Index Per Article: 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
|
36
|
Nyvang A, Olsen J. A relativistic configuration interaction method with general expansions and initial applications to electronic g-factors. J Chem Phys 2023; 159:044102. [PMID: 37486047 DOI: 10.1063/5.0152655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 06/28/2023] [Indexed: 07/25/2023] Open
Abstract
A new implementation of the orbital-based two-component relativistic configuration interaction approach is reported and applied to calculations of the electronic g-shifts of three diatomic radicals: AlO, HgF, and PdH. The new implementation augments efficient routines for the calculation of nonrelativistic Hamiltonians with new vectorized routines for the calculation of the action of the one-electron spin-orbit operator and allows efficient calculations for the expansion of generalized active space type. The program makes full use of double group as well as time-reversal symmetry. Particle-hole reorganization of the operators is used to improve the efficiency for expansions with nearly fully occupied orbital spaces. The flexibility of the algorithm and program is used to investigate the convergence of electronic g-shifts for the three diatomic radicals as functions of the active space, states included in the orbital optimization, and excitation levels. It was possible to converge to the valence limits within a few percent using expansions containing up to quadruple excitations. However, when excitations from the core orbitals were added, it was not possible to demonstrate convergence to within a few percent with expansions containing at most 10 × 109 determinants.
Collapse
Affiliation(s)
- Andreas Nyvang
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK 8000 Aarhus C, Denmark
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK 8000 Aarhus C, Denmark
| |
Collapse
|
37
|
Xu Y, Cheng Y, Song Y, Ma H. New Density Matrix Renormalization Group Approaches for Strongly Correlated Systems Coupled with Large Environments. J Chem Theory Comput 2023. [PMID: 37471519 DOI: 10.1021/acs.jctc.2c01316] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2023]
Abstract
Thanks to the high compression of the matrix product state (MPS) form of the wave function and the efficient site-by-site iterative sweeping optimization algorithm, the density matrix normalization group (DMRG) and its time-dependent variant (TD-DMRG) have been established as powerful computational tools in accurately simulating the electronic structure and quantum dynamics of strongly correlated molecules with a large number (101-2) of quantum degrees of freedom (active orbitals or vibrational modes). However, the quantitative characterization of the quantum many-body behaviors of realistic strongly correlated systems requires a further consideration of the interaction between the embedded active subsystem and the remaining correlated environment, e.g., a larger number (102-3) of external orbitals in electronic structure or infinite condensed-phase phononic modes in nucleus dynamics. To this end, we introduced three new post-DMRG and TD-DMRG approaches, namely (1) DMRG2sCI-MRCI and DMRG2sCI-ENPT by the reconstruction of selected configuration interaction (sCI) type of compact reference function from DMRG coefficients and the use of externally contracted MRCI (multireference configuration interaction) and Epstein-Nesbet perturbation theory (ENPT), without recourse to the expensive high order n-electron reduced density matrices (n-RDMs). (2) DMRG combined with RR-MRCI (renormalized residue-based MRCI), which improves the computational accuracy and efficiency of internally contracted (ic) MRCI by renormalizing the contracted bases with small-sized buffer environment(s) of a few external orbitals as probes based on quantum information theory. (3) HM (hierarchical mapping)-TD-DMRG in which a large environment is reduced to a small number of renormalized environmental modes (which accounts for the most vital system-environment interactions) through stepwise mapping transformation. These advances extend the efficacy of highly accurate DMRG/TD-DMRG computations to the quantitative characterization of the electronic structure and quantum dynamics in realistic strongly correlated systems coupled with large environments and are reviewed in this paper.
Collapse
Affiliation(s)
- Yihe Xu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yifan Cheng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yinxuan Song
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, China
| |
Collapse
|
38
|
Motta LC, Autschbach J. Actinide inverse trans influence versus cooperative pushing from below and multi-center bonding. Nat Commun 2023; 14:4307. [PMID: 37463900 DOI: 10.1038/s41467-023-39626-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Accepted: 06/22/2023] [Indexed: 07/20/2023] Open
Abstract
Actinide-ligand bonds with high multiplicities remain poorly understood. Decades ago, an effect known as 6p pushing from below (PFB) was proposed to enhance actinide covalency. A related effect-also poorly understood-is inverse trans influence (ITI). The present computational study of actinide-ligand covalent interactions with high bond multiplicities quantifies the energetic contributions from PFB and identifies a hitherto overlooked fourth bonding interaction for 2nd-row ligands in the studied organometallic systems. The latter are best described by a terminal O/N ligand exhibiting quadruple bonding interactions with the actinide. The 4th interaction may be characterized as a multi-center or charge-shift bond involving the trans ligand. It is shown in this work that the 4th bonding interaction is a manifestation of ITI, assisted by PFB, and provides a long-sought missing piece in the understanding of actinide chemistry.
Collapse
Affiliation(s)
- Laura C Motta
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, NY, 14260-3000, USA
- Department of Marine Chemistry & Geochemistry, Woods Hole Oceanographic Institution, Woods Hole, MA, 02543-1050, USA
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, NY, 14260-3000, USA.
| |
Collapse
|
39
|
Han R, Luber S, Li Manni G. Magnetic Interactions in a [Co(II) 3Er(III)(OR) 4] Model Cubane through Forefront Multiconfigurational Methods. J Chem Theory Comput 2023; 19:2811-2826. [PMID: 37126736 DOI: 10.1021/acs.jctc.2c01318] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Strong electron correlation effects are one of the major challenges in modern quantum chemistry. Polynuclear transition metal clusters are peculiar examples of systems featuring such forms of electron correlation. Multireference strategies, often based on but not limited to the concept of complete active space, are adopted to accurately account for strong electron correlation and to resolve their complex electronic structures. However, transition metal clusters already containing four magnetic centers with multiple unpaired electrons make conventional active space based strategies prohibitively expensive, due to their unfavorable scaling with the size of the active space. In this work, forefront techniques, such as density matrix renormalization group (DMRG), full configuration interaction quantum Monte Carlo (FCIQMC), and multiconfiguration pair-density functional theory (MCPDFT), are employed to overcome the computational limitation of conventional multireference approaches and to accurately investigate the magnetic interactions taking place in a [Co(II)3Er(III)(OR)4] (chemical formula [Co(II)3Er(III)(hmp)4(μ2-OAc)2(OH)3(H2O)], hmp = 2-(hydroxymethyl)-pyridine) model cubane water oxidation catalyst. Complete active spaces with up to 56 electrons in 56 orbitals have been constructed for the seven energetically lowest different spin states. Relative energies, local spin, and spin-spin correlation values are reported and provide crucial insights on the spin interactions for this model system, pivotal in the rationalization of the catalytic activity of this system in the water-splitting reaction. A ferromagnetic ground state is found with a very small, ∼50 cm-1, highest-to-lowest spin gap. Moreover, for the energetically lowest states, S = 3-6, the three Co(II) sites exhibit parallel aligned spins, and for the lower states, S = 0-2, two Co(II) sites retain strong parallel spin alignment.
Collapse
Affiliation(s)
- Ruocheng Han
- Department of Chemistry A, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Sandra Luber
- Department of Chemistry A, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Giovanni Li Manni
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| |
Collapse
|
40
|
Manawadu D, Valentine DJ, Barford W. Dynamical Simulations of Carotenoid Photoexcited States Using Density Matrix Renormalization Group Techniques. J Phys Chem A 2023; 127:3714-3727. [PMID: 37054397 PMCID: PMC10150368 DOI: 10.1021/acs.jpca.3c00988] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
We present a dynamical simulation scheme to model the highly correlated excited state dynamics of linear polyenes. We apply it to investigate the internal conversion processes of carotenoids following their photoexcitation. We use the extended Hubbard-Peierls model, H^UVP, to describe the π-electronic system coupled to nuclear degrees of freedom. This is supplemented by a Hamiltonian, H^ϵ, that explicitly breaks both the particle-hole and two-fold rotation symmetries of idealized carotenoid structures. The electronic degrees of freedom are treated quantum mechanically by solving the time-dependent Schrödinger equation using the adaptive time-dependent DMRG (tDMRG) method, while nuclear dynamics are treated via the Ehrenfest equations of motion. By defining adiabatic excited states as the eigenstates of the full Hamiltonian, H^=H^UVP+H^ϵ, and diabatic excited states as eigenstates of H^UVP, we present a computational framework to monitor the internal conversion process from the initial photoexcited 11Bu+ state to the singlet triplet-pair states of carotenoids. We further incorporate Lanczos-DMRG to the tDMRG-Ehrenfest method to calculate transient absorption spectra from the evolving photoexcited state. We describe in detail the accuracy and convergence criteria for DMRG, and show that this method accurately describes the dynamical processes of carotenoid excited states. We also discuss the effect of the symmetry-breaking term, H^ϵ, on the internal conversion process, and show that its effect on the extent of internal conversion can be described by a Landau-Zener-type transition. This methodological paper is a companion to our more explanatory discussion of carotenoid excited state dynamics in Manawadu, D.; Georges, T. N.; Barford, W. Photoexcited State Dynamics and Singlet Fission in Carotenoids. J. Phys. Chem. A 2023, 127, 1342.
Collapse
Affiliation(s)
- Dilhan Manawadu
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
- Linacre College, University of Oxford, Oxford OX1 3JA, United Kingdom
| | - Darren J Valentine
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
- Balliol College, University of Oxford, Oxford OX1 3BJ, United Kingdom
| | - William Barford
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
41
|
Morgante P, Peverati R. Comparison of the Performance of Density Functional Methods for the Description of Spin States and Binding Energies of Porphyrins. Molecules 2023; 28:molecules28083487. [PMID: 37110720 PMCID: PMC10146789 DOI: 10.3390/molecules28083487] [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: 03/28/2023] [Revised: 04/10/2023] [Accepted: 04/13/2023] [Indexed: 04/29/2023] Open
Abstract
This work analyzes the performance of 250 electronic structure theory methods (including 240 density functional approximations) for the description of spin states and the binding properties of iron, manganese, and cobalt porphyrins. The assessment employs the Por21 database of high-level computational data (CASPT2 reference energies taken from the literature). Results show that current approximations fail to achieve the "chemical accuracy" target of 1.0 kcal/mol by a long margin. The best-performing methods achieve a mean unsigned error (MUE) <15.0 kcal/mol, but the errors are at least twice as large for most methods. Semilocal functionals and global hybrid functionals with a low percentage of exact exchange are found to be the least problematic for spin states and binding energies, in agreement with the general knowledge in transition metal computational chemistry. Approximations with high percentages of exact exchange (including range-separated and double-hybrid functionals) can lead to catastrophic failures. More modern approximations usually perform better than older functionals. An accurate statistical analysis of the results also casts doubts on some of the reference energies calculated using multireference methods. Suggestions and general guidelines for users are provided in the conclusions. These results hopefully stimulate advances for both the wave function and the density functional side of electronic structure calculations.
Collapse
Affiliation(s)
- Pierpaolo Morgante
- Department of Chemistry and Chemical Engineering, Florida Institute of Technology, 150 W. University Blvd., Melbourne, FL 32901, USA
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, NY 14260, USA
| | - Roberto Peverati
- Department of Chemistry and Chemical Engineering, Florida Institute of Technology, 150 W. University Blvd., Melbourne, FL 32901, USA
| |
Collapse
|
42
|
Fang Y, Zhang F, Ye Z, Zhang H, Lu WC, Wu S, Yao YX, Wang CZ, Ho KM. Ground and excited states of even-numbered Hubbard ring at half-filling: comparison of the extended Gutzwiller approach with exact diagonalization. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 35:265602. [PMID: 36972616 DOI: 10.1088/1361-648x/acc7ed] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 03/27/2023] [Indexed: 06/18/2023]
Abstract
It remains a great challenge in condensed matter physics to develop a method to treat strongly correlated many-body systems with balanced accuracy and efficiency. We introduce an extended Gutzwiller (EG) method incorporating a manifold technique, which builds an effective manifold of the many-body Hilbert space, to describe the ground-state (GS) and excited-state (ES) properties of strongly correlated electrons. We systematically apply an EG projector onto the GS and ES of a non-interacting system. Diagonalization of the true Hamiltonian within the manifold formed by the resulting EG wavefunctions gives the approximate GS and ES of the correlated system. To validate this technique, we implement it on even-numbered fermionic Hubbard rings at half-filling with periodic boundary conditions, and compare the results with the exact diagonalization (ED) method. The EG method is capable of generating high-quality GS and low-lying ES wavefunctions, as evidenced by the high overlaps of wavefunctions between the EG and ED methods. Favorable comparisons are also achieved for other quantities including the total energy, the double occupancy, the total spin and the staggered magnetization. With the capability of accessing the ESs, the EG method can capture the essential features of the one-electron removal spectral function that contains contributions from states deep in the excited spectrum. Finally, we provide an outlook on the application of this method on large extended systems.
Collapse
Affiliation(s)
- Yimei Fang
- Department of Physics, OSED, Key Laboratory of Low Dimensional Condensed Matter Physics (Department of Education of Fujian Province), Jiujiang Research institute, Xiamen University, Xiamen 361005, People's Republic of China
| | - Feng Zhang
- Division of Materials Science and Engineering, Ames National Laboratory, Ames, IA 50011, United States of America
- Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, United States of America
| | - Zhuo Ye
- Division of Materials Science and Engineering, Ames National Laboratory, Ames, IA 50011, United States of America
| | - Han Zhang
- College of Physics, Qingdao University, Qingdao, Shandong 266071, People's Republic of China
| | - Wen-Cai Lu
- College of Physics, Qingdao University, Qingdao, Shandong 266071, People's Republic of China
| | - Shunqing Wu
- Department of Physics, OSED, Key Laboratory of Low Dimensional Condensed Matter Physics (Department of Education of Fujian Province), Jiujiang Research institute, Xiamen University, Xiamen 361005, People's Republic of China
| | - Yong-Xin Yao
- Division of Materials Science and Engineering, Ames National Laboratory, Ames, IA 50011, United States of America
- Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, United States of America
| | - Cai-Zhuang Wang
- Division of Materials Science and Engineering, Ames National Laboratory, Ames, IA 50011, United States of America
- Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, United States of America
| | - Kai-Ming Ho
- Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, United States of America
| |
Collapse
|
43
|
Lee S, Lee J, Zhai H, Tong Y, Dalzell AM, Kumar A, Helms P, Gray J, Cui ZH, Liu W, Kastoryano M, Babbush R, Preskill J, Reichman DR, Campbell ET, Valeev EF, Lin L, Chan GKL. Evaluating the evidence for exponential quantum advantage in ground-state quantum chemistry. Nat Commun 2023; 14:1952. [PMID: 37029105 PMCID: PMC10082187 DOI: 10.1038/s41467-023-37587-6] [Citation(s) in RCA: 25] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Accepted: 03/22/2023] [Indexed: 04/09/2023] Open
Abstract
Due to intense interest in the potential applications of quantum computing, it is critical to understand the basis for potential exponential quantum advantage in quantum chemistry. Here we gather the evidence for this case in the most common task in quantum chemistry, namely, ground-state energy estimation, for generic chemical problems where heuristic quantum state preparation might be assumed to be efficient. The availability of exponential quantum advantage then centers on whether features of the physical problem that enable efficient heuristic quantum state preparation also enable efficient solution by classical heuristics. Through numerical studies of quantum state preparation and empirical complexity analysis (including the error scaling) of classical heuristics, in both ab initio and model Hamiltonian settings, we conclude that evidence for such an exponential advantage across chemical space has yet to be found. While quantum computers may still prove useful for ground-state quantum chemistry through polynomial speedups, it may be prudent to assume exponential speedups are not generically available for this problem.
Collapse
Affiliation(s)
- Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, NY, 10027, USA
| | - Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Yu Tong
- Department of Mathematics, University of California, Berkeley, CA, 94720, USA
| | | | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, VA, 24061, USA
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Phillip Helms
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Johnnie Gray
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Zhi-Hao Cui
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Wenyuan Liu
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Michael Kastoryano
- AWS Center for Quantum Computing, Pasadena, CA, 91125, USA
- Amazon Quantum Solutions Lab, Seattle, WA, 98170, USA
| | - Ryan Babbush
- Google Quantum AI, 340 Main Street, Venice, CA, 90291, USA
| | - John Preskill
- AWS Center for Quantum Computing, Pasadena, CA, 91125, USA
- Institute for Quantum Information and Matter, California Institute of Technology, Pasadena, CA, 91125, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, NY, 10027, USA
| | | | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, VA, 24061, USA
| | - Lin Lin
- Department of Mathematics, University of California, Berkeley, CA, 94720, USA
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
| |
Collapse
|
44
|
Yehorova D, Kretchmer JS. A multi-fragment real-time extension of projected density matrix embedding theory: Non-equilibrium electron dynamics in extended systems. J Chem Phys 2023; 158:131102. [PMID: 37031109 DOI: 10.1063/5.0146973] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/07/2023] Open
Abstract
In this work, we derive a multi-fragment real-time extension of the projected density matrix embedding theory (pDMET) designed to treat non-equilibrium electron dynamics in strongly correlated systems. As in the previously developed static pDMET, the real time pDMET partitions the total system into many fragments; the coupling between each fragment and the rest of the system is treated through a compact representation of the environment in terms of a quantum bath. The real-time pDMET involves simultaneously propagating the wavefunctions for each separate fragment–bath embedding system along with an auxiliary mean-field wavefunction of the total system. The equations of motion are derived by (i) projecting the time-dependent Schrödinger equation in the fragment and bath space associated with each separate fragment and by (ii) enforcing the pDMET matching conditions between the global 1-particle reduced density matrix (1-RDM) obtained from the fragment calculations and the mean-field 1-RDM at all points in time. The accuracy of the method is benchmarked through comparisons to time-dependent density-matrix renormalization group and time-dependent Hartree–Fock (TDHF) theory; the methods were applied to a one- and two-dimensional single-impurity Anderson model and multi-impurity Anderson models with ordered and disordered distributions of the impurities. The results demonstrate a large improvement over TDHF and rapid convergence to the exact dynamics with an increase in fragment size. Our results demonstrate that the real-time pDMET is a promising and flexible method that balances accuracy and efficiency to simulate the non-equilibrium electron dynamics in heterogeneous systems of large size.
Collapse
Affiliation(s)
- Dariia Yehorova
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Joshua S. Kretchmer
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
45
|
Segatta F, Ruiz DA, Aleotti F, Yaghoubi M, Mukamel S, Garavelli M, Santoro F, Nenov A. Nonlinear Molecular Electronic Spectroscopy via MCTDH Quantum Dynamics: From Exact to Approximate Expressions. J Chem Theory Comput 2023; 19:2075-2091. [PMID: 36961952 PMCID: PMC10100531 DOI: 10.1021/acs.jctc.2c01059] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/26/2023]
Abstract
We present an accurate and efficient approach to computing the linear and nonlinear optical spectroscopy of a closed quantum system subject to impulsive interactions with an incident electromagnetic field. It incorporates the effect of ultrafast nonadiabatic dynamics by means of explicit numerical propagation of the nuclear wave packet. The fundamental expressions for the evaluation of first- and higher-order response functions are recast in a general form that can be used with any quantum dynamics code capable of computing the overlap of nuclear wave packets evolving in different states. Here we present the evaluation of these expressions with the multiconfiguration time-dependent Hartree (MCTDH) method. Application is made to pyrene, excited to its lowest bright excited state S2 which exhibits a sub-100-fs nonadiabatic decay to a dark state S1. The system is described by a linear vibronic coupling Hamiltonian, parametrized with multiconfiguration electronic structure methods. We show that the ultrafast nonadiabatic dynamics can have a remarkable effect on the spectral line shapes that goes beyond simple lifetime broadening. Furthermore, a widely employed approximate expression based on the time scale separation of dephasing and population relaxation is recast in the same theoretical framework. Application to pyrene shows the range of validity of such approximations.
Collapse
Affiliation(s)
- Francesco Segatta
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Daniel Aranda Ruiz
- ICMol, Universidad de Valencia, Catedrático José Beltrán Martínez, 2, 46980 Paterna, Spain
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Flavia Aleotti
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Martha Yaghoubi
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Shaul Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, United States
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| |
Collapse
|
46
|
Liao K, Zhai H, Christlmaier EM, Schraivogel T, Ríos PL, Kats D, Alavi A. Density Matrix Renormalization Group for Transcorrelated Hamiltonians: Ground and Excited States in Molecules. J Chem Theory Comput 2023; 19:1734-1743. [PMID: 36912635 DOI: 10.1021/acs.jctc.2c01207] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/14/2023]
Abstract
We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians and show applications in molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem, but this is currently limited to treating the ground state and has lower efficiency than the time-independent DMRG (TI-DMRG) due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely, replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. An accelerated basis set convergence rate is also observed, as expected, within the TC framework.
Collapse
Affiliation(s)
- Ke Liao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | | | - Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Pablo López Ríos
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany.,Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
47
|
Nowak A, Boguslawski K. A configuration interaction correction on top of pair coupled cluster doubles. Phys Chem Chem Phys 2023; 25:7289-7301. [PMID: 36810525 DOI: 10.1039/d2cp05171k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Numerous numerical studies have shown that geminal-based methods are a promising direction to model strongly correlated systems with low computational costs. Several strategies have been introduced to capture the missing dynamical correlation effects, which typically exploit a posteriori corrections to account for correlation effects associated with broken-pair states or inter-geminal correlations. In this article, we scrutinize the accuracy of the pair coupled cluster doubles (pCCD) method extended by configuration interaction (CI) theory. Specifically, we benchmark various CI models, including, at most double excitations against selected CC corrections as well as conventional single-reference CC methods. A simple Davidson correction is also tested. The accuracy of the proposed pCCD-CI approaches is assessed for challenging small model systems such as the N2 and F2 dimers and various di- and triatomic actinide-containing compounds. In general, the proposed CI methods considerably improve spectroscopic constants compared to the conventional CCSD approach, provided a Davidson correction is included in the theoretical model. At the same time, their accuracy lies between those of the linearized frozen pCCD and frozen pCCD variants.
Collapse
Affiliation(s)
- Artur Nowak
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| |
Collapse
|
48
|
Bensberg M, Reiher M. Corresponding Active Orbital Spaces along Chemical Reaction Paths. J Phys Chem Lett 2023; 14:2112-2118. [PMID: 36802629 PMCID: PMC9986954 DOI: 10.1021/acs.jpclett.2c03905] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2022] [Accepted: 02/10/2023] [Indexed: 06/18/2023]
Abstract
The accuracy of reaction energy profiles calculated with multiconfigurational electronic structure methods and corrected by multireference perturbation theory depends crucially on consistent active orbital spaces selected along the reaction path. However, it has been challenging to choose molecular orbitals that can be considered corresponding in different molecular structures. Here, we demonstrate how active orbital spaces can be selected consistently along reaction coordinates in a fully automatized way. The approach requires no structure interpolation between reactants and products. Instead, it emerges from a synergy of the Direct Orbital Selection orbital mapping ansatz combined with our fully automated active space selection algorithm autoCAS. We demonstrate our algorithm for the potential energy profile of the homolytic carbon-carbon bond dissociation and rotation around the double bond of 1-pentene in the electronic ground state. However, our algorithm also applies to electronically excited Born-Oppenheimer surfaces.
Collapse
|
49
|
Unsleber JP, Liu H, Talirz L, Weymuth T, Mörchen M, Grofe A, Wecker D, Stein CJ, Panyala A, Peng B, Kowalski K, Troyer M, Reiher M. High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation. J Chem Phys 2023; 158:084803. [PMID: 36859110 DOI: 10.1063/5.0136526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.
Collapse
Affiliation(s)
- Jan P Unsleber
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Hongbin Liu
- Microsoft Quantum, Redmond, Washington 98052, USA
| | | | - Thomas Weymuth
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Maximilian Mörchen
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Adam Grofe
- Microsoft Quantum, Redmond, Washington 98052, USA
| | - Dave Wecker
- Microsoft Quantum, Redmond, Washington 98052, USA
| | - Christopher J Stein
- Department of Chemistry, TUM School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, D-85748 Garching, Germany
| | - Ajay Panyala
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Bo Peng
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Karol Kowalski
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | | | - Markus Reiher
- Laboratory of Physical Chemistry and NCCR Catalysis, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
50
|
Xie Z, Song Y, Peng F, Li J, Cheng Y, Zhang L, Ma Y, Tian Y, Luo Z, Ma H. Kylin 1.0: An ab-initio density matrix renormalization group quantum chemistry program. J Comput Chem 2023; 44:1316-1328. [PMID: 36809661 DOI: 10.1002/jcc.27085] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Revised: 01/15/2023] [Accepted: 01/19/2023] [Indexed: 02/23/2023]
Abstract
The accurate evaluation of electron correlations is highly necessary for the proper descriptions of the electronic structures in strongly correlated molecules, ranging from bond-dissociating molecules, polyradicals, to large conjugated molecules and transition metal complexes. For this purpose, in this paper, a new ab-initio quantum chemistry program Kylin 1.0 for electron correlation calculations at various quantum many-body levels, including configuration interaction (CI), perturbation theory (PT), and density matrix renormalization group (DMRG), is presented. Furthermore, fundamental quantum chemistry methods such as Hartree-Fock self-consistent field (HF-SCF) and the complete active space SCF (CASSCF) are also implemented. The Kylin 1.0 program possesses the following features: (1) a matrix product operator (MPO) formulation-based efficient DMRG implementation for describing static electron correlation within a large active space composed of more than 100 orbitals, supporting both U 1 n × U 1 S z $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{U}{(1)}_{{\mathrm{S}}_{\mathrm{z}}} $$ and U 1 n × SU 2 S $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{SU}{(2)}_{\mathrm{S}} $$ symmetries; (2) an efficient second-order DMRG-self-consistent field (SCF) implementation; (3) an externally contracted multi-reference CI (MRCI) and Epstein-Nesbet PT with DMRG reference wave functions for including the remaining dynamic electron correlation outside the large active spaces. In this paper, we introduce the capabilities and numerical benchmark examples of the Kylin 1.0 program.
Collapse
Affiliation(s)
- Zhaoxuan Xie
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Yinxuan Song
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Fangwen Peng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Jianhao Li
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Yifan Cheng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Lingzhi Zhang
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Yingqi Tian
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Zhen Luo
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Haibo Ma
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China.,Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, China
| |
Collapse
|