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

by   Jason Kaye, et al.

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.


page 12

page 29

page 32

page 34


Eliminating artificial boundary conditions in time-dependent density functional theory using Fourier contour deformation

We present an efficient method for propagating the time-dependent Kohn-S...

Pseudospectral computational methods for the time-dependent Dirac equation in static curved spaces

Pseudospectral numerical schemes for solving the Dirac equation in gener...

A space-time discretization of a nonlinear peridynamic model on a 2D lamina

Peridynamics is a nonlocal theory for dynamic fracture analysis consisti...

A fast spectral method for electrostatics in doubly-periodic slit channels

We develop a fast method for computing the electrostatic energy and forc...

Multiple-scattering frequency-time hybrid solver for the wave equation in interior domains

This paper proposes a frequency-time hybrid solver for the time-dependen...

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

We introduce a Fourier-Bessel-based spectral solver for Cauchy problems ...

Fourier smoothed pre-corrected trapezoidal rule for solution of Lippmann-Schwinger integral equation

For the numerical solution of the Lippmann-Schwinger equation, while the...

Please sign up or login with your details

Forgot password? Click here to reset