Abstract:We introduce a numerical method for the solution of the time-dependent Schrodinger equation with a smooth potential, based on its reformulation as a Volterra integral equation. We present versions of the method both for periodic boundary conditions, and for free space problems with compactly supported initial data and potential. A spatially uniform electric field may be included, making the solver applicable to simulations of light-matter interaction. The primary computational challenge in using the Volterra formulation is the application of a space-time history dependent integral operator. This may be accomplished by projecting the solution onto a set of Fourier modes, and updating their coefficients from one time step to the next by a simple recurrence. In the periodic case, the modes are those of the usual Fourier series, and the fast Fourier transform (FFT) is used to alternate between physical and frequency domain grids. In the free space case, the oscillatory behavior of the spectral Green's function leads us to use a set of complex-frequency Fourier modes obtained by discretizing a contour deformation of the inverse Fourier transform, and we develop a corresponding fast transform based on the FFT. Our approach is related to pseudo-spectral methods, but applied to an integral rather than the usual differential formulation. This has several advantages: it avoids the need for artificial boundary conditions, admits simple, inexpensive high-order implicit time marching schemes, and naturally includes time-dependent potentials. We present examples in one and two dimensions showing spectral accuracy in space and eighth-order accuracy in time for both periodic and free space problems.

Combined finite difference and spectral methods for the numerical solution of hyperbolic equation with an integral condition

Efficient Chebyshev Spectral Method for Solving Linear Elliptic Pdes Using Quasi-Inverse Technique

A numerical method for solving hyperbolic partial differential equations with piecewise constant arguments and variable coefficients

An integral equation method for the advection-diffusion equation on time-dependent domains in the plane

Estimation of a heat source in a parabolic equation with nonlocal Wentzell boundary condition using a spectral technique

Application of a hybrid pseudospectral method to a new two-dimensional multi-term mixed sub-diffusion and wave-diffusion equation of fractional order

A meshless method for numerical solution of a linear hyperbolic equation with variable coefficients in two space dimensions

METHOD OF LINES AND FINITE DIFFERENCE SCHEMES WITH THE EXACT SPECTRUM FOR SOLUTION THE HYPERBOLIC HEAT CONDUCTION EQUATION

Application of spectral element method for solving Sobolev equations with error estimation

Numerical Solution of Partial Integro-Differential Equations by using Projection Method

The Fourier Spectral Method for the Cahn-Hilliard Equation

Second-Order Finite Difference/Spectral Element Formulation for Solving the Fractional Advection-Diffusion Equation

Solving Diffusion and Wave Equations Meshlessly via Helmholtz Equations

A Tailored Finite Point Method for the Helmholtz Equation with High Wave Numbers in Heterogeneous Medium

Approximating Continuous Spectra of Hyperbolic Systems with Summation-by-Parts Finite Difference Operators

Numerical solution of fractional Fredholm integro-differential equations by spectral method with fractional basis functions

A high-order integral equation-based solver for the time-dependent Schrodinger equation

A new spectral element method for numerical solution of partial differential equations on annular-type regions

On spectral numerical method for variable-order partial differential equations

Spectral solver for Cauchy problems in polar coordinates using discrete Hankel transforms

The application of integral equation methods to the numerical solution of some exterior boundary-value problems