Advances in High Energy Physics
Volume 2020 (2020), Article ID 7010957, 6 pages
https://doi.org/10.1155/2020/7010957

Research Article

Analytical Solution for the Gross-Pitaevskii Equation in Phase Space and Wigner Function

A. X. Martins,¹ R. A. S. Paiva,¹ G. Petronilo,¹ R. R. Luz,¹ R. G. G. Amorim,1,2,3 S. C. Ulhoa,1,3 and T. M. R. Filho¹

¹International Center of Physics, Instituto de Física, Universidade de Brasília, 70910-900 Brasília, DF, Brazil
²Faculdade Gama, Universidade de Brasília, 70910-900 Brasília, DF, Brazil
³Canadian Quantum Research Center, 204-3002 32 Ave Vernon, BC, Canada V1T 2L7

Correspondence should be addressed to R. G. G. Amorim; ronniamorim@gmail.com

Received 30 January 2020; Accepted 6 April 2020; Published 30 April 2020

Academic Editor: Salvatore Mignemi

Copyright © 2020 A. X. Martins et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In this work, we study symplectic unitary representations for the Galilei group. As a consequence a nonlinear Schrödinger equation is derived in phase space. The formalism is based on the noncommutative structure of the star product, and using the group theory approach as a guide a physically consistent theory is constructed in phase space. The state is described by a quasi-probability amplitude that is in association with the Wigner function. With these results, we solve the Gross-Pitaevskii equation in phase space and obtained the Wigner function for the system considered.

1. Introduction

A relevant equation that describes a variety physical phenomena, as a Bose-Einstein condensed, is the Gross-Pitaevskii equation [1]. The Gross-Pitaevskii model is an extension of the Schrödinger equation, and it is given by [2, 3]

where represents mass, is the interaction potential, and is the intensity of atomic interaction. The study of solutions for this equation is relevant in both ways, theoretical and applied viewpoints. In addition, an important case for the Gross-Pitaevskii system is their approach in quantum phase space, particularly the calculation of Wigner function for this system, in which it is not known in the literature.

In this context, the first formalism to quantum mechanics is phase space which was introduced by Wigner notion of phase space in 1932 [4]. He was motivated by the problem of finding a way to improve the quantum statistical mechanics. Wigner introduced his formalism by using a kind of Fourier transform of the density matrix, giving rise to what is nowadays called the Wigner function, where are the coordinates of a phase space manifold () [4–7]. The Wigner function has the same content of usual wave function obtained by the Schrödinger equation. However, the Wigner function is identified as a quasiprobability density in the sense that is real but not positively definite and as such cannot be interpreted as a probability. However, the integrals and are (true) distribution functions. The calculation of Wigner function is based in the following steps: (i) first, the Schrödinger equation for a specific potential must be solved; (ii) in sequence, using the solutions founded, the matrix density elements are calculated; (iii) finally, a kind of Fourier transform of matrix elements must be performed. We notice that in the Wigner approach it is complicated to treat nonlinear potentials such as the Gross-Pitaevskii system. For this reason, other methods to quantum mechanics in phase space are developed in the literature. There is an alternative method based on the following property of Wigner formalism: in the Wigner function approach, each operator, , defined in the Hilbert space, , is associated with a function, , in . This procedure is precisely specified by a mapping , such that, the associative algebra of operators defined in turns out to be an algebra in given by where the star product, , is defined by

The result is a noncommutative structure in that has been explored in different ways [5–27].

Using star operators defined in Wigner formalism, unitary representations of symmetry Lie groups have been developed on a symplectic [28–32]. The unitary representation of Galilei group leads to the Schrödinger equation in phase space. In the analog procedure, the scalar Lorentz group for spin 0 and spin 1/2 leads to the Klein-Gordon and Dirac equations in phase space. In both cases, relativistic and nonrelativistic, the wave functions are closely associated with the Wigner function [28, 29]. In terms of nonrelativistic quantum mechanics, the proposed formalism has been used to treat a nonlinear oscillator perturbatively, to study the notion of coherent states and to introduce a nonlinear Schrödinger equation from the point of view of phase space. In this context, there are a few examples of analytical solutions such as the harmonic oscillator [33], the Hydrogen atom [34], and some spin systems [35–37]. In the present work, we apply this symplectic formalism to find the Wigner function for the Gross-Pitaevskii model. We find an analytical solution for the wave function but the Wigner function is calculated up to a given order of approximation of the star product.

