MathSciDoc: An Archive for Mathematician ∫

We present a uniformly first order accurate numerical method for solving the Klein-Gordon-Zakharov (KGZ) system with two dimensionless parameters 0 < ε ≤ 1 and 0 <γ ≤ 1, which are inversely proportional to the plasma frequency and the acoustic speed, respectively. In the simultaneous high-plasma-frequency and subsonic limit regime, i.e. ε < γ → 0+, the KGZ system collapses to a cubic Schrödinger equation, and the solution propagates waves with O ( ε2)-wavelength in time and meanwhile contains rapid outgoing initial layers with speed O (1/ γ ) in space due to the incompatibility of the initial data. By presenting a multiscale decomposition of the KGZ system, we propose a multiscale time integrator Fourier pseudospectral method which is explicit, efficient and uniformly accurate for solving the KGZ system for all 0 < ε < γ ≤ 1. Numerical results are reported to show the efficiency and accuracy of scheme. Finally, the method is applied to investigate the convergence rates of the KGZ system to its limiting models when ε < γ → 0+.

A Deuflhard-type exponential integrator sine pseudospectral (DEI-SP) method is proposed and analyzed for solving the generalized improved Boussinesq (GIBq) equation. The numerical scheme is based on a second-order exponential integrator for time integration and a sine pseudospectral discretization in space. Rigorous analysis and abundant experiments show that the method converges quadratically and spectrally in time and space, respectively. Finally the DEI-SP method is applied to investigate the complicated and interesting long-time dynamics of the GIBq equation.

This paper is devoted to developing and analysing a highly accurate conservative method for solving the quantum Zakharov system. The scheme is based on a linearly-implicit compact finite difference discretization and conserve the mass as well as energy in discrete level. Detailed numerical analysis is presented which shows the method is fourth-order accurate in space and second-order accurate in time. Several numerical examples are reported to confirm the conservation properties and high accuracy of the proposed scheme. Finally the compact scheme is applied to study the convergence rate of the quantum Zakharov system to its limiting model in the semi-classical limit.

The logarithmic nonlinearity has been used in many partial differential equations (PDEs) for modeling problems in various applications. Due to the singularity of the logarithmic function, it introduces tremendous difficulties in establishing mathematical theories, as well as in designing and analyzing numerical methods for PDEs with such nonlinearity. Here, we take the logarithmic Schrödinger equation (LogSE) as a prototype model. Instead of regularizing f(ρ)=lnρ in the LogSE directly and globally as being done in the literature, we propose a local energy regularization (LER) for the LogSE by first regularizing F(ρ)=ρlnρ−ρ locally near ρ=0+ with a polynomial approximation in the energy functional of the LogSE and then obtaining an energy regularized logarithmic Schrödinger equation (ERLogSE) via energy variation. Linear convergence is established between the solutions of ERLogSE and LogSE in terms of a small regularization parameter 0<ε≪1. Moreover, the conserved energy of the ERLogSE converges to that of LogSE quadratically, which significantly improves the linear convergence rate of the regularization method in the literature. Error estimates are also presented for solving the ERLogSE by using Lie–Trotter splitting integrators. Numerical results are reported to confirm our error estimates of the LER and of the time-splitting integrators for the ERLogSE. Finally, our results suggest that the LER performs better than regularizing the logarithmic nonlinearity in the LogSE directly.

We establish uniform error bounds of time-splitting Fourier pseudospectral (TSFP) methods for the nonlinear Klein–Gordon equation (NKGE) with weak power-type nonlinearity and O(1) initial data, while the nonlinearity strength is characterized by εp with a constant p ∈ N+ and a dimensionless parameter ε ∈ (0, 1], for the long-time dynamics up to the time at O(ε^{−β}) with 0 ≤ β ≤ p. In fact, when 0 < ε ≪ 1, the problem is equivalent to the long-time dynamics of NKGE with small initial data and O(1) nonlinearity strength, while the amplitude of the initial data (and the solution) is at O(ε). By reformulating the NKGE into a relativistic nonlinear Schr ̈odinger equation, we adapt the TSFP method to discretize it numerically. By using the method of mathematical induction to bound the numerical solution, we prove uniform error bounds at O(h^m + ε^{p−β} τ^2) of the TSFP method with h mesh size, τ time step and m ≥ 2 depending on the regularity of the solution. The error bounds are uniformly accurate for the long-time simulation up to the time at O(ε^{−β}) and uniformly valid for ε ∈ (0, 1]. Especially, the error bounds are uniformly at the second order rate for the large time step τ = O(ε^{−(p−β)/2}) in the parameter regime 0 ≤ β < p. Numerical results are reported to confirm our error bounds in the long-time regime. Finally, the TSFP method and its error bounds are extended to a highly oscillatory complex NKGE which propagates waves with wavelength at O(1) in space and O(ε^β) in time and wave velocity at O(ε^{−β}).

