One of the most powerful nonlinear equations that has been discussed and used in different field of physics is the Non-linear Schrödinger equation (NLSE) in (2 + 1) dimension. The structure and integrability of (2 + 1) dimensional equations have received considerable attention in the last few years   . Non-linear Schrödinger type equations are a particular case of interest. These equations were discovered by Calogero  and then discussed by Zakharov  . The geometrical properties of these equations have been studied  .
Historically, more and more problems involved nonlinearity; much attention had been paid to the integrability of nonlinear models. Close relationships were observed between Painlevé test and the integrability of non-linear Partial Differential Equations (PDEs). Therefore it carries significant importance to investigate whether the non-linear PDEs have the Painlevé characteristics. In order to verify this proposition various methods have been historically used by researchers like Ablowitz-Ramani-Segur (ARS) method  , the Weiss-Tabor-Carnevale (WTC) method  , the Kruskal simplification method  , Conte invariant method  and Pickering approach method  .
Various methods have verified Painlevé property but the Kruskal’s simplification and the Weiss-Tabor-Carnevale (WTC) methods are frequently employed for this verification. WTC method is significantly helpful in Painlevé test where non-linear PDEs are involved as tests of symmetry, Hirota’s method of bilinear forms and special solutions. However, this is less significant when applied to the coupled equation system or equation where high resonance is observable. This study proposes another prospective approach i.e. to apply WTC method with Kruskal’s simplification for verification of Painlevé test.
By observing chaotic behavior of certain systems and their non-integrability, significant advancement had been observed with respect to research on the integrability of non-linear PDEs thus leading to analysis of manifold classes of equations employing Painlevé test as a base    . Painlevé test had being used in various situations, although was neither observed as a necessity nor a condition for integrability of PDEs. The current study proposes a Painleve test analysis using (2 + 1) dimensions, for testing the integrability of Non-linear Schrödinger equation (NLSE). It was observed that such approach could satisfy requirements of Cauchy-Kovalevskaya theorem while providing satisfactory results.
The (2 + 1) dimensional NLSE under taken in this study is:
In Equation (1), S is called a forcing term which shows the real quantity and or prescribes the elliptic or hyperbolic nature respectively.
2. Painlevé Analysis
Painlevé test is frequently used to determine integrability of non-linear equations. A PDE which lacks any movable dimension (Algebraic or Logarithmic) is called Painlevé-type (P-type)  . An ordinary differential equation (ODE) carries movable singularities without having movable branch points which place the same ODE into P-type equations. P-test to analyze the singularities structure for
ODE was explained by Ablowitz et al.  , further extended by Weiss-Tabor- Carnevale (WTC) to PDE  the method provides a great basis for investigation of integrability of various non-linear equations. A brief procedure for WTC is explained below.
Firstly, perform the leading order analysis, secondly, identify resonance position of arbitrary functions as per Laurent series and finally, existence of a significant number of arbitrary functions are verified without the introduction of movable critical manifold.
3. (2 + 1) Dimensional Non-Linear Schrodinger Equation
In this paper the Painlevé test is performed for (2 + 1) dimensional non-linear Schrödinger equation Equation (1) as described in three steps constituting the Kruskal’s simplification algorithms such as determination of leading order, determination of resonance position and verification of the resonance conditions.
3.1. Leading Order Analysis
For estimation of leading singularity let’s suppose that:
whereas will be all negative integers, and all are functions of . The same assumption were used by Jimbo et al.  and Goldstein & Infeld  the same were previously proposed in private communication by M. Kruskal  . As per original explanation of Weiss et al.,  is a function of . Substituting Equation (2) in Equation (1) followed by balancing the leading term of , and considering that these exponents will be all negative integers, get the only possible values i.e. and .
According to the condition
3.2. Determine the Resonance Position
As per WTC method, the system of Equation (1) which could satisfy the Painlevé test with the general solution is presented as:
This equation contains six arbitrary functions with in addition to the increasing function , where the final solution is a single value regarding the arbitrary movable non-characteristic singular manifold as shown in the Laurent series.
Now substituting Equation (2), Equation (3) and Equation (4) into Equation (1) and equating the coefficients of we get
That value at which the determinant of the system matrix vanishes is called resonance position. It is establish that
Therefore we get the resonances at
3.3. Coefficient of Expansion at the Resonance Position
Resonance at first position i = −1 is observed to relate the arbitrary singularity manifold . For the Equation (1) to pass the Painlevé test, it is required that second resonance position at i = 0, third at i = 2, fourth & fifth at i = 3, 3 and sixth at i = 4 should be identically satisfied. We utilize the Kruskal’s simplification for verification of resonance conditions, such as
With the resolution of Equation (7), Equation (3) is decreased to
Case a: For i = −1.
This point to the arbitrariness of
Case b: For i = 0.
This point to the arbitrariness of or (See Equation (3))
Case c: For i = 1,
Equating the coefficients of we get,
Case d: For i = 2,
Eliminating the coefficients of we get,
From Equation (8) it is observed that Equation (15) is identically satisfied viva Equation (8) and (12).
And hence there are only two Equations (13) and (14) for three unknown variables , and which clarifies that one of them is arbitrary. Therefore the resonance condition at i = 2 is satisfied.
Case e: For i = 3 (double resonance).
By collecting the coefficients of and according to Equation (8) we get,
Obviously due to Equation (12), Equation (18) is identically satisfied. Now we solve from Equation (17) we get,
Now putting Equation (19) into Equation (17) and using Equations (10)-(14), we find that Equation (17) is satisfied identically, which shows that and must be arbitrary if is determined by Equation (19). Hence the resonance condition at i = 3 is satisfied.
Case f: for i = 4.
Vanishing the coefficients of and according to Equation (8) and Equation (11) we get,
Firstly we simplify Equations (20) and (21) for and and then putting into Equation (22) according to Equations (8) (10) (11) (13) (14), and (19) we can find that Equation (22) is identically satisfied.
So any one of , and is arbitrary. Therefore the resonance condition at i = 4 is also satisfied.
Hence we proved that the Equation (1) satisfies the Painlevé test.
4. Associated Bäcklund Transformation
To construct the Bäcklund transformation of Equation (1), let us truncate the Laurent series at the constant level term to give
Here , and , .
where the pair of functions ; and satisfies Equation (1), and hence Equation (23) may be treated as the associated Bäcklund transformation of Equation (1).
5. Bilinear Form
In order to derive the Hirota’s bilinear form we consider the vacuum solutions in Equation (23) then we have
This suggest that we take the Hirota’s bilinear transformation in the form
In Equation (25) U, V and W are functions of variables x, y and z where U is a real function and V, W are complex functions. Using Equation (25) and the Hirota’s bilinear operator Equation (24), Equation (1) can be transforming into bilinear forms as,
Here D is called the bilinear operator defined by  .
The power-series can be further used with the bilinear form obtained in this paper to construct a soliton solution by expanding the dependent variables.
A demonstration on the system which passes the Painlevé test is presented in this paper and further explained that the system has multi-linear variable solutions having arbitrary function, also the associated Bäcklund transformation and bilinear form are obtained directly from the Painlevé test.
The analysis points that our equation having six positions of resonance at i = −1, 0, 2, 3, 3, 4, there are six arbitrary functions at these resonance positions. Therefore our equation conform the Painlevé criterion of integrability as explained in Cauchy-Kovalevskaya theorem. However, there is one further proposition of this paper that Painlevé analysis cannot be used to derive the Lax pair because in the case of coupled system there is still no specialized method for derivation of the Lax pair for coupled non-linear equations.
This work was supported by the National Natural Science Foundation of China (grant No. 11371361), the Innovation Team of Jiangsu Province hosted by China University of Mining and Technology (2014), and the Key Discipline Construction by China University of Mining and Technology (Grant No. XZD 201602). The authors thank to the reviewer’s for their helpful advices.