The paper is organized as follows: in Section 2, we write the nonlinear equation in phase space and we present the relation between phase space amplitude and Wigner function. In Section 3, we solve the Gross-Pitaevskii equation in phase space and calculate the Wigner function. In Section 4, we plot graphs of Wigner function and calculate nonnegativity parameter associated to the system. Finally, some closing comments are given in Section 5.

2. Nonlinear Schrödinger Equation in Phase Space

Using the star operators, , we define the momentum and position operators, respectively, by

Then, we introduce the following operators:

These operators, given in Equations (3), (4), (5), and (6), are defined in the Hilbert space, , constructed with complex functions in the phase space [28], and satisfy the set of commutation relations for the Galilei-Lie algebra, that is,

with all other commutation relations being null. This is the Galilei-Lie algebra with a central extension characterized by . The operators defining the Galilei symmetry , , , and are the generators of translations, boost, rotations, and time translations, respectively. and can be taken to be the physical observable of position and momentum. To be consistent, generators are interpreted as the angular momentum operator, and is taken as the Hamiltonian operator. The Casimir invariants of the Lie algebra are given by

where describes the Hamiltonian of a free particle and is associated with the spin degrees of freedom. First, we study the scalar representation, i.e. spin zero.

Using the time-translation generator, , we derive the time-evolution equation for , i.e.,

which is the Schrödinger equation in phase space [28]. The function is defined in a Hilbert space associated to phase space [28].

The association of with the Wigner function is given by [28]

This function satisfies the Liouville-von Neumann equation [28]. This provides a complete set of physical rules to interpret representations and opens the way to study other improvements. In this sense, the nonlinear Schrödinger equation in phase space is given by

where is intensity of atomic interaction. This equation describes several physical systems; in particular, it is used to study the Bose-Einstein condensation. Equation (20) is derived from the Lagrangian density

In the next section, we solve this equation and then calculate the Wigner function for this system by expanding the star product up to the second order of approximation.

3. Solution of Gross-Pitaevskii Equation and Wigner Function

In this section, we present a solution for the nonlinear Schrödinger equation in phase space and the associated Wigner function.

In order, the Schrödinger equation in phase space is written by

In this work, we address the stationary equation without external potential, i.e., . In this way, the equation above becomes

We now consider the solution of Equation (15) in the regions of constant potential, which may be taken to be without loss of generality. We note first that if vanishes anywhere in an interval, as for example at the edges of the box, then may be taken to be purely real throughout that interval. This can be done only if

Thus, we may remove the absolute value symbol in Equation (15). So letting , the equation becomes an ordinary nonlinear equation for a real function:

Letting , the solution of Equation (17) is

with , is the Jacobian elliptic function, and and will be determined by the boundary conditions, while and will be determined by substituting (18) into (17) and normalization. The boundary conditions

The boundary condition at the origin can be satisfied by taking . The function is periodic in with a period of , where is the elliptic integral of the first kind. Thus, the boundary equations are satisfied if , where . The number of nodes in the solution is . We then solve Equation (17) substituting (18), and using the Jacobian elliptic identities, this results in the equation for the amplitude and energy ,

Equation (18) becomes

The wave-function and the energy are determined up to factor . This result is similar to what is obtained in configuration space [1]. Using Equation (22), we calculate Wigner function for the Gross-Pitaevskii system by Equation (11). There are in fact some different approaches from our proposal to calculate the Wigner function, such as the use of the parity operator for this purpose [38, 39]. In this article, the Wigner function is calculated up to the second order in the expansion of the star product. In the next section, this function under these terms is analyzed.

4. Analysis of Solution

In this section, we plot the graphics to Wigner function founded in the section above and calculate the negativity parameter for this system.

4.1. Particle in a Box Limit

Both the zero density linear limit and the highly excited-state limit give the particle in a box limit type solution. Mathematically, and . Physically, . In this limit. and , so Equation (21) becomes

which converges to the linear quantum mechanics particle in a box as . One may obtain these results from the first order perturbation theory [31]. The behavior of the Wigner function given in Equation (22) for the three first energy levels can be visualized in Figures 1–3.

Figure 1: Wigner function for Gross-Pitaevskii model; .

Figure 2: Wigner function for Gross-Pitaevskii model; .

Figure 3: Wigner function for Gross-Pitaevskii model; .

Using the Wigner function, the negative parameter for the system is calculated. The results are presented in Table 1.As the negativity parameter is correlated with the nonclassicality of the system, it seems that in the limit where the BCE behaviors are in accordance with classical mechanics which is a consistent result. It is important to note that such classic behavior may be a consequence of the second order expansion of the star product. The investigation of a general solution is necessary to understand the role of the negativity parameter. Although the expansion of the star product is a legitimate procedure given the order of magnitude of the Planck constant, in a future work, we hope to apply a more precise method to understand which order of expansion is ideal.

