1
|
de Moura CEV, Sokolov AY. Efficient Spin-Adapted Implementation of Multireference Algebraic Diagrammatic Construction Theory. I. Core-Ionized States and X-ray Photoelectron Spectra. J Phys Chem A 2024; 128:5816-5831. [PMID: 38962857 DOI: 10.1021/acs.jpca.4c03161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2024]
Abstract
We present an efficient implementation of multireference algebraic diagrammatic construction theory (MR-ADC) for simulating core-ionized states and X-ray photoelectron spectra (XPS). Taking advantage of spin adaptation, automatic code generation, and density fitting, our implementation can perform calculations for molecules with more than 1500 molecular orbitals, incorporating static and dynamic correlation in the ground and excited electronic states. We demonstrate the capabilities of MR-ADC methods by simulating the XPS spectra of substituted ferrocene complexes and azobenzene isomers. For the ground electronic states of these molecules, the XPS spectra computed using the extended second-order MR-ADC method (MR-ADC(2)-X) are in a very good agreement with available experimental results. We further show that MR-ADC can be used as a tool for interpreting or predicting the results of time-resolved XPS measurements by simulating the core ionization spectra of azobenzene along its photoisomerization, including the XPS signatures of excited states and the minimum energy conical intersection. This work is the first in a series of publications reporting the efficient implementations of MR-ADC methods.
Collapse
Affiliation(s)
- Carlos E V de Moura
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
2
|
Giussani A, Segarra-Martí J. XMS-CASPT2//XMS-CASPT2 and XMS-CASPT2//CASSCF at comparison: The impact of dynamic correlation in the excited state optimization of nitronaphthalene. J Chem Phys 2024; 161:011103. [PMID: 38953440 DOI: 10.1063/5.0203550] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2024] [Accepted: 06/12/2024] [Indexed: 07/04/2024] Open
Abstract
Here, analytical extended multi-state complete active space second-order perturbation method (XMS-CASPT2) gradients are used to rationalize the decreasing triplet quantum yield trend in 2-nitronaphthalene, 1-nitronaphthalene, and 2-methyl-1-nitronaphthalene, a series of nitro-substituted aromatic compounds. Comparison with the XMS-CASPT2//CASSCF (where CASSCF stands for complete active space self-consistent field method) results highlights the importance of dynamic correlation in geometry optimization and challenges the validity of an XMS-CASPT2//CASSCF approach: XMS-CASPT2 S1 minima leads to planar structures, while CASSCF optimizations trigger a pyramidalization of the nitro group. The XMS-CASPT2 results correlate the reported decreasing triplet quantum yield trend in these species to a decrease in S1 to T2 population transfer and an increase in S1-S0 decay, while no such correlation is observed when using XMS-CASPT2//CASSCF data.
Collapse
Affiliation(s)
- Angelo Giussani
- Instituto de Ciencia Molecular, Universitat de València, Apartado22085, ES-46071 Valencia, Spain
| | - Javier Segarra-Martí
- Instituto de Ciencia Molecular, Universitat de València, Apartado22085, ES-46071 Valencia, Spain
| |
Collapse
|
3
|
Lechner MH, Papadopoulos A, Sivalingam K, Auer AA, Koslowski A, Becker U, Wennmohs F, Neese F. Code generation in ORCA: progress, efficiency and tight integration. Phys Chem Chem Phys 2024; 26:15205-15220. [PMID: 38767596 DOI: 10.1039/d4cp00444b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
An improved version of ORCA's automated generator environment (ORCA-AGE II) is presented. The algorithmic improvements and the move to C++ as the programming language lead to a performance gain of up to two orders of magnitude compared to the previously developed PYTHON toolchain. Additionally, the restructured modular design allows for far more complex code engines to be implemented readily. Importantly, we have realised an extremely tight integration with the ORCA host program. This allows for a workflow in which only the wavefunction Ansatz is part of the source code repository while all actual high-level code is generated automatically, inserted at the appropriate place in the host program before it is compiled and linked together with the hand written code parts. This construction ensures longevity and uniform code quality. Furthermore the new developments allow ORCA-AGE II to generate parallelised production-level code for highly complex theories, such as fully internally contracted multireference coupled-cluster theory (fic-MRCC) with an enormous number of contributing tensor contractions. We also discuss the automated implementation of nuclear gradients for arbitrary theories. All these improvements enable the implementation of theories that are too complex for the human mind and also reduce development times by orders of magnitude. We hope that this work enables researchers to concentrate on the intellectual content of the theories they develop rather than be concerned with technical details of the implementation.
Collapse
Affiliation(s)
- Marvin H Lechner
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Anastasios Papadopoulos
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Kantharuban Sivalingam
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Alexander A Auer
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Axel Koslowski
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Ute Becker
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Frank Wennmohs
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Frank Neese
- Department of Molecular Theory and Spectroscopy, Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| |
Collapse
|
4
|
Karak P, Moitra T, Banerjee A, Ruud K, Chakrabarti S. Accidental triplet harvesting in donor-acceptor dyads with low spin-orbit coupling. Phys Chem Chem Phys 2024; 26:5344-5355. [PMID: 38268441 DOI: 10.1039/d3cp04904c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2024]
Abstract
We present an accidental mechanism for efficient intersystem crossing (ISC) between singlet and triplet states with low spin-orbit coupling (SOC) in molecules having donor-acceptor (D-A) moieties separated by a Sigma bond. Our study shows that SOC between the lowest singlet excited state and the higher-lying triplet states, together with nuclear motion-driven coupling of this triplet state with lower-lying triplet states during the free rotation about a Sigma bond, is one of the possible ways to achieve the experimentally observed ISC rate for a class of D-A type photoredox catalysts. This mechanism is found to be the dominant contributor to the ISC process with the corresponding rate reaching a maximum at a dihedral angle in the range of 72°-78° between the D-A moieties of 10-(naphthalen-1-yl)-3,7-diphenyl-10H-phenoxazine and other molecules included in the study. We have further demonstrated that the same mechanism is operative in a specific spirobis[anthracene]dione molecule, where the D and A moieties are interlocked near to the optimal dihedral angle, indicating the plausible effectiveness of the proposed mechanism. The present finding is expected to have implications in strategies for the synthesis of new generations of triplet-harvesting organic molecules.
Collapse
Affiliation(s)
- Pijush Karak
- Department of Chemistry, University of Calcutta, 92 A.P.C Road, Kolkata - 700009, West Bengal, India.
| | - Torsha Moitra
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, 9037 Tromsø, Norway.
| | - Ambar Banerjee
- Department of Physics and Astronomy, X-ray Photon Science, Uppsala University, Ångströmlaboratoriet, Lägerhyddsvägen 1, 75120, Uppsala, Sweden.
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, 9037 Tromsø, Norway.
- Norwegian Defence Research Establishment, P.O.Box 25, 2027 Kjeller, Norway
| | - Swapan Chakrabarti
- Department of Chemistry, University of Calcutta, 92 A.P.C Road, Kolkata - 700009, West Bengal, India.
| |
Collapse
|
5
|
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: 9] [Impact Index Per Article: 9.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
|
6
|
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: 56] [Impact Index Per Article: 56.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
|
7
|
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
|
8
|
Bhowmick R, Roy Chowdhury S, Vlaisavljevich B. Molecular Geometry and Electronic Structure of Copper Corroles. Inorg Chem 2023; 62:13877-13891. [PMID: 37590888 DOI: 10.1021/acs.inorgchem.3c01779] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/19/2023]
Abstract
Copper corroles are known for their unique multiconfigurational electronic structures in the ground state, which arise from the transfer of electrons from the π orbitals of the corrole to the d-orbital of copper. While density functional theory (DFT) provides reasonably good molecular geometries, the determination of the ground spin state and the associated energetics is heavily influenced by functional choice, particularly the percentage of the Hartree-Fock exchange. Using extended multireference perturbation theory methods (XMS-CASPT2), the functional choice can be assessed. The molecular geometries and electronic structures of both the unsubstituted and the meso-triphenyl copper corroles were investigated. A minimal active space was employed for structural characterization, while larger active spaces are required to examine the electronic structure. The XMS-CASPT2 investigations conclusively identify the ground electronic state as a multiconfigurational singlet (S0) with three dominant electronic configurations in its lowest energy and characteristic saddled structure. In contrast, the planar geometry corresponds to the triplet state (T0), which is approximately 5 kcal/mol higher in energy compared to the S0 state for both the bare and substituted copper corroles. Notably, the planarity of the T0 geometry is reduced in the substituted corrole compared with that in the unsubstituted one. By analyzing the potential energy surface (PES) between the S0 and T0 geometries using XMS-CASPT2, the multiconfigurational electronic structure is shown to transition toward a single electron configuration as the saddling angle decreases (i.e., as one approaches the planar geometry). Despite the ability of the functionals to reproduce the minimum energy structures, only the TPSSh-D3 PES is reasonably close to the XMS-CASPT2 surface. Significant deviations along the PES are observed with other functionals.
Collapse
Affiliation(s)
- Rina Bhowmick
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Sabyasachi Roy Chowdhury
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| |
Collapse
|
9
|
Nishimoto Y. Analytic first-order derivatives of CASPT2 with IPEA shift. J Chem Phys 2023; 158:2888841. [PMID: 37144712 DOI: 10.1063/5.0147611] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2023] [Accepted: 04/18/2023] [Indexed: 05/06/2023] Open
Abstract
Complete active space second-order perturbation theory (CASPT2) is useful for accurately predicting properties of complex electronic structures, but it is well known that it systematically underestimates excitation energies. The underestimation can be corrected using the ionization potential-electron affinity (IPEA) shift. In this study, analytic first-order derivatives of CASPT2 with the IPEA shift are developed. CASPT2-IPEA is not invariant with respect to rotations among active molecular orbitals, and two additional constraint conditions are necessary in the CASPT2 Lagrangian to formulate analytic derivatives. The method developed here is applied to methylpyrimidine derivatives and cytosine, and minimum energy structures and conical intersections are located. By comparing energies relative to the closed-shell ground state, we find that the agreement with experiments and high-level calculations is indeed improved by the inclusion of the IPEA shift. The agreement of geometrical parameters with high-level calculations may also be improved in some cases.
Collapse
Affiliation(s)
- Yoshio Nishimoto
- Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
10
|
Karak P, Moitra T, Ruud K, Chakrabarti S. Photophysics of uracil: an explicit time-dependent generating function-based method combining both nonadiabatic and spin-orbit coupling effects. Phys Chem Chem Phys 2023; 25:8209-8219. [PMID: 36881024 DOI: 10.1039/d2cp05955j] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
Abstract
We present a composite framework for calculating the rates of non-radiative deactivation processes, namely internal conversion (IC) and intersystem crossing (ISC), on an equal footing by explicitly computing the non-adiabatic coupling (NAC) and spin-orbit coupling (SOC) constants, respectively. The stationary-state approach uses a time-dependent generating function based on Fermi's golden rule. We validate the applicability of the framework by computing the rate of IC for azulene, obtaining comparable rates to experimental and previous theoretical results. Next, we investigate the photophysics associated with the complex photodynamics of the uracil molecule. Interestingly, our simulated rates corroborate experimental observations. Detailed analyses using Duschinsky rotation matrices, displacement vectors and NAC matrix elements are presented to interpret the findings alongside testing the suitability of the approach for such molecular systems. The suitability of the Fermi's golden rule based method is explained qualitatively in terms of single-mode potential energy surfaces.
Collapse
Affiliation(s)
- Pijush Karak
- Department of Chemistry, University of Calcutta, 92 A.P.C Road, Kolkata-700009, West Bengal, India.
| | - Torsha Moitra
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø-The Arctic University of Norway, 9037 Tromsø, Norway.
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø-The Arctic University of Norway, 9037 Tromsø, Norway. .,Norwegian Defence Research Establishment, P.O.Box 25, 2027 Kjeller, Norway
| | - Swapan Chakrabarti
- Department of Chemistry, University of Calcutta, 92 A.P.C Road, Kolkata-700009, West Bengal, India.
| |
Collapse
|
11
|
Iino T, Shiozaki T, Yanai T. Algorithm for analytic nuclear energy gradients of state averaged DMRG-CASSCF theory with newly derived coupled-perturbed equations. J Chem Phys 2023; 158:054107. [PMID: 36754810 DOI: 10.1063/5.0130636] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
We present an algorithm for evaluating analytic nuclear energy gradients of the state-averaged density matrix renormalization group complete-active-space self-consistent field (SA-DMRG-CASSCF) theory based on the newly derived coupled-perturbed (CP) DMRG-CASSCF equations. The Lagrangian for the conventional SA-CASSCF analytic gradient theory is extended to the SA-DMRG-CASSCF variant that can fully consider a whole set of constraints on the parameters of multi-root canonical matrix product states formed at all the DMRG block configurations. An efficient algorithm to solve the CP-DMRG-CASSCF equations for determining the multipliers was developed. The complexity of the resultant analytic gradient algorithm is overall the same as that of the unperturbed SA-DMRG-CASSCF algorithm. In addition, a reduced-scaling approach was developed to directly compute the SA reduced density matrices (SA-RDMs) and their perturbed ones without calculating separate state-specific RDMs. As part of our implementation scheme, we neglect the term associated with the constraint on the active orbitals in terms of the active-active rotation in the Lagrangian. Thus, errors from the true analytic gradients may be caused in this scheme. The proposed gradient algorithm was tested with the spin-adapted implementation by checking how accurately the computed analytic energy gradients reproduce numerical gradients of the SA-DMRG-CASSCF energies using a common number of renormalized bases. The illustrative applications show that the errors are sufficiently small when using a typical number of the renormalized bases, which is required to attain adequate accuracy in DMRG's total energies.
Collapse
Affiliation(s)
- Tsubasa Iino
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| | - Toru Shiozaki
- Quantum Simulation Technologies, Inc., Boston, Massachusetts 02135, USA
| | - Takeshi Yanai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| |
Collapse
|
12
|
Abstract
We develop analytic gradients for selected configuration interaction wave functions. Despite all pairs of molecular orbitals now potentially having to be considered for the coupled perturbed Hartree-Fock equations, we show that degenerate orbital pairs belonging to different irreducible representations in the largest abelian subgroup do not need to be included and instabilities due to degeneracies are avoided. We introduce seminumerical gradients and use them to validate the analytic approach even when near degeneracies are present due to high-symmetry geometries being slightly distorted to break symmetry. The method is applied to carbon monoxide, ammonia, square planar H4, hexagonal planar H6, and methane for a range of bond lengths where we demonstrate that analytic gradients for selected configuration interaction can approach the quality of full configuration interaction yet only use a very small fraction of its determinants.
Collapse
|
13
|
Evangelista FA. Automatic derivation of many-body theories based on general Fermi vacua. J Chem Phys 2022; 157:064111. [PMID: 35963725 DOI: 10.1063/5.0097858] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper describes Wick&d, an implementation of the algebra of second-quantized operators normal ordered with respect to general correlated references and the corresponding Wick theorem [D. Mukherjee, Chem. Phys. Lett. 274, 561 (1997) and W. Kutzelnigg and D. Mukherjee, J. Chem. Phys. 107, 432 (1997)]. Wick&d employs a compact representation of operators and a backtracking algorithm to efficiently evaluate Wick contractions. Since Wick&d can handle both fully and partially contracted terms, it can be applied to both projective and Fock-space many-body formalisms. To demonstrate the usefulness of Wick&d, we use it to evaluate the single-reference coupled cluster equations up to octuple excitations and report an automated derivation and implementation of the second-order driven similarity renormalization group multireference perturbation theory.
Collapse
Affiliation(s)
- Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
14
|
Battaglia S, Fransén L, Fdez Galván I, Lindh R. Regularized CASPT2: an Intruder-State-Free Approach. J Chem Theory Comput 2022; 18:4814-4825. [PMID: 35876618 PMCID: PMC9367007 DOI: 10.1021/acs.jctc.2c00368] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
In this work we present a new approach to fix the intruder-state
problem (ISP) in CASPT2 based on σp regularization. The resulting σp-CASPT2 method is compared to previous techniques, namely, the real
and imaginary level shifts, on a theoretical basis and by performing
a series of systematic calculations. The analysis is focused on two
aspects, the effectiveness of σp-CASPT2 in removing the ISP and the sensitivity of the approach with
respect to the input parameter. We found that σp-CASPT2 compares favorably with respect to previous
approaches and that different versions, σ1-CASPT2
and σ2-CASPT2, have different potential application
domains. This analysis also reveals the unsuitability of the real
level shift technique as a general way to avoid the intruder-state
problem.
Collapse
Affiliation(s)
- Stefano Battaglia
- Department of Chemistry─BMC, Uppsala University, P. O. Box 576, SE-75123 Uppsala, Sweden
| | - Lina Fransén
- Department of Chemistry─BMC, Uppsala University, P. O. Box 576, SE-75123 Uppsala, Sweden
| | - Ignacio Fdez Galván
- Department of Chemistry─BMC, Uppsala University, P. O. Box 576, SE-75123 Uppsala, Sweden
| | - Roland Lindh
- Department of Chemistry─BMC, Uppsala University, P. O. Box 576, SE-75123 Uppsala, Sweden
| |
Collapse
|
15
|
Nishimoto Y, Battaglia S, Lindh R. Analytic First-Order Derivatives of (X)MS, XDW, and RMS Variants of the CASPT2 and RASPT2 Methods. J Chem Theory Comput 2022; 18:4269-4281. [PMID: 35699280 DOI: 10.1021/acs.jctc.2c00301] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Crossings between states involve complex electronic structures, making the accurate characterization of the crossing point difficult. In this study, the analytic derivatives of three complete active space second-order perturbation theory (CASPT2) variants as well as an extension of the restricted active space (RASPT2) are developed. These variants are applied to locating minimum energy conical intersections. Our results demonstrate that the three CASPT2 variants predict qualitatively similar results, but a recently developed variant, the rotated multistate CASPT2 (RMS-CASPT2), is least sensitive to the number of states considered in the calculation. We demonstrate that CASPT2 and the reference self-consistent field calculations predict qualitatively different energetics and bond lengths.
Collapse
Affiliation(s)
- Yoshio Nishimoto
- Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Stefano Battaglia
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Roland Lindh
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| |
Collapse
|
16
|
Affiliation(s)
- Milica Feldt
- Leibniz Institute for Catalysis: Leibniz-Institut fur Katalyse eV Theory & Catalysis Albert-Einstein-Str 29A 18059 Rostock GERMANY
| | - Quan Manh Phung
- Nagoya University: Nagoya Daigaku Department of Chemistry JAPAN
| |
Collapse
|
17
|
Park JW. Analytical Gradient Theory for Spin-Free State-Averaged Second-Order Driven Similarity Renormalization Group Perturbation Theory (SA-DSRG-MRPT2) and Its Applications for Conical Intersection Optimizations. J Chem Theory Comput 2022; 18:2233-2245. [PMID: 35229599 DOI: 10.1021/acs.jctc.1c01150] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Second-order multireference-driven similarity renormalization group perturbation theory (DSRG-MRPT2) provides an efficient means of correcting the dynamical correlation with the multiconfiguration reference function. The state-averaged DSRG-MRPT2 (SA-DSRG-MRPT2) method is the simplest means of treating the excited states with DSRG-MRPT2. In this method, the Hamiltonian dressed with dynamical correlation is diagonalized in the CASCI state subspace (SA-DSRG-MRPT2c) or the configuration subspace (SA-DSRG-MRPT2). This work develops analytical gradient theory for spin-free SA-DSRG-MRPT2(c) with the density-fitting approximation. We check the accuracy of the analytical gradients against the numerical gradients. We present applications for optimizing minimum energy conical intersections (MECI) of ethylene and retinal model chromophores (PSB3 and RPSB6). We investigate the dependence of the optimized geometries and energies on the flow parameters and reference relaxations. The smoothness of the SA-DSRG-MRPT2(c) potential energy surfaces near the reference (complete active space self-consistent field) MECI is comparable to the XMCQDPT2 one. These results render SA-DSRG-MRPT2(c) theory a promising approach for studies of conical intersections.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
18
|
Vitillo JG, Cramer CJ, Gagliardi L. Multireference Methods are Realistic and Useful Tools for Modeling Catalysis. Isr J Chem 2022. [DOI: 10.1002/ijch.202100136] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Jenny G. Vitillo
- Department of Science and High Technology and INSTM Università degli Studi dell'Insubria Via Valleggio 9 I-22100 Como Italy
| | - Christopher J. Cramer
- Underwriters Laboratories Inc. 333 Pfingsten Road Northbrook Illinois 60602 United States
| | - Laura Gagliardi
- Department of Chemistry Pritzker School of Molecular Engineering James Franck Institute University of Chicago Chicago Illinois 60637 United States
| |
Collapse
|
19
|
Finney BA, Chowdhury SR, Kirkvold C, Vlaisavljevich B. CASPT2 molecular geometries of Fe(II) spin-crossover complexes. Phys Chem Chem Phys 2022; 24:1390-1398. [PMID: 34981806 DOI: 10.1039/d1cp04885f] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. We recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe-ligand bond distances.
Collapse
Affiliation(s)
- Brian A Finney
- University of South Dakota, Department of Chemistry, 414 E Clark St., Vermillion SD, 57069, USA.
| | - Sabyasachi Roy Chowdhury
- University of South Dakota, Department of Chemistry, 414 E Clark St., Vermillion SD, 57069, USA.
| | - Clara Kirkvold
- University of South Dakota, Department of Chemistry, 414 E Clark St., Vermillion SD, 57069, USA.
| | - Bess Vlaisavljevich
- University of South Dakota, Department of Chemistry, 414 E Clark St., Vermillion SD, 57069, USA.
| |
Collapse
|
20
|
Alías-Rodríguez M, Huix-Rotllant M, de Graaf C. Quantum dynamics simulations of the thermal and light-induced high-spin to low-spin relaxation in Fe(bpy)3 and Fe(mtz)6. Faraday Discuss 2022; 237:93-107. [DOI: 10.1039/d2fd00027j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
21
|
Hocky GM, White AD. Natural language processing models that automate programming will transform chemistry research and teaching. DIGITAL DISCOVERY 2022; 1:79-83. [PMID: 35515080 PMCID: PMC8996826 DOI: 10.1039/d1dd00009h] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Accepted: 01/25/2022] [Indexed: 11/21/2022]
Abstract
Natural language processing models have emerged that can generate useable software and automate a number of programming tasks with high fidelity. These tools have yet to have an impact on the chemistry community. Yet, our initial testing demonstrates that this form of artificial intelligence is poised to transform chemistry and chemical engineering research. Here, we review developments that brought us to this point, examine applications in chemistry, and give our perspective on how this may fundamentally alter research and teaching. Natural language processing models have emerged that can generate useable software and automate a number of programming tasks with high fidelity.![]()
Collapse
Affiliation(s)
| | - Andrew D. White
- Department of Chemical Engineering, University of Rochester, USA
| |
Collapse
|
22
|
Shiozaki T, Vlaisavljevich B. Computational Spectroscopy of the Cr-Cr Bond in Coordination Complexes. Inorg Chem 2021; 60:19219-19225. [PMID: 34883014 DOI: 10.1021/acs.inorgchem.1c03005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report the accurate computational vibrational analysis of the Cr-Cr bond in dichromium complexes using second-order multireference complete active space methods (CASPT2), allowing direct comparison with experimental spectroscopic data both to facilitate interpreting the low-energy region of the spectra and to provide insights into the nature of the bonds themselves. Recent technological development by the authors has realized such computation for the first time. Accurate simulation of the vibrational structure of these compounds has been hampered by their notorious multiconfigurational electronic structure that yields bond distances that do not correlate with bond order. Some measured Cr-Cr vibrational stretching modes, ν(Cr2), have suggested weaker bonding, even for so-called ultrashort Cr-Cr bonds, while others are in line with the bond distance. Here, we optimize geometries and compute ν(Cr2) with CASPT2 for three well-characterized complexes, Cr2(O2CCH3)4(H2O)2, Cr2(mhp)4, and Cr2(dmp)4. We obtain CASPT2 harmonic ν(Cr2) modes in good agreement with experiment at 282 cm-1 for Cr2(mhp)4 and 353 cm-1 for Cr2(dmp)4, compute 50Cr and 54Cr isotope shifts, and demonstrate that the use of the so-called IPEA shift leads to improved Cr-Cr distances. Additionally, normal mode sampling was used to estimate anharmonicity along ν(Cr2), leading to an anharmonic mode of 272 cm-1 for Cr2(mhp)4 and 333 cm-1 for Cr2(dmp)4.
Collapse
Affiliation(s)
- Toru Shiozaki
- Quantum Simulation Technologies, Inc., Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| |
Collapse
|
23
|
Jara-Cortés J, Leal-Sánchez E, Francisco E, Pérez-Pimienta JA, Martín Pendás Á, Hernández-Trujillo J. Implementation of the interacting quantum atom energy decomposition using the CASPT2 method. Phys Chem Chem Phys 2021; 23:27508-27519. [PMID: 34874377 DOI: 10.1039/d1cp02837e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We present an implementation of the interacting quantum atom (IQA) energy decomposition scheme using the complete active space second-order perturbation theory (CASPT2). This combination yields a real-space interpretation tool with a proper account of the static and dynamic correlation that is particularly relevant for the description of processes in electronic excited states. The IQA/CASPT2 approach allows determination of the energy redistribution that takes place along a photophysical/photochemical deactivation path in terms of self- and interatomic contributions. The applicability of the method is illustrated by the description of representative processes spanning different bonding regimes: noble gas excimer and exciplex formation, the reaction of ozone with a chlorine atom, and the photodissociations of formaldehyde and cyclobutane. These examples show the versatility of using CASPT2 with the significant information provided by the IQA partition to describe chemical processes with a large multiconfigurational character.
Collapse
Affiliation(s)
- Jesús Jara-Cortés
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic 63155, Mexico.
| | - Edith Leal-Sánchez
- Departamento de Física y Química Teórica, Facultad de Química, UNAM, México City 04510, Mexico
| | - Evelio Francisco
- Departamento de Química Física y Analítica, Faculta de Química, Universidad de Oviedo, Oviedo 33006, Spain
| | - José A Pérez-Pimienta
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic 63155, Mexico.
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Faculta de Química, Universidad de Oviedo, Oviedo 33006, Spain
| | | |
Collapse
|
24
|
Wang S, Li C, Evangelista FA. Analytic Energy Gradients for the Driven Similarity Renormalization Group Multireference Second-Order Perturbation Theory. J Chem Theory Comput 2021; 17:7666-7681. [PMID: 34839660 DOI: 10.1021/acs.jctc.1c00980] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We derive analytic energy gradients of the driven similarity renormalization group (DSRG) multireference second-order perturbation theory (MRPT2) using the method of Lagrange multipliers. In the Lagrangian, we impose constraints for a complete-active-space self-consistent-field reference wave function and the semicanonical orthonormal molecular orbitals. Solving for the associated Lagrange multipliers is found to share the same asymptotic scaling of a single DSRG-MRPT2 energy computation. A pilot implementation of the DSRG-MRPT2 analytic gradients is used to optimize the geometry of the singlet and triplet states of p-benzyne. The equilibrium bond lengths and angles are similar to those computed via other MRPT2s and Mukherjee's multireference coupled cluster theory. An approximate DSRG-MRPT2 method that neglects the contributions of the three-body density cumulant is found to introduce negligible errors in the geometry of p-benzyne, lending itself to a promising low-cost approach for molecular geometry optimizations using large active spaces.
Collapse
Affiliation(s)
- Shuhe Wang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
25
|
Segarra‐Martí J, Bearpark MJ. Modelling Photoionisation in Isocytosine: Potential Formation of Longer-Lived Excited State Cations in its Keto Form. Chemphyschem 2021; 22:2172-2181. [PMID: 34370368 PMCID: PMC8597144 DOI: 10.1002/cphc.202100402] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2021] [Revised: 08/09/2021] [Indexed: 11/24/2022]
Abstract
Studying the effects of UV and VUV radiation on non-canonical DNA/RNA nucleobases allows us to compare how they release excess energy following absorption with respect to their canonical counterparts. This has attracted much research attention in recent years because of its likely influence on the origin of our genetic lexicon in prebiotic times. Here we present a CASSCF and XMS-CASPT2 theoretical study of the photoionisation of non-canonical pyrimidine nucleobase isocytosine in both its keto and enol tautomeric forms. We analyse their lowest energy cationic excited states including 2 π + , 2 n O + and 2 n N + and compare these to the corresponding electronic states in cytosine. Investigating lower-energy decay pathways we find - unexpectedly - that keto-isocytosine+ presents a sizeable energy barrier potentially inhibiting decay to its cationic ground state, whereas enol-isocytosine+ features a barrierless and consequently ultrafast pathway analogous to the one previously found for the canonical (keto) form of cytosine+ . Dynamic electron correlation reduces the energy barrier in the keto form substantially (by ∼1 eV) but it is nevertheless still present. We additionally compute the UV/Vis absorption signals of the structures encountered along these decay channels to provide spectroscopic fingerprints to assist future experiments in monitoring these intricate photo-processes.
Collapse
Affiliation(s)
- Javier Segarra‐Martí
- Department of ChemistryMolecular Sciences Research HubImperial College LondonWhite City Campus, 82 Wood LaneLondonW12 0BZUK
- Present address: Instituto de Ciencia MolecularUniversitat de ValenciaP.O. Box 22085ES-46071ValenciaSpain
| | - Michael J. Bearpark
- Department of ChemistryMolecular Sciences Research HubImperial College LondonWhite City Campus, 82 Wood LaneLondonW12 0BZUK
| |
Collapse
|
26
|
Park JW. Analytical Gradient Theory for Resolvent-Fitted Second-Order Extended Multiconfiguration Perturbation Theory (XMCQDPT2). J Chem Theory Comput 2021; 17:6122-6133. [PMID: 34582217 DOI: 10.1021/acs.jctc.1c00613] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
We present the formulation and implementation of an analytical gradient algorithm for extended multiconfiguration quasidegenerate perturbation theory (XMCQDPT2) with the resolvent-fitting approximation by Granovsky. This algorithm is powerful when optimizing molecular configurations with a moderate-sized active space and many electronic states. First, we present the powerfulness and accuracy of resolvent-fitting approximations compared to canonical XMCQDPT2 theory. Then, we demonstrate the utility of the current algorithm in frequency analyses, optimizing the minimum energy conical intersection geometries of the retinal chromophore model RPSB6 and evaluating nuclear gradients when there are many electronic states. Furthermore, we parallelize the algorithm using the OpenMP/MPI hybrid approach. Additionally, we report the computational cost and parallel efficiency of the program.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
27
|
Matsika S. Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections. Chem Rev 2021; 121:9407-9449. [PMID: 34156838 DOI: 10.1021/acs.chemrev.1c00074] [Citation(s) in RCA: 50] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.
Collapse
Affiliation(s)
- Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
28
|
Affiliation(s)
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL, USA
| |
Collapse
|
29
|
Guo Y, Sivalingam K, Neese F. Approximations of density matrices in N-electron valence state second-order perturbation theory (NEVPT2). I. Revisiting the NEVPT2 construction. J Chem Phys 2021; 154:214111. [PMID: 34240991 DOI: 10.1063/5.0051211] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Over the last decade, the second-order N-electron valence state perturbation theory (NEVPT2) has developed into a widely used multireference perturbation method. To apply NEVPT2 to systems with large active spaces, the computational bottleneck is the construction of the fourth-order reduced density matrix. Both its generation and storage become quickly problematic beyond the usual maximum active space of about 15 active orbitals. To reduce the computational cost of handling fourth-order density matrices, the cumulant approximation (CU) has been proposed in several studies. A more conventional strategy to address the higher-order density matrices is the pre-screening approximation (PS), which is the default one in the ORCA program package since 2010. In the present work, the performance of the CU, PS, and extended PS (EPS) approximations for the fourth-order density matrices is compared. Following a pedagogical introduction to NEVPT2, contraction schemes, as well as the approximations to density matrices, and the intruder state problem are discussed. The CU approximation, while potentially leading to large computational savings, virtually always leads to intruder states. With the PS approximation, the computational savings are more modest. However, in conjunction with conservative cutoffs, it produces stable results. The EPS approximation to the fourth-order density matrices can reproduce very accurate NEVPT2 results without any intruder states. However, its computational cost is not much lower than that of the canonical algorithm. Moreover, we found that a good indicator of intrude states problems in any approximation to high order density matrices is the eigenspectra of the Koopmans matrices.
Collapse
Affiliation(s)
- Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, China
| | - Kantharuban Sivalingam
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
30
|
Nishimoto Y. Analytic gradients for restricted active space second-order perturbation theory (RASPT2). J Chem Phys 2021; 154:194103. [PMID: 34240887 DOI: 10.1063/5.0050074] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
The computational cost of analytic derivatives in multireference perturbation theory is strongly affected by the size of the active space employed in the reference self-consistent field calculation. To overcome previous limits on the active space size, the analytic gradients of single-state restricted active space second-order perturbation theory (RASPT2) and its complete active space second-order perturbation theory (CASPT2) have been developed and implemented in a local version of OpenMolcas. Similar to previous implementations of CASPT2, the RASPT2 implementation employs the Lagrangian or Z-vector method. The numerical results show that restricted active spaces with up to 20 electrons in 20 orbitals can now be employed for geometry optimizations.
Collapse
Affiliation(s)
- Yoshio Nishimoto
- Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
31
|
Inai N, Yokogawa D, Yanai T. Investigating the Nonradiative Decay Pathway in the Excited State of Silepin Derivatives: A Study with Second-Order Multireference Perturbation Wavefunction Theory. J Phys Chem A 2021; 125:559-569. [PMID: 33416309 DOI: 10.1021/acs.jpca.0c08738] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
The fluorescence quantum yield for fluorescent organic molecules is an important molecular property, and tuning it up is desired for various applications. For the computational estimation of the fluorescence quantum yield, the theoretical prediction of the nonradiative decay rate constant has become an attractive subject of study. The rate constant of thermally activated nonradiative decay is related to the activation energy in the photoreaction; thus, the accuracy and reliability of the excited-state potential energies in the quantum chemical computation are critical. In this study, we employed a second-order multireference perturbation wavefunction theory for studying the thermally activated decay via conical intersection (CI) of 1,1-dimethyldibenzo[b,f]silepin derivatives. The correlation between the computed activation energy to reach the CI geometry in the S1 state and the experimentally determined fluorescence quantum yield implied that silepins nonradiatively decay via the CI triggered by the twisting of the central C-C bond. Geometry optimization of the transition state using multireference perturbation theory drastically reduced the estimated activation energy. Our computation gave reasonable predictions of the activation free energies of photoexcited 1,1-dimethyldibenzo[b,f]silepin. The energy profiles and geometry optimizations using proper quantum chemical methods played a critical role in reliable estimation of the rate constant and fluorescence quantum yield.
Collapse
Affiliation(s)
- Naoto Inai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
| | - Daisuke Yokogawa
- Department of Basic Science, Graduate School of Arts and Science, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Takeshi Yanai
- 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-8602, Japan
| |
Collapse
|
32
|
Song C, Neaton JB, Martínez TJ. Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method. J Chem Phys 2021; 154:014103. [PMID: 33412861 DOI: 10.1063/5.0035233] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a reduced scaling and exact reformulation of state specific complete active space second-order perturbation (CASPT2) analytical gradients in terms of the MP2 and Fock derivatives using the supporting subspace method. This work follows naturally from the supporting subspace formulation of the CASPT2 energy in terms of the MP2 energy using dressed orbitals and Fock builds. For a given active space configuration, the terms corresponding to the MP2-gradient can be evaluated with O(N5) operations, while the rest of the calculations can be computed with O(N3) operations using Fock builds, Fock gradients, and linear algebra. When tensor-hyper-contraction is applied simultaneously, the computational cost can be further reduced to O(N4) for a fixed active space size. The new formulation enables efficient implementation of CASPT2 analytical gradients by leveraging the existing graphical processing unit (GPU)-based MP2 and Fock routines. We present benchmark results that demonstrate the accuracy and performance of the new method. Example applications of the new method in ab initio molecular dynamics simulation and constrained geometry optimization are given.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Physics, University of California Berkeley, Berkeley, California 94720, USA
| | - Jeffrey B Neaton
- Department of Physics, University of California Berkeley, Berkeley, California 94720, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
33
|
Calvin JA, Peng C, Rishi V, Kumar A, Valeev EF. Many-Body Quantum Chemistry on Massively Parallel Computers. Chem Rev 2020; 121:1203-1231. [DOI: 10.1021/acs.chemrev.0c00006] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Justus A. Calvin
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Varun Rishi
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
34
|
Conti I, Cerullo G, Nenov A, Garavelli M. Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going. J Am Chem Soc 2020; 142:16117-16139. [PMID: 32841559 PMCID: PMC7901644 DOI: 10.1021/jacs.0c04952] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Computational spectroscopy is becoming a mandatory tool for the interpretation of the
complex, and often congested, spectral maps delivered by modern non-linear multi-pulse
techniques. The fields of Electronic Structure Methods,
Non-Adiabatic Molecular Dynamics, and Theoretical
Spectroscopy represent the three pillars of the virtual ultrafast
optical spectrometer, able to deliver transient spectra in
silico from first principles. A successful simulation strategy requires a
synergistic approach that balances between the three fields, each one having its very
own challenges and bottlenecks. The aim of this Perspective is to demonstrate that,
despite these challenges, an impressive agreement between theory and experiment is
achievable now regarding the modeling of ultrafast photoinduced processes in complex
molecular architectures. Beyond that, some key recent developments in the three fields
are presented that we believe will have major impacts on spectroscopic simulations in
the very near future. Potential directions of development, pending challenges, and
rising opportunities are illustrated.
Collapse
Affiliation(s)
- Irene Conti
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Giulio Cerullo
- Dipartimento di Fisica, Politecnico di Milano, IFN-CNR, Piazza Leonardo da Vinci 32, I-20133 Milano, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| |
Collapse
|
35
|
Park JW. Analytical First-Order Derivatives of Second-Order Extended Multiconfiguration Quasi-Degenerate Perturbation Theory (XMCQDPT2): Implementation and Application. J Chem Theory Comput 2020; 16:5562-5571. [PMID: 32786905 DOI: 10.1021/acs.jctc.0c00389] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Analytical gradient theory for the second-order extended multiconfiguration quasi-degenerate perturbation theory (XMCQDPT2), which can be regarded as the multistate version of the multireference second-order Møller-Plesset perturbation theory (MRMP2), is formulated and implemented. The theory is similar to the previous analytical gradient theory for MCQDPT2, but we take into account the intruder state avoidance (ISA) technique and the "extension" of the MCQDPT2 theory by Granovsky. Although the (X)MCQDPT2 theory is not invariant with respect to rotations among the active orbitals, the resulting analytical gradients are accurate. We demonstrate the utility of the current algorithm in optimizing the minimum energy conical intersections (MECIs) of ethylene, butadiene, benzene, the retinal model chromophore PSB3, and the green fluorescent protein model chromophore pHBI. The XMCQDPT2 MECIs are very similar to the XMS-CASPT2 MECIs in terms of molecular conformation and the computed energies. We also discuss possible improvements of the current algorithm.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
36
|
Song C, Martínez TJ. Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction. J Chem Phys 2020; 152:234113. [PMID: 32571032 DOI: 10.1063/5.0007417] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N3) operations and O(N2) memory, where N is the number of basis functions. This limits the overall computational scaling to O(N4) operations and O(N2) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2. In addition, we also investigate how the energy denominators can be factorized with the Laplace quadrature when some of the denominators are negative, which is critical for excited state calculations. An efficient implementation of the method has been developed using graphical processing units while also exploiting spatial sparsity in tensor operations. We benchmark the accuracy of the new method by comparison to non-THC formulated XMS-CASPT2 for the excited states of various molecules. In our tests, the THC approximation introduces negligible errors (≈0.01 eV) compared to the non-THC reference method. Scaling behavior and computational timings are presented to demonstrate performance. The new method is also interfaced with quantum mechanics/molecular mechanics (QM/MM). In an example study of green fluorescent protein, we show how the XMS-CASPT2 potential energy surfaces and excitation energies are affected by increasing the size of the QM region up to 278 QM atoms with more than 2300 basis functions.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
37
|
Aprà E, Bylaska EJ, de Jong WA, Govind N, Kowalski K, Straatsma TP, Valiev M, van Dam HJJ, Alexeev Y, Anchell J, Anisimov V, Aquino FW, Atta-Fynn R, Autschbach J, Bauman NP, Becca JC, Bernholdt DE, Bhaskaran-Nair K, Bogatko S, Borowski P, Boschen J, Brabec J, Bruner A, Cauët E, Chen Y, Chuev GN, Cramer CJ, Daily J, Deegan MJO, Dunning TH, Dupuis M, Dyall KG, Fann GI, Fischer SA, Fonari A, Früchtl H, Gagliardi L, Garza J, Gawande N, Ghosh S, Glaesemann K, Götz AW, Hammond J, Helms V, Hermes ED, Hirao K, Hirata S, Jacquelin M, Jensen L, Johnson BG, Jónsson H, Kendall RA, Klemm M, Kobayashi R, Konkov V, Krishnamoorthy S, Krishnan M, Lin Z, Lins RD, Littlefield RJ, Logsdail AJ, Lopata K, Ma W, Marenich AV, Martin Del Campo J, Mejia-Rodriguez D, Moore JE, Mullin JM, Nakajima T, Nascimento DR, Nichols JA, Nichols PJ, Nieplocha J, Otero-de-la-Roza A, Palmer B, Panyala A, Pirojsirikul T, Peng B, Peverati R, Pittner J, Pollack L, Richard RM, Sadayappan P, Schatz GC, Shelton WA, Silverstein DW, Smith DMA, Soares TA, Song D, Swart M, Taylor HL, Thomas GS, Tipparaju V, Truhlar DG, Tsemekhman K, Van Voorhis T, Vázquez-Mayagoitia Á, Verma P, Villa O, Vishnu A, Vogiatzis KD, Wang D, Weare JH, Williamson MJ, Windus TL, Woliński K, Wong AT, Wu Q, Yang C, Yu Q, Zacharias M, Zhang Z, Zhao Y, Harrison RJ. NWChem: Past, present, and future. J Chem Phys 2020; 152:184102. [PMID: 32414274 DOI: 10.1063/5.0004997] [Citation(s) in RCA: 293] [Impact Index Per Article: 73.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.
Collapse
Affiliation(s)
- E Aprà
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - E J Bylaska
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - W A de Jong
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - N Govind
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - K Kowalski
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - T P Straatsma
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - M Valiev
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - H J J van Dam
- Brookhaven National Laboratory, Upton, New York 11973, USA
| | - Y Alexeev
- Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - J Anchell
- Intel Corporation, Santa Clara, California 95054, USA
| | - V Anisimov
- Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - F W Aquino
- QSimulate, Cambridge, Massachusetts 02139, USA
| | - R Atta-Fynn
- Department of Physics, The University of Texas at Arlington, Arlington, Texas 76019, USA
| | - J Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - N P Bauman
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - J C Becca
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - D E Bernholdt
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | | | - S Bogatko
- 4G Clinical, Wellesley, Massachusetts 02481, USA
| | - P Borowski
- Faculty of Chemistry, Maria Curie-Skłodowska University in Lublin, 20-031 Lublin, Poland
| | - J Boschen
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| | - J Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, 18223 Prague 8, Czech Republic
| | - A Bruner
- Department of Chemistry and Physics, University of Tennessee at Martin, Martin, Tennessee 38238, USA
| | - E Cauët
- Service de Chimie Quantique et Photophysique (CP 160/09), Université libre de Bruxelles, B-1050 Brussels, Belgium
| | - Y Chen
- Facebook, Menlo Park, California 94025, USA
| | - G N Chuev
- Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region 142290, Russia
| | - C J Cramer
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - J Daily
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M J O Deegan
- SKAO, Jodrell Bank Observatory, Macclesfield SK11 9DL, United Kingdom
| | - T H Dunning
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - M Dupuis
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - K G Dyall
- Dirac Solutions, Portland, Oregon 97229, USA
| | - G I Fann
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - S A Fischer
- Chemistry Division, U. S. Naval Research Laboratory, Washington, DC 20375, USA
| | - A Fonari
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - H Früchtl
- EaStCHEM and School of Chemistry, University of St. Andrews, St. Andrews KY16 9ST, United Kingdom
| | - L Gagliardi
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - J Garza
- Departamento de Química, División de Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana-Iztapalapa, Col. Vicentina, Iztapalapa, C.P. 09340 Ciudad de México, Mexico
| | - N Gawande
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - S Ghosh
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 5545, USA
| | - K Glaesemann
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - A W Götz
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - J Hammond
- Intel Corporation, Santa Clara, California 95054, USA
| | - V Helms
- Center for Bioinformatics, Saarland University, D-66041 Saarbrücken, Germany
| | - E D Hermes
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, USA
| | - K Hirao
- Next-generation Molecular Theory Unit, Advanced Science Institute, RIKEN, Saitama 351-0198, Japan
| | - S Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - M Jacquelin
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - L Jensen
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - B G Johnson
- Acrobatiq, Pittsburgh, Pennsylvania 15206, USA
| | - H Jónsson
- Faculty of Physical Sciences, University of Iceland, Reykjavík, Iceland and Department of Applied Physics, Aalto University, FI-00076 Aalto, Espoo, Finland
| | - R A Kendall
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - M Klemm
- Intel Corporation, Santa Clara, California 95054, USA
| | - R Kobayashi
- ANU Supercomputer Facility, Australian National University, Canberra, Australia
| | - V Konkov
- Chemistry Program, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - S Krishnamoorthy
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M Krishnan
- Facebook, Menlo Park, California 94025, USA
| | - Z Lin
- Department of Physics, University of Science and Technology of China, Hefei, China
| | - R D Lins
- Aggeu Magalhaes Institute, Oswaldo Cruz Foundation, Recife, Brazil
| | | | - A J Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Cardiff, Wales CF10 3AT, United Kingdom
| | - K Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - W Ma
- Institute of Software, Chinese Academy of Sciences, Beijing, China
| | - A V Marenich
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - J Martin Del Campo
- Departamento de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, México City, Mexico
| | - D Mejia-Rodriguez
- Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida 32611, USA
| | - J E Moore
- Intel Corporation, Santa Clara, California 95054, USA
| | - J M Mullin
- DCI-Solutions, Aberdeen Proving Ground, Maryland 21005, USA
| | - T Nakajima
- Computational Molecular Science Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - D R Nascimento
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - J A Nichols
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - P J Nichols
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - J Nieplocha
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - A Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, 33006 Oviedo, Spain
| | - B Palmer
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - A Panyala
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - T Pirojsirikul
- Department of Chemistry, Prince of Songkla University, Hat Yai, Songkhla 90112, Thailand
| | - B Peng
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - R Peverati
- Chemistry Program, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - J Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., 18223 Prague 8, Czech Republic
| | - L Pollack
- StudyPoint, Boston, Massachusetts 02114, USA
| | | | - P Sadayappan
- School of Computing, University of Utah, Salt Lake City, Utah 84112, USA
| | - G C Schatz
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - W A Shelton
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | | | - D M A Smith
- Intel Corporation, Santa Clara, California 95054, USA
| | - T A Soares
- Dept. of Fundamental Chemistry, Universidade Federal de Pernambuco, Recife, Brazil
| | - D Song
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M Swart
- ICREA, 08010 Barcelona, Spain and Universitat Girona, Institut de Química Computacional i Catàlisi, Campus Montilivi, 17003 Girona, Spain
| | - H L Taylor
- CD-adapco/Siemens, Melville, New York 11747, USA
| | - G S Thomas
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - V Tipparaju
- Cray Inc., Bloomington, Minnesota 55425, USA
| | - D G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | | | - T Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Á Vázquez-Mayagoitia
- Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - P Verma
- 1QBit, Vancouver, British Columbia V6E 4B1, Canada
| | - O Villa
- NVIDIA, Santa Clara, California 95051, USA
| | - A Vishnu
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - K D Vogiatzis
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - D Wang
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong 250014, China
| | - J H Weare
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - M J Williamson
- Department of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - T L Windus
- Department of Chemistry, Iowa State University and Ames Laboratory, Ames, Iowa 50011, USA
| | - K Woliński
- Faculty of Chemistry, Maria Curie-Skłodowska University in Lublin, 20-031 Lublin, Poland
| | - A T Wong
- Qwil, San Francisco, California 94107, USA
| | - Q Wu
- Brookhaven National Laboratory, Upton, New York 11973, USA
| | - C Yang
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Q Yu
- AMD, Santa Clara, California 95054, USA
| | - M Zacharias
- Department of Physics, Technical University of Munich, 85748 Garching, Germany
| | - Z Zhang
- Stanford Research Computing Center, Stanford University, Stanford, California 94305, USA
| | - Y Zhao
- State Key Laboratory of Silicate Materials for Architectures, International School of Materials Science and Engineering, Wuhan University of Technology, Wuhan 430070, China
| | - R J Harrison
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
38
|
Winslow M, Cross WB, Robinson D. Comparison of Spin-Flip TDDFT-Based Conical Intersection Approaches with XMS-CASPT2. J Chem Theory Comput 2020; 16:3253-3263. [PMID: 32302484 PMCID: PMC8279405 DOI: 10.1021/acs.jctc.9b00917] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Determining conical intersection
geometries is of key importance
to understanding the photochemical reactivity of molecules. While
many small- to medium-sized molecules can be treated accurately using
multireference approaches, larger molecules require a less computationally
demanding approach. In this work, minimum energy crossing point conical
intersection geometries for a series of molecules have been studied
using spin-flip TDDFT (SF-TDDFT), within the Tamm-Dancoff Approximation,
both with and without explicit calculation of nonadiabatic coupling
terms, and compared with both XMS-CASPT2 and CASSCF calculated geometries.
The less computationally demanding algorithms, which do not require
explicit calculation of the nonadiabatic coupling terms, generally
fare well with the XMS-CASPT2 reference structures, while the relative
energetics are only reasonably replicated with the MECP structure
as
calculated with the BHHLYP functional and full nonadiabatic coupling
terms. We also demonstrate that, occasionally, CASSCF structures deviate
quantitatively from the XMS-CASPT2 structures, showing the importance
of including dynamical correlation.
Collapse
Affiliation(s)
- Max Winslow
- Department of Chemistry and Forensics, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, United Kingdom
| | - Warren B Cross
- Department of Chemistry and Forensics, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, United Kingdom
| | - David Robinson
- Department of Chemistry and Forensics, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, United Kingdom
| |
Collapse
|
39
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
40
|
List NH, Dempwolff AL, Dreuw A, Norman P, Martínez TJ. Probing competing relaxation pathways in malonaldehyde with transient X-ray absorption spectroscopy. Chem Sci 2020; 11:4180-4193. [PMID: 34122881 PMCID: PMC8152795 DOI: 10.1039/d0sc00840k] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Excited-state intramolecular hydrogen transfer (ESIHT) is a fundamental reaction relevant to chemistry and biology. Malonaldehyde is the simplest example of ESIHT, yet only little is known experimentally about its excited-state dynamics. Several competing relaxation pathways have been proposed, including internal conversion mediated by ESIHT and C[double bond, length as m-dash]C torsional motion as well as intersystem crossing. We perform an in silico transient X-ray absorption spectroscopy (TRXAS) experiment at the oxygen K-edge to investigate its potential to monitor the proposed ultrafast decay pathways in malonaldehyde upon photoexcitation to its bright S2(ππ*) state. We employ both restricted active space perturbation theory and algebraic-diagrammatic construction for the polarization propagator along interpolated reaction coordinates as well as representative trajectories from ab initio multiple spawning simulations to compute the TRXAS signals from the lowest valence states. Our study suggests that oxygen K-edge TRXAS can distinctly fingerprint the passage through the H-transfer intersection and the concomitant population transfer to the S1(nπ*) state. Potential intersystem crossing to T1(ππ*) is detectable from reappearance of the double pre-edge signature and reversed intensities. Moreover, the torsional deactivation pathway induces transient charge redistribution from the enol side towards the central C-atom and manifests itself as substantial shifts of the pre-edge features. Given the continuous advances in X-ray light sources, our study proposes an experimental route to disentangle ultrafast excited-state decay channels in this prototypical ESIHT system and provides a pathway-specific mapping of the TRXAS signal to facilitate the interpretation of future experiments.
Collapse
Affiliation(s)
- Nanna H List
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Adrian L Dempwolff
- Interdisciplinary Center for Scientific Computing, Heidelberg University Im Neuenheimer Feld 205 D-69120 Heidelberg Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University Im Neuenheimer Feld 205 D-69120 Heidelberg Germany
| | - Patrick Norman
- School of Engineering Sciences in Chemistry, Biotechnology and Health (CBH), Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology Sweden
| | - Todd J Martínez
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| |
Collapse
|
41
|
Saitow M, Yanai T. A multireference coupled-electron pair approximation combined with complete-active space perturbation theory in local pair-natural orbital framework. J Chem Phys 2020; 152:114111. [PMID: 32199413 DOI: 10.1063/1.5142622] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Complete-Active Space Second-order Perturbation Theory (CASPT2) has been one of the most widely-used methods for reliably calculating electronic structures of multireference systems. Because of its lowest level treatment of dynamic correlation, it has a high computational feasibility; however, its accuracy in some cases falls short of needs. Here, as a simple yet higher-order alternative, we introduce a hybrid theory of the CASPT2 and a multireference variant of the Coupled-Electron Pair Approximation (CEPA), which is a class of high level correlation theory. A central feature of our theory (CEPT2) is to use the two underlying theories for describing different divisions of correlation components based on the full internal contraction framework. The external components, which usually give a major contribution to the dynamic correlation, are intensively described using the CEPA Ansatz, while the rests are treated at the CASPT2 level. Furthermore, to drastically reduce the computational demands, we have incorporated the pair-natural orbital (PNO) method into our multireference implementations. This development, thus, requires highly complex derivations and coding, while it has been largely facilitated with an automatic expression and code generation technique. To highlight the accuracy of the CEPT2 approach and to assess the errors caused by the PNO truncation, benchmark calculations are shown on small- to medium-size molecules, illustrating the high accuracy of the present CEPT2 model. By tightening the truncation thresholds, the PNO-CEPT2 energy converges toward the canonical counterpart and is more accurate than that of PNO-CASPT2 as long as the same truncation thresholds are used.
Collapse
Affiliation(s)
- Masaaki Saitow
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| | - Takeshi Yanai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| |
Collapse
|
42
|
Park JW. Analytical Gradient Theory for Quasidegenerate N-Electron Valence State Perturbation Theory (QD-NEVPT2). J Chem Theory Comput 2019; 16:326-339. [DOI: 10.1021/acs.jctc.9b00919] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
43
|
Nishimoto Y. Analytic first-order derivatives of partially contracted n-electron valence state second-order perturbation theory (PC-NEVPT2). J Chem Phys 2019; 151:114103. [PMID: 31542000 DOI: 10.1063/1.5115819] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
A balanced treatment of dynamic and static electron correlation is important in computational chemistry, and multireference perturbation theory (MRPT) is able to do this at a reasonable computational cost. In this paper, analytic first-order derivatives, specifically gradients and dipole moments, are developed for a particular MRPT method, state-specific partially contracted n-electron valence state second-order perturbation theory (PC-NEVPT2). Only one linear equation needs to be solved for the derivative calculation if the Z-vector method is employed, which facilitates the practical application of this approach. A comparison of the calculated results with experimental geometrical parameters of O3 indicates excellent agreement although the calculated results for O3 - are slightly outside the experimental error bars. The 0-0 transition energies of various methylpyrimidines and trans-polyacetylene are calculated by performing geometry optimizations and seminumerical second-order geometrical derivative calculations. In particular, the deviations of 0-0 transition energies of trans-polyacetylene from experimental values are consistently less than 0.1 eV with PC-NEVPT2, indicating the reliability of the method. These results demonstrate the importance of adding dynamic electron correlation on top of methods dominated by static electron correlation and of developing analytic derivatives for highly accurate methods.
Collapse
Affiliation(s)
- Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
44
|
Chatterjee K, Sokolov AY. Second-Order Multireference Algebraic Diagrammatic Construction Theory for Photoelectron Spectra of Strongly Correlated Systems. J Chem Theory Comput 2019; 15:5908-5924. [DOI: 10.1021/acs.jctc.9b00528] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Koushik Chatterjee
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
45
|
Park JW. Analytical Gradient Theory for Strongly Contracted (SC) and Partially Contracted (PC) N-Electron Valence State Perturbation Theory (NEVPT2). J Chem Theory Comput 2019; 15:5417-5425. [DOI: 10.1021/acs.jctc.9b00762] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
46
|
Wang S, Li C, Evangelista FA. Analytic gradients for the single-reference driven similarity renormalization group second-order perturbation theory. J Chem Phys 2019; 151:044118. [PMID: 31370522 DOI: 10.1063/1.5100175] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
We derive and implement analytic energy gradients for the single-reference driven similarity renormalization group second-order perturbation theory (DSRG-PT2). The resulting equations possess an asymptotic scaling that is identical to that of the second-order Møller-Plesset perturbation theory (MP2), indicating that the exponential regularizer in the DSRG equations does not introduce formal difficulties in the gradient theory. We apply the DSRG-PT2 method to optimizing the geometries of 15 small molecules. The equilibrium bond lengths computed with DSRG-PT2 are found similar to those of MP2, yielding a mean absolute error of 0.0033 Å and a standard deviation of 0.0045 Å when compared with coupled cluster with singles, doubles, and perturbative triples.
Collapse
Affiliation(s)
- Shuhe Wang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Chenyang Li
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
47
|
Feng X, Epifanovsky E, Gauss J, Krylov AI. Implementation of analytic gradients for CCSD and EOM-CCSD using Cholesky decomposition of the electron-repulsion integrals and their derivatives: Theory and benchmarks. J Chem Phys 2019; 151:014110. [DOI: 10.1063/1.5100022] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Xintian Feng
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| |
Collapse
|
48
|
Segarra‐Martí J, Tran T, Bearpark MJ. Computing the Ultrafast and Radiationless Electronic Excited State Decay of Cytosine and 5‐methyl‐cytosine Cations: Uncovering the Role of Dynamic Electron Correlation. CHEMPHOTOCHEM 2019. [DOI: 10.1002/cptc.201900105] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Javier Segarra‐Martí
- Department of Chemistry, Molecular Sciences Research HubImperial College London White City Campus, 80 Wood Lane W12 0BZ London UK
| | - Thierry Tran
- Department of Chemistry, Molecular Sciences Research HubImperial College London White City Campus, 80 Wood Lane W12 0BZ London UK
| | - Michael J. Bearpark
- Department of Chemistry, Molecular Sciences Research HubImperial College London White City Campus, 80 Wood Lane W12 0BZ London UK
| |
Collapse
|
49
|
Kats D, Werner HJ. Multi-state local complete active space second-order perturbation theory using pair natural orbitals (PNO-MS-CASPT2). J Chem Phys 2019; 150:214107. [DOI: 10.1063/1.5097644] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Daniel Kats
- Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
50
|
Park JW, Al-Saadon R, Strand NE, Shiozaki T. Imaginary Shift in CASPT2 Nuclear Gradient and Derivative Coupling Theory. J Chem Theory Comput 2019; 15:4088-4098. [DOI: 10.1021/acs.jctc.9b00368] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, United States
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju, Chungbuk 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, United States
| | - Nils E. Strand
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, United States
| |
Collapse
|