Discrete solitons in nonlinear Schrödinger lattices with a power-law nonlinearity

We study the discrete nonlinear Schrödinger lattice model with the onsite nonlinearity of the general form, |u|^2σu. We systematically verify the conditions for the existence and stability of discrete solitons in the one-dimensional version of the model predicted by means of the variational approximation (VA), and demonstrate the following: monostability of fundamental solitons (FSs) in the case of the weak nonlinearity, 2σ+1<3.68; bistability, in a finite range of values of the soliton’s power, for 3.68<2σ+1<5; and the presence of a threshold (minimum norm of the FS), for 2σ+1≥5. We also perform systematic numerical simulations to study higher-order solitons in the same general model, i.e., bound states of the FSs. While all in-phase bound states are unstable, stability regions are identified for antisymmetric double solitons and their triple counterparts. These numerical findings are supplemented by an analytical treatment of the stability problem, which allows quantitively accurate predictions for the stability features of such multipulses. When these waveforms are found to be unstable, we show, by means of direct simulations, that they self-trap into a persistent lattice breather, or relax into a stable FS, or sometimes decay completely.     

An iterative starting method to control parasitism for the Leapfrog method

The Leapfrog method is a time-symmetric multistep method, widely used to solve the Euler equations and other Hamiltonian systems, due to its low cost and geometric properties. A drawback with Leapfrog is that it suffers from parasitism. This paper describes an iterative starting method, which may be used to reduce to machine precision the size of the parasitic components in the numerical solution at the start of the computation. The severity of parasitic growth is also a function of the differential equation, the main method and the time-step. When the tendency to parasitic growth is relatively mild, computational results indicate that using this iterative starting method may significantly increase the time-scale over which parasitic effects remain acceptably small. Using an iterative starting method, Leapfrog is applied to the cubic Schrödinger equation. The computational results show that the Hamiltonian and soliton behaviour are well-preserved over long time-scales.     

求解Schrdinger方程的高阶紧致差分格式

基于Richardson外推法提出了一种求解Schrdinger方程的高阶紧致差分方法.该方法首先利用二阶微商的四阶精度紧致差分逼近公式对原方程进行求解,然后利用Richardson外推技术外推一次,得到了Schrdinger方程具有O(r~4+h~4)精度的数值解.通过Fourier分析方法证明了该格式是无条件稳定的.数值实验验证了该方法的高阶精度及有效性.     

Eigenvalue estimates for Schrödinger operators on metric trees

We consider Schrödinger operators on radial metric trees and prove Lieb–Thirring and Cwikel–Lieb–Rozenblum inequalities for their negative eigenvalues. The validity of these inequalities depends on the volume growth of the tree. We show that the bounds are valid in the endpoint case and reflect the correct order in the weak or strong coupling limit.     

Short-range interaction in three-dimensional quantum mechanics

We show that it is possible to define shape-independent three-dimensional short-range quantum interactions in two-parameter form for non-spherical angular momentum channels through double rescaling of potential strength. Unlike the special case of l=0, where the zero-range limit of the system is renormalizable, the effective ranges diverge for l≠0 channels, and the system becomes trivial at zero-size limit. It is also shown that the two-parameter representation with finite interaction range is useful in analyzing phase shifts and describing resonances with accuracy in non-spherical scatterings.     

Fast Gaussian wavepacket transforms and Gaussian beams for the Schrödinger equation

This paper introduces a wavepacket-transform-based Gaussian beam method for solving the Schrödinger equation. We focus on addressing two computational issues of the Gaussian beam method: how to generate a Gaussian beam representation for general initial conditions and how to perform long time propagation for any finite period of time. To address the first question, we introduce fast Gaussian wavepacket transforms and develop on top of them an efficient initialization algorithm for general initial conditions. Based on this new initialization algorithm, we address the second question by reinitializing the beam representation when the beams become too wide. Numerical examples in one, two, and three dimensions demonstrate the efficiency and accuracy of the proposed algorithms. The methodology can be readily generalized to deal with other semi-classical quantum mechanical problems.     

The backward phase flow and FBI-transform-based Eulerian Gaussian beams for the Schrödinger equation

We propose the backward phase flow method to implement the Fourier–Bros–Iagolnitzer (FBI)-transform-based Eulerian Gaussian beam method for solving the Schrödinger equation in the semi-classical regime. The idea of Eulerian Gaussian beams has been first proposed in [12]. In this paper we aim at two crucial computational issues of the Eulerian Gaussian beam method: how to carry out long-time beam propagation and how to compute beam ingredients rapidly in phase space. By virtue of the FBI transform, we address the first issue by introducing the reinitialization strategy into the Eulerian Gaussian beam framework. Essentially we reinitialize beam propagation by applying the FBI transform to wavefields at intermediate time steps when the beams become too wide. To address the second issue, inspired by the original phase flow method, we propose the backward phase flow method which allows us to compute beam ingredients rapidly. Numerical examples demonstrate the efficiency and accuracy of the proposed algorithms.     

Entanglement conditions for four-mode states via the uncertainty relations in conjunction with partial transposition

From the Heisenberg uncertainty relation in conjunction with partial transposition, we derive a class of inequalities for detecting entanglements in four-mode states. The sufficient conditions for bipartite entangled states are presented. We also discuss the generalization of the entanglement conditions via the Schrödinger-Robertson indeterminacy relation, which are in general stronger than those based on the Heisenberg uncertainty relation.     

Propagation of light in(2+1)-dimensional nonlinear optical media with spatially inhomogeneous nonlinearities

The(2 1)-dimensional nonlinear Schr(o)dinger(NLS)equation with spatially inhomogeneous nonlinearities is investigated,which describes propagation of light in(2 1)-dimensional nonlinear optical media with inhomogeneous nonlinearities.New types of optical modes and nonlinear effects in optical media are presented numerically.The results reveal that the regular split of beam can be obtained in (2 1)-dimensional nonlinear optical media with inhomogeneous nonlinearities,by adjusting the guiding parameter.Furthermore,the stability of beam regular split is discussed numerically,and the results reveal that the beam regular split is stable to the finite initial perturbations.