Table 1: Parameter of negativity for , 200, and 300.

100
0

200
0

300
0

5. Concluding Remarks

In this work, we studied the nontrivial problem of the Gross-Pitaevskii equation in phase space, in which is a case of nonlinear Schrödinger. The Wigner function for this system was obtained. In our knowledge, this is the first time that an analytical solution of the Wigner function for the Gross-Pitaevskii system appears in the literature. The particle in a box solution is studied for the physical meaningful limit of the solution, . We studied the parameter of negativity of the system and concluded that in these limits it appears to have a purely classical behavior.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was partially supported by CAPES and CNPq of Brazil.

References

L. D. Carr, C. W. Clark, and W. P. Reinhardt, “Stationary solutions of the one-dimensional nonlinear Schrödinger equation. I. Case of repulsive nonlinearity,” Physical Review A, vol. 62, no. 6, 2000.

J. Kasprzak, M. Richard, S. Kundermann et al., “Bose-Einstein condensation of exciton polaritons,” Nature, vol. 443, no. 7110, pp. 409–414, 2006.

M. R. Andrews, C. G. Townsend, H. J. Miesner, D. S. Durfee, D. M. Kurn, and W. Ketterle, “Observation of interference between two Bose condensates,” Science, vol. 275, no. 5300, pp. 637–641, 1997.

E. Wigner, “On the Quantum corrections to thermodynamics equilibrium,” Physical Review, vol. 40, no. 5, pp. 749–759, 1932.

M. Hillery, R. F. O. 'Connell, M. O. Scully, and E. P. Wigner, “Distribution functions in physics: fundamentals,” Physics Reports, vol. 106, no. 3, pp. 121–167, 1984.

Y. S. Kim and M. E. Noz, Phase Space Picture and Quantum Mechanics-Group Theoretical Approach, W. Scientific, London, 1991.

T. Curtright, D. Fairlie, and C. Zachos, “Features of time-independent Wigner functions,” Physical Review D, vol. 58, no. 2, 1998.

C. Zachos, “Deformation quantization: quantum mechanics lives and works in phase-space,” International Journal of Modern Physics A: Particles and Fields; Gravitation; Cosmology; Nuclear Physics, vol. 17, no. 3, pp. 297–316, 2002.

F. C. Khanna, A. P. C. Malbouisson, J. M. C. Malbouisson, and A. E. Santana, Thermal Quantum Field Theory: Algebraic Aspects and Applications, W. Scientific, Singapore, 2009.

J. D. Vianna, M. C. B. Fernandes, and A. E. Santana, “Galilean-covariant Clifford algebras in the phase-space representation,” Foundations of Physics, vol. 35, no. 1, pp. 109–129, 2005.

H. Weyl, “Quantenmechanik und Gruppentheorie,” Zeitschrift für Physik, vol. 46, no. 1-2, pp. 1–46, 1927.

J. E. Moyal, “Quantum mechanics as a statistical theory,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 45, no. 1, pp. 99–124, 1949.

S. A. Smolyansky, A. V. Prozorkevich, G. Maino, and S. G. Mashnik, “A covariant generalization of the real-time Green's functions method in the theory of kinetic equations,” Annals of Physics, vol. 277, no. 2, pp. 193–218, 1999.

T. Curtright and C. Zachos, “Wigner trajectory characteristics in phase space and field theory,” Journal of Physics A, vol. 32, no. 5, pp. 771–779, 1999.

I. Galaviz, H. Garca-Compeán, M. Przanowski, and F. J. Turrubiates, Weyl-Wigner-Moyal for Fermi Classical Systems, http://arxiv.org/abs/0612245v1.

J. Dito, “Star‐products and nonstandard quantization for Klein–Gordon equation,” Journal of Mathematical Physics, vol. 33, no. 2, pp. 791–801, 1992.

G. Torres‐Vega and J. H. Frederick, “Quantum mechanics in phase space: new approaches to the correspondence principle,” The Journal of Chemical Physics, vol. 93, no. 12, pp. 8862–8874, 1990.

M. A. De Gosson, “Symplectically covariant Schrödinger equation in phase space,” Journal of Physics A: Mathematical and General, vol. 38, no. 42, pp. 9263–9287, 2005.

M. A. De Gosson, “Semi-classical propagation of wavepackets for the phase space Schrödinger equation: interpretation in terms of the Feichtinger algebra,” Journal of Physics A: Mathematical and Theoretical, vol. 41, no. 9, p. 095202, 2008.

