Multisymplectic implicit and explicit methods for Klein–Gordon–Schrödinger equations

We propose multisymplectic implicit and explicit Fourier pseudospectral methods for the Klein-Gordon-Schrödinger equations. We prove that the implicit method satisfies the charge conservation law exactly. Both methods provide accurate solutions in long-time computations and simulate the soliton collision well. Numerical results show the abilities of the two methods in preserving charge, energy, and momentum conservation laws.     

A conservative Fourier pseudospectral algorithm for a coupled nonlinear Schroedinger system

We derive a new method for a coupled nonlinear Schr/Sdinger system by using the square of first-order Fourier spectral differentiation matrix D1 instead of traditional second-order Fourier spectral differentiation matrix D2 to approximate the second derivative. We prove the proposed method preserves the charge and energy conservation laws exactly. In numerical tests, we display the accuracy of numerical solution and the role of the nonlinear coupling parameter in cases of soliton collisions. Numerical experiments also exhibit the excellent performance of the method in preserving the charge and energy conservation laws. These numerical results verify that the proposed method is both a charge-preserving and an energy-preserving algorithm.     

A conservative Fourier pseudospectral algorithm for a coupled nonlinear Schrdinger system

We derive a new method for a coupled nonlinear Schrdinger system by using the square of first-order Fourier spectral differentiation matrix D1 instead of traditional second-order Fourier spectral differentiation matrix D2 to approximate the second derivative.We prove the proposed method preserves the charge and energy conservation laws exactly.In numerical tests,we display the accuracy of numerical solution and the role of the nonlinear coupling parameter in cases of soliton collisions.Numerical experiments also exhibit the excellent performance of the method in preserving the charge and energy conservation laws.These numerical results verify that the proposed method is both a charge-preserving and an energy-preserving algorithm.     

A new finite difference scheme for a dissipative cubic nonlinear Schr6dinger equation

This paper considers the one-dimensional dissipative cubic nonlinear SchrSdinger equation with zero Dirichlet boundary conditions on a bounded domain. The equation is discretized in time by a linear implicit three-level central difference scheme, which has analogous discrete conservation laws of charge and energy. The convergence with two orders and the stability of the scheme are analysed using a priori estimates. Numerical tests show that the three-level scheme is more efficient.     

A conservative Fourier pseudospectral algorithm for the nonlinear Schrdinger equation

In this paper, we derive a new method for a nonlinear Schrodinger system by using the square of the first-order Fourier spectral differentiation matrix D1 instead of the traditional second-order Fourier spectral differentiation matrix D2 to approximate the second derivative. We prove that the proposed method preserves the charge and energy conservation laws exactly. A deduction argument is used to prove that the numerical solution is second-order convergent to the exact solutions in ｜｜·｜｜2 norm. Some numerical results are reported to illustrate the efficiency of the new scheme in preserving the charge and energy conservation laws.     

A new finite difference scheme for a dissipative cubic nonlinear Schrdinger equation

This paper considers the one-dimensional dissipative cubic nonlinear Schrdinger equation with zero Dirichlet boundary conditions on a bounded domain.The equation is discretized in time by a linear implicit three-level central difference scheme,which has analogous discrete conservation laws of charge and energy.The convergence with two orders and the stability of the scheme are analysed using a priori estimates.Numerical tests show that the three-level scheme is more efficient.     

Estimating topological charge of propagating vortex from single-shot non-imaged speckle

Encoding information using the topological charge of vortex beams has been proposed for optical communications. The conservation of the topological charge on propagation and the detection of the topological charge by a receiver are significant in these applications and have been well established in free-space. However, when vortex beams enter a diffuser,the wavefront is distorted, leading to a challenge in the conservation and detection of the topological charge. Here, we present a technique to measure the value of the topological charge of a vortex beam obscured in the randomly scattered light. The results of the numerical simulations and experiments are presented and are in good agreement. In particular, only a single-shot measurement is required to detect the topological charge of vortex beams, indicating that the method is applicable to a dynamic diffuser.     

On the energy conservation electrostatic particle-in-cell/Monte Carlo simulation: Benchmark and application to the radio frequency discharges

