Related papers: Fourth Order Gradient Symplectic Integrator Method…
We review and further develop the recently introduced numerical approach for scattering calculations based on a so called pseudo-time Schroedinger equation, which is in turn a modification of the damped Chebyshev polynomial expansion…
A novel symplectic algorithm is proposed to solve the Maxwell-Schr\"odinger (M-S) system for investigating light-matter interaction. Using the fourth-order symplectic integration and fourth-order collocated differences,…
We review a new iterative procedure to solve the low-lying states of the Schroedinger equation, done in collaboration with Richard Friedberg. For the groundstate energy, the $n^{th}$ order iterative energy is bounded by a finite limit,…
We show that a pseudospectral representation of the wavefunction using multiple spatial domains of variable size yields a highly accurate, yet efficient method to solve the time-dependent Schr\"odinger equation. The overall spatial domain…
We present a novel approach for solving the time-dependent Schr\"{o}dinger equation (TDSE). The method we propose converts the TDSE to an equivalent Volterra integral equation; introducing a global Lagrange interpolation of the integrand…
We consider the linear time-dependent Schr\"odinger equation with a time-dependent smooth potential on an unbounded domain. A Galerkin spectral method with a tensor-product Hermite basis is used as a discretization in space. Discretizing…
In this paper, we propose a first-order Fourier integrator for solving the cubic nonlinear Schr\"odinger equation in one dimension. The scheme is explicit and can be implemented using the fast Fourier transform. By a rigorous analysis, we…
This paper introduces an efficient high-order numerical method for solving the 1D stationary Schr\"odinger equation in the highly oscillatory regime. Building upon the ideas from [Arnold, Ben Abdallah, Negulescu, SIAM J. Numer. Anal.,…
We have proposed new algorithms for the numerical integration of the equations of motion for classical spin systems. In close analogy to symplectic integrators for Hamiltonian equations of motion used in Molecular Dynamics these algorithms…
The Schr\"odingerization method converts linear partial and ordinary differential equations with non-unitary dynamics into systems of Schr\"odinger-type equations with unitary evolution. It does so via the so-called warped phase…
A propagation method for time-dependent Schr\"odinger equations with an explicitly time-dependent Hamiltonian is developed where time ordering is achieved iteratively. The explicit time-dependence of the time-dependent Schr\"odinger…
A parareal algorithm based on an exponential $\theta$-scheme is proposed for the stochastic Schr\"odinger equation with weak damping and additive noise. It proceeds as a two-level temporal parallelizable integrator with the exponential…
Exponential integrators are time stepping schemes which exactly solve the linear part of a semilinear ODE system. This class of schemes requires the approxima- tion of a matrix exponential in every step, and one successful modern method is…
High-order methods gain increased attention in computational fluid dynamics. However, due to the time step restrictions arising from the semi-implicit time stepping for the incompressible case, the potential advantage of these methods…
The formulation of the time-dependent Schrodinger equation in terms of coupled-cluster theory is outlined, with emphasis on the bivariational framework and its classical Hamiltonian structure. An indefinite inner product is introduced,…
We introduce a new iterative method for computing solutions of elliptic equations with random rapidly oscillating coefficients. Similarly to a multigrid method, each step of the iteration involves different computations meant to address…
We present an approach to numerically solving the time-dependent Schroedinger equation and other parabolic equations by the split-step technique with fast Fourier transform, which suppresses the backreflection of waves from the grid…
Many force-gradient explicit symplectic integration algorithms have been designed for the Hamiltonian $H=T (\mathbf{p})+V(\mathbf{q})$ with kinetic energy $T(\mathbf{p})=\mathbf{p}^2/2$ in the existing references. When the force-gradient…
Numeric modeling of electromagnetics and acoustics frequently entails matrix-vector multiplication with block Toeplitz structure. When the corresponding block Toeplitz matrix is not highly sparse, e.g. when considering the electromagnetic…
We discuss a new approach to solve the low lying states of the Schroedinger equation. For a fairly large class of problems, this new approach leads to convergent iterative solutions, in contrast to perturbative series expansions. These…