Pruning the Hamiltonian Matrix in MULTIMODE: Test for C2H4 and Application to CH3NO2 Using a New Ab Initio Potential Energy Surface

Resource Optimization for Quantum Dynamics with Tensor Networks: Quantum and Classical Algorithms

The exponential scaling of the quantum degrees of freedom with the size of the system is one of the biggest challenges in computational chemistry and particularly in quantum dynamics. We present a tensor network approach for the time-evolution of the nuclear degrees of freedom of multiconfigurational chemical systems at a reduced storage and computational complexity. We also present quantum algorithms for the resultant dynamics. To preserve the compression advantage achieved via tensor network decompositions, we present an adaptive algorithm for the regularization of nonphysical bond dimensions, preventing the potentially exponential growth of these with time. While applicable to any quantum dynamical problem, our method is particularly valuable for dynamical simulations of nuclear chemical systems. Our algorithm is demonstrated using ab initio potentials obtained for a symmetric hydrogen-bonded system, namely, the protonated 2,2'-bipyridine, and compared to exact diagonalization numerical results.

A New A Priori Method to Avoid Calculation of Negligible Hamiltonian Matrix Elements in CI Calculation

Hamiltonian matrices typically contain many elements that are negligibly small compared to the diagonal elements, even with methods to prune the underlying basis. Because for general potentials the calculation of H-matrix elements is a major part of the computational effort to obtain eigenvalues and eigenfunctions of the Hamiltonian, there is strong motivation to investigate locating these negligible elements without calculating them or at least avoid calculating them. We recently demonstrated an effective means to "learn" negligible elements using machine learning classification (J. Chem. Phys. 2023, 159, 071101). Here we present a simple, new method to avoid calculating them by using a cut-off value for the absolute difference in the quantum numbers for the bra and ket. This method is demonstrated for many of the same case studies as were used in the paper above, namely for realistic H-matrices of H2O, the vinyl radical, C2H3, and glycine, C2H5NO2. The new method is compared to the recently reported machine learning approach. In addition, we point out an important synergy between the two methods.

Exact tunneling splittings from symmetrized path integrals

We develop a new simulation technique based on path-integral molecular dynamics for calculating ground-state tunneling splitting patterns from ratios of symmetrized partition functions. In particular, molecular systems are rigorously projected onto their J = 0 rotational state by an "Eckart spring" that connects two adjacent beads in a ring polymer. Using this procedure, the tunneling splitting can be obtained from thermodynamic integration at just one (sufficiently low) temperature. Converged results are formally identical to the values that would have been obtained by solving the full rovibrational Schrödinger equation on a given Born-Oppenheimer potential energy surface. The new approach is showcased with simulations of hydronium and methanol, which are in good agreement with wavefunction-based calculations and experimental measurements. The method will be of particular use for the study of low-barrier methyl rotations and other floppy modes, where instanton theory is not valid.

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.

Graph-Theoretic Molecular Fragmentation for Potential Surfaces Leads Naturally to a Tensor Network Form and Allows Accurate and Efficient Quantum Nuclear Dynamics

Molecular fragmentation methods have revolutionized quantum chemistry. Here, we use a graph-theoretically generated molecular fragmentation method, to obtain accurate and efficient representations for multidimensional potential energy surfaces and the quantum time-evolution operator, which plays a critical role in quantum chemical dynamics. In doing so, we find that the graph-theoretic fragmentation approach naturally reduces the potential portion of the time-evolution operator into a tensor network that contains a stream of coupled lower-dimensional propagation steps to potentially achieve quantum dynamics with reduced complexity. Furthermore, the fragmentation approach used here has previously been shown to allow accurate and efficient computation of post-Hartree-Fock electronic potential energy surfaces, which in many cases has been shown to be at density functional theory cost. Thus, by combining the advantages of molecular fragmentation with the tensor network formalism, the approach yields an on-the-fly quantum dynamics scheme where both the electronic potential calculation and nuclear propagation portion are enormously simplified through a single stroke. The method is demonstrated by computing approximations to the propagator and to potential surfaces for a set of coupled nuclear dimensions within a protonated water wire problem exhibiting the Grotthuss mechanism of proton transport. In all cases, our approach has been shown to reduce the complexity of representing the quantum propagator, and by extension action of the propagator on an initial wavepacket, by several orders, with minimal loss in accuracy.

Semiclassical and VSCF/VCI Calculations of the Vibrational Energies of trans- and gauche-Ethanol Using a CCSD(T) Potential Energy Surface