We benchmark and analyze the error of energy conservation （EC） scheme in particle-in-cell/Monte Carlo （PIC/MC） algorithms by simulating the radio frequency discharge. The plasma heating behaviors and electron distributing functions obtained by one-dimensional （1D） simulation are analyzed. Both explicit and implicit algorithms are checked. The results showed that the EC scheme can eliminated the self-heating with wide grid spacing in both cases with a small reduction of the accuracies. In typical parameters, the EC implicit scheme has higher precision than EC explicit scheme. Some ＂numerical cooling＂ behaviors are observed and analyzed. Some other errors are also analyzed. The analysis showed that the EC implicit scheme can be used to qualitative estimation of some discharge problems with much less computational resource cost without much loss of accuracies.     

Construction of Third-Order Diagonal Implicit Runge-Kutta Methods for Stiff Problems

We presents a new third-order diagonally implicit Runge~Kutta integration formula for stiff initial value problems, designed to be A-stable methods. The stability of the methods is analyzed and numerical results are shown to verify the conclusions.     

Explicit multi-symplectic method for the Zakharov—Kuznetsov equation

We propose an explicit multi-symplectic method to solve the two-dimensional Zakharov-Kuznetsov equation. The method combines the multi-symplectic Fourier pseudospectral method for spatial discretization and the Euler method for temporal discretization. It is verified that the proposed method has corresponding discrete multi-symplectic conservation laws. Numerical simulations indicate that the proposed scheme is characterized by excellent conservation.     

Simulation of shock-induced instability using an essentially conservative adaptive CE/SE method

An essentially conservative adaptive space time conservation element and solution element(CE/SE)method is proposed for the effective simulation of shock-induced instability with low computational cost.Its implementation is based on redefined conservation elements(CEs)and solution elements(SEs),optimized interpolations and a Courant number insensitive CE/SE scheme.This approach is used in two applications,the Woodward double Mach reflection and a twocomponent Richtmyer–Meshkov instability experiment.This scheme reveals the essential features of the investigated cases,captures small unstable structures,and yields a solution that is consistent with the results from experiments or other high order methods.     

The Squeezing Effect in a Mesoscopic RLC Circuit

Using the path integral method we derive quantum wave function and quantum fluctuation of charge and current in the mesoscopic RLC circuit.We find that the quantum fluctuation of charge decreases with time,oppositely,the quantum fluctuation of current increases with time monotonously.Therefore there is a squeezing effect in the circuit.If some more charge devices are used in the mesoscopic-damped circuit,the quantum noise can be reduced.We also find that uncertainty relation of charge and current periodically varies with the period π/2 in the under-damped case.     

Current-driven transformations of a skyrmion tube and a bobber in stepped nanostructures of chiral magnets

It has been over 30 years since the concept of optical vortices was first proposed by Coullet et al.in 1989,and the field of structured beams has grown extremely.In the last two decades,partially coherent vortex beams(PCVBs)have received increasing interest in the fields of optical manipulation,optical communication,optical imaging,etc.,and great progress has been made in the area of the coherence singularities,generation methods,topological charge measurements,and promising applications of PCVBs.In this review,we firstly outline the basic concepts of PCVBs.We explicate the relationship between the coherence vortices and optical vortices,and the evolution behavior of optical vortices to coherence vortices is summarized in detail.We discuss a special form of coherence singularity,ring dislocation,mathematically and physically.The ring dislocation in the correlation functions under low coherence is dependent on the mode indices,which provide a feasible approach to measure mode indices of PCVBs.Subsequently,we summarize the various methods for measuring the topological charge of PCVBs,highlight the measurement method based on the cross-correlation function,and a physical explanation on the relation between ring dislocation and topological charge is given.After that,we review the recent advances on experimental generation of several kinds of PCVBs.Lastly,we give an overview on the potential applications of PCVBs.     

ISSDE: A Monte Carlo implicit simulation code based on Stratonovich SDE approach of Coulomb collision

