The quadrature discretization method (QDM) in the solution of the Schrödinger equation with nonclassical basis functions

Klein-Gordon equation on a Lagrange mesh

The Lagrange-mesh method is an approximate variational method which provides accurate solutions of the Schrödinger equation for bound-state and scattering few-body problems. The stationary Klein-Gordon equation depends quadratically on the energy. For a central potential, it is solved on a Lagrange-Laguerre mesh by iteration. Results are tested with the Coulomb potential for which exact solutions are available. A high accuracy is obtained with a rather small number of mesh points. For various potentials and levels, few iterations provide accurate energies and mean values in short computer times. Analytical expressions of the wave functions are available.

Algebraic DVR Approaches Applied to Piecewise Potentials: Symmetry and Degeneracy

Algebraic discrete variable representation (DVR) methods that have been recently proposed are applied to describe 1D and 2D piecewise potentials. First, it is shown that it is possible to use a DVR approach to describe 1D square well potentials testing the wave functions with exact results. Thereafter, Morse and Pöschl-Teller (PT) potentials are described with multistep piecewise potentials in order to explore the sensibility of the potential to reproduce the transition from a pure square well energy pattern to an anharmonic energy spectrum. Once the properties of the different algebraic DVR approaches are identified, the 2D square potential as a function of the potential depth is studied. We show that this system displays natural degeneracy, accidental degeneracy and systematic accidental degeneracy. The latter appears only for a confined potential, where the symmetry group is identified and irreducible representations are constructed. One particle confined in a rectangular well potential with commensurate sides is also analyzed. It is proved that the systematic accidental degeneracy appearing in this system is removed for finite potential depth.

Algebraic DVR Approaches Applied to Describe the Stark Effect

Two algebraic approaches based on a discrete variable representation are introduced and applied to describe the Stark effect in the non-relativistic Hydrogen atom. One approach consists of considering an algebraic representation of a cutoff 3D harmonic oscillator where the matrix representation of the operators r2 and p2 are diagonalized to define useful bases to obtain the matrix representation of the Hamiltonian in a simple form in terms of diagonal matrices. The second approach is based on the U(4) dynamical algebra which consists of the addition of a scalar boson to the 3D harmonic oscillator space keeping constant the total number of bosons. This allows the kets associated with the different subgroup chains to be linked to energy, coordinate and momentum representations, whose involved branching rules define the discrete variable representation. Both methods, although originating from the harmonic oscillator basis, provide different convergence tests due to the fact that the associated discrete bases turn out to be different. These approaches provide powerful tools to obtain the matrix representation of 3D general Hamiltonians in a simple form. In particular, the Hydrogen atom interacting with a static electric field is described. To accomplish this task, the diagonalization of the exact matrix representation of the Hamiltonian is carried out. Particular attention is paid to the subspaces associated with the quantum numbers n=2,3 with m=0.

Diffusion in a bistable system: The eigenvalue spectrum of the Fokker-Planck operator and Kramers' reaction rate theory

The time-dependent solution of the Fokker-Planck equation with bistable potentials is considered in terms of the eigenfunctions and eigenvalues of the linear Fokker-Planck operator. The Fokker-Planck equation is the high friction limit of the corresponding Kramers' equation. Two different potentials are considered defined with a constant diffusion coefficient, ε, and position-dependent drift coefficients. The smallest nonzero eigenvalue of the Fokker-Planck operator, λ_{1}, provides the long-time rate coefficient for the transformation of the different species in the two stable states. A novel pseudospectral method with nonclassical polynomials is applied to this class of systems. The convergence of the eigenvalues and eigenfunctions of the Fokker-Planck operator versus the number of basis functions is studied and compared with previous results. The results are consistent with Kramers' theory, and a linear relationship between lnλ_{1} and 1/ε for sufficiently small ε values is verified. A comparison with analytic approximations to λ_{1} is provided.

Variational properties of the discrete variable representation: discrete variable representation via effective operators

