1
|
Wang YS, Zhong Manis JX, Rohan MC, Orlando TM, Kretchmer JS. Modeling Intermolecular Coulombic Decay with Non-Hermitian Real-Time Time-Dependent Density Functional Theory. J Phys Chem Lett 2024:7806-7813. [PMID: 39052307 DOI: 10.1021/acs.jpclett.4c01146] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
In this work, we investigate the capability of using real-time time-dependent density functional theory (RT-TDDFT) in conjunction with a complex absorbing potential (CAP) to simulate the intermolecular Coulombic decay (ICD) processes following the ionization of an inner-valence electron. We examine the ICD dynamics in a series of noncovalent bonded dimer systems, including hydrogen-bonded and purely van der Waals (VdW)-bonded systems. In comparison to previous work, we show that RT-TDDFT simulations with a CAP correctly capture the ICD phenomenon in systems exhibiting a stronger binding energy. The calculated time scales for ICD of the studied systems are in the range of 5-50 fs, in agreement with previous studies. However, there is a breakdown in the accuracy of the methodology for the pure VdW-bonded systems. Overall, the presented RT-TDDFT/CAP methodology provides a powerful tool for differentiating between competing electronic relaxation pathways following inner-valence or core ionization without necessitating any a priori assumptions.
Collapse
Affiliation(s)
- Yi-Siang Wang
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - James X Zhong Manis
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Matthew C Rohan
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Thomas M Orlando
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Joshua S Kretchmer
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| |
Collapse
|
2
|
Liu Z, Lyu N, Hu Z, Zeng H, Batista VS, Sun X. Benchmarking various nonadiabatic semiclassical mapping dynamics methods with tensor-train thermo-field dynamics. J Chem Phys 2024; 161:024102. [PMID: 38980091 DOI: 10.1063/5.0208708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2024] [Accepted: 06/20/2024] [Indexed: 07/10/2024] Open
Abstract
Accurate quantum dynamics simulations of nonadiabatic processes are important for studies of electron transfer, energy transfer, and photochemical reactions in complex systems. In this comparative study, we benchmark various approximate nonadiabatic dynamics methods with mapping variables against numerically exact calculations based on the tensor-train (TT) representation of high-dimensional arrays, including TT-KSL for zero-temperature dynamics and TT-thermofield dynamics for finite-temperature dynamics. The approximate nonadiabatic dynamics methods investigated include mixed quantum-classical Ehrenfest mean-field and fewest-switches surface hopping, linearized semiclassical mapping dynamics, symmetrized quasiclassical dynamics, the spin-mapping method, and extended classical mapping models. Different model systems were evaluated, including the spin-boson model for nonadiabatic dynamics in the condensed phase, the linear vibronic coupling model for electronic transition through conical intersections, the photoisomerization model of retinal, and Tully's one-dimensional scattering models. Our calculations show that the optimal choice of approximate dynamical method is system-specific, and the accuracy is sensitively dependent on the zero-point-energy parameter and the initial sampling strategy for the mapping variables.
Collapse
Affiliation(s)
- Zengkui Liu
- Division of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, USA
| | - Ningyi Lyu
- Division of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Zhubin Hu
- Division of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China
| | - Hao Zeng
- Division of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai 200124, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, USA
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China
| |
Collapse
|
3
|
Weymuth T, Unsleber JP, Türtscher PL, Steiner M, Sobez JG, Müller CH, Mörchen M, Klasovita V, Grimmel SA, Eckhoff M, Csizi KS, Bosia F, Bensberg M, Reiher M. SCINE-Software for chemical interaction networks. J Chem Phys 2024; 160:222501. [PMID: 38857173 DOI: 10.1063/5.0206974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 05/09/2024] [Indexed: 06/12/2024] Open
Abstract
The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge.
Collapse
Affiliation(s)
- Thomas Weymuth
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan P Unsleber
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Paul L Türtscher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Miguel Steiner
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan-Grimo Sobez
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Charlotte H Müller
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Maximilian Mörchen
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Veronika Klasovita
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stephanie A Grimmel
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Marco Eckhoff
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Katja-Sophia Csizi
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Francesco Bosia
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Moritz Bensberg
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
4
|
Hino K, Kurashige Y. Encoding a Many-Body Potential Energy Surface into a Grid-Based Matrix Product Operator. J Chem Theory Comput 2024; 20:3839-3849. [PMID: 38647101 DOI: 10.1021/acs.jctc.4c00046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
An efficient algorithm for compressing a given many-body potential energy surface (PES) of molecular systems into a grid-based matrix product operator (MPO) is proposed. The PES is once represented by a full-dimensional or truncated many-body expansion form, which is obtained by ab initio calculations at each grid mesh point, and then all terms in the expansion are compressed and merged into a single MPO while maintaining the bond dimension of the MPO as small as possible. It was shown that the ab initio PES of the H2CO was compressed by more than 2 orders of magnitude in the size of the site operators without loss of accuracy. By the use of grid basis, the tensor rank of the site operators of the MPO is reduced from four to three due to the diagonal nature of the position-dependent operators on grid basis, which significantly reduces the computational cost of the tensor contractions required in the real and imaginary time evolution of the matrix product state (MPS) wave functions with the grid-based MPO (Grid-MPO) Hamiltonian. Similar to other grid-based methods, Grid-MPO is easily applicable to any kinds of potentials of molecular systems, such as analytical empirical model potentials expressed by position operators and ab initio potentials, if the values at the grid points are available. Using the Grid-MPO combined with the MPS, we calculated the time correlation function of the Eigen cation H 3 O + ( H 2 O ) 3 to predict the infrared spectrum and compared with the experimental and the previous theoretical studies. The actual scaling with the size of systems was examined for the multidimensional Henon-Heiles Hamiltonian. It was shown that the method is considerably accelerated by the graphic processing unit (GPU) because the sizes of site operators were kept small and all tensors were able to be stored on the VRAM of a GPU.
Collapse
Affiliation(s)
- Kentaro Hino
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
- FOREST, JST, Honcho 4-1-8, Kawaguchi, Saitama 332-0012, Japan
- CREST, JST, Honcho 4-1-8, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
5
|
Glaser N, Baiardi A, Reiher M. Flexible DMRG-Based Framework for Anharmonic Vibrational Calculations. J Chem Theory Comput 2023; 19:9329-9343. [PMID: 38060309 PMCID: PMC10753801 DOI: 10.1021/acs.jctc.3c00902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 10/16/2023] [Accepted: 10/23/2023] [Indexed: 12/08/2023]
Abstract
We present a novel formulation of the vibrational density matrix renormalization group (vDMRG) algorithm tailored to strongly anharmonic molecules described by general, high-dimensional model representations of potential energy surfaces. For this purpose, we extend the vDMRG framework to support vibrational Hamiltonians expressed in the so-called n-mode second-quantization formalism. The resulting n-mode vDMRG method offers full flexibility with respect to both the functional form of the PES and the choice of the single-particle basis set. We leverage this framework to apply, for the first time, vDMRG based on an anharmonic modal basis set optimized with the vibrational self-consistent field algorithm on an on-the-fly constructed PES. We also extend the n-mode vDMRG framework to include excited-state-targeting algorithms in order to efficiently calculate anharmonic transition frequencies. We demonstrate the capabilities of our novel n-mode vDMRG framework for methyloxirane, a challenging molecule with 24 coupled vibrational modes.
Collapse
Affiliation(s)
- Nina Glaser
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Alberto Baiardi
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
6
|
Jensen AB, Højlund MG, Zoccante A, Madsen NK, Christiansen O. Efficient time-dependent vibrational coupled cluster computations with time-dependent basis sets at the two-mode coupling level: Full and hybrid TDMVCC[2]. J Chem Phys 2023; 159:204106. [PMID: 38010335 DOI: 10.1063/5.0175506] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Accepted: 11/05/2023] [Indexed: 11/29/2023] Open
Abstract
The computation of the nuclear quantum dynamics of molecules is challenging, requiring both accuracy and efficiency to be applicable to systems of interest. Recently, theories have been developed for employing time-dependent basis functions (denoted modals) with vibrational coupled cluster theory (TDMVCC). The TDMVCC method was introduced along with a pilot implementation, which illustrated good accuracy in benchmark computations. In this paper, we report an efficient implementation of TDMVCC, covering the case where the wave function and Hamiltonian contain up to two-mode couplings. After a careful regrouping of terms, the wave function can be propagated with a cubic computational scaling with respect to the number of degrees of freedom. We discuss the use of a restricted set of active one-mode basis functions for each mode, as well as two interesting limits: (i) the use of a full active basis where the variational modal determination amounts essentially to the variational determination of a time-dependent reference state for the cluster expansion; and (ii) the use of a single function as an active basis for some degrees of freedom. The latter case defines a hybrid TDMVCC/TDH (time-dependent Hartree) approach that can obtain even lower computational scaling. The resulting computational scaling for hybrid and full TDMVCC[2] is illustrated for polyaromatic hydrocarbons with up to 264 modes. Finally, computations on the internal vibrational redistribution of benzoic acid (39 modes) are used to show the faster convergence of TDMVCC/TDH hybrid computations towards TDMVCC compared to simple neglect of some degrees of freedom.
Collapse
Affiliation(s)
| | - Mads Greisen Højlund
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Alberto Zoccante
- Dipartimento di Scienze e Innovazione Tecnologica, Università del Piemonte Orientale (UPO), Via T. Michel 11, 15100 Alessandria, Italy
| | - Niels Kristian Madsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
7
|
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
|
8
|
Liu Z, Hu H, Sun X. Multistate Reaction Coordinate Model for Charge and Energy Transfer Dynamics in the Condensed Phase. J Chem Theory Comput 2023; 19:7151-7170. [PMID: 37815937 PMCID: PMC10601487 DOI: 10.1021/acs.jctc.3c00770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Indexed: 10/12/2023]
Abstract
Constructing multistate model Hamiltonians from all-atom electronic structure calculations and molecular dynamics simulations is crucial for understanding charge and energy transfer dynamics in complex condensed phases. The most popular two-level system model is the spin-boson Hamiltonian, where the nuclear degrees of freedom are represented as shifted normal modes. Recently, we proposed the general multistate nontrivial extension of the spin-boson model, i.e., the multistate harmonic (MSH) model, which is constructed by extending the spatial dimensions of each nuclear mode so as to satisfy the all-atom reorganization energy restrictions for all pairs of electronic states. In this work, we propose the multistate reaction coordinate (MRC) model with a primary reaction coordinate and secondary bath modes as in the Caldeira-Leggett form but in extended spatial dimensions. The MRC model is proven to be equivalent to the MSH model and offers an intuitive physical picture of the nuclear-electronic feedback in nonadiabatic processes such as the inherent trajectory of the reaction coordinate. The reaction coordinate is represented in extended dimensions, carrying the entire reorganization energies and bilinearly coupled to the secondary bath modes. We demonstrate the MRC model construction for photoinduced charge transfer in an organic photovoltaic caroteniod-porphyrin-C60 molecular triad dissolved in tetrahydrofuran as well as excitation energy transfer in a photosynthetic light-harvesting Fenna-Matthews-Olson complex. The MRC model provides an effective and robust platform for investigating quantum dissipative dynamics in complex condensed-phase systems since it allows a consistent description of realistic spectral density, state-dependent system-bath couplings, and heterogeneous environments due to static disorder in reorganization energies.
Collapse
Affiliation(s)
- Zengkui Liu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
| | - Haorui Hu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| | - Xiang Sun
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
- Shanghai
Frontiers Science Center of Artificial Intelligence and Deep Learning, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| |
Collapse
|
9
|
Iyengar SS, Kumar A, Saha D, Sabry A. Synthesis of Hidden Subgroup Quantum Algorithms and Quantum Chemical Dynamics. J Chem Theory Comput 2023; 19:6082-6092. [PMID: 37703187 DOI: 10.1021/acs.jctc.3c00404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2023]
Abstract
We describe a general formalism for quantum dynamics and show how this formalism subsumes several quantum algorithms, including the Deutsch, Deutsch-Jozsa, Bernstein-Vazirani, Simon, and Shor algorithms as well as the conventional approach to quantum dynamics based on tensor networks. The common framework exposes similarities among quantum algorithms and natural quantum phenomena: we illustrate this connection by showing how the correlated behavior of protons in water wire systems that are common in many biological and materials systems parallels the structure of Shor's algorithm.
Collapse
Affiliation(s)
- Srinivasan S Iyengar
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405-7102, United States
- Quantum Science and Engineering Center (QSEc), Indiana University, Bloomington, Indiana 47405-7102, United States
| | - Anup Kumar
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405-7102, United States
| | - Debadrita Saha
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405-7102, United States
| | - Amr Sabry
- Quantum Science and Engineering Center (QSEc), Indiana University, Bloomington, Indiana 47405-7102, United States
- Department of Computer Science, Luddy School of Informatics, Computing, and Engineering, Indiana University, Bloomington, Indiana 47405-7102, United States
| |
Collapse
|
10
|
Wang Y, Ren J, Shuai Z. Minimizing non-radiative decay in molecular aggregates through control of excitonic coupling. Nat Commun 2023; 14:5056. [PMID: 37598183 PMCID: PMC10439946 DOI: 10.1038/s41467-023-40716-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2022] [Accepted: 08/03/2023] [Indexed: 08/21/2023] Open
Abstract
The widely known "Energy Gap Law" (EGL) predicts a monotonically exponential increase in the non-radiative decay rate (knr) as the energy gap narrows, which hinders the development of near-infrared (NIR) emissive molecular materials. Recently, several experiments proposed that the exciton delocalization in molecular aggregates could counteract EGL to facilitate NIR emission. In this work, the nearly exact time-dependent density matrix renormalization group (TD-DMRG) method is developed to evaluate the non-radiative decay rate for exciton-phonon coupled molecular aggregates. Systematical numerical simulations show, by increasing the excitonic coupling, knr will first decrease, then reach a minimum, and finally start to increase to follow EGL, which is an overall result of two opposite effects of a smaller energy gap and a smaller effective electron-phonon coupling. This anomalous non-monotonic behavior is found robust in a number of models, including dimer, one-dimensional chain, and two-dimensional square lattice. The optimal excitonic coupling strength that gives the minimum knr is about half of the monomer reorganization energy and is also influenced by system size, dimensionality, and temperature.
Collapse
Affiliation(s)
- Yuanheng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084, Beijing, People's Republic of China
| | - Jiajun Ren
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 100875, Beijing, People's Republic of China.
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084, Beijing, People's Republic of China.
- School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen, 518172, People's Republic of China.
| |
Collapse
|
11
|
Yehorova D, Kretchmer JS. A multi-fragment real-time extension of projected density matrix embedding theory: Non-equilibrium electron dynamics in extended systems. J Chem Phys 2023; 158:131102. [PMID: 37031109 DOI: 10.1063/5.0146973] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/07/2023] Open
Abstract
In this work, we derive a multi-fragment real-time extension of the projected density matrix embedding theory (pDMET) designed to treat non-equilibrium electron dynamics in strongly correlated systems. As in the previously developed static pDMET, the real time pDMET partitions the total system into many fragments; the coupling between each fragment and the rest of the system is treated through a compact representation of the environment in terms of a quantum bath. The real-time pDMET involves simultaneously propagating the wavefunctions for each separate fragment–bath embedding system along with an auxiliary mean-field wavefunction of the total system. The equations of motion are derived by (i) projecting the time-dependent Schrödinger equation in the fragment and bath space associated with each separate fragment and by (ii) enforcing the pDMET matching conditions between the global 1-particle reduced density matrix (1-RDM) obtained from the fragment calculations and the mean-field 1-RDM at all points in time. The accuracy of the method is benchmarked through comparisons to time-dependent density-matrix renormalization group and time-dependent Hartree–Fock (TDHF) theory; the methods were applied to a one- and two-dimensional single-impurity Anderson model and multi-impurity Anderson models with ordered and disordered distributions of the impurities. The results demonstrate a large improvement over TDHF and rapid convergence to the exact dynamics with an increase in fragment size. Our results demonstrate that the real-time pDMET is a promising and flexible method that balances accuracy and efficiency to simulate the non-equilibrium electron dynamics in heterogeneous systems of large size.
Collapse
Affiliation(s)
- Dariia Yehorova
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Joshua S. Kretchmer
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
12
|
Abstract
We present a nonadiabatic classical-trajectory approach that offers the best of both worlds between fewest-switches surface hopping (FSSH) and quasiclassical mapping dynamics. This mapping approach to surface hopping (MASH) propagates the nuclei on the active adiabatic potential-energy surface, such as in FSSH. However, unlike in FSSH, transitions between active surfaces are deterministic and occur when the electronic mapping variables evolve between specified regions of the electronic phase space. This guarantees internal consistency between the active surface and the electronic degrees of freedom throughout the dynamics. MASH is rigorously derivable from exact quantum mechanics as a limit of the quantum-classical Liouville equation (QCLE), leading to a unique prescription for momentum rescaling and frustrated hops. Hence, a quantum-jump procedure can, in principle, be used to systematically converge the accuracy of the results to that of the QCLE. This jump procedure also provides a rigorous framework for deriving approximate decoherence corrections similar to those proposed for FSSH. We apply MASH to simulate the nonadiabatic dynamics in various model systems and show that it consistently produces more accurate results than FSSH at a comparable computational cost.
Collapse
|
13
|
Patil P, Heyl M, Alet F. Anomalous relaxation of density waves in a ring-exchange system. Phys Rev E 2023; 107:034119. [PMID: 37072977 DOI: 10.1103/physreve.107.034119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2022] [Accepted: 02/27/2023] [Indexed: 04/20/2023]
Abstract
We present the analysis of the slowing down exhibited by stochastic dynamics of a ring-exchange model on a square lattice, by means of numerical simulations. We find the preservation of coarse-grained memory of initial state of density-wave types for unexpectedly long times. This behavior is inconsistent with the prediction from a low frequency continuum theory developed by assuming a mean-field solution. Through a detailed analysis of correlation functions of the dynamically active regions, we exhibit an unconventional transient long ranged structure formation in a direction which is featureless for the initial condition, and argue that its slow melting plays a crucial role in the slowing-down mechanism. We expect our results to be relevant also for the dynamics of quantum ring-exchange dynamics of hard-core bosons and more generally for dipole moment conserving models.
Collapse
Affiliation(s)
- Pranay Patil
- Max-Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
| | - Markus Heyl
- Max-Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany
- Theoretical Physics III, Center for Electronic Correlations and Magnetism, Institute of Physics, University of Augsburg, 86135 Augsburg, Germany
| | - Fabien Alet
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
| |
Collapse
|
14
|
Xu Y, Liu C, Ma H. Hierarchical Mapping for Efficient Simulation of Strong System-Environment Interactions. J Chem Theory Comput 2023; 19:426-435. [PMID: 36626721 DOI: 10.1021/acs.jctc.2c00851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Quantum dynamics (QD) simulation is a powerful tool for interpreting ultrafast spectroscopy experiments and unraveling their microscopic mechanism in out-of-equilibrium excited state behaviors in various chemical, biological, and material systems. Although state-of-the-art numerical QD approaches such as the time-dependent density matrix renormalization group (TD-DMRG) already greatly extended the solvable system size of general linearly coupled exciton-phonon models with up to a few hundred phonon modes, the accurate simulation of larger system sizes or strong system-environment interactions is still computationally highly challenging. Based on quantum information theory (QIT), in this work, we realize that only a small number of effective phonon modes couple to the excitonic system directly regardless of a large or even infinite number of modes in the condensed phase environment. On top of the identified small number of direct effective modes, we propose a hierarchical mapping (HM) approach through performing block Lanczos transformations on the remaining indirect modes, which transforms the Hamiltonian matrix to a nearly block-tridiagonal form and eliminates the long-range interactions. Numerical tests on model spin-boson systems and realistic singlet fission models in a rubrene crystal environment with up to 7000 modes and strong system-environment interactions indicate HM can reduce the system size by 1-2 orders of magnitude and accelerate the calculation by ∼80% without losing accuracy.
Collapse
Affiliation(s)
- Yihe Xu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Chungen Liu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, China
| |
Collapse
|
15
|
Miessen A, Ollitrault PJ, Tacchino F, Tavernelli I. Quantum algorithms for quantum dynamics. NATURE COMPUTATIONAL SCIENCE 2023; 3:25-37. [PMID: 38177956 DOI: 10.1038/s43588-022-00374-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2021] [Accepted: 11/12/2022] [Indexed: 01/06/2024]
Abstract
Among the many computational challenges faced across different disciplines, quantum-mechanical systems pose some of the hardest ones and offer a natural playground for the growing field of quantum technologies. In this Perspective, we discuss quantum algorithmic solutions for quantum dynamics, reporting on the latest developments and offering a viewpoint on their potential and current limitations. We present some of the most promising areas of application and identify possible research directions for the coming years.
Collapse
Affiliation(s)
| | - Pauline J Ollitrault
- IBM Quantum, IBM Research - Zurich, Rüschlikon, Switzerland
- QC Ware, Palo Alto, CA, USA
| | | | | |
Collapse
|
16
|
Gelin MF, Chen L, Domcke W. Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals. Chem Rev 2022; 122:17339-17396. [PMID: 36278801 DOI: 10.1021/acs.chemrev.2c00329] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Lipeng Chen
- Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching,Germany
| |
Collapse
|
17
|
Højlund MG, Jensen AB, Zoccante A, Christiansen O. Bivariational time-dependent wave functions with biorthogonal adaptive basis sets: General formulation and regularization of equations of motion through polar decomposition. J Chem Phys 2022; 157:234104. [PMID: 36550053 DOI: 10.1063/5.0127431] [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/30/2022] Open
Abstract
We derive general bivariational equations of motion (EOMs) for time-dependent wave functions with biorthogonal time-dependent basis sets. The time-dependent basis functions are linearly parameterized and their fully variational time evolution is ensured by solving a set of so-called constraint equations, which we derive for arbitrary wave function expansions. The formalism allows division of the basis set into an active basis and a secondary basis, ensuring a flexible and compact wave function. We show how the EOMs specialize to a few common wave function forms, including coupled cluster and linearly expanded wave functions. It is demonstrated, for the first time, that the propagation of such wave functions is not unconditionally stable when a secondary basis is employed. The main signature of the instability is a strong increase in non-orthogonality, which eventually causes the calculation to fail; specifically, the biorthogonal active bra and ket bases tend toward spanning different spaces. Although formally allowed, this causes severe numerical issues. We identify the source of this problem by reparametrizing the time-dependent basis set through polar decomposition. Subsequent analysis allows us to remove the instability by setting appropriate matrix elements to zero. Although this solution is not fully variational, we find essentially no deviation in terms of autocorrelation functions relative to the variational formulation. We expect that the results presented here will be useful for the formal analysis of bivariational time-dependent wave functions for electronic and nuclear dynamics in general and for the practical implementation of time-dependent CC wave functions in particular.
Collapse
Affiliation(s)
- Mads Greisen Højlund
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | | | - Alberto Zoccante
- Dipartimento di Scienze e Innovazione Tecnologica, Universitá del Piemonte Orientale (UPO), Via T. Michel 11, 15100 Alessandria, Italy
| | - Ove Christiansen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
18
|
Kumar A, DeGregorio N, Ricard T, Iyengar SS. Graph-Theoretic Molecular Fragmentation for Potential Surfaces Leads Naturally to a Tensor Network Form and Allows Accurate and Efficient Quantum Nuclear Dynamics. J Chem Theory Comput 2022; 18:7243-7259. [PMID: 36332133 DOI: 10.1021/acs.jctc.2c00484] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Molecular fragmentation methods have revolutionized quantum chemistry. Here, we use a graph-theoretically generated molecular fragmentation method, to obtain accurate and efficient representations for multidimensional potential energy surfaces and the quantum time-evolution operator, which plays a critical role in quantum chemical dynamics. In doing so, we find that the graph-theoretic fragmentation approach naturally reduces the potential portion of the time-evolution operator into a tensor network that contains a stream of coupled lower-dimensional propagation steps to potentially achieve quantum dynamics with reduced complexity. Furthermore, the fragmentation approach used here has previously been shown to allow accurate and efficient computation of post-Hartree-Fock electronic potential energy surfaces, which in many cases has been shown to be at density functional theory cost. Thus, by combining the advantages of molecular fragmentation with the tensor network formalism, the approach yields an on-the-fly quantum dynamics scheme where both the electronic potential calculation and nuclear propagation portion are enormously simplified through a single stroke. The method is demonstrated by computing approximations to the propagator and to potential surfaces for a set of coupled nuclear dimensions within a protonated water wire problem exhibiting the Grotthuss mechanism of proton transport. In all cases, our approach has been shown to reduce the complexity of representing the quantum propagator, and by extension action of the propagator on an initial wavepacket, by several orders, with minimal loss in accuracy.
Collapse
Affiliation(s)
- Anup Kumar
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Nicole DeGregorio
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Timothy Ricard
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
19
|
Rajpurohit S, Simoni J, Tan LZ. Photo-induced phase-transitions in complex solids. NANOSCALE ADVANCES 2022; 4:4997-5008. [PMID: 36504738 PMCID: PMC9680828 DOI: 10.1039/d2na00481j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Accepted: 11/01/2022] [Indexed: 06/17/2023]
Abstract
Photo-induced phase-transitions (PIPTs) driven by highly cooperative interactions are of fundamental interest as they offer a way to tune and control material properties on ultrafast timescales. Due to strong correlations and interactions, complex quantum materials host several fascinating PIPTs such as light-induced charge density waves and ferroelectricity and have become a desirable setting for studying these PIPTs. A central issue in this field is the proper understanding of the underlying mechanisms driving the PIPTs. As these PIPTs are highly nonlinear processes and often involve multiple time and length scales, different theoretical approaches are often needed to understand the underlying mechanisms. In this review, we present a brief overview of PIPTs realized in complex materials, followed by a discussion of the available theoretical methods with selected examples of recent progress in understanding of the nonequilibrium pathways of PIPTs.
Collapse
Affiliation(s)
| | - Jacopo Simoni
- Molecular Foundry, Lawrence Berkeley National Laboratory USA
| | - Liang Z Tan
- Molecular Foundry, Lawrence Berkeley National Laboratory USA
| |
Collapse
|
20
|
Blavier M, Gelfand N, Levine R, Remacle F. Entanglement of electrons and nuclei: A most compact representation of the molecular wave function. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.139885] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
21
|
Ge Y, Li W, Ren J, Shuai Z. Computational Method for Evaluating the Thermoelectric Power Factor for Organic Materials Modeled by the Holstein Model: A Time-Dependent Density Matrix Renormalization Group Formalism. J Chem Theory Comput 2022; 18:6437-6446. [PMID: 36174220 DOI: 10.1021/acs.jctc.2c00651] [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/28/2022]
Abstract
Organic/polymeric materials are of emerging importance for thermoelectric conversion. The soft nature of these materials implies strong electron-phonon coupling, often leading to carrier localization. This poses great challenges for the conventional Boltzmann transport description based on relaxation time approximation and band structure calculations. In this work, combining the Kubo formula with the finite-temperature time-dependent density matrix renormalization group (FT-TD-DMRG) in the grand canonical ensemble, we developed a nearly exact algorithm to calculate the thermoelectric power factor PF = α2 σ, where α is the Seebeck coefficient and σ is the electrical conductivity, and apply the algorithm to Holstein Hamiltonian with electron-phonon coupling to model organic materials. Our algorithm can provide a unified description covering the weak coupling limit described by the bandlike Boltzmann transport to the strong coupling hopping limit.
Collapse
Affiliation(s)
- Yufei Ge
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084Beijing, P. R. China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084Beijing, P. R. China
| | - Jiajun Ren
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, 100875Beijing, P. R. China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084Beijing, P. R. China.,School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen518172, Guangdong, P. R. China
| |
Collapse
|
22
|
Baiardi A, Lesiuk M, Reiher M. Explicitly Correlated Electronic Structure Calculations with Transcorrelated Matrix Product Operators. J Chem Theory Comput 2022; 18:4203-4217. [PMID: 35666238 DOI: 10.1021/acs.jctc.2c00167] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, we present the first implementation of the transcorrelated electronic Hamiltonian in an optimization procedure for matrix product states by the density matrix renormalization group (DMRG) algorithm. In the transcorrelation ansatz, the electronic Hamiltonian is similarity-transformed with a Jastrow factor to describe the cusp in the wave function at electron-electron coalescence. As a result, the wave function is easier to approximate accurately with the conventional expansion in terms of one-particle basis functions and Slater determinants. The transcorrelated Hamiltonian in first quantization comprises up to three-body interactions, which we deal with in the standard way by applying robust density fitting to two- and three-body integrals entering the second-quantized representation of this Hamiltonian. The lack of hermiticity of the transcorrelated Hamiltonian is taken care of along the lines of the first work on transcorrelated DMRG [ J. Chem. Phys. 2020, 153, 164115] by encoding it as a matrix product operator and optimizing the corresponding ground state wave function with imaginary-time time-dependent DMRG. We demonstrate our quantum chemical transcorrelated DMRG approach at the example of several atoms and first-row diatomic molecules. We show that transcorrelation improves the convergence rate to the complete basis set limit in comparison to conventional DMRG. Moreover, we study extensions of our approach that aim at reducing the cost of handling the matrix product operator representation of the transcorrelated Hamiltonian.
Collapse
Affiliation(s)
- Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Michał Lesiuk
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland.,Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
23
|
Lyu N, Soley MB, Batista VS. Tensor-Train Split-Operator KSL (TT-SOKSL) Method for Quantum Dynamics Simulations. J Chem Theory Comput 2022; 18:3327-3346. [PMID: 35649210 DOI: 10.1021/acs.jctc.2c00209] [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/30/2022]
Abstract
Numerically exact simulations of quantum reaction dynamics, including nonadiabatic effects in excited electronic states, are essential to gain fundamental insights into ultrafast chemical reactivity and rigorous interpretations of molecular spectroscopy. Here, we introduce the tensor-train split-operator KSL (TT-SOKSL) method for quantum simulations in tensor-train (TT)/matrix product state (MPS) representations. TT-SOKSL propagates the quantum state as a tensor train using the Trotter expansion of the time-evolution operator, as in the tensor-train split-operator Fourier transform (TT-SOFT) method. However, the exponential operators of the Trotter expansion are applied using a rank-adaptive TT-KSL scheme instead of using the scaling and squaring approach as in TT-SOFT. We demonstrate the accuracy and efficiency of TT-SOKSL as applied to simulations of the photoisomerization of the retinal chromophore in rhodopsin, including nonadiabatic dynamics at a conical intersection of potential energy surfaces. The quantum evolution is described in full dimensionality by a time-dependent wavepacket evolving according to a two-state 25-dimensional model Hamiltonian. We find that TT-SOKSL converges faster than TT-SOFT with respect to the maximally allowed memory requirement of the tensor-train representation and better preserves the norm of the time-evolving state. When compared to the corresponding simulations based on the TT-KSL method, TT-SOKSL has the advantage of avoiding the need to construct the matrix product state Laplacian by exploiting the linear scaling of multidimensional tensor-train Fourier transforms.
Collapse
Affiliation(s)
- Ningyi Lyu
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
| | - Micheline B Soley
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States.,Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States.,Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States
| |
Collapse
|
24
|
Hino K, Kurashige Y. Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials. J Chem Theory Comput 2022; 18:3347-3356. [PMID: 35606892 DOI: 10.1021/acs.jctc.2c00243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH2O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C2nH2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.
Collapse
Affiliation(s)
- Kentaro Hino
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
25
|
Li W, Ren J, Yang H, Shuai Z. On the fly swapping algorithm for ordering of degrees of freedom in density matrix renormalization group. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:254003. [PMID: 35378514 DOI: 10.1088/1361-648x/ac640e] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/23/2022] [Accepted: 04/04/2022] [Indexed: 06/14/2023]
Abstract
Density matrix renormalization group (DMRG) and its time-dependent variants have found widespread applications in quantum chemistry, includingab initioelectronic structure of complex bio-molecules, spectroscopy for molecular aggregates, and charge transport in bulk organic semiconductors. The underlying wavefunction ansatz for DMRG, matrix product state (MPS), requires mapping degrees of freedom (DOF) into a one-dimensional topology. DOF ordering becomes a crucial factor for DMRG accuracy. In this work, we propose swapping neighboring DOFs during the DMRG sweeps for DOF ordering, which we term 'on the fly swapping' (OFS) algorithm. We show that OFS is universal for both static and time-dependent DMRG with minimum computational overhead. Examples are given for one dimensional antiferromagnetic Heisenberg model,ab initioelectronic structure of N2molecule, and the S1/S2internal conversion dynamics of pyrazine molecule. It is found that OFS can indeed improve accuracy by finding better DOF ordering in all cases.
Collapse
Affiliation(s)
- Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| | - Hengrui Yang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| |
Collapse
|
26
|
Ren J, Li W, Jiang T, Wang Y, Shuai Z. Time‐dependent density matrix renormalization group method for quantum dynamics in complex systems. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1614] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Yuanheng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| |
Collapse
|
27
|
Xu Y, Xie Z, Xie X, Schollwöck U, Ma H. Stochastic Adaptive Single-Site Time-Dependent Variational Principle. JACS AU 2022; 2:335-340. [PMID: 35252984 PMCID: PMC8889605 DOI: 10.1021/jacsau.1c00474] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/25/2021] [Indexed: 06/14/2023]
Abstract
In recent years, the time-dependent variational principle (TDVP) method based on the matrix product state (MPS) wave function formulation has shown its great power in performing large-scale quantum dynamics simulations for realistic chemical systems with strong electron-vibration interactions. In this work, we propose a stochastic adaptive single-site TDVP (SA-1TDVP) scheme to evolve the bond-dimension adaptively, which can integrate the traditional advantages of both the high efficiency of the single-site TDVP (1TDVP) variant and the high accuracy of the two-site TDVP (2TDVP) variant. Based on the assumption that the level statistics of entanglement Hamiltonians, which originate from the reduced density matrices of the MPS method, follows a Poisson or Wigner distribution, as generically predicted by random-matrix theory, additional random singular values are generated to expand the bond-dimension automatically. Tests on simulating the vibrationally resolved quantum dynamics and absorption spectra in the pyrazine molecule and perylene bisimide (PBI) J-aggregate trimer as well as a spin-1/2 Heisenberg chain show that it can be automatic and as accurate as 2TDVP but reduce the computational time remarkably.
Collapse
Affiliation(s)
- Yihe Xu
- School
of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle
Emissions Control, Nanjing University, Nanjing 210023, China
| | - Zhaoxuan Xie
- School
of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle
Emissions Control, Nanjing University, Nanjing 210023, China
| | - Xiaoyu Xie
- Department
of Chemistry, University of Liverpool, Liverpool L69 3BX, United Kingdom
| | - Ulrich Schollwöck
- Arnold
Sommerfeld Center of Theoretical Physics, Department of Physics, University of Munich, Theresienstrasse 37, 80333 Munich, Germany
- Munich
Center for Quantum Science and Technology (MCQST), 80799 Munich, Germany
| | - Haibo Ma
- School
of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle
Emissions Control, Nanjing University, Nanjing 210023, China
| |
Collapse
|
28
|
Feldmann R, Muolo A, Baiardi A, Reiher M. Quantum Proton Effects from Density Matrix Renormalization Group Calculations. J Chem Theory Comput 2022; 18:234-250. [PMID: 34978441 DOI: 10.1021/acs.jctc.1c00913] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
We recently introduced [J. Chem. Phys. 2020, 152, 204103] the nuclear-electronic all-particle density matrix renormalization group (NEAP-DMRG) method to solve the molecular Schrödinger equation, based on a stochastically optimized orbital basis, without invoking the Born-Oppenheimer approximation. In this work, we combine the DMRG method with the nuclear-electronic Hartree-Fock (NEHF-DMRG) approach, treating nuclei and electrons on the same footing. Inter- and intraspecies correlations are described within the DMRG method without truncating the excitation degree of the full configuration interaction wave function. We extend the concept of orbital entanglement and mutual information to nuclear-electronic wave functions and demonstrate that they are reliable metrics to detect strong correlation effects. We apply the NEHF-DMRG method to the HeHHe+ molecular ion, to obtain accurate proton densities, ground-state total energies, and vibrational transition frequencies by comparison with state-of-the-art data obtained with grid-based approaches and modern configuration interaction methods. For HCN, we improve on the accuracy of the latter approaches with respect to both the ground-state absolute energy and proton density, which is a major challenge for multireference nuclear-electronic state-of-the-art methods.
Collapse
Affiliation(s)
- Robin Feldmann
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Andrea Muolo
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, 91904 Jerusalem, Israel
| | - Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
29
|
Baiardi A, Grimmel SA, Steiner M, Türtscher PL, Unsleber JP, Weymuth T, Reiher M. Expansive Quantum Mechanical Exploration of Chemical Reaction Paths. Acc Chem Res 2022; 55:35-43. [PMID: 34918903 DOI: 10.1021/acs.accounts.1c00472] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Quantum mechanical methods have been well-established for the elucidation of reaction paths of chemical processes and for the explicit dynamics of molecular systems. While they are usually deployed in routine manual calculations on reactions for which some insights are already available (typically from experiment), new algorithms and continuously increasing capabilities of modern computer hardware allow for exploratory open-ended computational campaigns that are unbiased and therefore enable unexpected discoveries. Highly efficient and even automated procedures facilitate systematic approaches toward the exploration of uncharted territory in molecular transformations and dynamics. In this work, we elaborate on such explorative approaches that range from reaction network explorations with (stationary) quantum chemical methods to explorative molecular dynamics and migrant wave packet dynamics. The focus is on recent developments that cover the following strategies. (i) Pruning search options for elementary reaction steps by heuristic rules based on the first-principles of quantum mechanics: Rules are required for reducing the combinatorial explosion of potentially reactive atom pairings, and rooting them in concepts derived from the electronic wave function makes them applicable to any molecular system. (ii) Enforcing reactive events by external biases: Inducing a reaction requires constraints that steer and direct elementary-step searches, which can be formulated in terms of forces, velocities, or supplementary potentials. (iii) Manual steering facilitated by interactive quantum mechanics: As ultrafast quantum chemical methods allow for real-time manual interactions with molecular systems, human-intuition-guided paths can be easily explored with suitable human-machine interfaces. (iv) New approaches for transition-state optimization with continuous curve representations can provide stable schemes to be driven in an automated way by allowing for an efficient tuning of the curve's parameters (instead of a manipulation of a collection of structures along the path), and (v) reactive molecular dynamics and direct wave packet propagation exploit the equations of motion of an underlying mechanical theory (usually, classical Newtonian mechanics or Schrödinger quantum mechanics). Explorative approaches are likely to replace the current state of the art in computational chemistry, because they reduce the human effort to be invested in reaction path elucidations, they are less prone to errors and bias-free, and they cover more extensive regions of the relevant configuration space. As a result, computational investigations that rely on these techniques are more likely to deliver surprising discoveries.
Collapse
Affiliation(s)
- Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stephanie A. Grimmel
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Miguel Steiner
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Paul L. Türtscher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan P. Unsleber
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Thomas Weymuth
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
30
|
Kochetov V, Bokarev SI. RhoDyn: A ρ-TD-RASCI Framework to Study Ultrafast Electron Dynamics in Molecules. J Chem Theory Comput 2021; 18:46-58. [PMID: 34965135 DOI: 10.1021/acs.jctc.1c01097] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
This article presents the program module RhoDyn as part of the OpenMOLCAS project intended to study ultrafast electron dynamics within the density-matrix-based time-dependent restricted active space configuration interaction framework (ρ-TD-RASCI). The formalism allows for the treatment of spin-orbit coupling effects, accounts for nuclear vibrations in the form of a vibrational heat bath, and naturally incorporates (auto)ionization effects. Apart from describing the theory behind and the program workflow, the paper also contains examples of its application to the simulations of the linear L2,3 absorption spectra of a titanium complex, high harmonic generation in the hydrogen molecule, ultrafast charge migration in benzene and iodoacetylene, and spin-flip dynamics in the core excited states of iron complexes.
Collapse
Affiliation(s)
- Vladislav Kochetov
- Institut für Physik, Universität Rostock, A.-Einstein-Strasse 23-24, 18059 Rostock, Germany
| | - Sergey I Bokarev
- Institut für Physik, Universität Rostock, A.-Einstein-Strasse 23-24, 18059 Rostock, Germany
| |
Collapse
|
31
|
Abstract
We introduce DMRG[FEAST], a new method for optimizing excited-state many-body wave functions with the density matrix renormalization group (DMRG) algorithm. Our approach applies the FEAST algorithm, originally designed for large-scale diagonalization problems, to matrix product state wave functions. We show that DMRG[FEAST] enables the stable optimization of both low- and high-energy eigenstates, therefore overcoming the limitations of state-of-the-art excited-state DMRG algorithms. We demonstrate the reliability of DMRG[FEAST] by calculating anharmonic vibrational excitation energies of molecules with up to 30 fully coupled degrees of freedom.
Collapse
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Anna Klára Kelemen
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
32
|
Soley MB, Bergold P, Gorodetsky AA, Batista VS. Functional Tensor-Train Chebyshev Method for Multidimensional Quantum Dynamics Simulations. J Chem Theory Comput 2021; 18:25-36. [PMID: 34898201 DOI: 10.1021/acs.jctc.1c00941] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Methods for efficient simulations of multidimensional quantum dynamics are essential for theoretical studies of chemical systems where quantum effects are important, such as those involving rearrangements of protons or electronic configurations. Here, we introduce the functional tensor-train Chebyshev (FTTC) method for rigorous nuclear quantum dynamics simulations. FTTC is essentially the Chebyshev propagation scheme applied to the initial state represented in a continuous analogue tensor-train format. We demonstrate the capabilities of FTTC as applied to simulations of proton quantum dynamics in a 50-dimensional model of hydrogen-bonded DNA base pairs.
Collapse
Affiliation(s)
- Micheline B Soley
- Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States.,Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520, United States
| | - Paul Bergold
- Zentrum Mathematik, Technical University of Munich, Boltzmannstr. 3, 85748 Garching, Germany
| | - Alex A Gorodetsky
- Department of Aerospace Engineering, University of Michigan, 1320 Beal Avenue, Ann Arbor, Michigan 48109-2140, United States
| | - Victor S Batista
- Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States.,Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520, United States.,Energy Sciences Institute, Yale University, P.O. Box 27394, West Haven, Connecticut 06516-7394, United States
| |
Collapse
|
33
|
Li W, Ma H, Li S, Ma J. Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning. Chem Sci 2021; 12:14987-15006. [PMID: 34909141 PMCID: PMC8612375 DOI: 10.1039/d1sc02574k] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2021] [Accepted: 10/12/2021] [Indexed: 12/11/2022] Open
Abstract
Electronic structure methods based on quantum mechanics (QM) are widely employed in the computational predictions of the molecular properties and optoelectronic properties of molecular materials. The computational costs of these QM methods, ranging from density functional theory (DFT) or time-dependent DFT (TDDFT) to wave-function theory (WFT), usually increase sharply with the system size, causing the curse of dimensionality and hindering the QM calculations for large sized systems such as long polymer oligomers and complex molecular aggregates. In such cases, in recent years low scaling QM methods and machine learning (ML) techniques have been adopted to reduce the computational costs and thus assist computational and data driven molecular material design. In this review, we illustrated low scaling ground-state and excited-state QM approaches and their applications to long oligomers, self-assembled supramolecular complexes, stimuli-responsive materials, mechanically interlocked molecules, and excited state processes in molecular aggregates. Variable electrostatic parameters were also introduced in the modified force fields with the polarization model. On the basis of QM computational or experimental datasets, several ML algorithms, including explainable models, deep learning, and on-line learning methods, have been employed to predict the molecular energies, forces, electronic structure properties, and optical or electrical properties of materials. It can be conceived that low scaling algorithms with periodic boundary conditions are expected to be further applicable to functional materials, perhaps in combination with machine learning to fast predict the lattice energy, crystal structures, and spectroscopic properties of periodic functional materials.
Collapse
Affiliation(s)
- Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| | - Haibo Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
- Jiangsu Key Laboratory of Advanced Organic Materials, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University Nanjing 210023 China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| | - Jing Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
- Jiangsu Key Laboratory of Advanced Organic Materials, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University Nanjing 210023 China
| |
Collapse
|
34
|
Gelin MF, Velardo A, Borrelli R. Efficient quantum dynamics simulations of complex molecular systems: A unified treatment of dynamic and static disorder. J Chem Phys 2021; 155:134102. [PMID: 34624969 DOI: 10.1063/5.0065896] [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/05/2023] Open
Abstract
We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su-Schrieffer-Heeger (Holstein-Peierls) model Hamiltonian.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | | | | |
Collapse
|
35
|
Jiang T, Ren J, Shuai Z. Chebyshev Matrix Product States with Canonical Orthogonalization for Spectral Functions of Many-Body Systems. J Phys Chem Lett 2021; 12:9344-9352. [PMID: 34549961 DOI: 10.1021/acs.jpclett.1c02688] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We propose a method to calculate the spectral functions of many-body systems by Chebyshev expansion in the framework of matrix product states coupled with canonical orthogonalization (coCheMPS). The canonical orthogonalization can improve the accuracy and efficiency significantly because the orthogonalized Chebyshev vectors can provide an ideal basis for constructing the effective Hamiltonian in which the exact recurrence relation can be retained. In addition, not only the spectral function but also the excited states and eigenenergies can be directly calculated, which is usually impossible for other MPS-based methods such as time-dependent formalism or correction vector. The remarkable accuracy and efficiency of coCheMPS over other methods are demonstrated by calculating the spectral functions of spin chain and ab initio hydrogen chain. For the first time we demonstrate that Chebyshev MPS can be used to deal with ab initio electronic Hamiltonian effectively. We emphasize the strength of coCheMPS to calculate the low excited states of systems with sparse discrete spectrum. We also caution the application for electron-phonon systems with dense density of states.
Collapse
Affiliation(s)
- Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
36
|
Gelin MF, Borrelli R. Simulation of Nonlinear Femtosecond Signals at Finite Temperature via a Thermo Field Dynamics-Tensor Train Method: General Theory and Application to Time- and Frequency-Resolved Fluorescence of the Fenna-Matthews-Olson Complex. J Chem Theory Comput 2021; 17:4316-4331. [PMID: 34076412 DOI: 10.1021/acs.jctc.1c00158] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | | |
Collapse
|
37
|
Baiardi A. Electron Dynamics with the Time-Dependent Density Matrix Renormalization Group. J Chem Theory Comput 2021; 17:3320-3334. [PMID: 34043347 DOI: 10.1021/acs.jctc.0c01048] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
In this work, we simulate the electron dynamics in molecular systems with the time-dependent density matrix renormalization group (TD-DMRG) algorithm. We leverage the generality of the so-called tangent-space TD-DMRG formulation and design a computational framework in which the dynamics is driven by the exact nonrelativistic electronic Hamiltonian. We show that by parametrizing the wave function as a matrix product state, we can accurately simulate the dynamics of systems including up to 20 electrons and 32 orbitals. We apply the TD-DMRG algorithm to three problems that are hardly targeted by time-independent methods: the calculation of molecular (hyper)polarizabilities, the simulation of electronic absorption spectra, and the study of ultrafast ionization dynamics.
Collapse
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, Zürich 8093, Switzerland
| |
Collapse
|
38
|
Wang Y, Ren J, Shuai Z. Evaluating the anharmonicity contributions to the molecular excited state internal conversion rates with finite temperature TD-DMRG. J Chem Phys 2021; 154:214109. [PMID: 34240969 DOI: 10.1063/5.0052804] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we propose a new method to calculate molecular nonradiative electronic relaxation rates based on the numerically exact time-dependent density matrix renormalization group theory. This method could go beyond the existing frameworks under the harmonic approximation (HA) of the potential energy surface (PES) so that the anharmonic effect could be considered, which is of vital importance when the electronic energy gap is much larger than the vibrational frequency. We calculate the internal conversion (IC) rates in a two-mode model with Morse potential to investigate the validity of HA. We find that HA is unsatisfactory unless only the lowest several vibrational states of the lower electronic state are involved in the transition process when the adiabatic excitation energy is relatively low. As the excitation energy increases, HA first underestimates and then overestimates the IC rates when the excited state PES shifts toward the dissociative side of the ground state PES. On the contrary, HA slightly overestimates the IC rates when the excited state PES shifts toward the repulsive side. In both cases, a higher temperature enlarges the error of HA. As a real example to demonstrate the effectiveness and scalability of the method, we calculate the IC rates of azulene from S1 to S0 on the ab initio anharmonic PES approximated by the one-mode representation. The calculated IC rates of azulene under HA are consistent with the analytically exact results. The rates on the anharmonic PES are 30%-40% higher than the rates under HA.
Collapse
Affiliation(s)
- Yuanheng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
39
|
Yan Y, Xu M, Li T, Shi Q. Efficient propagation of the hierarchical equations of motion using the Tucker and hierarchical Tucker tensors. J Chem Phys 2021; 154:194104. [PMID: 34240893 DOI: 10.1063/5.0050720] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We develop new methods to efficiently propagate the hierarchical equations of motion (HEOM) by using the Tucker and hierarchical Tucker (HT) tensors to represent the reduced density operator and auxiliary density operators. We first show that by employing the split operator method, the specific structure of the HEOM allows a simple propagation scheme using the Tucker tensor. When the number of effective modes in the HEOM increases and the Tucker representation becomes intractable, the split operator method is extended to the binary tree structure of the HT representation. It is found that to update the binary tree nodes related to a specific effective mode, we only need to propagate a short matrix product state constructed from these nodes. Numerical results show that by further employing the mode combination technique commonly used in the multi-configuration time-dependent Hartree approaches, the binary tree representation can be applied to study excitation energy transfer dynamics in a fairly large system including over 104 effective modes. The new methods may thus provide a promising tool in simulating quantum dynamics in condensed phases.
Collapse
Affiliation(s)
- Yaming Yan
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Meng Xu
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Tianchu Li
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| |
Collapse
|
40
|
Peng J, Xie Y, Hu D, Lan Z. Analysis of bath motion in MM-SQC dynamics via dimensionality reduction approach: Principal component analysis. J Chem Phys 2021; 154:094122. [PMID: 33685149 DOI: 10.1063/5.0039743] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The system-plus-bath model is an important tool to understand the nonadiabatic dynamics of large molecular systems. Understanding the collective motion of a large number of bath modes is essential for revealing their key roles in the overall dynamics. Here, we applied principal component analysis (PCA) to investigate the bath motion in the basis of a large dataset generated from the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian nonadiabatic dynamics for the excited-state energy transfer in the Frenkel-exciton model. The PCA method clearly elucidated that two types of bath modes, which either display strong vibronic coupling or have frequencies close to that of the electronic transition, are important to the nonadiabatic dynamics. These observations were fully consistent with the physical insights. The conclusions were based on the PCA of the trajectory data and did not involve significant pre-defined physical knowledge. The results show that the PCA approach, which is one of the simplest unsupervised machine learning dimensionality reduction methods, is a powerful one for analyzing complicated nonadiabatic dynamics in the condensed phase with many degrees of freedom.
Collapse
Affiliation(s)
- Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
41
|
Madsen NK, Jensen RB, Christiansen O. Calculating vibrational excitation energies using tensor-decomposed vibrational coupled-cluster response theory. J Chem Phys 2021; 154:054113. [DOI: 10.1063/5.0037240] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Niels Kristian Madsen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Rasmus Berg Jensen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| |
Collapse
|
42
|
Yan Y, Xing T, Shi Q. A new method to improve the numerical stability of the hierarchical equations of motion for discrete harmonic oscillator modes. J Chem Phys 2020; 153:204109. [DOI: 10.1063/5.0027962] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yaming Yan
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Tao Xing
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| |
Collapse
|
43
|
Madsen NK, Hansen MB, Christiansen O, Zoccante A. Time-dependent vibrational coupled cluster with variationally optimized time-dependent basis sets. J Chem Phys 2020; 153:174108. [DOI: 10.1063/5.0024428] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Niels Kristian Madsen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Mads Bøttger Hansen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Alberto Zoccante
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
44
|
Abstract
We introduce the transcorrelated Density Matrix Renormalization Group (tcDMRG) theory for the efficient approximation of the energy for strongly correlated systems. tcDMRG encodes the wave function as a product of a fixed Jastrow or Gutzwiller correlator and a matrix product state. The latter is optimized by applying the imaginary-time variant of time-dependent (TD) DMRG to the non-Hermitian transcorrelated Hamiltonian. We demonstrate the efficiency of tcDMRG with the example of the two-dimensional Fermi-Hubbard Hamiltonian, a notoriously difficult target for the DMRG algorithm, for different sizes, occupation numbers, and interaction strengths. We demonstrate fast energy convergence of tcDMRG, which indicates that tcDMRG could increase the efficiency of standard DMRG beyond quasi-monodimensional systems and provides a generally powerful approach toward the dynamic correlation problem of DMRG.
Collapse
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
45
|
Sasmal S, Vendrell O. Non-adiabatic quantum dynamics without potential energy surfaces based on second-quantized electrons: Application within the framework of the MCTDH method. J Chem Phys 2020; 153:154110. [PMID: 33092359 DOI: 10.1063/5.0028116] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
A first principles quantum formalism to describe the non-adiabatic dynamics of electrons and nuclei based on a second quantization representation (SQR) of the electronic motion combined with the usual representation of the nuclear coordinates is introduced. This procedure circumvents the introduction of potential energy surfaces and non-adiabatic couplings, providing an alternative to the Born-Oppenheimer approximation. An important feature of the molecular Hamiltonian in the mixed first quantized representation for the nuclei and the SQR representation for the electrons is that all degrees of freedom, nuclear positions and electronic occupations, are distinguishable. This makes the approach compatible with various tensor decomposition Ansätze for the propagation of the nuclear-electronic wavefunction. Here, we describe the application of this formalism within the multi-configuration time-dependent Hartree framework and its multilayer generalization, corresponding to Tucker and hierarchical Tucker tensor decompositions of the wavefunction, respectively. The approach is applied to the calculation of the photodissociation cross section of the HeH+ molecule under extreme ultraviolet irradiation, which features non-adiabatic effects and quantum interferences between the two possible fragmentation channels, He + H+ and He+ + H. These calculations are compared with the usual description based on ab initio potential energy surfaces and non-adiabatic coupling matrix elements, which fully agree. The proof-of-principle calculations serve to illustrate the advantages and drawbacks of this formalism, which are discussed in detail, as well as possible ways to overcome them. We close with an outlook of possible application domains where the formalism might outperform the usual approach, for example, in situations that combine a strong static correlation of the electrons with non-adiabatic electronic-nuclear effects.
Collapse
Affiliation(s)
- Sudip Sasmal
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuneheimer Feld 229, 69120 Heidelberg, Germany
| | - Oriol Vendrell
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuneheimer Feld 229, 69120 Heidelberg, Germany
| |
Collapse
|
46
|
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
|
47
|
Ren J, Li W, Jiang T, Shuai Z. A general automatic method for optimal construction of matrix product operators using bipartite graph theory. J Chem Phys 2020; 153:084118. [PMID: 32872857 DOI: 10.1063/5.0018149] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Constructing matrix product operators (MPOs) is at the core of the modern density matrix renormalization group (DMRG) and its time dependent formulation. For the DMRG to be conveniently used in different problems described by different Hamiltonians, in this work, we propose a new generic algorithm to construct the MPO of an arbitrary operator with a sum-of-products form based on the bipartite graph theory. We show that the method has the following advantages: (i) it is automatic in that only the definition of the operator is required; (ii) it is symbolic thus free of any numerical error; (iii) the complementary operator technique can be fully employed so that the resulting MPO is globally optimal for any given order of degrees of freedom; and (iv) the symmetry of the system could be fully employed to reduce the dimension of MPO. To demonstrate the effectiveness of the new algorithm, the MPOs of Hamiltonians ranging from the prototypical spin-boson model and the Holstein model to the more complicated ab initio electronic Hamiltonian and the anharmonic vibrational Hamiltonian with the sextic force field are constructed. It is found that for the former three cases, our automatic algorithm can reproduce exactly the same MPOs as the optimally hand-crafted ones already known in the literature.
Collapse
Affiliation(s)
- Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
48
|
Li W, Ren J, Shuai Z. Finite-Temperature TD-DMRG for the Carrier Mobility of Organic Semiconductors. J Phys Chem Lett 2020; 11:4930-4936. [PMID: 32492339 DOI: 10.1021/acs.jpclett.0c01072] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
A large number of nonadiabatic dynamical studies have been applied to reveal the nature of carrier transport in organic semiconductors with different approximations. We present here a "nearly exact" graphical-process-unit-based finite-temperature time-dependent density matrix renormalization group (TD-DMRG) method to evaluate the carrier mobility in organic semiconductors, as described by the electron-phonon model, in particular, in rubrene crystal, one of the prototypical organic semiconductors, with parameters derived from first-principles. We find that (i) TD-DMRG is a general and robust method that can bridge the gap between hopping and band pictures, covering a wide range of electronic coupling strengths and (ii) with realistic parameters, TD-DMRG is able to account for the experimentally observed "band-like" transport behavior (∂μ/∂T < 0) in rubrene. We further study the long-standing puzzle of the isotope effect for charge transport and unambiguously demonstrate that the negative isotope effect (∂μ/∂m < 0 where m is the atomic mass) should be universal.
Collapse
Affiliation(s)
- Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
49
|
Kowalski K, Bauman NP. Sub-system quantum dynamics using coupled cluster downfolding techniques. J Chem Phys 2020; 152:244127. [DOI: 10.1063/5.0008436] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Nicholas P. Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| |
Collapse
|
50
|
Aquilante F, Autschbach J, Baiardi A, Battaglia S, Borin VA, Chibotaru LF, Conti I, De Vico L, Delcey M, Fdez Galván I, Ferré N, Freitag L, Garavelli M, Gong X, Knecht S, Larsson ED, Lindh R, Lundberg M, Malmqvist PÅ, Nenov A, Norell J, Odelius M, Olivucci M, Pedersen TB, Pedraza-González L, Phung QM, Pierloot K, Reiher M, Schapiro I, Segarra-Martí J, Segatta F, Seijo L, Sen S, Sergentu DC, Stein CJ, Ungur L, Vacher M, Valentini A, Veryazov V. Modern quantum chemistry with [Open]Molcas. J Chem Phys 2020; 152:214117. [PMID: 32505150 DOI: 10.1063/5.0004835] [Citation(s) in RCA: 226] [Impact Index Per Article: 56.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.
Collapse
Affiliation(s)
- 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, Buffalo, New York 14260-3000, USA
| | - Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stefano Battaglia
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Veniamin A Borin
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Liviu F Chibotaru
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Irene Conti
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100 Siena, Italy
| | - Mickaël Delcey
- Department of Chemistry - Ångström Laboratory, Uppsala University, SE-751 21 Uppsala, Sweden
| | - Ignacio Fdez Galván
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Nicolas Ferré
- Aix-Marseille University, CNRS, Institut Chimie Radicalaire, Marseille, France
| | - Leon Freitag
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Xuejun Gong
- Department of Chemistry, University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Stefan Knecht
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Ernst D Larsson
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, Lund 22100, Sweden
| | - Roland Lindh
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Marcus Lundberg
- Department of Chemistry - Ångström Laboratory, Uppsala University, SE-751 21 Uppsala, Sweden
| | - Per Åke Malmqvist
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, Lund 22100, Sweden
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Jesper Norell
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Michael Odelius
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100 Siena, Italy
| | - Thomas B Pedersen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Laura Pedraza-González
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100 Siena, Italy
| | - Quan M Phung
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Chikusa, Nagoya 464-8602, Japan
| | - Kristine Pierloot
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Javier Segarra-Martí
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, London W12 0BZ, United Kingdom
| | - Francesco Segatta
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Luis Seijo
- Departamento de Química, Instituto Universitario de Ciencia de Materiales Nicolás Cabrera, and Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Saumik Sen
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | | | - Christopher J Stein
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Liviu Ungur
- Department of Chemistry, University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Morgane Vacher
- Laboratoire CEISAM - UMR CNRS 6230, Université de Nantes, 44300 Nantes, France
| | - Alessio Valentini
- Theoretical Physical Chemistry, Research Unit MolSys, Université de Liège, Allée du 6 Août, 11, 4000 Liège, Belgium
| | - Valera Veryazov
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, Lund 22100, Sweden
| |
Collapse
|