A recent full-dimensional Δ-Machine learning potential energy surface (PES) for ethanol is employed in semiclassical and vibrational self-consistent field (VSCF) and virtual-state configuration interaction (VCI) calculations, using MULTIMODE, to determine the anharmonic vibrational frequencies of vibration for both the trans and gauche conformers of ethanol. Both semiclassical and VSCF/VCI energies agree well with the experimental data. We find significant mixing between the VSCF basis states due to Fermi resonances between bending and stretching modes. The same effects are also accurately described by the full-dimensional semiclassical calculations. These are the first high-level anharmonic calculations using a PES, in particular a "gold-standard" CCSD(T) one.

Performance of Vibrational Self-Consistent Field Theory for Accurate Potential Energy Surfaces: Fundamentals, Excited States, and Intensities

The performance of vibrational structure calculations beyond harmonic approximation in the framework of the vibrational self-consistent field method with second-order perturbation corrections (VSCF-PT2) is investigated in conjunction with very accurate potential energy surfaces (PESs) given by various coupled-cluster electronic structure theories. The quality of anharmonic calculations depends on the accuracy of the underlying multidimensional PES obtained from its functional form, which is given by the level of electronic structure theory. Two such highest levels of typical coupled-cluster electronic structure methods, CCSD and the ″gold standard″ CCSD(T), along with their variants such as CCD, CR-CCL (completely renormalized CR-CC(2,3) approach), and CCSD(TQ) are tested for the construction of accurate anharmonic potentials without any fitting or ad hoc scaling and using cc-pVTZ basis sets. The accuracy of VSCF-PT2 theory in comparison to experimental values is tested for a series of 16 molecules with 135 fundamental bands, 64 overtones, and combination bands and also for 39 intensities. It is found that CCD and CCSD bind the potential tighter than CCSD(T) and the computed VSCF-PT2 transitions are more blue-shifted showing higher deviation from the experiment. In general, VSCF-PT2 results computed at the CCSD(T) potential offer a good cost/accuracy ratio, with the mean absolute deviation and the mean absolute percentage error with the experiment being ∼16 cm^-1 and 1.38, respectively, for fundamentals. Additionally, while the CR-CCL and CCSD(TQ) methods offer similar levels of accuracies as compared to CCSD(T), the former offers a better accuracy/cost ratio than the latter and is a suitable alternative to CCSD(T).

MULTIMODE Calculations of Vibrational Spectroscopy and 1d Interconformer Tunneling Dynamics in Glycine Using a Full-Dimensional Potential Energy Surface

A full-dimensional, permutationally invariant polynomial potential energy surface for glycine recently reported (R. Conte et al., J. Chem. Phys. 2020, 153, 244301) is used with the code MULTIMODE to determine the IR absorption spectra for Conformers I and II using a new separable dipole moment function. The calculated spectra agree well with the experimental ones. The full-dimensional nature of the potential allows us also to examine dynamical results, such as tunneling rates. Remarkably, using a one-dimensional path based on the potential energy surface to estimate the tunneling rate from Conformer VI to Conformer I, good agreement is found with the recent experimental measurement. Finally a brief comparison of our potential energy surface with a recently reported sGDML one is made.

A-VCI: A flexible method to efficiently compute vibrational spectra

The adaptive vibrational configuration interaction algorithm has been introduced as a new method to efficiently reduce the dimension of the set of basis functions used in a vibrational configuration interaction process. It is based on the construction of nested bases for the discretization of the Hamiltonian operator according to a theoretical criterion that ensures the convergence of the method. In the present work, the Hamiltonian is written as a sum of products of operators. The purpose of this paper is to study the properties and outline the performance details of the main steps of the algorithm. New parameters have been incorporated to increase flexibility, and their influence has been thoroughly investigated. The robustness and reliability of the method are demonstrated for the computation of the vibrational spectrum up to 3000 cm^-1 of a widely studied 6-atom molecule (acetonitrile). Our results are compared to the most accurate up to date computation; we also give a new reference calculation for future work on this system. The algorithm has also been applied to a more challenging 7-atom molecule (ethylene oxide). The computed spectrum up to 3200 cm^-1 is the most accurate computation that exists today on such systems.

Quantum and classical IR spectra of (HCOOH)2, (DCOOH)2 and (DCOOD)2 using ab initio potential energy and dipole moment surfaces

Full-dimensional (24 modes) quantum calculation of the IR spectrum of (DCOOD)₂, and comparison with classical MD one.

An ab initio potential energy surface for the formic acid dimer: zero-point energy, selected anharmonic fundamental energies, and ground-state tunneling splitting calculated in relaxed 1–4-mode subspaces

We report a full-dimensional, permutationally invariant potential energy surface (PES) for the cyclic formic acid dimer.