Search

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. [ABSTRACT FROM AUTHOR]

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 CH3 and CH4. [ABSTRACT FROM AUTHOR]

Although very useful, the original multi-configuration time-dependent Hartree (MCTDH) method has two weaknesses: (1) its cost scales exponentially with the number of atoms in the system; (2) the standard MCTDH implementation requires that the potential energy surface (PES) be in the sum-of-product (SOP) form in order to reduce the cost of computing integrals in the MCTDH basis. One way to deal with (1) is to lump coordinates into groups. This is mode combination (MC). One way to deal with (2) is to reformulate MCTDH using collocation so that there are no integrals. In this paper, we combine MC and collocation to formulate a MC collocation multi-configuration time-dependent Hartree (MC-C-MCTDH) method. In practice, its cost does not scale exponentially with the number of atoms, and it can be used with any general PES; the PES need not be an SOP and need not have a special form. No integrals and, hence, no quadratures are necessary. We demonstrate the accuracy and efficiency of the new method by computing vibrational energy eigenstates of methyl radical, methane, and acetonitrile. To do this, we use MC-C-MCTDH with a variant of improved relaxation, derived by evaluating a residual at points. Because the MC basis functions are multivariate, collocation points in multi-dimensional spaces are required. We use two types of collocation points: (1) discrete variable representation-like points obtained from (approximate) simultaneous diagonalization of matrices and (2) Leja points, which are known to be good interpolation points, determined from a generalized recipe suitable for any basis. [ABSTRACT FROM AUTHOR]

Standard multiconfiguration time-dependent Hartree (MCTDH) calculations use a direct product basis and rely on the potential being a sum of products (SOPs). The size of the direct product MCTDH basis scales exponentially with the number of atoms. Accurate potentials may not be SOPs. We introduce an MCTDH approach that uses a pruned basis and a collocation grid. Pruning the basis significantly reduces its size. Collocation makes it possible to do calculations using a potential that is not a SOP. The collocation point set is a Smolyak grid. Strategies using pruned MCTDH bases already exist, but they work only if the potential is a SOP. Strategies for using MCTDH with collocation also exist, but they work only if the MCTDH basis is a direct product. In this paper, we combine a pruned basis with collocation. This makes it possible to mitigate the direct-product basis size problem and do calculations when the potential is not a SOP. Because collocation is used, there are no integrals and no need for quadrature. All required matrix-vector products can be evaluated sequentially. We use nested sets of collocation points and hierarchical basis functions. They permit efficient inversion of the (large) matrix whose elements are basis functions evaluated at points, which is necessary to transform values of functions at points to basis coefficients. The inversion technique could be used outside of chemical physics. We confirm the validity of this new pruned, collocation-based (PC-)MCTDH approach by calculating the first 50 vibrational eigenenergies of CH2NH. [ABSTRACT FROM AUTHOR]

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, CH3CN, 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. [ABSTRACT FROM AUTHOR]

In this paper, we propose a pruned, nondirect product multi-configuration time dependent Hartree (MCTDH) method for solving the Schrödinger equation. MCTDH uses optimized 1D basis functions, called single particle functions, but the size of the standard direct product MCTDH basis scales exponentially with D, the number of coordinates. We compare the pruned approach to standard MCTDH calculations for basis sizes small enough that the latter are possible and demonstrate that pruning the basis reduces the CPU cost of computing vibrational energy levels of acetonitrile (D = 12) by more than two orders of magnitude. Using the pruned method, it is possible to do calculations with larger bases, for which the cost of standard MCTDH calculations is prohibitive. Pruning the basis complicates the evaluation of matrix-vector products. In this paper, they are done term by term for a sum-of-products Hamiltonian. When no attempt is made to exploit the fact that matrices representing some of the factors of a term are identity matrices, one needs only to carefully constrain indices. In this paper, we develop new ideas that make it possible to further reduce the CPU time by exploiting identity matrices. [ABSTRACT FROM AUTHOR]