We consider the Schrödinger equation with a nondispersive logarithmic nonlinearity and a repulsive harmonic potential. For a suitable range of the coefficients, there exist two positive stationary solutions, each one generating a continuous family of solitary waves. These solutions are Gaussian, and turn out to be orbitally unstable. We also discuss the notion of ground state in this setting: for any natural definition, the set of ground states is empty.

We investigate variational principles for the approximation of quantum dynamics that apply for approximation manifolds that do not have complex linear tangent spaces. The first one, dating back to McLachlan (1964) minimizes the residuum of the time-dependent Schrödinger equation, while the second one, originating from the lecture notes of Kramer--Saraceno (1981), imposes the stationarity of an action functional. We characterize both principles in terms of metric and a symplectic orthogonality conditions, consider their conservation properties, and derive an elementary a-posteriori error estimate. As an application, we revisit the time-dependent Hartree approximation and frozen Gaussian wave packets.

We consider the nonlinear Schr{ö}dinger equation with a defocusing nonlinearity which is mass-(super)critical and energy-subcritical. We prove uniform in time error estimates for the Lie-Trotter time splitting discretization. This uniformity in time is obtained thanks to a vectorfield which provides time decay estimates for the exact and numerical solutions. This vectorfield is classical in scattering theory, and requires several technical modifications compared to previous error estimates for splitting methods.

We consider the Schrodinger equation with a logarithmic nonlinearity and a repulsive harmonic potential. Depending on the parameters of the equation, the solution may or may not be dispersive. When dispersion occurs, it does with an exponential rate in time. To control this, we change the unknown function through a generalized lens transform. This approach neutralizes the possible boundary effects, and could be used in the case of the Schrodinger equation without potential. We then employ standard splitting methods on the new equation via a nonuniform grid, after the logarithmic nonlinearity has been regularized. We also discuss the case of a power nonlinearity and give some results concerning the error estimates of the first-order Lie-Trotter splitting method for both cases of nonlinearities. Finally extensive numerical experiments are reported to investigate the dynamics of the equations.

We present a uniformly accurate finite difference method and establish rigorously its uniform error bounds for the Zakharov system (ZS) with a dimensionless parameter 0 < ε ≤ 1, which is inversely proportional to the speed of sound. In the subsonic limit regime, i.e., 0 < ε << 1, the solution propagates highly oscillatory waves and/or rapid outgoing initial layers due to the perturbation of the wave operator in the ZS and/or the incompatibility of the initial data, which is characterized by two parameters α ≥ 0 and β ≥ −1. Specifically, the solution propagates waves with wavelength of O(ε) and O(1) in time and space, respectively, and amplitude at O(ε^{min{2,α,1+β}}). This high oscillation of the solution in time brings significant difficulties in designing numerical methods and establishing their error bounds, especially in the subsonic limit regime. A uniformly accurate finite difference method is proposed by reformulating the ZS into an asymptotic consistent formulation and adopting an integral approximation of the oscillatory term. By applying the energy method and using the limiting equation via a nonlinear Schr ̈odinger equation with an oscillatory potential, we rigorously establish two independent error bounds at O(h^2 + τ^2/ε) and O(h^2 + τ^2 + τε^{α^∗} + ε^{1+α^∗}), respectively, with h the mesh size, τ the time step and α^∗ = min{1, α, 1 + β}. Thus we obtain error bounds at O(h^2 + τ^{4/3}) and O(h^2 + τ^{1+α^∗/(2+α^∗)}) for well-prepared (α^∗ = 1) and ill-prepared (0 ≤ α^∗ < 1) initial data, respectively, which are uniform in both space and time for 0 < ε ≤ 1 and optimal at the second order in space. Other techniques in the analysis include the cut-off technique for treating the nonlinearity and inverse estimates to bound the numerical solution. Numerical results are reported to demonstrate our error bounds.