A multi-dimensional Smolyak collocation method in curvilinear coordinates for computing vibrational spectra

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

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.

Nonadiabatic Conical Intersection Dynamics in the Local Diabatic Representation with Strang Splitting and Fourier Basis

We develop and implement an exact conical intersection nonadiabatic wave packet dynamics method that combines the local diabatic representation, Strang splitting for the total molecular propagator, and discrete variable representation with uniform grids. By employing the local diabatic representation, this method captures all nonadiabatic effects, including nonadiabatic transitions, electronic coherences, and geometric phase. Moreover, it is free of singularities in the first and second derivative couplings and does not require the electronic wave function to be continuous with respect to the nuclear coordinates. We further show that in contrast to the adiabatic representation, the split-operator method can be directly applied to the full molecular propagator with the locally diabatic ansatz. The Fourier series, employed as the primitive nuclear basis functions, is universal and can be applied to all types of reactive coordinates. The combination of local diabatic representation, Strang splitting, and Fourier basis allows numerically exact modeling of conical intersection quantum dynamics directly with adiabatic electronic states that can be obtained from standard electronic structure computations.

Computing vibrational spectra using a new collocation method with a pruned basis and more points than basis functions: Avoiding quadrature

We present a new collocation method for computing the vibrational spectrum of a polyatomic molecule. Some form of quadrature or collocation is necessary when the potential energy surface does not have a simple form that simplifies the calculation of the potential matrix elements required to do a variational calculation. With quadrature, better accuracy is obtained by using more points than basis functions. To achieve the same advantage with collocation, we introduce a collocation method with more points than basis functions. Critically important, the method can be used with a large basis because it is incorporated into an iterative eigensolver. Previous collocation methods with more points than functions were incompatible with iterative eigensolvers. We test the new ideas by computing energy levels of molecules with as many as six atoms. We use pruned bases but expect the new method to be advantageous whenever one uses a basis for which it is not possible to find an accurate quadrature with about as many points as there are basis functions. For our test molecules, accurate energy levels are obtained even using non-optimal, simple, equally spaced points.

Using Collocation to Solve the Schrödinger Equation

We review the collocation approach to the solution of the Schrödinger equation and its uses in applications. Interrelations between collocation and other methods are highlighted. We also stress advantages and disadvantages of the rectangular collocation formulation. Using collocation makes it possible to use any, e.g. optimized, coordinates and basis functions, including nonintegrable basis functions, and provides a straightforward way of dealing with singularities in the potential. In addition, we stress that using collocation facilitates tuning the shape of basis functions and the placement of points, both of which can be done with machine-learning methods. Applications to electronic and vibrational problems are reviewed focusing on calculations for molecules on surfaces for which there are few variational calculations. Collocation has advantages when potential energy surfaces are unavailable, in particular, for molecule-surface systems, and for systems for which standard direct product quadrature grids, often used with variational methods, are costly.

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.

Smolyak Algorithm Adapted to a System-Bath Separation: Application to an Encapsulated Molecule with Large-Amplitude Motions

A Smolyak algorithm adapted to system-bath separation is proposed for rigorous quantum simulations. This technique combines a sparse grid method with the system-bath concept in a specific configuration without limitations on the form of the Hamiltonian, thus achieving a highly efficient convergence of the excitation transitions for the "system" part. Our approach provides a general way to overcome the perennial convergence problem for the standard Smolyak algorithm and enables the simulation of floppy molecules with more than a hundred degrees of freedom. The efficiency of the present method is illustrated on the simulation of H₂ caged in an sII clathrate hydrate including two kinds of cage modes. The transition energies are converged by increasing the number of normal modes of water molecules. Our results confirm the triplet splittings of both translational and rotational (j = 1) transitions of the H₂ molecule. Furthermore, they show a slight increase of the translational transitions with respect to the ones in a rigid cage.

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).

A rectangular collocation multi-configuration time-dependent Hartree (MCTDH) approach with time-independent points for calculations on general potential energy surfaces

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 CH₃ and CH₄.

Using collocation to study the vibrational dynamics of molecules

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.

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.

Using collocation and a hierarchical basis to solve the vibrational Schrödinger equation

We show that it is possible to compute vibrational energy levels of polyatomic molecules with a collocation method and a basis of products of one-dimensional harmonic oscillator functions pruned so that it does not include functions for which the indices of many of the one-dimensional functions are nonzero. Functions with many nonzero indices are coupled only by terms that depend simultaneously on many coordinates, and they are typically small. The collocation equation is derived without invoking differences of interpolation operators, which simplifies implementation of the method. This, however, requires inverting a matrix whose elements are values of the pruned basis functions at the collocation points. The collocation points are the points on a Smolyak grid whose size is equal to the size of the pruned basis set. The Smolyak grid is built from symmetrized Leja points. Because both the basis and the grid are not tensor products, the inverse is not straightforward. It can be done by using so-called hierarchical 1-D basis functions. They are defined so that the matrix whose elements are the 1-D hierarchical basis functions evaluated at points is lower triangular. We test the method by applying it to compute 100 energy levels of CH₂NH with an iterative eigensolver.

Systematically expanding nondirect product bases within the pruned multi-configuration time-dependent Hartree (MCTDH) method: A comparison with multi-layer MCTDH

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, CH₃CN, 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.