1
|
Malpathak S, Ananth N. Semiclassical dynamics in Wigner phase space I: Adiabatic hybrid Wigner dynamics. J Chem Phys 2024; 161:094109. [PMID: 39234962 DOI: 10.1063/5.0223185] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2024] [Accepted: 08/12/2024] [Indexed: 09/06/2024] Open
Abstract
The Wigner phase space formulation of quantum mechanics is a complete framework for quantum dynamic calculations that elegantly highlights connections with classical dynamics. In this series of two articles, building upon previous efforts, we derive the full hierarchy of approximate semiclassical (SC) dynamic methods for adiabatic and non-adiabatic problems in Wigner phase space. In Paper I, focusing on adiabatic single surface processes, we derive the well-known double Herman-Kluk (DHK) approximation for real-time correlation functions in Wigner phase space and connect it to the linearized SC (LSC) approximation through a stationary phase approximation. We exploit this relationship to introduce a new hybrid SC method, termed Adiabatic Hybrid Wigner Dynamics (AHWD) that allows for a few important "system" degrees of freedom (dofs) to be treated at the DHK level, while treating the rest of the dofs (the "bath") at the LSC level. AHWD is shown to accurately capture quantum interference effects in models of coupled oscillators and the decoherence of vibrational probability density of a model I2 Morse oscillator coupled to an Ohmic thermal bath. We show that AHWD significantly mitigates the sign problem and employs reduced dimensional prefactors bringing calculations of complex system-bath problems within the reach of SC methods. Paper II focuses on extending this hybrid SC dynamics to nonadiabatic processes.
Collapse
Affiliation(s)
- Shreyas Malpathak
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
2
|
Glaser N, Baiardi A, Lieberherr AZ, Reiher M. Vibrational Entanglement through the Lens of Quantum Information Measures. J Phys Chem Lett 2024; 15:6958-6965. [PMID: 38940632 PMCID: PMC11247497 DOI: 10.1021/acs.jpclett.4c01298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2024]
Abstract
We introduce a quantum information analysis of vibrational wave functions to understand complex vibrational spectra of molecules with strong anharmonic couplings and vibrational resonances. For this purpose, we define one- and two-modal entropies to guide the identification of strongly coupled vibrational modes and to characterize correlations within modal basis sets. We evaluate these descriptors for multiconfigurational vibrational wave functions which we calculate with the n-mode vibrational density matrix renormalization group algorithm. Based on the quantum information measures, we present a vibrational entanglement analysis of the vibrational ground and excited states of CO2, which display strong anharmonic effects due to the symmetry-induced and accidental (near-) degeneracies. We investigate the entanglement signature of the Fermi resonance and discuss the maximally entangled state arising from the two degenerate bending modes.
Collapse
Affiliation(s)
- Nina Glaser
- ETH Zürich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Alberto Baiardi
- ETH Zürich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Annina Z Lieberherr
- ETH Zürich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
3
|
Huang QR, Yano K, Yang Y, Fujii A, Kuo JL. Near-infrared spectroscopy of H 3O +⋯X n (X = Ar, N 2, and CO, n = 1-3). Phys Chem Chem Phys 2024; 26:10757-10768. [PMID: 38516880 DOI: 10.1039/d4cp00458b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2024]
Abstract
Near-infrared (NIR) spectra of H3O+⋯Xn (X = Ar, N2, and CO, n = 1-3) in the first overtone region of OH-stretching vibrations (4800-7000 cm-1) were measured. Not only OH-stretching overtones but also several combination bands are major features in this region, and assignments of these observed bands are not obvious at a glance. High-precision anharmonic vibrational simulations based on the discrete variable representation approach were performed. The simulated spectra show good agreement with the observed ones and provide firm assignments of the observed bands, except in the case of X = CO, in which higher order vibrational mode couplings seem significant. This agreement demonstrates that the present system can be a benchmark for high precision anharmonic vibrational computations of NIR spectra. Band broadening in the observed spectra becomes remarkable with an increase of the interaction with the solvent molecule (X). The origin of the band broadening is explored by rare gas tagging experiments and anharmonic vibrational simulations of hot bands.
Collapse
Affiliation(s)
- Qian-Rui Huang
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, Taiwan.
| | - Kazuyoshi Yano
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai, Japan.
| | - Yaodi Yang
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai, Japan.
| | - Asuka Fujii
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai, Japan.
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, Taiwan.
| |
Collapse
|
4
|
Ellerbrock R, Johnson KG, Seritan S, Hoppe H, Zhang JH, Lenzen T, Weike T, Manthe U, Martínez TJ. QuTree: A tree tensor network package. J Chem Phys 2024; 160:112501. [PMID: 38497471 DOI: 10.1063/5.0180233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2023] [Accepted: 02/20/2024] [Indexed: 03/19/2024] Open
Abstract
We present QuTree, a C++ library for tree tensor network approaches. QuTree provides class structures for tensors, tensor trees, and related linear algebra functions that facilitate the fast development of tree tensor network approaches such as the multilayer multiconfigurational time-dependent Hartree approach or the density matrix renormalization group approach and its various extensions. We investigate the efficiency of relevant tensor and tensor network operations and show that the overhead for managing the network structure is negligible, even in cases with a million leaves and small tensors. QuTree focuses on providing simple, high-level routines while retaining easy access to the backend to facilitate novel developments. We demonstrate the capabilities of the package by computing the eigenstates of coupled harmonic oscillator Hamiltonians and performing random circuit simulations on a virtual quantum computer.
Collapse
Affiliation(s)
- Roman Ellerbrock
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - K Grace Johnson
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Stefan Seritan
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Hannes Hoppe
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - J H Zhang
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Tim Lenzen
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - Thomas Weike
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - Uwe Manthe
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| |
Collapse
|
5
|
Finney JM, McCoy AB. Correlations between the Structures and Spectra of Protonated Water Clusters. J Phys Chem A 2024; 128:868-879. [PMID: 38265889 DOI: 10.1021/acs.jpca.3c07338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2024]
Abstract
Badger's rule-like correlations between OH stretching frequencies and intensities and the OH bond length are used to develop a spectral mapping procedure for studies of pure and protonated water clusters. This approach utilizes the vibrationally averaged OH bond lengths, which were obtained from diffusion Monte Carlo simulations that were performed using the general potential developed by Yu and Bowman. Good agreement is achieved between the spectra obtained using this approach and previously reported spectra for H+(H2O)n clusters, with n = 3, 4, and 5, as well as their perdeuterated analogues. The analysis of the spectra obtained by this spectral mapping approach supports previous work that assigned the spectrum of H+(H2O)6 to a mixture of Eigen and Zundel-like structures. Analysis of the calculated spectra also suggests a reassignment of the frequency of one of the transitions that involves the OH stretching vibration of the OH bonds in the hydronium core in the Eigen-like structure of H+(H2O)6 from 1917 cm-1 to roughly 2100 cm-1. For D+(D2O)6, comparison of the measured spectrum to those obtained by using the spectral mapping approach suggests that the carrier of the measured spectrum is one or more of the isomers of D+(D2O)6 that contain a four-membered ring and two flanking water molecules. While there are several candidate structures, the two flanking water molecules most likely form a chain that is bound to the hydronium core.
Collapse
Affiliation(s)
- Jacob M Finney
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
6
|
Hoppe H, Manthe U. Eigenstate calculation in the state-averaged (multi-layer) multi-configurational time-dependent Hartree approach. J Chem Phys 2024; 160:034104. [PMID: 38230812 DOI: 10.1063/5.0188748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 12/15/2023] [Indexed: 01/18/2024] Open
Abstract
A new approach for the calculation of eigenstates with the state-averaged (multi-layer) multi-configurational time-dependent Hartree (MCTDH) approach is presented. The approach is inspired by the recent work of Larsson [J. Chem. Phys. 151, 204102 (2019)]. It employs local optimization of the basis sets at each node of the multi-layer MCTDH tree and successive downward and upward sweeps to obtain a globally converged result. At the top node, the Hamiltonian represented in the basis of the single-particle functions (SPFs) of the first layer is diagonalized. Here p wavefunctions corresponding to the p lowest eigenvalues are computed by a block Lanczos approach. At all other nodes, a non-linear operator consisting of the respective mean-field Hamiltonian matrix and a projector onto the space spanned by the respective SPFs is considered. Here, the eigenstate corresponding to the lowest eigenvalue is computed using a short iterative Lanczos scheme. Two different examples are studied to illustrate the new approach: the calculation of the vibrational states of methyl and acetonitrile. The calculations for methyl employ the single-layer MCTDH approach, a general potential energy surface, and the correlation discrete variable representation. A five-layer MCTDH representation and a sum of product-type Hamiltonian are used in the acetonitrile calculations. Very fast convergence and order of magnitude reductions in the numerical effort compared to the previously used block relaxation scheme are found. Furthermore, a detailed comparison with the results of Avila and Carrington [J. Chem. Phys. 134, 054126 (2011)] for acetonitrile highlights the potential problems of convergence tests for high-dimensional systems.
Collapse
Affiliation(s)
- Hannes Hoppe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
7
|
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
|
8
|
Tan JA, Barbarona RF, Kuo JL. Approaching the "Zundel" Limit: Tuning the Vibrational Coupling in N 2H +Ng, Ng = {He, Ne, Ar, Kr, Xe, and Rn}. J Phys Chem A 2023. [PMID: 37418837 DOI: 10.1021/acs.jpca.3c01956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2023]
Abstract
The diazenylium ion (N2H+) is a ubiquitous ion in dense molecular clouds. This ion is often used as a dense gas tracer in outer space. Most of the previous works on diazenylium ion have focused on the shared-proton stretch band, νH+. In this work, we have performed reduced-dimensional calculations to investigate the vibrational structure of N2H+Ng, Ng = {He, Ne, Ar, Kr, Xe, and Rn}. We demonstrate a few interesting things about this system. First, the vibrational coupling in N2H+ can be tuned to switch on interesting anharmonic effects such as Fermi resonance or combination bands by tagging it with different noble gases. Second, a comparison of the vibrational spectrum from N2H+He to N2H+Rn shows that the νH+ can be swept from an "Eigen-like" to a "Zundel-like" limiting case. Anharmonic calculations were performed using a multilevel approach, which utilized the MP2 and CCSD(T) levels of theories. Binding energies for the elimination of Ng in N2H+Ng are also reported.
Collapse
Affiliation(s)
- Jake A Tan
- Institute of Atomic and Molecular Science, Academia Sinica, Taipei 10617, Taiwan
| | - Rona F Barbarona
- Institute of Atomic and Molecular Science, Academia Sinica, Taipei 10617, Taiwan
- Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei 11529, Taiwan
- Department of Physics, National Central University, Taoyuan 32001, Taiwan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Science, Academia Sinica, Taipei 10617, Taiwan
- Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei 11529, Taiwan
| |
Collapse
|
9
|
Wang Y, Chen Z, Yang Y. Calculating Vibrational Excited State Absorptions with Excited State Constrained Minimized Energy Surfaces. J Phys Chem A 2023. [PMID: 37335973 DOI: 10.1021/acs.jpca.3c01420] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
The modeling and interpretation of vibrational spectra are crucial for studying reaction dynamics using vibrational spectroscopy. Most prior theoretical developments focused on describing fundamental vibrational transitions while fewer developments focused on vibrational excited state absorptions. In this study, we present a new method that uses excited state constrained minimized energy surfaces (CMESs) to describe vibrational excited state absorptions. The excited state CMESs are obtained similarly to the previous ground state CMES development in our group but with additional wave function orthogonality constraints. Using a series of model systems, including the harmonic oscillator, Morse potential, double-well potential, quartic potential, and two-dimensional anharmonic potential, we demonstrate that this new procedure provides good estimations of the transition frequencies for vibrational excited state absorptions. These results are significantly better than those obtained from harmonic approximations using conventional potential energy surfaces, demonstrating the promise of excited state CMES-based methods for calculating vibrational excited state absorptions in real systems.
Collapse
Affiliation(s)
- Yiwen Wang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
10
|
Penfold TJ, Eng J. Mind the GAP: quantifying the breakdown of the linear vibronic coupling Hamiltonian. Phys Chem Chem Phys 2023; 25:7195-7204. [PMID: 36820783 DOI: 10.1039/d2cp05576g] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
Excited state dynamics play a critical role across a broad range of scientific fields. Importantly, the highly non-equilibrium nature of the states generated by photoexcitation means that excited state simulations should usually include an accurate description of the coupled electronic-nuclear motion, which often requires solving the time-dependent Schrödinger equation (TDSE). One of the biggest challenges for these simulations is the requirement to calculate the PES over which the nuclei evolve. An effective approach for addressing this challenge is to use the approximate linear vibronic coupling (LVC) Hamiltonian, which enables a model potential to be parameterised using relatively few quantum chemistry calculations. However, this approach is only valid provided there are no large amplitude motions in the excited state dynamics. In this paper we introduce and deploy a metric, the global anharmonicity parameter (GAP), which can be used to assess the accuracy of an LVC potential. Following its derivation, we illustrate its utility by applying it to three molecules exhibiting different rigidity in their excited states.
Collapse
Affiliation(s)
- Thomas J Penfold
- Chemistry, School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK.
| | - Julien Eng
- Chemistry, School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK.
| |
Collapse
|
11
|
Mendive-Tapia D, Meyer HD, Vendrell O. Optimal Mode Combination in the Multiconfiguration Time-Dependent Hartree Method through Multivariate Statistics: Factor Analysis and Hierarchical Clustering. J Chem Theory Comput 2023; 19:1144-1156. [PMID: 36716214 DOI: 10.1021/acs.jctc.2c01089] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The multiconfiguration time-dependent Hartree (MCTDH) method and its multilayer extension (ML-MCTDH) are powerful algorithms for the efficient computation of nuclear quantum dynamics in high-dimensional systems. By providing time-dependent variational orbitals and an optimal choice of layered effective degrees of freedom, one is able to reduce the computational cost to an amenable number of configurations. However, choices related to selecting properly the mode grouping and tensor tree are strongly system dependent and, thus far, subjectively based on intuition and/or experience. Therefore, herein we detail a new protocol based on multivariate statistics─more specifically, factor analysis and hierarchical clustering─for a reliable and convenient guiding in the optimal design of such complex "system-of-systems" tensor-network decompositions. The advantages of employing the new algorithm and its applicability are tested on water and two floppy protonated water clusters with large amplitude motions.
Collapse
Affiliation(s)
- David Mendive-Tapia
- Theoretische Chemie, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120Heidelberg, Germany
| | - Hans-Dieter Meyer
- Theoretische Chemie, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120Heidelberg, Germany
| | - Oriol Vendrell
- Theoretische Chemie, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120Heidelberg, Germany
| |
Collapse
|
12
|
Tan JA, Takahashi K. Vibrational spectrum of a 1D oscillator: The quantum, the Wigner, and the classical ways. J CHIN CHEM SOC-TAIP 2022. [DOI: 10.1002/jccs.202200398] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Jake A. Tan
- Institute of Atomic and Molecular Sciences Academia Sinica Taipei Taiwan
| | - Kaito Takahashi
- Institute of Atomic and Molecular Sciences Academia Sinica Taipei Taiwan
| |
Collapse
|
13
|
Schröder M, Gatti F, Lauvergnat D, Meyer HD, Vendrell O. The coupling of the hydrated proton to its first solvation shell. Nat Commun 2022; 13:6170. [PMID: 36257946 PMCID: PMC9579203 DOI: 10.1038/s41467-022-33650-w] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2022] [Accepted: 09/28/2022] [Indexed: 11/30/2022] Open
Abstract
The Zundel ([Formula: see text]) and Eigen ([Formula: see text]) cations play an important role as intermediate structures for proton transfer processes in liquid water. In the gas phase they exhibit radically different infrared (IR) spectra. The question arises: is there a least common denominator structure that explains the IR spectra of both, the Zundel and Eigen cations, and hence of the solvated proton? Full dimensional quantum simulations of these protonated cations demonstrate that two dynamical water molecules and an excess proton constitute this fundamental subunit. Embedded in the static environment of the parent Eigen cation, this subunit reproduces the positions and broadenings of its main excess-proton bands. In isolation, its spectrum reverts to the well-known Zundel ion. Hence, the dynamics of this subunit polarized by an environment suffice to explain the spectral signatures and anharmonic couplings of the solvated proton in its first solvation shell.
Collapse
Affiliation(s)
- Markus Schröder
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120, Heidelberg, Germany.
| | - Fabien Gatti
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay UMR 8214, 91405, Orsay, France
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR 8000, 91405, Orsay, France
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120, Heidelberg, Germany
| | - Oriol Vendrell
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120, Heidelberg, Germany.
| |
Collapse
|
14
|
Ellerbrock R, Manthe U. A non-hierarchical correlation discrete variable representation. J Chem Phys 2022; 156:134107. [PMID: 35395891 DOI: 10.1063/5.0088509] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. The non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH3 demonstrate the accuracy of the non-hierarchical CDVR approach.
Collapse
Affiliation(s)
- Roman Ellerbrock
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
15
|
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
|
16
|
Lindoy LP, Kloss B, Reichman DR. Time evolution of ML-MCTDH wavefunctions. II. Application of the projector splitting integrator. J Chem Phys 2021; 155:174109. [PMID: 34742222 DOI: 10.1063/5.0070043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 106 bath modes and find that for these problems, the PSI requires roughly 3-4 orders of magnitude fewer Hamiltonian evaluations and 2-3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2-3/1-2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ∼1.3×109 parameters to obtain accurate dynamics.
Collapse
Affiliation(s)
- Lachlan P Lindoy
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Benedikt Kloss
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| |
Collapse
|
17
|
Lindoy LP, Kloss B, Reichman DR. Time evolution of ML-MCTDH wavefunctions. I. Gauge conditions, basis functions, and singularities. J Chem Phys 2021; 155:174108. [PMID: 34742180 DOI: 10.1063/5.0070042] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
We derive a family of equations-of-motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than those used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOM approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs.
Collapse
Affiliation(s)
- Lachlan P Lindoy
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Benedikt Kloss
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| |
Collapse
|
18
|
Weike T, Manthe U. The multi-configurational time-dependent Hartree approach in optimized second quantization: thermal ensembles and statistical sampling. Chem Phys 2021. [DOI: 10.1016/j.chemphys.2021.111413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
19
|
Hanson MD, Readnour JA, Hassanali AA, Corcelli SA. Coupled Local-Mode Approach for the Calculation of Vibrational Spectra: Application to Protonated Water Clusters. J Phys Chem Lett 2021; 12:9226-9232. [PMID: 34529914 DOI: 10.1021/acs.jpclett.1c02254] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Spectroscopic studies of protonated water clusters (PWCs) have yielded enormous insights into the fundamental nature of the hydrated proton. Here, we introduce a new coupled local-mode (CLM) approach to calculate PWC OH stretch vibrational spectra. The CLM method combines a sampling of representative configurations from density functional theory (DFT)-based ab initio molecular dynamics (AIMD) simulations with DFT calculations of local-mode vibrational frequencies and couplings. Calculations of inhomogeneous OH stretch vibrational spectra for H+(H2O)4 and H+(H2O)21 agree well with experiment and higher-level calculations, and decompositions of the calculated spectra in terms of the coupled modes aids in the interpretation of the spectra. This observation is consistent with the idea that capturing anharmonicity and coupling is as important to accuracy as the underlying level of electronic structure theory. The CLM calculations can easily discern the configuration that dominates the experimental measurement for H+(H2O)5, which can adopt several low-energy conformations.
Collapse
Affiliation(s)
- Matthew D Hanson
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Janel A Readnour
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Ali A Hassanali
- Condensed Matter and Statistical Physics, International Centre for Theoretical Physics, Strada Costiera, 11 I - 34151 Trieste, Italy
| | - Steven A Corcelli
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| |
Collapse
|
20
|
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: 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
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
|
21
|
Weike T, Manthe U. Symmetries in the multi-configurational time-dependent Hartree wavefunction representation and propagation. J Chem Phys 2021; 154:194108. [PMID: 34240912 DOI: 10.1063/5.0054105] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In multi-configurational time-dependent Hartree (MCTDH) approaches, different multi-layered wavefunction representations can be used to represent the same physical wavefunction. Transformations between different equivalent representations of a physical wavefunction that alter the tree structure used in the multi-layer MCTDH wavefunction representation interchange the role of single-particle functions (SPFs) and single-hole functions (SHFs) in the MCTDH formalism. While the physical wavefunction is invariant under these transformations, this invariance does not hold for the standard multi-layer MCTDH equations of motion. Introducing transformed SPFs, which obey normalization conditions typically associated with SHFs, revised equations of motion are derived. These equations do not show the singularities resulting from the inverse single-particle density matrix and are invariant under tree transformations. Based on the revised equations of motion, a new integration scheme is introduced. The scheme combines the advantages of the constant mean-field approach of Beck and Meyer [Z. Phys. D 42, 113 (1997)] and the singularity-free integrator suggested by Lubich [Appl. Math. Res. Express 2015, 311]. Numerical calculations studying the spin boson model in high dimensionality confirm the favorable properties of the new integration scheme.
Collapse
Affiliation(s)
- Thomas Weike
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
22
|
Mainali S, Gatti F, Iouchtchenko D, Roy PN, Meyer HD. Comparison of the multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) method and the density matrix renormalization group (DMRG) for ground state properties of linear rotor chains. J Chem Phys 2021; 154:174106. [PMID: 34241072 DOI: 10.1063/5.0047090] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We demonstrate the applicability of the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method to the problem of computing ground states of one-dimensional chains of linear rotors with dipolar interactions. Specifically, we successfully obtain energies, entanglement entropies, and orientational correlations that are in agreement with the Density Matrix Renormalization Group (DMRG), which has been previously used for this system. We find that the entropies calculated by ML-MCTDH for larger system sizes contain nonmonotonicity, as expected in the vicinity of a second-order quantum phase transition between ordered and disordered rotor states. We observe that this effect remains when all couplings besides nearest-neighbor are omitted from the Hamiltonian, which suggests that it is not sensitive to the rate of decay of the interactions. In contrast to DMRG, which is tailored to the one-dimensional case, ML-MCTDH (as implemented in the Heidelberg MCTDH package) requires more computational time and memory, although the requirements are still within reach of commodity hardware. The numerical convergence and computational demand of two practical implementations of ML-MCTDH and DMRG are presented in detail for various combinations of system parameters.
Collapse
Affiliation(s)
- Samrit Mainali
- Université Paris-Saclay, Institut des Sciences Moléculaires d'Orsay ISMO, UMR CNRS 8214, F-91405 Orsay, France
| | - Fabien Gatti
- Université Paris-Saclay, Institut des Sciences Moléculaires d'Orsay ISMO, UMR CNRS 8214, F-91405 Orsay, France
| | - Dmitri Iouchtchenko
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| |
Collapse
|
23
|
Zhao B, Manthe U. Direct product-type grid representations for angular coordinates in extended space and their application in the MCTDH approach. J Chem Phys 2021; 154:104115. [PMID: 33722051 DOI: 10.1063/5.0045054] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Multi-configurational time-dependent Hartree (MCTDH) calculations using time-dependent grid representations can be used to accurately simulate high-dimensional quantum dynamics on general ab initio potential energy surfaces. Employing the correlation discrete variable representation, sets of direct product type grids are employed in the calculation of the required potential energy matrix elements. This direct product structure can be a problem if the coordinate system includes polar and azimuthal angles that result in singularities in the kinetic energy operator. In the present work, a new direct product-type discrete variable representation (DVR) for arbitrary sets of polar and azimuthal angles is introduced. It employs an extended coordinate space where the range of the polar angles is taken to be [-π, π]. The resulting extended space DVR resolves problems caused by the singularities in the kinetic energy operator without generating a very large spectral width. MCTDH calculations studying the F·CH4 complex are used to investigate important properties of the new scheme. The scheme is found to allow for more efficient integration of the equations of motion compared to the previously employed cot-DVR approach [G. Schiffel and U. Manthe, Chem. Phys. 374, 118 (2010)] and decreases the required central processing unit times by about an order of magnitude.
Collapse
Affiliation(s)
- Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
24
|
Wang H, Meyer HD. Importance of Appropriately Regularizing the ML-MCTDH Equations of Motion. J Phys Chem A 2021; 125:3077-3087. [DOI: 10.1021/acs.jpca.0c11221] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Haobin Wang
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217-3364, United States
| | - Hans-Dieter Meyer
- Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| |
Collapse
|
25
|
Carpenter WB, Yu Q, Hack JH, Dereka B, Bowman JM, Tokmakoff A. Decoding the 2D IR spectrum of the aqueous proton with high-level VSCF/VCI calculations. J Chem Phys 2020; 153:124506. [PMID: 33003749 DOI: 10.1063/5.0020279] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
The aqueous proton is a common and long-studied species in chemistry, yet there is currently intense interest devoted to understanding its hydration structure and transport dynamics. Typically described in terms of two limiting structures observed in gas-phase clusters, the Zundel H5O2 + and Eigen H9O4 + ions, the aqueous structure is less clear due to the heterogeneity of hydrogen bonding environments and room-temperature structural fluctuations in water. The linear infrared (IR) spectrum, which reports on structural configurations, is challenging to interpret because it appears as a continuum of absorption, and the underlying vibrational modes are strongly anharmonically coupled to each other. Recent two-dimensional IR (2D IR) experiments presented strong evidence for asymmetric Zundel-like motifs in solution, but true structure-spectrum correlations are missing and complicated by the anharmonicity of the system. In this study, we employ high-level vibrational self-consistent field/virtual state configuration interaction calculations to demonstrate that the 2D IR spectrum reports on a broad distribution of geometric configurations of the aqueous proton. We find that the diagonal 2D IR spectrum around 1200 cm-1 is dominated by the proton stretch vibrations of Zundel-like and intermediate geometries, broadened by the heterogeneity of aqueous configurations. There is a wide distribution of multidimensional potential shapes for the proton stretching vibration with varying degrees of potential asymmetry and confinement. Finally, we find specific cross peak patterns due to aqueous Zundel-like species. These studies provide clarity on highly debated spectral assignments and stringent spectroscopic benchmarks for future simulations.
Collapse
Affiliation(s)
- William B Carpenter
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - John H Hack
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Bogdan Dereka
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
26
|
Boutwell D, Okere O, Omodemi O, Toledo A, Barrios A, Olocha M, Kaledin M. Analysis of the Proton Transfer Bands in the Infrared Spectra of Linear N 2H +···OC and N 2D +···OC Complexes Using Electric Field-Driven Classical Trajectories. J Phys Chem A 2020; 124:7549-7558. [PMID: 32808782 DOI: 10.1021/acs.jpca.0c06756] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this work, we describe ab initio calculations and assignment of infrared (IR) spectra of hydrogen-bonded ion-molecular complexes that involve a fluxional proton: the linear N2H+···OC and N2D+···OC complexes. Given the challenges of describing fluxional proton dynamics and especially its IR activity, we use electric field-driven classical trajectories, i.e., the driven molecular dynamics (DMD) method that was developed by us in recent years and for similar applications, in conjunction with high-level electronic structure theory. Namely, we present a modified and a numerically efficient implementation of DMD specifically for direct (or "on the fly") calculations, which we carry out at the MP2-F12/AVDZ level of theory for the potential energy surface (PES) and MP2/AVDZ for the dipole moment surfaces (DMSs). Detailed analysis of the PES, DMS, and the time-dependence of the first derivative of the DMS, referred to as the driving force, for the highly fluxional vibrations involving H+/D+ revealed that the strongly non-harmonic PES and non-linear DMS yield remarkably complex vibrational spectra. Interestingly, the classical trajectories reveal a doublet in the proton transfer part of the spectrum with the two peaks at 1800 and 1980 cm-1. We find that their shared intensity is due to a Fermi-like resonance interaction, within the classical limit, of the H+ parallel stretch fundamental and an H+ perpendicular bending overtone. This doublet is also observed in the deuterated species at 1360 and 1460 cm-1.
Collapse
Affiliation(s)
- Dalton Boutwell
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Onyinye Okere
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Oluwaseun Omodemi
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Alexander Toledo
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Antonio Barrios
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Monique Olocha
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| |
Collapse
|
27
|
Schmitz G, Klinting EL, Christiansen O. A Gaussian process regression adaptive density guided approach for potential energy surface construction. J Chem Phys 2020; 153:064105. [DOI: 10.1063/5.0015344] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
28
|
Klinting EL, Lauvergnat D, Christiansen O. Vibrational Coupled Cluster Computations in Polyspherical Coordinates with the Exact Analytical Kinetic Energy Operator. J Chem Theory Comput 2020; 16:4505-4520. [DOI: 10.1021/acs.jctc.0c00261] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
29
|
Weike T, Manthe U. The multi-configurational time-dependent Hartree approach in optimized second quantization: Imaginary time propagation and particle number conservation. J Chem Phys 2020; 152:034101. [DOI: 10.1063/1.5140984] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Thomas Weike
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
30
|
Schröder M. Transforming high-dimensional potential energy surfaces into a canonical polyadic decomposition using Monte Carlo methods. J Chem Phys 2020; 152:024108. [DOI: 10.1063/1.5140085] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Markus Schröder
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| |
Collapse
|
31
|
Bertaina G, Di Liberto G, Ceotto M. Reduced rovibrational coupling Cartesian dynamics for semiclassical calculations: Application to the spectrum of the Zundel cation. J Chem Phys 2019; 151:114307. [PMID: 31542046 DOI: 10.1063/1.5114616] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study the vibrational spectrum of the protonated water dimer, by means of a divide-and-conquer semiclassical initial value representation of the quantum propagator, as a first step in the study of larger protonated water clusters. We use the potential energy surface from the work of Huang et al. [J. Chem. Phys. 122, 044308 (2005)]. To tackle such an anharmonic and floppy molecule, we employ fully Cartesian dynamics and carefully reduce the coupling to global rotations in the definition of normal modes. We apply the time-averaging filter and obtain clean power spectra relative to suitable reference states that highlight the spectral peaks corresponding to the fundamental excitations of the system. Our trajectory-based approach allows for the physical interpretation of the very challenging proton transfer modes. We find that it is important, for such a floppy molecule, to selectively avoid initially exciting lower energy modes, in order to obtain cleaner spectra. The estimated vibrational energies display a mean absolute error (MAE) of ∼29 cm-1 with respect to available multiconfiguration time-dependent Hartree calculations and MAE ∼ 14 cm-1 when compared to the optically active experimental excitations of the Ne-tagged Zundel cation. The reasonable scaling in the number of trajectories for Monte Carlo convergence is promising for applications to higher dimensional protonated cluster systems.
Collapse
Affiliation(s)
- G Bertaina
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - G Di Liberto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - M Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
32
|
|
33
|
Tan JA, Kuo JL. A theoretical study on the infrared signatures of proton-bound rare gas dimers (Rg-H +-Rg), Rg = {Ne, Ar, Kr, and Xe}. J Chem Phys 2019; 150:124305. [PMID: 30927880 DOI: 10.1063/1.5090031] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The infrared spectrum of proton-bound rare gas dimers has been extensively studied via matrix isolation spectroscopy. However, little attention has been paid on their spectrum in the gas phase. Most of the Rg2H+ has not been detected outside the matrix environment. Recently, ArnH+ (n = 3-7) has been first detected in the gas-phase [D. C. McDonald et al., J. Chem. Phys. 145, 231101 (2016)]. In that work, anharmonic theory can reproduce the observed vibrational structure. In this paper, we extend the existing theory to examine the vibrational signatures of Rg2H+, Rg = {Ne, Ar, Kr, and Xe}. The successive binding of Rg to H+ was investigated through the calculation of stepwise formation energies. It was found that this binding is anti-cooperative. High-level full-dimensional potential energy surfaces at the CCSD(T)/aug-cc-pVQZ//MP2/aug-cc-pVQZ were constructed and used in the anharmonic calculation via discrete variable representation. We found that the potential coupling between the symmetric and asymmetric Rg-H+ stretch (ν1 and ν3 respectively) causes a series of bright n1ν1 + ν3 progressions. From Ne2H+ to Xe2H+, an enhancement of intensities for these bands was observed.
Collapse
Affiliation(s)
- Jake A Tan
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Taipei 10617, Taiwan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Taipei 10617, Taiwan
| |
Collapse
|
34
|
Tan JA, Kuo JL. Multilevel Approach for Direct VSCF/VCI MULTIMODE Calculations with Applications to Large “Zundel” Cations. J Chem Theory Comput 2018; 14:6405-6416. [DOI: 10.1021/acs.jctc.8b00679] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jake A. Tan
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Taipei 10617, Taiwan (ROC)
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Taipei 10617, Taiwan (ROC)
| |
Collapse
|
35
|
McDonald DC, Wagner JP, McCoy AB, Duncan MA. Near-Infrared Spectroscopy and Anharmonic Theory of Protonated Water Clusters: Higher Elevations in the Hydrogen Bonding Landscape. J Phys Chem Lett 2018; 9:5664-5671. [PMID: 30205006 DOI: 10.1021/acs.jpclett.8b02499] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Near-infrared spectroscopy measurements are presented for protonated water clusters, H+(H2O) n, in the size range of n = 1-8. Clusters are produced in a pulsed-discharge supersonic expansion, mass selected, and studied with infrared laser photodissociation spectroscopy in the regions of 3600-4550 and 4850-7350 cm-1. Although there is some variation with cluster size, the main features of these spectra are a broad absorption near 5300 cm-1, a sharp doublet near 7200 cm-1, as well as a structured absorption near 4100 cm-1 for n ≥ 2. The vibrational patterns measured for the hydronium, Zundel, and Eigen ions are compared to those predicted by different forms of anharmonic theory. Second-order vibrational perturbation theory (VPT2) and a local mode treatment of the OH stretches both capture key aspects of the spectra but suffer understandable deficiencies in the quantitative description of band positions and intensities.
Collapse
Affiliation(s)
- D C McDonald
- Department of Chemistry , University of Georgia , Athens , Georgia 30602 , United States
| | - J P Wagner
- Department of Chemistry , University of Georgia , Athens , Georgia 30602 , United States
| | - A B McCoy
- Department of Chemistry , University of Washington , Seattle , Washington 98195 , United States
| | - M A Duncan
- Department of Chemistry , University of Georgia , Athens , Georgia 30602 , United States
| |
Collapse
|
36
|
Affiliation(s)
- Haobin Wang
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217-3364, USA and Beijing Computational Science Research Center, No. 10 East Xibeiwang Road, Haidian District, Beijing 100193, China
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| |
Collapse
|
37
|
Di Liberto G, Conte R, Ceotto M. "Divide and conquer" semiclassical molecular dynamics: A practical method for spectroscopic calculations of high dimensional molecular systems. J Chem Phys 2018; 148:014307. [PMID: 29306274 DOI: 10.1063/1.5010388] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
We extensively describe our recently established "divide-and-conquer" semiclassical method [M. Ceotto, G. Di Liberto, and R. Conte, Phys. Rev. Lett. 119, 010401 (2017)] and propose a new implementation of it to increase the accuracy of results. The technique permits us to perform spectroscopic calculations of high-dimensional systems by dividing the full-dimensional problem into a set of smaller dimensional ones. The partition procedure, originally based on a dynamical analysis of the Hessian matrix, is here more rigorously achieved through a hierarchical subspace-separation criterion based on Liouville's theorem. Comparisons of calculated vibrational frequencies to exact quantum ones for a set of molecules including benzene show that the new implementation performs better than the original one and that, on average, the loss in accuracy with respect to full-dimensional semiclassical calculations is reduced to only 10 wavenumbers. Furthermore, by investigating the challenging Zundel cation, we also demonstrate that the "divide-and-conquer" approach allows us to deal with complex strongly anharmonic molecular systems. Overall the method very much helps the assignment and physical interpretation of experimental IR spectra by providing accurate vibrational fundamentals and overtones decomposed into reduced dimensionality spectra.
Collapse
Affiliation(s)
- 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
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
38
|
Otto F, Chiang YC, Peláez D. Accuracy of Potfit-based potential representations and its impact on the performance of (ML-)MCTDH. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2017.11.013] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
39
|
Affiliation(s)
- Hans-Dieter Meyer
- Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| | - Haobin Wang
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217-3364, USA and Beijing Computational Science Research Center, No. 10 East Xibeiwang Road, Haidian District, Beijing 100193, China
| |
Collapse
|
40
|
Bulik IW, Frisch MJ, Vaccaro PH. Fixed-Node, Importance-Sampling Diffusion Monte Carlo for Vibrational Structure with Accurate and Compact Trial States. J Chem Theory Comput 2018; 14:1554-1563. [DOI: 10.1021/acs.jctc.8b00016] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ireneusz W. Bulik
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Michael J. Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, United States
| | - Patrick H. Vaccaro
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| |
Collapse
|
41
|
Avila G, Carrington T. Reducing the cost of using collocation to compute vibrational energy levels: Results for CH2NH. J Chem Phys 2017; 147:064103. [DOI: 10.1063/1.4994920] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gustavo Avila
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
42
|
Schröder M, Meyer HD. Transforming high-dimensional potential energy surfaces into sum-of-products form using Monte Carlo methods. J Chem Phys 2017; 147:064105. [DOI: 10.1063/1.4991851] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
43
|
Larsson HR, Tannor DJ. Dynamical pruning of the multiconfiguration time-dependent Hartree (DP-MCTDH) method: An efficient approach for multidimensional quantum dynamics. J Chem Phys 2017; 147:044103. [DOI: 10.1063/1.4993219] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- H. R. Larsson
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, 24098 Kiel, Germany
- Department of Chemical Physics, Weizmann Institute of Science, 76100 Rehovot, Israel
| | - D. J. Tannor
- Department of Chemical Physics, Weizmann Institute of Science, 76100 Rehovot, Israel
| |
Collapse
|
44
|
Bulik IW, Frisch MJ, Vaccaro PH. Vibrational self-consistent field theory using optimized curvilinear coordinates. J Chem Phys 2017; 147:044110. [DOI: 10.1063/1.4995440] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ireneusz W. Bulik
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| | - Michael J. Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | - Patrick H. Vaccaro
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| |
Collapse
|
45
|
Pitsevich G, Malevich A, Kozlovskaya E, Shalamberidze E, Doroshenko I, Pogorelov V, Mahnach E, Sapeshko V, Balevicius V. MP4 study of the multimode coupling in protonated water dimer. J Mol Struct 2017. [DOI: 10.1016/j.molstruc.2017.03.022] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
46
|
Manthe U. Wavepacket dynamics and the multi-configurational time-dependent Hartree approach. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:253001. [PMID: 28430111 DOI: 10.1088/1361-648x/aa6e96] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
Multi-configurational time-dependent Hartree (MCTDH) based approaches are efficient, accurate, and versatile methods for high-dimensional quantum dynamics simulations. Applications range from detailed investigations of polyatomic reaction processes in the gas phase to high-dimensional simulations studying the dynamics of condensed phase systems described by typical solid state physics model Hamiltonians. The present article presents an overview of the different areas of application and provides a comprehensive review of the underlying theory. The concepts and guiding ideas underlying the MCTDH approach and its multi-mode and multi-layer extensions are discussed in detail. The general structure of the equations of motion is highlighted. The representation of the Hamiltonian and the correlated discrete variable representation (CDVR), which provides an efficient multi-dimensional quadrature in MCTDH calculations, are discussed. Methods which facilitate the calculation of eigenstates, the evaluation of correlation functions, and the efficient representation of thermal ensembles in MCTDH calculations are described. Different schemes for the treatment of indistinguishable particles in MCTDH calculations and recent developments towards a unified multi-layer MCTDH theory for systems including bosons and fermions are discussed.
Collapse
Affiliation(s)
- Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
47
|
Yagi K, Thomsen B. Infrared Spectra of Protonated Water Clusters, H+(H2O)4, in Eigen and Zundel Forms Studied by Vibrational Quasi-Degenerate Perturbation Theory. J Phys Chem A 2017; 121:2386-2398. [DOI: 10.1021/acs.jpca.6b11189] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory and ‡iTHES, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Bo Thomsen
- Theoretical
Molecular Science Laboratory and ‡iTHES, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|
48
|
Pitsevich G, Malevich A, Kozlovskaya E, Mahnach E, Doroshenko I, Pogorelov V, Pettersson LGM, Sablinskas V, Balevicius V. MP4 Study of the Anharmonic Coupling of the Shared Proton Stretching Vibration of the Protonated Water Dimer in Equilibrium and Transition States. J Phys Chem A 2017; 121:2151-2165. [PMID: 28187260 DOI: 10.1021/acs.jpca.7b00536] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The structure and harmonic and anharmonic IR spectra of the protonated water dimer (PWD) were calculated in C1, C2, and Cs symmetry at the MP4/acc-pVTZ level of theory. We found that structure and IR spectra are practically identical in C2 and C1 symmetry, demonstrating that an equilibrium C1 configuration of the PWD is not realized. Anharmonic coupling of the shared proton stretching vibration with all other modes in the PWD in C2 and Cs symmetry was the focus of this investigation. For this purpose, 28 two-dimensional potential energy surfaces (2D PES) were built at the MP4/acc-pVTZ level of theory and the corresponding vibrational Schrödinger equations were solved using the DVR method. Differences in the coupling of the investigated mode with other modes in the C2 and Cs configurations, along with some factors that determine the red- or blue-shift of the stretching vibration frequency, were analyzed. We obtained a rather reasonable value of the stretching frequency of the bridging proton (1058.4 cm-1) unperturbed by Fermi resonance. The Fermi resonance between the fundamental vibration ν7 and the combined vibration ν2 + ν6 of the same symmetry was analyzed through anharmonic second-order perturbation theory calculations, as well as by 3D PES constructed using Q2, Q6, and Q7 as normal coordinates. A significant (up to 50%) transfer of intensity from the fundamental vibration to the combined one was found. We have estimated the frequency of the bridging proton stretching vibration in the Cs configuration of the PWD based on calculations of the intrinsic anharmonicity and anharmonic double modes interactions at the MP4/acc-pVTZ level of theory (1261 cm-1).
Collapse
Affiliation(s)
- G Pitsevich
- Department of Physical Optics, Belarusian State University , Nezavisimosti ave., 4, 220030 Minsk, Belarus
| | - A Malevich
- Department of Physical Optics, Belarusian State University , Nezavisimosti ave., 4, 220030 Minsk, Belarus
| | - E Kozlovskaya
- Department of Physical Optics, Belarusian State University , Nezavisimosti ave., 4, 220030 Minsk, Belarus
| | - E Mahnach
- Department of Physical Optics, Belarusian State University , Nezavisimosti ave., 4, 220030 Minsk, Belarus
| | - I Doroshenko
- Taras Shevchenko National University of Kyiv , Volodymyrska str., 64\13, 01601, Kyiv, Ukraine
| | - V Pogorelov
- Taras Shevchenko National University of Kyiv , Volodymyrska str., 64\13, 01601, Kyiv, Ukraine
| | - Lars G M Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University , S-106 91 Stockholm, Sweden
| | - V Sablinskas
- Faculty of Physics, Vilnius University , Sauletekio al. 9-3, LT-10222 Vilnius, Lithuania
| | - V Balevicius
- Faculty of Physics, Vilnius University , Sauletekio al. 9-3, LT-10222 Vilnius, Lithuania
| |
Collapse
|
49
|
Manthe U, Weike T. On the multi-layer multi-configurational time-dependent Hartree approach for bosons and fermions. J Chem Phys 2017; 146:064117. [DOI: 10.1063/1.4975662] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Sibaev M, Crittenden DL. An efficient and numerically stable procedure for generating sextic force fields in normal mode coordinates. J Chem Phys 2016; 144:214107. [DOI: 10.1063/1.4953080] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- M. Sibaev
- Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| | - D. L. Crittenden
- Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|