Discrete variable representations of complicated kinetic energy operators

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.

Treating linear molecules in calculations of rotation-vibration spectra

In this article, a numerical implementation of the exact kinetic energy operator (KEO) for triatomic molecules (symmetric of XY₂-type and asymmetric of YXZ-type) is presented. The implementation is based on the valence coordinates with the bisecting (XY₂-type molecules) and bond-vector (YXZ) embeddings and includes the treatment of the singularity at linear geometry. The KEO is represented in a sum-of-product form. The singularity caused by the undetermined angle at the linear configuration is resolved with the help of the associated Legendre and Laguerre polynomials used as parameterized bending basis functions in the finite basis set representation. The exact KEO implementation is combined with the variational solver theoretical rovibrational energies, equipped with a general automatic symmetry-adaptation procedure and efficient basis step contraction schemes, providing a powerful computational solver of triatomic molecules for accurate computations of highly excited ro-vibrational spectra. The advantages of different basis set choices are discussed. Examples of specific applications for computing hot spectra of linear molecules are given.

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.

Discrete variable representation of the Smoluchowski equation using a sinc basis set

We present a new general framework for solving the monodimensional Smoluchowski equation using a discrete variable representation (DVR) based on the so called sinc basis set. The reliability of our implementation is assessed by comparing the convergence of diffusive operator eigenvalues calculated using our method and using a simple finite difference scheme for some model diffusive problems. The results here presented open encouraging possibilities for dealing with more complicated systems, where additional coordinate dependent terms in the equation or multidimensional treatments are needed and traditional methods often become unfeasible.

Spectral convergence of the quadrature discretization method in the solution of the Schrödinger and Fokker-Planck equations: Comparison with sinc methods

Spectral methods based on nonclassical polynomials and Fourier basis functions or sinc interpolation techniques are compared for several eigenvalue problems for the Fokker-Planck and Schrodinger equations. A very rapid spectral convergence of the eigenvalues versus the number of quadrature points is obtained with the quadrature discretization method (QDM) and the appropriate choice of the weight function. The QDM is a pseudospectral method and the rate of convergence is compared with the sinc method reported by Wei [J. Chem. Phys., 110, 8930 (1999)]. In general, sinc methods based on Fourier basis functions with a uniform grid provide a much slower convergence. The paper considers Fokker-Planck equations (and analogous Schrodinger equations) for the thermalization of electrons in atomic moderators and for a quartic potential employed to model chemical reactions. The solution of the Schrodinger equation for the vibrational states of I2 with a Morse potential is also considered.

Excited States of Weakly Bound Bosonic Clusters: Discrete Variable Representation and Quantum Monte Carlo

We compare two approaches for the accurate calculation of the energy levels of weakly bound boson trimers. The first approach is based on correlation function Monte Carlo employing optimized trial functions, while the second approach is based on the discrete variable representation. A trimer with atoms of half the mass of neon is used as a test problem for benchmark calculations. The two approaches yield identical results, within error bars, for all the J = 0 energy levels below the dissociation threshold. The relative merits of the two techniques are discussed, and a perspective is given for extension to larger clusters.

Treating singularities present in the Sutcliffe-Tennyson vibrational Hamiltonian in orthogonal internal coordinates

Two methods are developed, when solving the related time-independent Schrodinger equation (TISE), to cope with the singular terms of the vibrational kinetic energy operator of a triatomic molecule given in orthogonal internal coordinates. The first method provides a mathematically correct treatment of all singular terms. The vibrational eigenfunctions are approximated by linear combinations of functions of a three-dimensional nondirect-product basis, where basis functions are formed by coupling Bessel-DVR functions, where DVR stands for discrete variable representation, depending on distance-type coordinates and Legendre polynomials depending on angle bending. In the second method one of the singular terms related to a distance-type coordinate, deemed to be unimportant for spectroscopic applications, is given no special treatment. Here the basis set is obtained by taking the direct product of a one-dimensional DVR basis with a two-dimensional nondirect-product basis, the latter formed by coupling Bessel-DVR functions and Legendre polynomials. With the basis functions defined, matrix representations of the TISE are set up and solved numerically to obtain the vibrational energy levels of H3+. The numerical calculations show that the first method treating all singularities is computationally inefficient, while the second method treating properly only the singularities having physical importance is quite efficient.

Energy levels and wave functions of weakly-bound 4Hex 20NeyH (x+y=2) systems using Pekeris coordinates and a symmetry-adapted Lanczos approach

Energy levels and wave functions of floppy triatomic rare gas hydrides are calculated using a Pekeris coordinate system and the importance of various triangular configurations is assessed through the calculation of reduced distribution functions and relative weights. The calculations are performed using a symmetry-adapted Lanczos recursion within the discrete variable representation. For the 4He2H- anion, the present results are compared with those obtained from calculations based on other methods, and the accuracy of the present method is discussed. Calculations are also performed for the case of 4He2H and 20Ne2H, as well as for the mixed 4He 20NeH neutrals. Our results show that no bound states are found for 4He2H while only one bound state is found for both the 20Ne2H and 4He 20NeH complexes. Interestingly, a very important and common property of these systems is that there is a significant contribution from linear configurations to their bound states.

First principles simulation of the UV absorption spectrum of ethylene using the vertical Franck-Condon approach

A new method which we refer to as vertical Franck-Condon is proposed to calculate electronic absorption spectra of polyatomic molecules. In accord with the short-time picture of spectroscopy, the excited-state potential energy surface is expanded at the ground-state equilibrium geometry and the focus of the approach is more on the overall shape of the spectrum and the positions of the band maxima, rather than the precise position of the 0-0 lines. The Born-Oppenheimer approximation and the separability of the excited-state potential energy surface along the excited-state normal mode coordinates are assumed. However, the potential surface is not necessarily approximated as harmonic oscillator potentials along the individual normal modes. Instead, depending upon the nature of the potential surface along a particular normal mode, it is treated either in the harmonic approximation or the full one-dimensional potential is considered along this mode. The vertical Franck-Condon approach is applicable therefore even in cases where the excited state potential energy surface is highly anharmonic and the conventional harmonic Franck-Condon approach is inadequate. As an application of the method, the ultraviolet spectrum of ethylene between 6.2 eV (50,000 cm(-1)) and 8.7 eV (70,000 cm(-1)) is simulated, using the Similarity Transformed Equation of Motion Coupled-Cluster method to describe the required features of the potential energy surfaces. The spectrum is shown to be a result of sharp doublet structures stemming from the pi --> 3s (Rydberg) state superimposed on top of a broad band resulting from the pi --> pi* (valence) state. For the Rydberg state, the symmetric C=C stretch and the torsion mode contribute to the spectrum, while the broad valence band results from excitation into the C=C stretch, CH2 scissors, and the torsion mode. For both states, the potential along the torsion mode is highly anharmonic and the full treatment of the potential along this mode in the vertical Franck-Condon method is required.