With the help of similarity transformation, we obtain analytical spatiotemporal self-similar solutions of the nonautonomous (3+1)-dimensional cubic-quintic Gross-Pitaevskii equation with time-dependent diffraction,...With the help of similarity transformation, we obtain analytical spatiotemporal self-similar solutions of the nonautonomous (3+1)-dimensional cubic-quintic Gross-Pitaevskii equation with time-dependent diffraction, nonlinearity, harmonic potential and gain or loss when two constraints are satisfied. These constraints between the system parameters hint that self-similar solutions form and transmit stably without the distortion of shape based on the exact balance between the diffraction, nonlinearity and the gain/loss. Based on these analytical results, we investigate the dynamic behaviours in a periodic distributed amplification system.展开更多
By the generalized sub-equation expansion method and symbolic computation, this paper investigates the (3+1)dimensional Gross-Pitaevskii equation with time-and space-dependent potential, time-dependent nonlinearity...By the generalized sub-equation expansion method and symbolic computation, this paper investigates the (3+1)dimensional Gross-Pitaevskii equation with time-and space-dependent potential, time-dependent nonlinearity, and gain or loss. As a result, rich exact analytical solutions are obtained, which include bright and dark solitons, Jacobi elliptic function solutions and Weierstrass elliptic function solutions. With computer simulation, the main evolution features of some of these solutions are shown by some figures. Nonlinear dynamics of a soliton pulse is also investigated under the different regimes of soliton management.展开更多
The Homotopy analysis method (HAM) is adopted to find the approximate analytical solutions of the Gross- Pitaevskii equation, a nonlinear Schrodinger equation is used in simulation of Bose-Einstein condensates trapp...The Homotopy analysis method (HAM) is adopted to find the approximate analytical solutions of the Gross- Pitaevskii equation, a nonlinear Schrodinger equation is used in simulation of Bose-Einstein condensates trapped in a harmonic potential. Comparisons between the analytical solutions and the numerical solutions have been made. The results indicate that they fit very well with each other when the atomic interaction is weak.展开更多
In this paper, applying the dependent and independent variables transformations as well as the Jacobi elliptic function expansion method, the envelope periodic solutions to one-dimensional Gross-Pitaevskii equation in...In this paper, applying the dependent and independent variables transformations as well as the Jacobi elliptic function expansion method, the envelope periodic solutions to one-dimensional Gross-Pitaevskii equation in Bose-Einstein condensates are obtained.展开更多
The Gross-Pitaevskii equation (GPE), that describes the wave function of a number of coherent Bose particles contained in a trap, contains the cube of the normalized wave function, times a factor proportional to the n...The Gross-Pitaevskii equation (GPE), that describes the wave function of a number of coherent Bose particles contained in a trap, contains the cube of the normalized wave function, times a factor proportional to the number of coherent atoms. The square of the wave function, times the above mentioned factor, is defined as the Hartree potential. A method implemented here for the numerical solution of the GPE consists in obtaining the Hartree potential iteratively, starting with the Thomas Fermi approximation to this potential. The energy eigenvalues and the corresponding wave functions for each successive potential are obtained by a spectral method described previously. After approximately 35 iterations a stability of eight significant figures for the energy eigenvalues is obtained. This method has the advantage of being physically intuitive, and could be extended to the calculation of a shell-model potential in nuclear physics, once the Pauli exclusion principle is allowed for.展开更多
Modulational instability conditions for the generation of localized structures in the context of matter waves in Bose-Einstein condensates are investigated analytically and numerically. The model is based on a modifie...Modulational instability conditions for the generation of localized structures in the context of matter waves in Bose-Einstein condensates are investigated analytically and numerically. The model is based on a modified Gross-Pitaevskii equation, which account for the energy dependence of the two-body scattering amplitude. It is shown that the modified term due to the quantum fluctuations modify significantly the modulational instability gain. Direct numerical simulations of the full modified Gross-Pitaevskii equation are performed, and it is found that the modulated plane wave evolves into a train of pulses, which is destroyed at longer times due to the effects of quantum fluctuations.展开更多
In this paper, we investigate the Gross-Pitaevskii (GP) equation which describes the propagation of an electron plasma wave packet with a large wavelength and small amplitude in a medium with a parabolic density and c...In this paper, we investigate the Gross-Pitaevskii (GP) equation which describes the propagation of an electron plasma wave packet with a large wavelength and small amplitude in a medium with a parabolic density and constant interactional damping by the Covariant Prolongation Structure Theory. As a result, we obtain general forms of Lax-Pair representations. In addition, some hidden structural symmetries that govern the dynamics of the GP equation such as SL(2,R), SL(2,C), Virasoro algebra, SU(1,1) and SU(2) are unearthed. Using the Riccati form of the linear eigenvalue problem, infinite number of conservation laws of the GP equation is explicitly constructed and the exact analytical soliton solutions are obtained by employing the simple and straightforward Hirota’s bilinear method.展开更多
We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation (BEC). This scheme combined with the sixth-or...We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation (BEC). This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.展开更多
The stability analysis of the Rosenbrock method for the numerical solutions of system of delay differential equations was studied. The stability behavior of Rosenbrock method was analyzed for the solutions of linear t...The stability analysis of the Rosenbrock method for the numerical solutions of system of delay differential equations was studied. The stability behavior of Rosenbrock method was analyzed for the solutions of linear test equation. The result that the Rosenbrock method is GP-stable if and only if it is A-stable is obtained.展开更多
A(3+1)-dimensional Gross-Pitaevskii(GP)equation with time variable coefficients is considered,andis transformed into a standard nonlinear Schr(o|¨)dinger(NLS)equation.Exact solutions of the(3+1)D GP equationare c...A(3+1)-dimensional Gross-Pitaevskii(GP)equation with time variable coefficients is considered,andis transformed into a standard nonlinear Schr(o|¨)dinger(NLS)equation.Exact solutions of the(3+1)D GP equationare constructed via those of the NLS equation.By applying specific time-modulated nonlinearities,dispersions,andpotentials,the dynamics of the solutions can be controlled.Solitary and periodic wave solutions with snaking andbreathing behavior are reported.展开更多
The orthometric height (OH) system plays a key role in geodesy, and it has broad applications in various fields and activities. Based on general relativity theory (GRT), on an arbitrary equi-geo- potential surface, th...The orthometric height (OH) system plays a key role in geodesy, and it has broad applications in various fields and activities. Based on general relativity theory (GRT), on an arbitrary equi-geo- potential surface, there does not exist the gravity frequency shift of an electromagnetic wave signal. However, between arbitrary two different equi-geopotential surfaces, there exists the gra- vity frequency shift of the signal. The relationship between the geopotential difference and the gravity frequency shift between arbitrary two points P and Q is referred to as the gravity frequency shift equation. Based on this equation, one can determine the geopotential difference as well as the OH difference between two separated points P and Q either by using electromagnetic wave signals propagated between P and Q, or by using the Global Positioning System (GPS) satellite signals received simultaneously by receivers at P and Q. Suppose an emitter at P emits a signal with frequency f towards a receiver at Q, and the received frequency of the signal at Q is , or suppose an emitter on board a flying GPS satellite emits signals with frequency f towards two receivers at P and Q on ground, and the received frequencies of the signals at P and Q are and , respectively, then, the geopoten-tial dif- ference between these two points can be determined based on the geopotential frequen- cy shift equation, using either the gravity frequency shift ? f or ? , and the corresponding OH difference is further determined based on the Bruns’ formula. Besides, using this approach a unified world height datum system might be realized, because P and Q could be chosen quite arbitrarily, e.g., they are located on two separated continents or islands.展开更多
A local refinement hybrid scheme(LRCSPH-FDM)is proposed to solve the two-dimensional(2D)time fractional nonlinear Schrodinger equation(TF-NLSE)in regularly or irregularly shaped domains,and extends the scheme to predi...A local refinement hybrid scheme(LRCSPH-FDM)is proposed to solve the two-dimensional(2D)time fractional nonlinear Schrodinger equation(TF-NLSE)in regularly or irregularly shaped domains,and extends the scheme to predict the quantum mechanical properties governed by the time fractional Gross-Pitaevskii equation(TF-GPE)with the rotating Bose-Einstein condensate.It is the first application of the purely meshless method to the TF-NLSE to the author’s knowledge.The proposed LRCSPH-FDM(which is based on a local refinement corrected SPH method combined with FDM)is derived by using the finite difference scheme(FDM)to discretize the Caputo TF term,followed by using a corrected smoothed particle hydrodynamics(CSPH)scheme continuously without using the kernel derivative to approximate the spatial derivatives.Meanwhile,the local refinement technique is adopted to reduce the numerical error.In numerical simulations,the complex irregular geometry is considered to show the flexibility of the purely meshless particle method and its advantages over the grid-based method.The numerical convergence rate and merits of the proposed LRCSPH-FDM are illustrated by solving several 1D/2D(where 1D stands for one-dimensional)analytical TF-NLSEs in a rectangular region(with regular or irregular particle distribution)or in a region with irregular geometry.The proposed method is then used to predict the complex nonlinear dynamic characters of 2D TF-NLSE/TF-GPE in a complex irregular domain,and the results from the posed method are compared with those from the FDM.All the numerical results show that the present method has a good accuracy and flexible application capacity for the TF-NLSE/GPE in regions of a complex shape.展开更多
We develop an exponential spline interpolation method to solve the nonlinear Schrödinger equation. The truncation error and stability analysis of the method are investigated and the method is shown to be uncon...We develop an exponential spline interpolation method to solve the nonlinear Schrödinger equation. The truncation error and stability analysis of the method are investigated and the method is shown to be unconditionally stable. The conservation quantities are computed to determine the conservation properties of the problem. We will describe the method and present numerical tests by two problems. The numerical simulations results demonstrate the well performance of the proposed method.展开更多
基金Project supported by the National Natural Science Foundations of China (Grant No. 11005092)the Program for Innovative Research Team of Young Teachers (Grant No. 2009RC01)the Scientific Research and Developed Fund of Zhejiang Agricultural and Forestry University,China (Grant No. 2009FK42)
文摘With the help of similarity transformation, we obtain analytical spatiotemporal self-similar solutions of the nonautonomous (3+1)-dimensional cubic-quintic Gross-Pitaevskii equation with time-dependent diffraction, nonlinearity, harmonic potential and gain or loss when two constraints are satisfied. These constraints between the system parameters hint that self-similar solutions form and transmit stably without the distortion of shape based on the exact balance between the diffraction, nonlinearity and the gain/loss. Based on these analytical results, we investigate the dynamic behaviours in a periodic distributed amplification system.
基金Project supported by Zhejiang Provincial Natural Science Foundations of China (Grant No. Y6090592)National Natural Science Foundation of China (Grant Nos. 11041003 and 10735030)+1 种基金Ningbo Natural Science Foundation (Grant Nos. 2010A610095,2010A610103,and 2009B21003)K.C. Wong Magna Fund in Ningbo University of China
文摘By the generalized sub-equation expansion method and symbolic computation, this paper investigates the (3+1)dimensional Gross-Pitaevskii equation with time-and space-dependent potential, time-dependent nonlinearity, and gain or loss. As a result, rich exact analytical solutions are obtained, which include bright and dark solitons, Jacobi elliptic function solutions and Weierstrass elliptic function solutions. With computer simulation, the main evolution features of some of these solutions are shown by some figures. Nonlinear dynamics of a soliton pulse is also investigated under the different regimes of soliton management.
基金Project supported by the National Natural Science Foundation of China(Grant No.11047010)the Key Project Foundation of the Education Ministry of China(Grant No.209128)
文摘The Homotopy analysis method (HAM) is adopted to find the approximate analytical solutions of the Gross- Pitaevskii equation, a nonlinear Schrodinger equation is used in simulation of Bose-Einstein condensates trapped in a harmonic potential. Comparisons between the analytical solutions and the numerical solutions have been made. The results indicate that they fit very well with each other when the atomic interaction is weak.
基金Supported by National Natural Science Foundation of China under Grant No. 90511009
文摘In this paper, applying the dependent and independent variables transformations as well as the Jacobi elliptic function expansion method, the envelope periodic solutions to one-dimensional Gross-Pitaevskii equation in Bose-Einstein condensates are obtained.
文摘The Gross-Pitaevskii equation (GPE), that describes the wave function of a number of coherent Bose particles contained in a trap, contains the cube of the normalized wave function, times a factor proportional to the number of coherent atoms. The square of the wave function, times the above mentioned factor, is defined as the Hartree potential. A method implemented here for the numerical solution of the GPE consists in obtaining the Hartree potential iteratively, starting with the Thomas Fermi approximation to this potential. The energy eigenvalues and the corresponding wave functions for each successive potential are obtained by a spectral method described previously. After approximately 35 iterations a stability of eight significant figures for the energy eigenvalues is obtained. This method has the advantage of being physically intuitive, and could be extended to the calculation of a shell-model potential in nuclear physics, once the Pauli exclusion principle is allowed for.
文摘Modulational instability conditions for the generation of localized structures in the context of matter waves in Bose-Einstein condensates are investigated analytically and numerically. The model is based on a modified Gross-Pitaevskii equation, which account for the energy dependence of the two-body scattering amplitude. It is shown that the modified term due to the quantum fluctuations modify significantly the modulational instability gain. Direct numerical simulations of the full modified Gross-Pitaevskii equation are performed, and it is found that the modulated plane wave evolves into a train of pulses, which is destroyed at longer times due to the effects of quantum fluctuations.
文摘In this paper, we investigate the Gross-Pitaevskii (GP) equation which describes the propagation of an electron plasma wave packet with a large wavelength and small amplitude in a medium with a parabolic density and constant interactional damping by the Covariant Prolongation Structure Theory. As a result, we obtain general forms of Lax-Pair representations. In addition, some hidden structural symmetries that govern the dynamics of the GP equation such as SL(2,R), SL(2,C), Virasoro algebra, SU(1,1) and SU(2) are unearthed. Using the Riccati form of the linear eigenvalue problem, infinite number of conservation laws of the GP equation is explicitly constructed and the exact analytical soliton solutions are obtained by employing the simple and straightforward Hirota’s bilinear method.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11571366 and 11501570the Open Foundation of State Key Laboratory of High Performance Computing of China+1 种基金the Research Fund of National University of Defense Technology under Grant No JC15-02-02the Fund from HPCL
文摘We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation (BEC). This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.
文摘The stability analysis of the Rosenbrock method for the numerical solutions of system of delay differential equations was studied. The stability behavior of Rosenbrock method was analyzed for the solutions of linear test equation. The result that the Rosenbrock method is GP-stable if and only if it is A-stable is obtained.
基金Supported by the National Natural Science Foundations of China under Grant Nos.10735030,10475055,10675065 and 90503006National Basic Research Program of China (973 Program 2007CB814800) and PCSIRT (IRT0734)+1 种基金the Research Fund of Postdoctoral of China under Grant No.20070410727Specialized Research Fund for the Doctoral Program of Higher Education under Grant No.20070248120
文摘A(3+1)-dimensional Gross-Pitaevskii(GP)equation with time variable coefficients is considered,andis transformed into a standard nonlinear Schr(o|¨)dinger(NLS)equation.Exact solutions of the(3+1)D GP equationare constructed via those of the NLS equation.By applying specific time-modulated nonlinearities,dispersions,andpotentials,the dynamics of the solutions can be controlled.Solitary and periodic wave solutions with snaking andbreathing behavior are reported.
文摘The orthometric height (OH) system plays a key role in geodesy, and it has broad applications in various fields and activities. Based on general relativity theory (GRT), on an arbitrary equi-geo- potential surface, there does not exist the gravity frequency shift of an electromagnetic wave signal. However, between arbitrary two different equi-geopotential surfaces, there exists the gra- vity frequency shift of the signal. The relationship between the geopotential difference and the gravity frequency shift between arbitrary two points P and Q is referred to as the gravity frequency shift equation. Based on this equation, one can determine the geopotential difference as well as the OH difference between two separated points P and Q either by using electromagnetic wave signals propagated between P and Q, or by using the Global Positioning System (GPS) satellite signals received simultaneously by receivers at P and Q. Suppose an emitter at P emits a signal with frequency f towards a receiver at Q, and the received frequency of the signal at Q is , or suppose an emitter on board a flying GPS satellite emits signals with frequency f towards two receivers at P and Q on ground, and the received frequencies of the signals at P and Q are and , respectively, then, the geopoten-tial dif- ference between these two points can be determined based on the geopotential frequen- cy shift equation, using either the gravity frequency shift ? f or ? , and the corresponding OH difference is further determined based on the Bruns’ formula. Besides, using this approach a unified world height datum system might be realized, because P and Q could be chosen quite arbitrarily, e.g., they are located on two separated continents or islands.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11501495,51779215,and 11672259)the Postdoctoral Science Foundation of China(Grant Nos.2015M581869 and 2015T80589)the Jiangsu Government Scholarship for Overseas Studies,China(Grant No.JS-2017-227)。
文摘A local refinement hybrid scheme(LRCSPH-FDM)is proposed to solve the two-dimensional(2D)time fractional nonlinear Schrodinger equation(TF-NLSE)in regularly or irregularly shaped domains,and extends the scheme to predict the quantum mechanical properties governed by the time fractional Gross-Pitaevskii equation(TF-GPE)with the rotating Bose-Einstein condensate.It is the first application of the purely meshless method to the TF-NLSE to the author’s knowledge.The proposed LRCSPH-FDM(which is based on a local refinement corrected SPH method combined with FDM)is derived by using the finite difference scheme(FDM)to discretize the Caputo TF term,followed by using a corrected smoothed particle hydrodynamics(CSPH)scheme continuously without using the kernel derivative to approximate the spatial derivatives.Meanwhile,the local refinement technique is adopted to reduce the numerical error.In numerical simulations,the complex irregular geometry is considered to show the flexibility of the purely meshless particle method and its advantages over the grid-based method.The numerical convergence rate and merits of the proposed LRCSPH-FDM are illustrated by solving several 1D/2D(where 1D stands for one-dimensional)analytical TF-NLSEs in a rectangular region(with regular or irregular particle distribution)or in a region with irregular geometry.The proposed method is then used to predict the complex nonlinear dynamic characters of 2D TF-NLSE/TF-GPE in a complex irregular domain,and the results from the posed method are compared with those from the FDM.All the numerical results show that the present method has a good accuracy and flexible application capacity for the TF-NLSE/GPE in regions of a complex shape.
文摘We develop an exponential spline interpolation method to solve the nonlinear Schrödinger equation. The truncation error and stability analysis of the method are investigated and the method is shown to be unconditionally stable. The conservation quantities are computed to determine the conservation properties of the problem. We will describe the method and present numerical tests by two problems. The numerical simulations results demonstrate the well performance of the proposed method.