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.

Solving the two-dimensional time-dependent Schrödinger equation using the Sinc collocation method and double exponential transformations

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

High-oder Symplectic FDTD Scheme for Solving Time-Dependent Schrdinger Equation

Exponential Collocation Methods Based on Continuous Finite Element Approximations for Efficiently Solving the Cubic Schrödinger Equation

On numerical solutions of the time-dependent Schrödinger equation

Exponential Collocation Methods for the Cubic Schrödinger Equation

A Time Two-Mesh Compact Difference Method for the One-Dimensional Nonlinear Schrödinger Equation

Several Numerical Simulation Methods for the Time-Dependent Schrödinger Equation

New optical soliton solutions of the popularized anti-cubic nonlinear Schrödinger equation versus its numerical treatment

Numerical Solutions of Schrödinger Equations in ℝ3

Computational Solution of the Time-Fractional Schrödinger Equation by Using Trigonometric B-Spline Collocation Method

Accurate Solution of the Nonlinear Schrödinger Equation via Conservative Multiple-Relaxation ImEx Methods

B-spline finite element studies of the non-linear Schrödinger equation

Retrieval of soliton solutions for 4th-order (2+1)-dimensional Schrödinger equation with higher-order odd and even terms by modified Sardar sub-equation method

Analytical and numerical aspects of certain nonlinear evolution equations. II. Numerical, nonlinear Schrödinger equation

Tenth-Order Accurate Numerical Method for Solving the Time-Dependent Schrödinger Equation

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

The Solutions for Nonlinear Schrödinger Equations

Efficiently and accurately simulating multi-dimensional M-coupled nonlinear Schrödinger equations with fourth-order time integrators and Fourier spectral method

Numerical solutions of Korteweg de Vries and Korteweg de Vries-Burger's equations using computer programming

Hermite Neural Network Simulation for Solving the 2D Schrodinger Equation