Related papers: Explicit schemes for time propagating many-body wa…
Efficient and accurate numerical propagation of the time dependent Schroedinger equation is a problem with applications across a wide range of physics. This paper develops an efficient, trivially parallelizeable method for relaxing a trial…
In this study, a variety of methods are tested and compared for the numerical solution of the Schr\"odinger equation for few-body systems with explicitely time-dependent Hamiltonians, with the aim to find the optimal one. The configuration…
We present quantum algorithms for solving two problems regarding stochastic processes. The first algorithm prepares the thermal Gibbs state of a quantum system and runs in time almost linear in $\sqrt{N \beta/{\cal Z}}$ and polynomial in…
This work is devoted to the numerical simulation of nonlinear Schr\"odinger and Klein-Gordon equations. We present a general strategy to construct numerical schemes which are uniformly accurate with respect to the oscillation frequency.…
We present a novel numerical method and algorithm for the solution of the 3D axially symmetric time-dependent Schr\"odinger equation in cylindrical coordinates, involving singular Coulomb potential terms besides a smooth time-dependent…
Large-scale linear, time-invariant (LTI) dynamical systems are widely used to characterize complicated physical phenomena. We propose a two-stage algorithm to reduce the order of a large-scale LTI system given samples of its transfer…
Study of far-from-equilibrium thermalization dynamics in quantum materials, including the dynamics of different types of quasiparticles, is becoming increasingly crucial. However, the inherent complexity of either the full quantum…
A new parallel-in-time iterative method is proposed for solving the homogeneous second-order wave equation. The new method involves a coarse scale propagator, allowing for larger time steps, and a fine scale propagator which fully resolves…
We consider the initial/boundary value problem for the fractional diffusion and diffusion-wave equations involving a Caputo fractional derivative in time. We develop two "simple" fully discrete schemes based on the Galerkin finite element…
We provide and analyze the high order algorithms for the model describing the functional distributions of particles performing anomalous motion with power-law jump length and tempered power-law waiting time. The model is derived in [Wu,…
Global instability analysis of flows is often performed via time-stepping methods, based on the Arnoldi algorithm. When setting up these methods, several computational parameters must be chosen, which affect intrinsic errors of the…
Exponential time differencing methods is a power tool for high-performance numerical simulation of computationally challenging problems in condensed matter physics, fluid dynamics, chemical and biological physics, where mathematical models…
We discuss double ionization of atoms in strong laser pulses using a reduced dimensionality model. Following the insights obtained from an analysis of the classical mechanics of the process, we confine each electron to move along the lines…
In this paper, a linearized fully discrete scheme is proposed to solve the two-dimensional nonlinear time fractional Schr\"odinger equation with weakly singular solutions, which is constructed by using L1 scheme for Caputo fractional…
Since the kinetic and the potential energy term of the real time nonlinear Schr\"odinger equation can each be solved exactly, the entire equation can be solved to any order via splitting algorithms. We verified the fourth-order convergence…
The threshold dynamics algorithm of Merriman, Bence, and Osher is only first order accurate in the two-phase setting. Its accuracy degrades further to half order in the multi-phase setting, a shortcoming it has in common with other related,…
This paper constructs a predictor-corrector technique with orthogonal spline collocation finite element method for simulating a FitzHugh-Nagumo system subject to suitable initial and boundary conditions. The developed computational…
Most numerical methods for time integration use real-valued time steps. Complex time steps, however, can provide an additional degree of freedom, as we can select the magnitude of the time step in both the real and imaginary directions. We…
In this paper, the trial function method is employed to find the exact solutions for high-order nonlinear Schr\"odinger equations with time-dependent coefficients. This system describes the propagation of ultrashort light pulses in…
The eigenvalue-function pair of the 3D Schr\"odinger equation can be efficiently computed by use of high order, imaginary time propagators. Due to the diffusion character of the kinetic energy operator in imaginary time, algorithms…