Covariant Prolongation Structure, Conservation Laws and Soliton Solutions of the Gross-Pitaevskii Equation in the Bose-Einstein Condensate ()
1. Introduction
Nonlinear evolution equations (NLEEs) have been studied in diverse areas in physics and applied mathematics such as plasma physics, nonlinear optical fibers, condensed matter etc [1] [2] [3] . The inhomogeneous nonlinear Schrödinger equation among those (NLEEs) 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 [4] [5] .
As it is known, the investigation of integrability of certain kinds of (NLEEs) by many researchers has generated a great deals of attention over the past years and now many methods to analyze the complete integrability of nonlinear evolution equations are developed. Among them, Wahlquist and Estabrook’s prolongation technique [6] - [23] is a powerful and systematic method to test the integrability of the physically important nonlinear evolution equations. By virtue of these techniques, Morris [23] has analyzed the multi-component nonlinear Schrödinger equations.
In 1980’s, based upon the nonlinear connection theory proposed by Lu et al. [24] Guo et al. [25] [26] [27] have proposed a covariant geometry theory for the prolongation structure of the nonlinear evolution equation.
Since the experimental realization of the Bose-Einstein condensate (BEC) for rubidium and sodium [28] [29] , study on the properties of the BEC systems has attracted peoples attention [30] [31] . The BECs have also been observed in certain atomic gases such as the lithium, hydrogen, helium and potassium atomic gases [32] [33] . For a cigar-shaped BEC at a relatively low density, when the energy of the two-body interaction is much less than the kinetic energy in the transverse direction, the BEC system can be treated as quasi-one-dimensional [34] [35] . In the following paper, we plan to investigate the following Gross-Pitaevskii equation in the Bose-Einstein condensate [36] [37] by the covariant prolongation structure theory:
(1)
where
represents the mean-field wave function of the Bose-Einstein condensate;
and
are the normalized distance and retarded time, res- pectively;
and
are all the real numbers;
is the nonlinearity para- meter;
is the gain
or loss
term;
represents the linear external potential, while
accounts for the harmonic external potential.
In soliton theory and other fields of science and engineering, the language of technical computing played a very important role in analytically dealing with large amounts of complicated and tedious algebraic calculations [38] [39] . In this paper, we will employ symbolic computation to study the integrability aspects and relevant soliton structures of Gross-Pitaevskii equation in the Bose-Einstein condensate [36] [37] .
However, to our knowledge, for Equation (1), Lax-Pair, Conservation laws, multiple soliton solutions via Hirota’s method and symbolic computation have not been discussed yet. Motivated by the above, a Lax-Pair based on the generators of some hidden structural symmetries governed the dynamics of the (GP) equation will be got in section 3. In section 4, an infinite sequence of conservation laws of Equation (1) are obtained. In addition, in section 5 we present the exact one and two soliton solutions of the Gross-Pitaevskii Equation in the Bose-Einstein condensate. Finally, the conclusion will be addressed in section 6.
2. Covariant Theory for Prolongation Structure of Nonlinear Evolution Equations
For a given
-dimensional nonlinear evolution equation with two independent variables
and
, we can introduce a set of new variables
and define a set of 2-forms
such that it constitutes a differential closed ideal, which lead to the
-dimensional nonlinear evolution equation if the ideal is restricted on the solution space
.
Now we take
as the base space,
named prolongation variables as the fiber space and
as the structure group generated by the prolongation algebra
. We can consider a principle bundle
and the bundle
associated with
. Define the local cross-section on
,
, and its covariant derivatives:
(2)
where
is the dimension of the representation space of the prolongation algebra,
are the coefficients of the connection on the principal bundle
and
are the coefficients of the generators of the prolongation algebra
.
Then we introduce the following connection 1-forms:
(3)
are the structure constants of the prolongation algebra
. Using the induced connection
, we can define the following covariant exterior derivative:
(4)
where
and
are the curvature coefficients on
and the torsion coefficients in the fiber space
, respectively, and given by
(5)
Requiring
is an extended closed ideal, we may derive the following equation from Equation (4)
(6)
where
and
are the zero and one forms on the base manifold
, respectively. Equation (6) can decomposed into the following equations:
(7)
Equation (7) is called the fundamental equation of the prolongation structure [26] [27] . In general, we may completely determine the prolongation structure of a given nonlinear system when the solution of the one fundamental equation can be found.
3. Covariant Prolongation Structure of Gross-Pitaevskii Equation in the Bose-Einstein Condensate
In order to express Equation (1) in differential forms, we add the conjugate equation of Equation (1) to Equation (1) and obtain the following system :
(8)
We define the independent variables as
. The Gross-Pitaevskii equation can then be expressed in the following set of two-forms given by [26] [27]
(9)
where the letter
denotes the exterior derivative and the symbol
represents the exterior product. In order to ensure complete equivalence between the forms (9) and the Gross-Pitaevskii Equation (8), the ideal
must be closed, i.e.,
. In this closed ideal any local surface element which annuls the
also annuls their exterior derivatives
. In order to establish the prolongation structure, we extend the above ideal by adding to it a connection 1-forms, defined by [26] [27]
(10)
where
, and
are the prolongation variable. For some suitably chosen prolongation variables and imposing the closed condition of the extended ideal
under covariant exterior derivative, it leads to the covariant fundamental equations.
Substituting the above two forms
into the fundamental equation Equation (7), we have [26] [27]
(11)
Then, substituting the first equation of Equation (5) into Equation (11), we have the following over-determined difference equations
(12)
Solving the over-determined difference equations Equation (12), we obtain the following solutions
(13)
with
as the hidden spectral parameter and the other components are zero.
Let us use the two dimensional linear representation of
[22] [23] given by,
(14)
Setting the transformation
, which leads to
(15)
we therefore derive the following Lax-Pairs, given by
(16)
where
represent the generators of the
-symmetry [22] [23] .
On the other hand, by selecting the matrix representation of a generators of a
symmetry, the Lax-representation associated to such an algebra is then given by
(17)
where (
,
) are the generators of
Lie algebra [8] .
Besides the previous symmetries, we select the generators of the
- symmetry [23] and we obtain the following Lax-representation
(18)
where
are the generators of
Lie algebra [23] .
Another Lax-representation can be derived in the form
(19)
and
(20)
where
and
are the generators of a
-symmetry [22] [23] and centreless Virasoro Lie algebra [9] [10] .
From the previous discussion, it appears that the dynamics of the the Gross- Pitaevskii (GP) equation modeled by Equation (1), are basically governed by internal structural symmetries, including the Virasoro algebra,
,
,
and
. Such symmetries have some physical implications. For example, the
-symmetries show that the system (1) possesses some conserved quantities that are rotationally and hyperbolically invariant, respectively. Thus, we have shown that Eqution (1) is Lax integrable by giving its corresponding Lax-Pair Equations (16)-(20).
4. Conservation Laws of the Gross-Pitaevskii Equation in the Bose-Einstein Condensate
In the following, we will prove the existence of infinitely-many conservation laws, which further verifies the integrability of Equation (1).
By means of the one dimensional linear representation of
, we derive the Riccati equations
(21)
Then setting [40] ,
(22)
and substituting it into Equation (21), then comparing the coefficient of
, we have
(23)
from which we obtain
(24)
From the compatibility condition, the infinitely-many conservation laws for Equation (1) can be expressed as [41]
(25)
where the conserved density
and the conserved flow
are the following
(26)
Using the vanishing boundary condition, we can give the three constants of motions from the obtained conservation laws,
(27)
5. Exact multisoliton solutions of the Gross-Pitaevskii Equation in the Bose-Einstein Condensate
In order to derive the analytical soliton solutions to Equation (1), we will employ the Hirota bilinear method [42] [43] [44] , which is an efficient and direct approach to construct soliton solutions to nonlinear evolution equations via the bilinear forms from the dependent variables transformation.
To get the bilinear forms for Equation (1) we introduce the dependent variable transformation
(28)
where
is the complex differentiable function, and
is a real one. Substituting relation (28) into Equation (1), the bilinear equations of Equation (1) turns out to be in the following forms
(29)
where
denotes the Hirota’s derivative [42] [43] [44] .
To construct the soliton solutions of Equation (1), we expand
and
with respect to a formal expansion parameter
as
(30)
where
is the complex differentiable function, and
is a real one
To derive the one-soliton solutions to Equation (1), we truncate expressions Equation (30) as
and
, setting
and substituting then into Bilinear forms Equation (29). We obtain the one-soliton solutions to Equation (1) as
(31)
where
(32)
Figure 1. The one-soliton solution via expression (31) with parameters
;
;
;
;
.
it is depicted in Figure 1.
Similarly, in order to derive the two-soliton solutions, we can choose
(33)
where
(34)
The two-soliton solutions to Equation (1) is written as
(35)
Figure 2. The two-soliton solutions via expression (35) with parameters (a)
;
;
;
;
;
;
, (b)
;
;
;
;
;
;
.
and the corresponding depiction is shown in Figure 2.
If one and two-soliton solutions are calculated, then it is possible to generate the multi-soliton solution in the same way.
6. Conclusion
Throughout the present paper, we investigated the prolongation structure of the Gross-Pitaevskii 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 from the viewpoint of covariant prolongation structure. As a result, we have unearthed some hidden structural symmetries governing the dynamics of the Gross-Pitaevskii equation such as
,
, Virasoro algebra,
and
. Such symmetries have some physical implications. For example, the
-symmetries show that the system (1) possesses some conserved quantities that are rotationally and hyperbolically invariant, respectively. Thus, we have shown that Equation (1) is Lax integrable by giving its corresponding Lax-Pair Equations (16)-(20). In addition, infinite number of conservation Laws, one and two soliton solutions using Hirota bilinear method have been constructed. The prolongation structure analysis performed in the present study to the system (1) has revealed an infinite number of conserved quantities which stand as strong proof of integrability of this equation.
Acknowledgements
The authors would like to express their sincere thanks to the editors, anonymous reviewers and all the members of our discussion group for their critical comments and appropriate suggestions which have made this paper more precise and readable.