An intertwined method for making low-rank, sum-of-product basis functions that makes it possible to compute vibrational spectra of molecules with more than 10 atoms

Kylin-V: An open-source package calculating the dynamic and spectroscopic properties of large systems

Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples.

Using nested tensor train contracted basis functions with group theoretical techniques to compute (ro)-vibrational spectra of molecules with non-Abelian groups

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.

A numerical-tensorial "hybrid" nuclear motion Hamiltonian and dipole moment operator for spectra calculation of polyatomic nonrigid molecules

The analysis and modeling of high-resolution spectra of nonrigid molecules require a specific Hamiltonian and group-theoretical formulation that differs significantly from that of more familiar rigid systems. Within the framework of Hougen-Bunker-Johns (HBJ) theory, this paper is devoted to the construction of a nonrigid Hamiltonian based on a suitable combination of numerical calculations for the nonrigid part in conjunction with the irreducible tensor operator method for the rigid part. For the first time, a variational calculation from ab initio potential energy surfaces is performed using the HBJ kinetic energy operator built from vibrational, large-amplitude motion, and rotational tensor operators expressed in terms of curvilinear and normal coordinates. Group theory for nonrigid molecules plays a central role in the characterization of the overall tunneling splittings and is discussed in the present approach. The construction of the dipole moment operator is also examined. Validation tests consisting of a careful convergence study of the energy levels as well as a comparison of results obtained from independent computer codes are given for the nonrigid molecules CH2, CH3, NH3, and H2O2. This work paves the way for the modeling of high-resolution spectra of larger nonrigid systems.

Exact quantum dynamics developments for floppy molecular systems and complexes

Molecular rotation, vibration, internal rotation, isomerization, tunneling, intermolecular dynamics of weakly and strongly interacting systems, intra-to-inter-molecular energy transfer, hindered rotation and hindered translation over surfaces are important types of molecular motions. Their fundamentally correct and detailed description can be obtained by solving the nuclear Schrödinger equation on a potential energy surface. Many of the chemically interesting processes involve quantum nuclear motions which are 'delocalized' over multiple potential energy wells. These 'large-amplitude' motions in addition to the high dimensionality of the vibrational problem represent challenges to the current (ro)vibrational methodology. A review of the quantum nuclear motion methodology is provided, current bottlenecks of solving the nuclear Schrödinger equation are identified, and solution strategies are reviewed. Technical details, computational results, and analysis of these results in terms of limiting models and spectroscopically relevant concepts are highlighted for selected numerical examples.

Novel methodology for systematically constructing global effective models from ab initio-based surfaces: A new insight into high-resolution molecular spectra analysis

In this paper, a novel methodology is presented for the construction of ab initio effective rotation-vibration spectroscopic models from potential energy and dipole moment surfaces. Non-empirical effective Hamiltonians are obtained via the block-diagonalization of selected variationally computed eigenvector matrices. For the first time, the derivation of an effective dipole moment is carried out in a systematic way. This general approach can be implemented quite easily in most of the variational computer codes and turns out to be a clear alternative to the rather involved Van Vleck perturbation method. Symmetry is exploited at all stages to translate first-principles calculations into a set of spectroscopic parameters to be further refined on experiment. We demonstrate on H₂CO, PH₃, CH₄, C₂H₄, and SF₆ that the proposed effective model can provide crucial information to spectroscopists within a very short time compared to empirical spectroscopic models. This approach brings a new insight into high-resolution spectrum analysis of polyatomic molecules and will be also of great help in the modeling of hot atmospheres where completeness is important.

Computing vibrational energy levels by solving linear equations using a tensor method with an imposed rank

Present day computers do not have enough memory to store the high-dimensional tensors required when using a direct product basis to compute vibrational energy levels of a polyatomic molecule with more than about five atoms. One way to deal with this problem is to represent tensors using a tensor format. In this paper, we use the canonical polyadic (CP) format. Energy levels are computed by building a basis from vectors obtained by solving linear equations. The method can be thought of as a CP realization of a block inverse iteration method with multiple shifts. The CP rank of the tensors is fixed, and the linear equations are solved with an method. There is no need for rank reduction and no need for orthogonalization, and tensors with a rank larger than the fixed rank used to solve the linear equations are never generated. The ideas are tested by computing vibrational energy levels of a 64-D bilinearly coupled model Hamiltonian and of acetonitrile (12-D).

Fingerprint region of the formic acid dimer: variational vibrational computations in curvilinear coordinates

Curvilinear kinetic energy models are developed for variational nuclear motion computations including the inter- and the low-frequency intra-molecular degrees of freedom of the formic acid dimer. The coupling of the inter- and intra-molecular modes is studied by solving the vibrational Schrödinger equation for a series of vibrational models, from two up to ten active vibrational degrees of freedom by selecting various combinations of active modes and constrained coordinate values. Vibrational states, nodal assignment, and infrared vibrational intensity information is computed using the full-dimensional potential energy surface (PES) and electric dipole moment surface developed by Qu and Bowman [Phys. Chem. Chem. Phys., 2016, 18, 24835; J. Chem. Phys., 2018, 148, 241713]. Good results are obtained for several fundamental and combination bands in comparison with jet-cooled vibrational spectroscopy experiments, but the description of the ν₈ and ν₉ fundamental vibrations, which are close in energy and have the same symmetry, appears to be problematic. For further progress in comparison with experiment, the potential energy surface, and in particular, its multi-dimensional couplings representation, requires further improvement.

