We study the gap solitons and their stability properties in a Bose-Einstein condensation (BEC) under three-body interaction loaded in a Jacobian elliptic sine potential, which can be described by a cubic-quintic Gross-Pitaevskii equation (GPE) in the mean-field approximation. Firstly, the GPE is transformed into a stationary cubic-quintic nonlinear Schrödinger equation (NLSE) by the multi-scale method. A class of analytical solution of the NLSE is presented to describe the gap solitons. It is shown analytically that the amplitude of the gap soliton decreases as the two-body or three-body interaction strength increases. Secondly, many kinds of gap solitons, including the fundamental soliton and the sub-fundamental soliton, are obtained numerically by the Newton-Conjugate-Gradient (NCG) method. There are two families of fundamental solitons: one is the on-site soliton and the other is the off-site soliton. All of them are bifurcated from the Bloch band. Both in-phase and out-phase dipole solitons for off-site solitons do exist in such a nonlinear system. The numerical results also indicate that the amplitude of the gap soliton decreases as the nonlinear interaction strength increases, which accords well with the analytical prediction. Finally, long-time dynamical evolution for the GPE is performed by the time-splitting Fourier spectrum method to investigate the dynamical stability of gap solitons. It is shown that the on-site solitons are always dynamically stable, while the off-site solitons are always unstable. However, both stable and unstable in-phase or out-phase dipole solitons, which are not bifurcated from the Bloch band, indeed exist. For a type of out-phase soliton, there is a critical value <inline-formula><tex-math id="M1">\begin{document}$ q_c$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="1-20191278_M1.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="1-20191278_M1.png"/></alternatives></inline-formula> when the chemical potential <i>μ</i> is fixed. The solitons are linearly stable as <inline-formula><tex-math id="M2">\begin{document}$ q>q_c$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="1-20191278_M2.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="1-20191278_M2.png"/></alternatives></inline-formula>, while they are linearly unstable for <inline-formula><tex-math id="M3">\begin{document}$ q<q_c$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="1-20191278_M3.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="1-20191278_M3.png"/></alternatives></inline-formula>. Therefore, the modulus <i>q</i> plays an important role in the stability of gap solitons. One can change the dynamical behavior of gap solitons by adjusting the modulus of external potential in experiment. We also find that there exists a kind of gap soliton, in which the soliton is dynamically unstable if only the two-body interaction is considered, but it becomes stable when the three-body interaction is taken into account. This indicates that the three-body interaction has influence on the stability of gap solitons.
Read full abstract