1
|
Malpathak S, Church MS, Ananth N. A Semiclassical Framework for Mixed Quantum Classical Dynamics. J Phys Chem A 2022; 126:6359-6375. [PMID: 36070472 DOI: 10.1021/acs.jpca.2c03467] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Semiclassical (SC) approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level quantum mechanical description is only required for a handful of modes, while the majority of environment modes that do not play a key role in the reactive event of interest are well served with a lower level of theory. In this feature, we introduce modified Filinov filtration as a powerful tool to construct mixed quantum-classical SC theories where different subsystems can be quantized to different extents without introducing ad hoc intersubsystem interaction terms. We demonstrate that these Filinov-based SC methods can systematically tune between quantum and classical limit SC behavior, offering a practical way forward to accurate and computationally efficient simulations of high-dimensional quantum processes.
Collapse
Affiliation(s)
- Shreyas Malpathak
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| | - Matthew S Church
- Department of Chemistry, Brown University, Providence, Rhode Island 02906, United States
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| |
Collapse
|
2
|
Buchholz M, Grossmann F, Ceotto M. Simplified approach to the mixed time-averaging semiclassical initial value representation for the calculation of dense vibrational spectra. J Chem Phys 2018; 148:114107. [DOI: 10.1063/1.5020144] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Max Buchholz
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Frank Grossmann
- Institut für Theoretische Physik, Technische Universität Dresden, 01062 Dresden, Germany
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
3
|
Karsten S, Ivanov SD, Bokarev SI, Kühn O. Quasi-classical approaches to vibronic spectra revisited. J Chem Phys 2018; 148:102337. [DOI: 10.1063/1.5011764] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Sven Karsten
- Institute of Physics, Rostock University, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Sergei D. Ivanov
- Institute of Physics, Rostock University, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Sergey I. Bokarev
- Institute of Physics, Rostock University, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Oliver Kühn
- Institute of Physics, Rostock University, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| |
Collapse
|
4
|
Ceotto M, Di Liberto G, Conte R. Semiclassical "Divide-and-Conquer" Method for Spectroscopic Calculations of High Dimensional Molecular Systems. PHYSICAL REVIEW LETTERS 2017; 119:010401. [PMID: 28731742 DOI: 10.1103/physrevlett.119.010401] [Citation(s) in RCA: 47] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/16/2017] [Indexed: 05/11/2023]
Abstract
A new semiclassical "divide-and-conquer" method is presented with the aim of demonstrating that quantum dynamics simulations of high dimensional molecular systems are doable. The method is first tested by calculating the quantum vibrational power spectra of water, methane, and benzene-three molecules of increasing dimensionality for which benchmark quantum results are available-and then applied to C_{60}, a system characterized by 174 vibrational degrees of freedom. Results show that the approach can accurately account for quantum anharmonicities, purely quantum features like overtones, and the removal of degeneracy when the molecular symmetry is broken.
Collapse
Affiliation(s)
- Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - Giovanni Di Liberto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
5
|
Gabas F, Conte R, Ceotto M. On-the-Fly ab Initio Semiclassical Calculation of Glycine Vibrational Spectrum. J Chem Theory Comput 2017; 13:2378-2388. [PMID: 28489368 PMCID: PMC5472367 DOI: 10.1021/acs.jctc.6b01018] [Citation(s) in RCA: 49] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We
present an on-the-fly ab initio semiclassical study of vibrational
energy levels of glycine, calculated by Fourier transform of the wavepacket
correlation function. It is based on a multiple coherent states approach
integrated with monodromy matrix regularization for chaotic dynamics.
All four lowest-energy glycine conformers are investigated by means
of single-trajectory semiclassical spectra obtained upon classical
evolution of on-the-fly trajectories with harmonic zero-point energy.
For the most stable conformer I, direct dynamics trajectories are
also run for each vibrational mode with energy equal to the first
harmonic excitation. An analysis of trajectories evolved up to 50 000
atomic time units demonstrates that, in this time span, conformers
II and III can be considered as isolated species, while conformers
I and IV show a pretty facile interconversion. Therefore, previous
perturbative studies based on the assumption of isolated conformers
are often reliable but might be not completely appropriate in the
case of conformer IV and conformer I for which interconversion occurs
promptly.
Collapse
Affiliation(s)
- Fabio Gabas
- Dipartimento di Chimica, Università degli Studi di Milano , via Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano , via Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano , via Golgi 19, 20133 Milano, Italy
| |
Collapse
|
6
|
Beutier J, Vuilleumier R, Bonella S, Ciccotti G. Gas phase infrared spectra from quasi-classical Kubo time correlation functions. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1064550] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- J. Beutier
- Département de Chimie, École Normale Supérieure-PSL Research University , Paris, France
- Département de Chimie, Sorbonne Universités , UPMC Univ Paris 06, Paris, France
- Département de Chimie, CNRS , UMR 8640 PASTEUR, Paris, France
- Dipartimento di Fisica, Università di Roma ‘La Sapienza’ , Roma, Italy
| | - R. Vuilleumier
- Département de Chimie, École Normale Supérieure-PSL Research University , Paris, France
- Département de Chimie, Sorbonne Universités , UPMC Univ Paris 06, Paris, France
- Département de Chimie, CNRS , UMR 8640 PASTEUR, Paris, France
| | - S. Bonella
- Department of Chemistry, CECAM Centre Européen de Calcul Atomique et Moléculaire, École Polytechnique Fédérale de Lausanne , Batochime, Lausanne, Switzerland
| | - G. Ciccotti
- Dipartimento di Fisica, Università di Roma ‘La Sapienza’ , Roma, Italy
- School of Physics, University College of Dublin UCD , Dublin 4, Ireland
| |
Collapse
|
7
|
Song L, Shi Q. A new approach to calculate charge carrier transport mobility in organic molecular crystals from imaginary time path integral simulations. J Chem Phys 2015; 142:174103. [DOI: 10.1063/1.4919061] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Linze Song
- 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
| | - 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
| |
Collapse
|
8
|
Beutier J, Borgis D, Vuilleumier R, Bonella S. Computing thermal Wigner densities with the phase integration method. J Chem Phys 2014; 141:084102. [DOI: 10.1063/1.4892597] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
9
|
|
10
|
Beutier J, Monteferrante M, Bonella S, Vuilleumier R, Ciccotti G. Gas phase infrared spectra via the phase integration quasi-classical method. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.843776] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- J. Beutier
- Département de Chimie, UMR 8640 CNRS-ENS-UPMC, École Normale Supérieure, 24 rue Lhomond, 75005, Paris, France
| | - M. Monteferrante
- Dipartimento di Fisica and CNISM Unità 1, Università di Roma ‘La Sapienza’, Ple A. Moro 2, 00185, Rome, Italy
| | - S. Bonella
- Dipartimento di Fisica and CNISM Unità 1, Università di Roma ‘La Sapienza’, Ple A. Moro 2, 00185, Rome, Italy
| | - R. Vuilleumier
- Département de Chimie, UMR 8640 CNRS-ENS-UPMC, École Normale Supérieure, 24 rue Lhomond, 75005, Paris, France
| | - G. Ciccotti
- Dipartimento di Fisica and CNISM Unità 1, Università di Roma ‘La Sapienza’, Ple A. Moro 2, 00185, Rome, Italy
| |
Collapse
|
11
|
Borgis D, Assaraf R, Rotenberg B, Vuilleumier R. Computation of pair distribution functions and three-dimensional densities with a reduced variance principle. Mol Phys 2013. [DOI: 10.1080/00268976.2013.838316] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|