Neural Network Potential Energy Surfaces for Small Molecules and Reactions

We review progress in neural network (NN)-based methods for the construction of interatomic potentials from discrete samples (such as ab initio energies) for applications in classical and quantum dynamics including reaction dynamics and computational spectroscopy. The main focus is on methods for building molecular potential energy surfaces (PES) in internal coordinates that explicitly include all many-body contributions, even though some of the methods we review limit the degree of coupling, due either to a desire to limit computational cost or to limited data. Explicit and direct treatment of all many-body contributions is only practical for sufficiently small molecules, which are therefore our primary focus. This includes small molecules on surfaces. We consider direct, single NN PES fitting as well as more complex methods that impose structure (such as a multibody representation) on the PES function, either through the architecture of one NN or by using multiple NNs. We show how NNs are effective in building representations with low-dimensional functions including dimensionality reduction. We consider NN-based approaches to build PESs in the sums-of-product form important for quantum dynamics, ways to treat symmetry, and issues related to sampling data distributions and the relation between PES errors and errors in observables. We highlight combinations of NNs with other ideas such as permutationally invariant polynomials or sums of environment-dependent atomic contributions, which have recently emerged as powerful tools for building highly accurate PESs for relatively large molecular and reactive systems.

Exact quantum dynamics background of dispersion interactions: case study for CH4·Ar in full (12) dimensions

A full-dimensional ab initio potential energy surface of spectroscopic quality is developed for the van-der-Waals complex of a methane molecule and an argon atom. Variational vibrational states are computed on this surface including all twelve (12) vibrational degrees of freedom of the methane-argon complex using the GENIUSH computer program and the Smolyak sparse grid method. The full-dimensional computations make it possible to study the fine details of the interaction and distortion effects and to make a direct assessment of the reduced-dimensionality models often used in the quantum dynamics study of weakly-bound complexes. A 12-dimensional (12D) vibrational computation including only a single harmonic oscillator basis function (9D) to describe the methane fragment (for which we use the ground-state effective structure as the reference structure) has a 0.40 cm^-1 root-mean-square error (rms) with respect to the converged 12D bound-state excitation energies, which is less than half of the rms of the 3D model set up with the 〈r〉₀ methane structure. Allowing 10 basis functions for the methane fragment in a 12D computation performs much better than the 3D models by reducing the rms of the bound state vibrational energies to 0.07 cm^-1. The full-dimensional potential energy surface correctly describes the dissociation of the system, which together with further development of the variational (ro)vibrational methodology opens a route to the study of the role of dispersion forces in the excited methane vibrations and the energy transfer from the intra- to the intermolecular vibrational modes.

Toward breaking the curse of dimensionality in (ro)vibrational computations of molecular systems with multiple large-amplitude motions

Methodological progress is reported in the challenging direction of a black-box-type variational solution of the (ro)vibrational Schrödinger equation applicable to floppy, polyatomic systems with multiple large-amplitude motions. This progress is achieved through the combination of (i) the numerical kinetic-energy operator (KEO) approach of Mátyus et al. [J. Chem. Phys. 130, 134112 (2009)] and (ii) the Smolyak nonproduct grid method of Avila and Carrington, Jr. [J. Chem. Phys. 131, 174103 (2009)]. The numerical representation of the KEO makes it possible to choose internal coordinates and a body-fixed frame best suited for the molecular system. The Smolyak scheme reduces the size of the direct-product grid representation by orders of magnitude, while retaining some of the useful features of it. As a result, multidimensional (ro)vibrational states are computed with system-adapted coordinates, a compact basis- and grid-representation, and an iterative eigensolver. Details of the methodological developments and the first numerical applications are presented for the CH₄·Ar complex treated in full (12D) vibrational dimensionality.

Anharmonic vibrational spectroscopy of polycyclic aromatic hydrocarbons (PAHs)

While powerful techniques exist to accurately account for anharmonicity in vibrational molecular spectroscopy, they are computationally very expensive and cannot be routinely employed for large species and/or at non-zero vibrational temperatures. Motivated by the study of Polycyclic Aromatic Hydrocarbon (PAH) emission in space, we developed a new code, which takes into account all modes and can describe all infrared transitions including bands becoming active due to resonances as well as overtone, combination, and difference bands. In this article, we describe the methodology that was implemented and discuss how the main difficulties were overcome, so as to keep the problem tractable. Benchmarking with high-level calculations was performed on a small molecule. We carried out specific convergence tests on two prototypical PAHs, pyrene (C₁₆H₁₀) and coronene (C₂₄H₁₂), aiming at optimising tunable parameters to achieve both acceptable accuracy and computational costs for this class of molecules. We then report the results obtained at 0 K for pyrene and coronene, comparing the calculated spectra with available experimental data. The theoretical band positions were found to be significantly improved compared to harmonic density functional theory calculations. The band intensities are in reasonable agreement with experiments, the main limitation being the accuracy of the underlying calculations of the quartic force field. This is a first step toward calculating moderately high-temperature spectra of PAHs and other similarly rigid molecules using Monte Carlo sampling.