1
|
Rey M, Carrington T. Using nested tensor train contracted basis functions with group theoretical techniques to compute (ro)-vibrational spectra of molecules with non-Abelian groups. J Chem Phys 2024; 161:044102. [PMID: 39037133 DOI: 10.1063/5.0219434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2024] [Accepted: 06/30/2024] [Indexed: 07/23/2024] Open
Abstract
In this paper, we use nested tensor-train contractions to compute vibrational and ro-vibrational energy levels of molecules with five and six atoms. At each step, we fully exploit symmetry by using symmetry adapted basis functions obtained from an irreducible tensor method. Contracted basis functions are determined by diagonalizing reduced dimensional Hamiltonian matrices. The size of matrices of eigenvectors, used to account for coupling between groups of coordinates, is reduced by discarding rows and columns. The size of the matrices that must be diagonalized is thus substantially reduced, making it possible to use direct eigensolvers, even for molecules with five and six atoms. The symmetry-adapted contracted vibrational basis functions have been used to compute J = 0 energy levels of the CH3CN (C3v) and J > 0 levels of CH4.
Collapse
Affiliation(s)
- Michaël Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
2
|
Wodraszka R, Carrington T. Using a pruned basis and a sparse collocation grid with more points than basis functions to do efficient and accurate MCTDH calculations with general potential energy surfaces. J Chem Phys 2024; 160:214121. [PMID: 38836450 DOI: 10.1063/5.0214557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Accepted: 05/14/2024] [Indexed: 06/06/2024] Open
Abstract
We propose a new collocation multi-configuration time-dependent Hartree (MCTDH) method. It reduces point-set error by using more points than basis functions. Collocation makes it possible to use MCTDH with a general potential energy surface without computing any integrals. The collocation points are associated with a basis larger than the basis used to represent wavefunctions. Both bases are obtained from a direct product basis built from single-particle functions by imposing a pruning condition. The collocation points are those on a sparse grid. Heretofore, collocation MCTDH calculations with more points than basis functions have only been possible if both the collocation grid and the basis set are direct products. In this paper, we exploit a new pseudo-inverse to use both more points than basis functions and a pruned basis and grid. We demonstrate that, for a calculation of the lowest 50 vibrational states (energy levels and wavefunctions) of CH2NH, errors can be reduced by two orders of magnitude by increasing the number of points, without increasing the basis size. This is true also when unrefined time-independent points are used.
Collapse
Affiliation(s)
- Robert Wodraszka
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
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
|
6
|
Richings GW, Habershon S. Predicting Molecular Photochemistry Using Machine-Learning-Enhanced Quantum Dynamics Simulations. Acc Chem Res 2022; 55:209-220. [PMID: 34982533 DOI: 10.1021/acs.accounts.1c00665] [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/30/2022]
Abstract
The processes which occur after molecules absorb light underpin an enormous range of fundamental technologies and applications, including photocatalysis to enable new chemical transformations, sunscreens to protect against the harmful effects of UV overexposure, efficient photovoltaics for energy generation from sunlight, and fluorescent probes to image the intricate details of complex biomolecular structures. Reflecting this broad range of applications, an enormously versatile set of experiments are now regularly used to interrogate light-driven chemical dynamics, ranging from the typical ultrafast transient absorption spectroscopy used in many university laboratories to the inspiring central facilities around the world, such as the next-generation of X-ray free-electron lasers.Computer simulations of light-driven molecular and material dynamics are an essential route to analyzing the enormous amount of transient electronic and structural data produced by these experimental sources. However, to date, the direct simulation of molecular photochemistry remains a frontier challenge in computational chemical science, simultaneously demanding the accurate treatment of molecular electronic structure, nuclear dynamics, and the impact of nonadiabatic couplings.To address these important challenges and to enable new computational methods which can be integrated with state-of-the-art experimental capabilities, the past few years have seen a burst of activity in the development of "direct" quantum dynamics methods, merging the machine learning of potential energy surfaces (PESs) and nonadiabatic couplings with accurate quantum propagation schemes such as the multiconfiguration time-dependent Hartree (MCTDH) method. The result of this approach is a new generation of direct quantum dynamics tools in which PESs are generated in tandem with wave function propagation, enabling accurate "on-the-fly" simulations of molecular photochemistry. These simulations offer an alternative route toward gaining quantum dynamics insights, circumventing the challenge of generating ab initio electronic structure data for PES fitting by instead only demanding expensive energy evaluations as and when they are needed.In this Account, we describe the chronological evolution of our own contributions to this field, focusing on describing the algorithmic developments that enable direct MCTDH simulations for complex molecular systems moving on multiple coupled electronic states. Specifically, we highlight active learning strategies for generating PESs during grid-based quantum chemical dynamics simulations, and we discuss the development and impact of novel diabatization schemes to enable direct grid-based simulations of photochemical dynamics; these developments are highlighted in a series of benchmark molecular simulations of systems containing multiple nuclear degrees of freedom moving on multiple coupled electronic states. We hope that the ongoing developments reported here represent a major step forward in tools for modeling excited-state chemistry such as photodissociation, proton and electron transfer, and ultrafast energy dissipation in complex molecular systems.
Collapse
Affiliation(s)
- Gareth W. Richings
- Department of Chemistry, University of Warwick, Coventry, United Kingdom CV4 7AL
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Coventry, United Kingdom CV4 7AL
| |
Collapse
|
7
|
Wodraszka R, Carrington T. A rectangular collocation multi-configuration time-dependent Hartree (MCTDH) approach with time-independent points for calculations on general potential energy surfaces. J Chem Phys 2021; 154:114107. [PMID: 33752363 DOI: 10.1063/5.0046425] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce a collocation-based multi-configuration time-dependent Hartree (MCTDH) method that uses more collocation points than basis functions. We call it the rectangular collocation MCTDH (RC-MCTDH) method. It does not require that the potential be a sum of products. RC-MCTDH has the important advantage that it makes it simple to use time-independent collocation points. When using time-independent points, it is necessary to evaluate the potential energy function only once and not repeatedly during an MCTDH calculation. It is inexpensive and straightforward to use RC-MCTDH with combined modes. Using more collocation points than basis functions enables one to reduce errors in energy levels without increasing the size of the single-particle function basis. On the contrary, whenever a discrete variable representation is used, the only way to reduce the quadrature error is to increase the basis size, which then also reduces the basis-set error. We demonstrate that with RC-MCTDH and time-independent points, it is possible to calculate accurate eigenenergies of CH3 and CH4.
Collapse
Affiliation(s)
- Robert Wodraszka
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
8
|
Carrington T. Using collocation to study the vibrational dynamics of molecules. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 248:119158. [PMID: 33218875 DOI: 10.1016/j.saa.2020.119158] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/23/2020] [Accepted: 10/26/2020] [Indexed: 06/11/2023]
Abstract
In this paper, I review collocation methods for solving the time-independent and the time-dependent Schroedinger equation. Unlike traditional variational methods, collocation methods do not require integrals and quadrature. Either collocation or quadrature is necessary if the potential does not have a special form. If the basis is a direct product of univariate bases and the quadrature grid is also a direct product, there exist variational methods that do not require quadrature approximations for potential energy matrix elements. These methods, however, do require storing, in computer memory, vectors with as many components as there are quadrature points. For this reason direct-product variational methods are poor for problems with more than five atoms. There are well established ideas for reducing the size of the basis in a variational calculation. Three such ideas are: 1) prune the direct product basis; 2) use basis functions that are products of multivariate functions; 3) optimise the basis functions (e.g. Multiconfiguration time-dependent Hartree). Reducing the basis size, however, is not enough to the make variational methods tractable because, for all three of these ideas, quadrature rears its ugly head. Collocation is an attractive alternative to variational methods.
Collapse
Affiliation(s)
- Tucker Carrington
- Department of Chemistry, Queen's University, Kingston, Ontario K7L 3N6, Canada.
| |
Collapse
|
9
|
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
|
10
|
Madsen NK, Hansen MB, Worth GA, Christiansen O. MR-MCTDH[n]: Flexible Configuration Spaces and Nonadiabatic Dynamics within the MCTDH[n] Framework. J Chem Theory Comput 2020; 16:4087-4097. [DOI: 10.1021/acs.jctc.0c00379] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
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
| | - Graham A. Worth
- Department of Chemistry, University College London, 20, Gordon St., WC1H 0AJ London, United Kingdom
| | - Ove Christiansen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK−8000 Aarhus C, Denmark
| |
Collapse
|
11
|
Wang XG, Carrington T. A variational calculation of vibrational levels of vinyl radical. J Chem Phys 2020; 152:204311. [PMID: 32486683 DOI: 10.1063/5.0007225] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report the vibrational energy levels of vinyl radical (VR) that are computed with a Lanczos eigensolver and a contracted basis. Many of the levels of the two previous VR variational calculations differ significantly and differ also from those reported in this paper. We identify the source of and correct symmetry errors on the potential energy surfaces used in the previous calculations. VR has two equivalent equilibrium structures. By plotting wavefunction cuts, we show that two tunneling paths play an important role. Using the computed wavefunctions, it is possible to assign many states and thereby to determine tunneling splittings that are compared with their experimental counterparts. Our computed red shift of the hot band at 2897.23 cm-1, observed by Dong et al. [J Chem. Phys. 128, 044305 (2008)], is 4.47 cm-1, which is close to the experimental value of 4.63 cm-1.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
12
|
Wodraszka R, Carrington T. A collocation-based multi-configuration time-dependent Hartree method using mode combination and improved relaxation. J Chem Phys 2020; 152:164117. [DOI: 10.1063/5.0006081] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Robert Wodraszka
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
13
|
Madsen NK, Hansen MB, Worth GA, Christiansen O. Systematic and variational truncation of the configuration space in the multiconfiguration time-dependent Hartree method: The MCTDH[n] hierarchy. J Chem Phys 2020; 152:084101. [DOI: 10.1063/1.5142459] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.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
| | - Graham A. Worth
- Department of Chemistry, University College London, 20, Gordon St., WC1H 0AJ London, United Kingdom
| | - Ove Christiansen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK–8000 Aarhus C, Denmark
| |
Collapse
|
14
|
Choi S, Vaníček J. A time-reversible integrator for the time-dependent Schrödinger equation on an adaptive grid. J Chem Phys 2019; 151:234102. [PMID: 31864241 DOI: 10.1063/1.5127856] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
One of the most accurate methods for solving the time-dependent Schrödinger equation uses a combination of the dynamic Fourier method with the split-operator algorithm on a tensor-product grid. To reduce the number of required grid points, we let the grid move together with the wavepacket but find that the naïve algorithm based on an alternate evolution of the wavefunction and grid destroys the time reversibility of the exact evolution. Yet, we show that the time reversibility is recovered if the wavefunction and grid are evolved simultaneously during each kinetic or potential step; this is achieved by using the Ehrenfest theorem together with the splitting method. The proposed algorithm is conditionally stable, symmetric, and time-reversible and conserves the norm of the wavefunction. The preservation of these geometric properties is shown analytically and demonstrated numerically on a three-dimensional harmonic model and collinear model of He-H2 scattering. We also show that the proposed algorithm can be symmetrically composed to obtain time-reversible integrators of an arbitrary even order. We observed 10 000-fold speedup by using the tenth-order instead of the second-order method to obtain a solution with a time discretization error below 10-9. Moreover, using the adaptive grid instead of the fixed grid resulted in a 64-fold reduction in the required number of grid points in the harmonic system and made it possible to simulate the He-H2 scattering for six times longer while maintaining reasonable accuracy. Applicability of the algorithm to high-dimensional quantum dynamics is demonstrated using the strongly anharmonic eight-dimensional Hénon-Heiles model.
Collapse
Affiliation(s)
- Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
15
|
Köhler F, Keiler K, Mistakidis SI, Meyer HD, Schmelcher P. Dynamical pruning of the non-equilibrium quantum dynamics of trapped ultracold bosons. J Chem Phys 2019. [DOI: 10.1063/1.5104344] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- F. Köhler
- Center for Optical Quantum Technologies, Department of Physics, University of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
- Hamburg Center for Ultrafast Imaging, University of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - K. Keiler
- Center for Optical Quantum Technologies, Department of Physics, University of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - S. I. Mistakidis
- Center for Optical Quantum Technologies, Department of Physics, University of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - H.-D. Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - P. Schmelcher
- Center for Optical Quantum Technologies, Department of Physics, University of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
- Hamburg Center for Ultrafast Imaging, University of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
16
|
Wodraszka R, Carrington T. A pruned collocation-based multiconfiguration time-dependent Hartree approach using a Smolyak grid for solving the Schrödinger equation with a general potential energy surface. J Chem Phys 2019; 150:154108. [DOI: 10.1063/1.5093317] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Robert Wodraszka
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
17
|
Viglaska D, Rey M, Nikitin AV, Tyuterev VG. Isotopic and symmetry breaking effects on phosphine spectra under H → D substitutions from ab initio variational calculations. J Chem Phys 2018; 149:174305. [PMID: 30409001 DOI: 10.1063/1.5045525] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Variationally computed infrared spectra in the range [0-5000] cm-1 are reported for the deuterated PH2D and PHD2 molecules from accurate potential energy and dipole moment surfaces initially derived for the major isotopologue PH3( C 3 v ). Energy level and line intensity calculations were performed by using a normal-mode model combined with isotopic and symmetry transformations for the H → D substitutions. Theoretical spectra were computed at 296 K up to J max = 30 and will be made available through the TheoReTS information system (http://theorets.univ-reims.fr, http://theorets.tsu.ru). For the very first time, ab initio intensity predictions of PH2D/PHD2 are in good qualitative agreement with the literature. This work will be useful for spectral intensity analysis for which accurate spectral intensity data are still missing.
Collapse
Affiliation(s)
- Dominika Viglaska
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| | - Michael Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| | - Andrei V Nikitin
- Laboratory of Theoretical Spectroscopy, Institute of Atmospheric Optics, SB RAS, 634055 Tomsk, Russia
| | - Vladimir G Tyuterev
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| |
Collapse
|
18
|
Thomas PS, Carrington T, Agarwal J, Schaefer HF. Using an iterative eigensolver and intertwined rank reduction to compute vibrational spectra of molecules with more than a dozen atoms: Uracil and naphthalene. J Chem Phys 2018; 149:064108. [DOI: 10.1063/1.5039147] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Affiliation(s)
- Phillip S. Thomas
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Jay Agarwal
- Center for Computational Chemistry, University of Georgia, Athens, Georgia 30602-0525, USA
| | - Henry F. Schaefer
- Center for Computational Chemistry, University of Georgia, Athens, Georgia 30602-0525, USA
| |
Collapse
|
19
|
Larsson HR, Riedel J, Wei J, Temps F, Hartke B. Resonance dynamics of DCO (X̃ A′2) simulated with the dynamically pruned discrete variable representation (DP-DVR). J Chem Phys 2018; 148:204309. [DOI: 10.1063/1.5026459] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Henrik R. Larsson
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Jens Riedel
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Jie Wei
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Friedrich Temps
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Bernd Hartke
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| |
Collapse
|
20
|
Jerke J, Poirier B. Two-body Schrödinger wave functions in a plane-wave basis via separation of dimensions. J Chem Phys 2018; 148:104101. [DOI: 10.1063/1.5017621] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jonathan Jerke
- Department of Chemistry and Biochemistry, and Department of Physics, Texas Tech University, P.O. Box 41061, Lubbock, Texas 79409-1061, USA
- Department of Physics, Texas Southern University, Houston, Texas 77004, USA
| | - Bill Poirier
- Department of Chemistry and Biochemistry, and Department of Physics, Texas Tech University, P.O. Box 41061, Lubbock, Texas 79409-1061, USA
| |
Collapse
|
21
|
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
|
22
|
Wodraszka R, Carrington T. A new collocation-based multi-configuration time-dependent Hartree (MCTDH) approach for solving the Schrödinger equation with a general potential energy surface. J Chem Phys 2018; 148:044115. [DOI: 10.1063/1.5018793] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Robert Wodraszka
- Department of Chemistry, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Department of Chemistry, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
23
|
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
|
24
|
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
|
25
|
Wodraszka R, Carrington T. Systematically expanding nondirect product bases within the pruned multi-configuration time-dependent Hartree (MCTDH) method: A comparison with multi-layer MCTDH. J Chem Phys 2017; 146:194105. [PMID: 28527461 DOI: 10.1063/1.4983281] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose a pruned multi-configuration time-dependent Hartree (MCTDH) method with systematically expanding nondirect product bases and use it to solve the time-independent Schrödinger equation. No pre-determined pruning condition is required to select the basis functions. Using about 65 000 basis functions, we calculate the first 69 vibrational eigenpairs of acetonitrile, CH3CN, to an accuracy better than that achieved in a previous pruned MCTDH calculation which required more than 100 000 basis functions. In addition, we compare the new pruned MCTDH method with the established multi-layer MCTDH (ML-MCTDH) scheme and determine that although ML-MCTDH is somewhat more efficient when low or intermediate accuracy is desired, pruned MCTDH is more efficient when high accuracy is required. In our largest calculation, the vast majority of the energies have errors smaller than 0.01 cm-1.
Collapse
Affiliation(s)
- Robert Wodraszka
- Department of Chemistry, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Department of Chemistry, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
26
|
Leclerc A, Thomas PS, Carrington T. Comparison of different eigensolvers for calculating vibrational spectra using low-rank, sum-of-product basis functions. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1249980] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Arnaud Leclerc
- Université de Lorraine , UMR CNRS 7565 SRSMC, Metz, France
| | - Phillip S. Thomas
- Chemistry Department, Queen's University , Kingston, Ontario, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University , Kingston, Ontario, Canada
| |
Collapse
|
27
|
Brown J, Carrington T. Using an expanding nondirect product harmonic basis with an iterative eigensolver to compute vibrational energy levels with as many as seven atoms. J Chem Phys 2016; 145:144104. [DOI: 10.1063/1.4963916] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- James Brown
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|