430 research outputs found
Multi-Symplectic Simulation on Soliton-Collision for Nonlinear Perturbed Schrödinger Equation
Seeking solitary wave solutions and revealing their interactional characteristics for nonlinear evolution equations help us lot to comprehend the motion laws of the microparticles. As a local nonlinear dynamic behavior, the soliton-collision is difficult to be reproduced numerically. In this paper, the soliton-collision process in the nonlinear perturbed Schrödinger equation is simulated employing the multi-symplectic method. The multi-symplectic formulations are derived including the multi-symplectic form and three local conservation laws of the nonlinear perturbed Schrödinger equation. Employing the implicit midpoint rule, we construct a multi-symplectic scheme, which is equivalent to the Preissmann box scheme, for the nonlinear perturbed Schrödinger equation. The elegant structure-preserving properties of the multi-symplectic scheme are illustrated by the tiny maximum absolute residual of the discrete multi-symplectic structure at each time step in the numerical simulations. The effects of the perturbation strength on the soliton-collision in the nonlinear perturbed Schrödinger equation are reported in the numerical results in detail
Symplectic structure-preserving integrators for the two-dimensional Gross–Pitaevskii equation for BEC
AbstractSymplectic integrators have been developed for solving the two-dimensional Gross–Pitaevskii equation. The equation is transformed into a Hamiltonian form with symplectic structure. Then, symplectic integrators, including the midpoint rule, and a splitting symplectic scheme are developed for treating this equation. It is shown that the proposed codes fulfill the discrete charge conservation law. Furthermore, the global error of the numerical solution is theoretically estimated. The theoretical analysis is supported by some numerical simulations
A time-splitting pseudospectral method for the solution of the Gross-Pitaevskii equations using spherical harmonics with generalised-Laguerre basis functions
We present a method for numerically solving a Gross-Pitaevskii system of equations with a harmonic and a toroidal external potential that governs the dynamics of one- and two-component Bose-Einstein condensates. The method we develop maintains spectral accuracy by employing Fourier or spherical harmonics in the angular coordinates combined with generalised-Laguerre basis functions in the radial direction. Using an error analysis, we show that the method presented leads to more accurate results than one based on a sine transform in the radial direction when combined with a time-splitting method for integrating the equations forward in time. In contrast to a number of previous studies, no assumptions of radial or cylindrical symmetry is assumed allowing the method to be applied to 2D and 3D time-dependent simulations. This is accomplished by developing an efficient algorithm that accurately performs the generalised-Laguerre transforms of rotating Bose-Einstein condensates for different orders of the Laguerre polynomials. Using this spatial discretisation together with a second order Strang time-splitting method, we illustrate the scheme on a number of 2D and 3D computations of the ground state of a non-rotating and rotating condensate. Comparisons between previously derived theoretical results for these ground state solutions and our numerical computations show excellent agreement for these benchmark problems. The method is further applied to simulate a number of time-dependent problems including the Kelvin-Helmholtz instability in a two-component rotating condensate and the motion of quantised vortices in a 3D condensate
Fourth order real space solver for the time-dependent Schr\"odinger equation with singular Coulomb potential
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 potential. We use fourth order finite difference real space
discretization, with special formulae for the arising Neumann and Robin
boundary conditions along the symmetry axis. Our propagation algorithm is based
on merging the method of the split-operator approximation of the exponential
operator with the implicit equations of second order cylindrical 2D
Crank-Nicolson scheme. We call this method hybrid splitting scheme because it
inherits both the speed of the split step finite difference schemes and the
robustness of the full Crank-Nicolson scheme. Based on a thorough error
analysis, we verified both the fourth order accuracy of the spatial
discretization in the optimal spatial step size range, and the fourth order
scaling with the time step in the case of proper high order expressions of the
split-operator. We demonstrate the performance and high accuracy of our hybrid
splitting scheme by simulating optical tunneling from a hydrogen atom due to a
few-cycle laser pulse with linear polarization
R-adaptive multisymplectic and variational integrators
Moving mesh methods (also called r-adaptive methods) are space-adaptive
strategies used for the numerical simulation of time-dependent partial
differential equations. These methods keep the total number of mesh points
fixed during the simulation, but redistribute them over time to follow the
areas where a higher mesh point density is required. There are a very limited
number of moving mesh methods designed for solving field-theoretic partial
differential equations, and the numerical analysis of the resulting schemes is
challenging. In this paper we present two ways to construct r-adaptive
variational and multisymplectic integrators for (1+1)-dimensional Lagrangian
field theories. The first method uses a variational discretization of the
physical equations and the mesh equations are then coupled in a way typical of
the existing r-adaptive schemes. The second method treats the mesh points as
pseudo-particles and incorporates their dynamics directly into the variational
principle. A user-specified adaptation strategy is then enforced through
Lagrange multipliers as a constraint on the dynamics of both the physical field
and the mesh points. We discuss the advantages and limitations of our methods.
Numerical results for the Sine-Gordon equation are also presented.Comment: 65 pages, 13 figure
- …