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    . 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   .
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  -  is a powerful and systematic method to test the integrability of the physically important nonlinear evolution equations. By virtue of these techniques, Morris  has analyzed the multi-component nonlinear Schrödinger equations.
In 1980’s, based upon the nonlinear connection theory proposed by Lu et al.  Guo et al.    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   , study on the properties of the BEC systems has attracted peoples attention   . The BECs have also been observed in certain atomic gases such as the lithium, hydrogen, helium and potassium atomic gases   . 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   . In the following paper, we plan to investigate the following Gross-Pitaevskii equation in the Bose-Einstein condensate   by the covariant prolongation structure theory:
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   . 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   .
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:
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:
are the structure constants of the prolongation algebra . Using the induced connection , we can define the following covariant exterior derivative:
where and are the curvature coefficients on and the torsion coefficients in the fiber space , respectively, and given by
Requiring is an extended closed ideal, we may derive the following equation from Equation (4)
where and are the zero and one forms on the base manifold , respectively. Equation (6) can decomposed into the following equations:
Equation (7) is called the fundamental equation of the prolongation structure   . 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 :
We define the independent variables as . The Gross-Pitaevskii equation can then be expressed in the following set of two-forms given by  
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  
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  
Then, substituting the first equation of Equation (5) into Equation (11), we have the following over-determined difference equations
Solving the over-determined difference equations Equation (12), we obtain the following solutions
with as the hidden spectral parameter and the other components are zero.
Let us use the two dimensional linear representation of   given by,
Setting the transformation , which leads to
we therefore derive the following Lax-Pairs, given by
where represent the generators of the -symmetry   .
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
where ( , ) are the generators of Lie algebra  .
Besides the previous symmetries, we select the generators of the - symmetry  and we obtain the following Lax-representation
where are the generators of Lie algebra  .
Another Lax-representation can be derived in the form
where and are the generators of a -symmetry   and centreless Virasoro Lie algebra   .
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
Then setting  ,
and substituting it into Equation (21), then comparing the coefficient of , we have
from which we obtain
From the compatibility condition, the infinitely-many conservation laws for Equation (1) can be expressed as 
where the conserved density and the conserved flow are the following
Using the vanishing boundary condition, we can give the three constants of motions from the obtained conservation laws,
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    , 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
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
where denotes the Hirota’s derivative    .
To construct the soliton solutions of Equation (1), we expand and with respect to a formal expansion parameter as
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
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
The two-soliton solutions to Equation (1) is written as
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.
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.
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.