D. Galetti and A. F. R. de Toledo Piza, “Symmetries and time evolution in discrete phase spaces: a soluble model calculation,” Physica A, vol. 214, no. 2, pp. 207–228, 1995.

L. P. Horwitz, S. Shashoua, and W. C. Schive, “A manifestly covariant relativistic Boltzmann equation for the evolution of a system of events,” Physica A, vol. 161, no. 2, pp. 300–338, 1989.

P. R. Holland, “Relativistic algebraic spinors and quantum motions in phase space,” Foundations of Physics, vol. 16, no. 8, pp. 701–719, 1986.

M. C. B. Fernandes, A. E. Santana, and J. D. M. Vianna, “Galilean Duffin Kemmer Petiau algebra and symplectic structure,” Journal of Physics A: Mathematical and General, vol. 36, no. 13, pp. 3841–3854, 2003.

A. E. D. Santana, A. M. Neto, J. D. M. Vianna, and F. C. Khanna, “Symmetry groups, density-matrix equations and covariant Wigner functions,” Physica A: Statistical Mechanics and its Applications, vol. 280, no. 3-4, pp. 405–436, 2000.

D. Bohm and B. J. Hiley, “On a quantum algebraic approach to a generalized phase space,” Foundations of Physics, vol. 11, no. 3-4, pp. 179–203, 1981.

M. C. B. Andrade, A. E. Santana, and J. D. M. Vianna, “Poincaré-Lie algebra and relativistic phase space,” Journal of Physics A: Mathematical and General, vol. 33, no. 22, pp. 4015–4024, 2000.

M. A. Alonso, G. S. Pogosyan, and K. B. Wolf, “Wigner functions for curved spaces. I. On hyperboloids,” Journal of Mathematical Physics, vol. 43, no. 12, pp. 5857–5871, 2002.

M. D. Oliveira, M. C. B. Fernandes, F. C. Khanna, A. E. D. Santana, and J. D. M. Vianna, “Symplectic quantum mechanics,” Annals of Physics, vol. 312, no. 2, pp. 492–510, 2004.

R. G. G. Amorim, M. C. B. Fernandes, F. C. Khanna, A. E. Santana, and J. D. M. Vianna, “Non-commutative geometry and symplectic field theory,” Physics Letters A, vol. 361, no. 6, pp. 464–471, 2007.

R. G. G. Amorim, F. C. Khanna, A. E. Santana, and J. D. M. Vianna, “Perturbative symplectic field theory and Wigner function,” Physica A, vol. 388, no. 18, pp. 3771–3778, 2009.

M. C. B. Fernandes, F. C. Khanna, M. G. R. Martins, A. E. Santana, and J. D. M. Vianna, “Non-linear Liouville and Shrödinger equations in phase space,” Physica A, vol. 389, no. 17, pp. 3409–3419, 2010.

L. M. Abreu, A. E. D. Santana, and A. Ribeiro Filho, “The Cangemi–Jackiw manifold in high dimensions and symplectic structure,” Annals of Physics, vol. 297, no. 2, pp. 396–408, 2002.

P. Campos, M. G. R. Martins, M. C. B. Fernandes, and J. D. M. Vianna, “Quantum mechanics on phase space: the hydrogen atom and its Wigner functions,” Ann. Phys., vol. 390, pp. 60–70, 2018.

J. P. Dahl and M. Springborg, “Wigner's phase space function and atomic structure,” Molecular Physics, vol. 47, no. 5, pp. 1001–1019, 1982.

J. C. Várilly and J. M. Gracia-Bondía, “The moyal representation for spin,” Annals of Physics, vol. 190, no. 1, pp. 107–148, 1989.

B. Koczor, R. Zeier, and S. J. Glaser, “Continuous phase spaces and the time evolution of spins: star products and spin-weighted spherical harmonics,” Journal of Physics A: Mathematical and Theoretical, vol. 52, no. 5, 2019.

B. Koczor, R. Zeier, and S. J. Glaser, “Time evolution of coupled spin systems in a generalized Wigner representation,” Annals of Physics, vol. 408, pp. 1–50, 2019.

A. Grossmann, “Parity operator and quantization of δ-functions,” Communications in Mathematical Physics, vol. 48, no. 3, pp. 191–194, 1976, http://arxiv.org/abs/1811.05872.

B. Koczor, F. vom Ende, M. de Gosson, S. J. Glaser, and R. Zeier, 2018, http://arxiv.org/abs/1811.05872.