A variational finite basis representation/discrete variable representation (FBR/DVR) Hamiltonian operator has been introduced. By calculating its matrix elements exactly one obtains, depending on the choice of the basis set, either a variational FBR or a variational DVR. The domain of grid points on which the FBR/DVR is variational has been shown to consist of the subsets of the set of grid points one obtains by diagonalizing commuting variational basis representations of the coordinate operators. The variational property implies that the optimal of the subsets of a fixed number of points, i.e., the subset which gives the possible highest accuracy eigenpairs, gives the DVR of the smallest trace. The symmetry properties of the variational FBR/DVR Hamiltonian operator are analyzed and methods to incorporate symmetry into FBR/DVR calculations are discussed. It is shown how the Fourier-basis FBR/DVR suitable to solving periodic systems arise within the theory presented. Numerical examples are given to illustrate the theoretical results. The use of variational effective Hamiltonian and coordinate operators has been instrumental in this study. They have been introduced in a novel way by exploiting quasi-Hermiticity.

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.

Application of Coulomb wave function discrete variable representation to atomic systems in strong laser fields

We present an efficient and accurate grid method for solving the time-dependent Schrodinger equation for an atomic system interacting with an intense laser pulse. Instead of the usual finite difference (FD) method, the radial coordinate is discretized using the discrete variable representation (DVR) constructed from Coulomb wave functions. For an accurate description of the ionization dynamics of atomic systems, the Coulomb wave function discrete variable representation (CWDVR) method needs three to ten times fewer grid points than the FD method. The resultant grid points of the CWDVR are distributed unevenly so that one has a finer grid near the origin and a coarser one at larger distances. The other important advantage of the CWDVR method is that it treats the Coulomb singularity accurately and gives a good representation of continuum wave functions. The time propagation of the wave function is implemented using the well-known Arnoldi method. As examples, the present method is applied to multiphoton ionization of both the H atom and the H(-) ion in intense laser fields. The short-time excitation and ionization dynamics of H by an abruptly introduced static electric field is also investigated. For a wide range of field parameters, ionization rates calculated using the present method are in excellent agreement with those from other accurate theoretical calculations.

Optimal grids for generalized finite basis and discrete variable representations: definition and method of calculation

The method of optimal generalized finite basis and discrete variable representations (FBR and DVR) generalizes the standard, Gaussian quadrature grid-classical orthonormal polynomial basis-based FBR/DVR method to general sets of grid points and to general, nondirect product, and/or nonpolynomial bases. Here, it is shown how an optimal set of grid points can be obtained for an optimal generalized FBR/DVR calculation with a given truncated basis. Basis set optimized and potential optimized grids are defined. The optimized grids are shown to minimize a function of grid points derived by relating the optimal generalized FBR of a Hamiltonian operator to a non-Hermitian effective Hamiltonian matrix. Locating the global minimum of this function can be reduced to finding the zeros of a function in the case of one dimensional problems and to solving a system of D nonlinear equations repeatedly in the case of D>1 dimensional problems when there is an equal number of grid points and basis functions. Gaussian quadrature grids are shown to be basis optimized grids. It is demonstrated by a numerical example that an optimal generalized FBR/DVR calculation of the eigenvalues of a Hamiltonian operator with potential optimized grids can have orders of magnitude higher accuracy than a variational calculation employing the same truncated basis. Nevertheless, for numerical integration with the optimal generalized FBR quadrature rule basis optimized grids are the best among grids of the same number of points. The notions of Gaussian quadrature and Gaussian quadrature accuracy are extended to general, multivariable basis functions.

Solving the bound-state Schrodinger equation by reproducing kernel interpolation

Based on reproducing kernel Hilbert space theory and radial basis approximation theory, a grid method is developed for numerically solving the N-dimensional bound-state Schrodinger equation. Central to the method is the construction of an appropriate bounded reproducing kernel (RK) Lambda(alpha)( ||r ||) from the linear operator -nabla(2)(r)+lambda(2) where nabla(2)(r) is the N-dimensional Laplacian, lambda>0 is a parameter related to the binding energy of the system under study, and the real number alpha>N. The proposed (Sobolev) RK Lambda(alpha)(r,r(')) is shown to be a positive-definite radial basis function, and it matches the asymptotic solutions of the bound-state Schrodinger equation. Numerical tests for the one-dimensional (1D) Morse potential and 2D Henon-Heiles potential reveal that the method can accurately and efficiently yield all the energy levels up to the dissociation limit. Comparisons are also made with the results based on the distributed Gaussian basis method in the 1D case as well as the distributed approximating functional method in both 1D and 2D cases.