Parameterized Bases for Calculating Vibrational Spectra Directly from ab Initio Data Using Rectangular Collocation

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.

The loss of the property of locality of the kernel in high-dimensional Gaussian process regression on the example of the fitting of molecular potential energy surfaces

Kernel-based methods, including Gaussian process regression (GPR) and generally kernel ridge regression, have been finding increasing use in computational chemistry, including the fitting of potential energy surfaces and density functionals in high-dimensional feature spaces. Kernels of the Matern family, such as Gaussian-like kernels (basis functions), are often used which allow imparting to them the meaning of covariance functions and formulating GPR as an estimator of the mean of a Gaussian distribution. The notion of locality of the kernel is critical for this interpretation. It is also critical to the formulation of multi-zeta type basis functions widely used in computational chemistry. We show, on the example of fitting of molecular potential energy surfaces of increasing dimensionality, the practical disappearance of the property of locality of a Gaussian-like kernel in high dimensionality. We also formulate a multi-zeta approach to the kernel and show that it significantly improves the quality of regression in low dimensionality but loses any advantage in high dimensionality, which is attributed to the loss of the property of locality.

Machine learning in computational chemistry: interplay between (non)linearity, basis sets, and dimensionality

Machine learning (ML) based methods and tools have now firmly established themselves in physical chemistry and in particular in theoretical and computational chemistry and in materials chemistry. The generality of popular ML techniques such as neural networks or kernel methods (Gaussian process and kernel ridge regression and their flavors) permitted their application to diverse problems from prediction of properties of functional materials (catalysts, solid state ionic conductors, etc.) from descriptors to the building of interatomic potentials (where ML is currently routinely used in applications) and electron density functionals. These ML techniques are assumed to have superior expressive power of nonlinear methods, and are often used "as is", with concepts such as "non-parametric" or "deep learning" used without a clear justification for their need or advantage over simpler and more robust alternatives. In this Perspective, we highlight some interrelations between popular ML techniques and traditional linear regressions and basis expansions and demonstrate that in certain regimes (such as a very high dimensionality) these approximations might collapse. We also discuss ways to recover the expressive power of a nonlinear approach and to help select hyperparameters with the help of high-dimensional model representation and to obtain elements of insight while preserving the generality of the method.

Computational vibrational spectroscopy of molecule-surface interactions: what is still difficult and what can be done about it

Interactions of molecules with solid surfaces are responsible for key functionalities for a range of currently actively pursued technologies, including heterogeneous catalysis for synthesis or decomposition of molecules, sensitization, surface...

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.

Inverse Multiquadratic Functions as the Basis for the Rectangular Collocation Method to Solve the Vibrational Schrödinger Equation

We explore the use of inverse multiquadratic (IMQ) functions as basis functions when solving the vibrational Schrödinger equation with the rectangular collocation method. The quality of the vibrational spectrum of formaldehyde (in six dimensions) is compared to that obtained using Gaussian basis functions when using different numbers of width-optimized IMQ functions. The effects of the ratio of the number of collocation points to the number of basis functions and of the choice of the IMQ exponent are studied. We show that the IMQ basis can be used with parameters where the IMQ function is not integrable. We find that the quality of the spectrum with IMQ basis functions is somewhat lower that that with a Gaussian basis when the basis size is large, and for a range of IMQ exponents. The IMQ functions are; however, advantageous when a small number of functions is used or with a small number of collocation points (e.g., when using square collocation).

Towards Accurate Spectroscopic Identification of Species at Catalytic Surfaces: Anharmonic Vibrations of Formate on AuPt

We present a calculation of vibrational frequencies of formate on the AuPt(111) surface alloy including full anharmonicity and coupling of all six intramolecular degrees of freedom. This species is a key intermediate in methanol oxidation on this material. We use a modified version of the method of Manzhos and Carrington to compute the spectrum directly from a small number (<10,000) of DFT single-point energies, bypassing the construction of a potential energy surface. This is the first such calculation for a 4-atomic species at a surface. The spectrum is obtained using rectangular collocation and a small basis set of parameterized Hermite functions. The achievable accuracy of the order of several cm-1 corresponds to the typical experimental resolution. Using normal coordinates makes the equations simple and general and easily applicable to other systems. This calculation is doable on a PC. We predict that anharmonicity and coupling lower the fundamental frequencies by dozens of cm-1, which could affect species assignment.