A Monte Carlo implicit simulation program, Implicit Stratonovich Stochastic Differential Equations(ISSDE), is developed for solving stochastic differential equations(SDEs) that describe plasmas with Coulomb collision. The basic idea of the program is the stochastic equivalence between the Fokker–Planck equation and the Stratonovich SDEs. The splitting method is used to increase the numerical stability of the algorithm for dynamics of charged particles with Coulomb collision. The cases of Lorentzian plasma, Maxwellian plasma and arbitrary distribution function of background plasma have been considered. The adoption of the implicit midpoint method guarantees exactly the energy conservation for the diffusion term and thus improves the numerical stability compared with conventional Runge–Kutta methods. ISSDE is built with C++ and has standard interfaces and extensible modules. The slowing down processes of electron beams in unmagnetized plasma and relaxation process in magnetized plasma are studied using the ISSDE, which shows its correctness and reliability.     

Lower Exciton Number Strong Light Matter Interaction in Plasmonic Tweezers

The plasmonic nanocavity is an excellent platform for the study of light matter interaction within a sub-diffraction volume under ambient conditions.We design a structure of plasmonic tweezers,which can trap molecular Jaggregates and also serve as a plasmonic cavity with which to investigate strong light matter interaction.The optical response of the cavity is calculated via finite-difference time-domain methods,and the optical force is evaluated based on the Maxwell stress tensor method.With the help of the coupled oscillator model and virtual exciton theory,we investigate the strong coupling progress at the lower level of excitons,finding that a Rabi splitting of 230 meV can be obtained in a single exciton system.We further analyze the relationship between optical force and model volume in the coupling system.The proposed method offers a way to locate molecular J-aggregates in plasmonic tweezers for investigating optical force performance and strong light matter interaction.     

Scattering properties of the ultracold 133Cs2 triplet state

The elastic scattering properties of ultracold 133Cs2 triplet state are investigated in detail.We construct a potential curve of the 133Cs2 triplet state,based on the latest ab initio molecular potential data and show how the scattering parameters are obtained by using three methods:the Numerov method,the semiclassical method and the variable phase method,where the scattering lengths of the 133Cs2 triplet state,i.e.301.79a0,300.67a0 and 310.81a0 are obtained respectively,with a0 being the Bohr radius.We also calculate the effective range and the number of bound states for the 133Cs2 triplet state.Our results are in agreement with the recent experimental data and the theoretical calculations.This confirms that the results of the scattering properties of the ultracold 133Cs2 triplet state,calculated by using these three methods,are reliable.     

Comparison of the PIC model and the Lie algebraic method in the simulation of intense continuous beam transport

Both the PIC (Particle-In-Cell) model and the Lie algebraic method can be used to simulate the transport of intense continuous beams. The PIC model is to calculate the space charge field, which is blended into the external field, and then simulate the trajectories of particles in the total field; the Lie algebraic method is to simulate the intense continuous beam transport with transport matrixes. Two simulation codes based on the two methods are developed respectively, and the simulated results of transport in a set of electrostatic lenses are compared. It is found that the results from the two codes are in agreement with each other, and both approaches have their own merits.     

The Squeezing Effect in a Mesoscopic RLC Circuit

Using the path integral method we derive quantum wave function and quantum fluctuations of charge andcurrent in the mesoscopic RLC circuit. We find that the quantum fluctuation of charge decreases with time, oppositely,the quantum fluctuation of current increases with time monotonously. Therefore there is a squeezing effect in the circuit.If some more charge devices are used in the mesoscopic-damped circuit, the quantum noise can be reduced. We also findthat uncertainty relation of charge and current periodically varies with the period π/2 in the under-damped case.     

Forward pion-nucleon charge exchange reaction and Regge constraints

We present our recent study of pion-nucleon charge exchange amplitudes above 2 GeV. We analyze the forward pion-nucleon charge exchange reaction data in a Regge model and compare the resulting amplitudes with those from the Karlsruhe-Helsinki and George-Washington-University partial-wave analyses. We explore possible high-energy constraints for theoretical baryon resonance analyses in the energy region above 2 GeV. Our results show that for the pion-nucleon charge exchange reaction, the appropriate energy region for matching meson-nucleon dynamics to diffractive scattering should be around 3 GeV for the helicity flip amplitude.