All articles published by MDPI are made immediately available worldwide under an open access license. No special
permission is required to reuse all or part of the article published by MDPI, including figures and tables. For
articles published under an open access Creative Common CC BY license, any part of the article may be reused without
permission provided that the original article is clearly cited. For more information, please refer to
https://www.mdpi.com/openaccess.
Feature papers represent the most advanced research with significant potential for high impact in the field. A Feature
Paper should be a substantial original Article that involves several techniques or approaches, provides an outlook for
future research directions and describes possible research applications.
Feature papers are submitted upon individual invitation or recommendation by the scientific editors and must receive
positive feedback from the reviewers.
Editor’s Choice articles are based on recommendations by the scientific editors of MDPI journals from around the world.
Editors select a small number of articles recently published in the journal that they believe will be particularly
interesting to readers, or important in the respective research area. The aim is to provide a snapshot of some of the
most exciting work published in the various research areas of the journal.
We propose the approach to constructing semiclassical spectral series for the generalized multidimensional stationary Gross–Pitaevskii equation with a nonlocal interaction term. The eigenvalues and eigenfunctions semiclassically concentrated on a curve are obtained. The curve is described by the dynamic system of moments of solutions to the nonlocal Gross–Pitaevskii equation. We solve the eigenvalue problem for the nonlocal stationary Gross–Pitaevskii equation basing on the semiclassical asymptotics found for the Cauchy problem of the parametric family of linear equations associated with the time-dependent Gross–Pitaevskii equation in the space of extended dimension. The approach proposed uses symmetries of equations in the space of extended dimension.
The collective modes of coherent quantum atom ensembles in the Bose–Einstein condensate (BEC) confined by a trap field in the mean-field approximation are described by the nonlinear Schrödinger equation, the Gross–Pitaevskii equation (GPE) (see, e.g., [1,2,3,4]):
Here, , is the Laplace operator, m is the particle mass, is the nonlinearity parameter, is the external field (trap) potential, the function is the macroscopic wavefunction of the condensate (the order parameter) such that is proportional to the BEC density and its phase gradient is proportional to the BEC velocity. For the sake of simplicity, we omit the argument ℏ of the function in formulas where it does not cause the confusion.The nonlinear interaction operator in (1) is usually chosen as
and it describes the local interaction with where a is the s-wave scattering length. However, to take into account long-range interaction, the nonlocal operators are considered here
The local operator of (2) can be treated as the limit of the nonlocal version (3) with a delta-like function that corresponds to the short-range interaction. Some applications in physics, specifically the BEC with the dipolar-to-dipolar interaction [5,6,7,8], require accounting the long-range interaction when one can notrequiring accounting the long-range interaction do not allow one to get rid of nonlocality for the sake of simplicity. In these cases, the combination of the local and nonlocal terms is often used. We call Equation (1) with the term the local GPE, and with the nonlocal GPE.
Substitution of into Equation (1) yields the time- independent GPE
which describes the stationary states of the BEC. If solutions of (4) are restricted to the unit norm
then eigenvalues of the nonlinear eigenvalue problem (4) and (5) correspond to the chemical potential and eigenfunctions minimize the energy functional
where E is the energy per the BEC particle. The ground state solution corresponds to the global minimum of the functional (6) while the exited states meet the local minima [9,10].
For the weak interaction (a small parameter ), the stationary solutions can be found using the perturbation theory as a correction to solutions of the so-called associated linear Schödinger equation that inherits the GPE linear part and has not the interaction term, i.e., it is the Equation (4) with . Solutions of the GPE obtained in this way are termed solutions with linear counterpart in [11]. The opposite limit of strong nonlinearity was also considered using the Thomas–Fermi approximation when the kinetic term can be ignored and the local GPE (4) takes the form of an algebraic equation. For the nonlocal case, the GPE (4) becomes an integral equation. Moreover, it was shown that there are stationary solutions of the local GPE that do not correspond to any eigenfunctions of the associated linear Schödinger equation for multiwell external potential in [12] (the GPE solutions without linear counterpart).
The minima of the functional (6) can be found numerically by iterative procedure (see, e.g., the review in [13] and references therein). The issue is that the iterative procedure needs a proper initial approximation for the respective stationary states (4). For the weak nonlinearity, the stationary states of the associated linear Schrödinger equation can be used as the initial data, and, for the strong nonlinearity, the initial data can be obtained from the Thomas–Fermi approximation [9,14,15]. Numerical methods are usually limited to calculating only the ground state [13] and do not allow one to construct spectral series in the general case.
The study of the GPE solutions with and without linear counterpart was continued in [16,17,18,19,20].
Note that the solutions with linear counterpart inherit the symmetries of the external field potential while these symmetries are broken for solutions without linear counterpart.
In order to avoid the limitations of the weak or strong nonlinearity, the semiclassical approximation can be used. For example, in [21], the semiclassical approximation was used to study the non-stationary local GPE with the special type of external field and different orders of the nonlinear term.
Here, we propose an approach to the spectral problem (4) that deals with the arbitrary strength of the nonlinearity and is applicable for a wide range of the external fields in n-dimensional nonlocal GPE. Moreover, we study the nonlocal GPE solutions localized on curves reflecting the BEC topological features. We apply the well-known WKB-Maslov theory of semiclassical approximation [22,23,24] to the problem under consideration. The approach to the eigenvalue problem (4) and (5) is based on [25], where the leading term of semiclassical asymptotics for the Cauchy problem has been constructed explicitly in a class of functions localized in a neighborhood of a curve in the phase space of a dynamic system of moments of the nonlocal GPE solution termed the Hamilton–Ehrenfest system. The last one for the n-dimensional manifolds in the phase space was studied in [26,27,28]. As opposed to the linear case ( in (4)), the variance matrix of the solutions constructed in [25] may be bounded due to the nonlinearity. Such solutions of the Cauchy problems are termed the semiclassical soliton-type solutions [29], and we use exactly such particular solutions for the construction of the spectral series of the nonlocal GPE.
Differ toDiffering from the associated linear Schödinger equation in the sense of [11], our method uses a parametric family of the associated linear Gross–Pitaevskii equations that include terms corresponding to both the external field and the interaction. Among the solutions of these equations, the solution of the initial nonlinear GPE (4) can be found.
Our consideration is not limited to the convolution nonlinearity with the kernel , which naturally arise in the BEC application. Although the convolution nonlinearity is usually simpler for calculations, it is unsubstantial for our approach. Moreover, the consideration of the wider class of kernels allows one to apply this approach to the GPE with the modulated nonlinearity that was studied in [18,19].
The WKB-Maslov method sometimes leads to quite complex classical equations (the Hamilton–Ehrenfest system) that need semi-analytical considerations. However, the development of numerical-analytical math packages, such as Wolfram Mathematica, breathed new life into the semiclassical method [30]. In [31], it was applied to the nonlocal generalization of the Fokker–Plank equation with solutions localized on a zero-dimensional manifold (a point). In [32,33], the WKB-Maslov method allowed to construct asymptotic solutions concentrated on curves for the nonlocal Fisher–Kolmogorov–Petrovskii–Piskunov equation. Semiclassical spectral series were constructed in [34,35] for the Hartree type equation with the special type of the external potential .
Our method focuses on formal mathematical structures, so it is applicable for constructing solutions of a wider class of equations. In particular, the approach can be generalized for nonlinear Schrödinger-like equations such as Lugiato–Lefever equation, which is an amplitude equation for an optical pulse in the pumped ring resonator [36]. The interest in this equation was encouraged by discovery of the frequency comb associated with the Kerr temporal solitons [37,38,39]. The results of our paper may be of value for the theoretical study of the Bose–Einstein condensation phenomena that occurs in a gas of atoms and of various quasiparticles [3,40,41,42,43,44]. Note that the conception of condensation is even involved in the study of dark matter [45,46].
The paper is structured as follows. In Section 2, we describe the asymptotic particular solutions of the non-stationary nonlocal GPE concentrated on curves. In Section 3, we construct the asymptotic solutions of the stationary nonlocal GPE based on the Section 2. In Section 3, the semiclassical spectral series are obtained. In Section 4, our formalism is illustrated by an example for the specific 2D stationary nonlocal GPE. The semiclassical solutions of the spectral problem for this equation are obtained in the analytical form. In Section 5, the concluding remarks are given. In Appendix, the auxiliary equations involved are given.
2. Nonlocal Gross–Pitaevskii Equation
The semiclassical approach to the spectral problem (4) can be developed for the nonlocal stationary GPE of a more general form
Here, , , , and is a linear partial differential operator depending on with a symbol smoothly depending on and (). Our approach to the problem (7) relies on the method proposed in [25] for solution of the Cauchy problem in the semiclassical approximation for the non-stationary nonlocal GPE:
This method allows one to construct the asymptotical solutions to the problem (8) in the class of functions semiclassically concentrated on curves of the phase space of the dynamical system of moments of the desired solution.
The core idea of the approach proposed is as follows. The method of [25] provides a practical possibility to select from the set of asymptotical solutions to the Cauchy problem (8), found for various initial functions , those that have the form
In the process of finding the functions (9), we come to and which are the eigenfunctions and eigenvalues of the problem (7).
Consider solutions of (8) concentrated in a vicinity of the one-dimensional manifold (the curve) in the -dimensional phase space , , :
Here, s is the curve parameter and:
For the functions semiclassically concentrated on the curve , the dynamics of the curve in the limits is determined by the -type Hamilton–Ehrenfest system (the term implies that the parameter s is a one-dimensional variable and the vector corresponds to the first order moments of the function , see [25] for details):
Here, is the -norm, which is conserved for exact solutions of the non-stationary nonlocal GPE (8). Hereinafter, for the sake of simplicity, we will use the unit normalization .
The class of functions semiclassically concentrated on curves is given by
where is a class of trajectory-concentrated functions [47]:
In this definition, the function belongs to a Schwarz space in the variable , smoothly depends on t and s, and regularly depends on as ; the function together with functions define the class and are to be determined; the family of hypersurfaces are given by relation [23,24]
Hereinafter, we assume that
This condition ensures the nondegeneracy of the localization curve (10) and the nonsingularity of the BEC density on the curve.
The construction of solutions of the form
to the non-stationary GPE (8) in the class of functions in the semiclassical approximation with the accuracy of can be reduced to the solution of the associated linear GPE:
with the additional algebraic constraint
Here, we denoted
and functions , , where , are the solutions of the second order Hamilton–Ehrenfest system with integration constants determined from the initial condition for the function . In general case, these functions are given by:
Here, is the mean value of the operator with a Weyl symbol in the class :
where is -dimensional vector of variables that complement the variable to form the coordinate system in , is the domain of the variables , and the operator acts subject to the condition .
The particular solution of (18) with auxiliary algebraic constraint (19) can be written in terms of solutions of the variational system, which is the linear Hamiltonian system:
where , is unit symplectic matrix, ( is the unit -matrix). The function is the particular solutions of (23). Let linear independent solutions of (23) be known and satisfy the conditions
Here, is the Kronecker symbol, is the skew-orthogonal product.
Denote . Then, the particular solutions of (18) reads [47]
where N is the normalization coefficient, and
The conditions (16) and (24) ensures the non-degeneracy of matrix . In order to make sure that function (40) satisfies the auxiliary constraint (19) it is enough to check the validity of the identity
which results from the property
following from the definition of the matrix .
Let us introduce first order operators of the form
Due to (24), the following commutation relations hold:
These operators are the symmetry operators of (18) that preserve the property (19) of the solutions. Thus, the operators , are creation operators of (18), and the operators , are annihilation operators. The function (40) is the “vacuum” solution. Let us construct the “occupation number representation” for (18) by determination of the following countable set of solutions:
Here, is a -dimensional multiindex, , , :
In the next Section, we will construct stationary solutions of the GPE based on the solutions (31) and (25).
3. Stationary Solutions of the GPE
In order to obtain solutions that do not depend on time, we construct time-invariant curves. Let us consider curves determined by the following form of functions :
Here, a parameter is the ratio between the periods by s and by t. The curve determined by such the functions is mapped into itself with change of t, i.e., does not depend on t. By the change of variables and , the system (12) can be written in the following form:
In this case, the equation of hypersurfaces has the form or
Let the solutions of the second order Hamilton–Ehrenfest satisfy the relations , , i.e.
The second order Hamilton–Ehrenfest system and the way of obtaining the functions satisfying the relations (36) are considered in the Appendix A.
For the solutions of the -type Hamilton–Ehrenfest system satisfying (33), the variational system (23) reads
where we put
and made the change of the variable . Since is the periodic function, we can pose the Floquet problem for the variational system (37). We know at least one solution of the Floquet problem for (37), namely (the T-periodic solution). Let there be more linear independent solutions of the Floquet problem for (37) that satisfy the following Floquet conditions:
and the skew-orthogonality condition (24). Since only linearly stable are of interest for us, we consider the solutions of (37) with the real Floquet exponent .
If conditions (33), (36) and (38) are met, the particular solutions (25), (31) of (18) read , where
Then, considering (9) and (35), the relation (17) is as follows:
The transition (42) from the extended -dimensional space to the n-dimensional configuration space can give the multi-valued function . However, the function is the single-valued one if the following periodicity condition holds:
In the next Section, we will study what constraints are caused by this condition.
4. Semiclassical Spectral Series
In view of the theorems given in Appendix B and relation (39), the condition (43) for solutions (40) and (41)The misprint in this reference was corrected reads [24]
Let us denote , where
and assume
where regularly depends on ℏ.
We consider the periodic trajectories of the -type Hamilton–Ehrenfest system with the period determined by the following integral of motion of (34):
Then, the quantization condition (44) is satisfied if
The first equation in (48) implicitly determines . The energy per particle for constructed semiclassical spectral series is given by
Note that the quantization condition (48) describes the spectrum well in the linear theory when .
5. Example
In this section, we will illustrate the formalism of our approach by a simple but nontrivial example that admits the analytical integration of the Hamilton–Ehrenfest system. We consider the following two-dimensional nonlocal Gross–Pitaevskii equation:
Here, the units of length and energy are given by such the scales that the term corresponds to in (1) in those notations. This equation can be treated as the nonlocal generalization of the model that describes the BEC in a radially symmetric trap. When , the nonlinear term in (50) corresponds to the repulsive interaction. For the unit normalization of the wave function , the coefficient is proportional to the number of condensate particles. Since the equation has the axial symmetry, the -type Hamilton–Ehrenfest system (34) admits solutions corresponding to the circumference in the configuration space. We look for a solution to the Equation (12) of the form
where the parameter s is the angle in polar coordinates. The -type Hamilton–Ehrenfest system (34) for the GPE (50) reads
From the first equation of the system (52) and (51), we have
For the function (51), the integrals in the Equations (52) are as follows:
where , are modified Bessel functions of the first kind, . Here, we have taken into consideration that and have used the property
As a result, the -type Hamilton–Ehrenfest system yields the following equation for R:
Consider the approximation . It meets the strongly local interaction. The following asymptotic estimate holds for the modified Bessel functions:
Then, we have
This equation has the unique solution. Thus, the right-hand side of (53) satisfies the solvability condition with respect to R for at least sufficiently small (the strongly local interaction). Note that the Equation (53) degenerates into for , i.e., the period T depends on the chemical potential only in the presence of the nonlinearity in (50). The first quantization condition in (48) for the solutions (53) reads
The quantized values of for and for various values of , are shown in Figure 1. Note that the eigenvalue linearly depends on l in the linear case (). The greater the positive value of , the greater a contribution of the nonlinearity to this dependence. The increase in leads to the decrease in the range of eigenvalues for the given range of quantum numbers l. For the large values of , there are not closed trajectories .
The energy per particle the zeroth order approximation with respect to ℏ is as follows:
The Equation (15) for the function for the curve (51) reads
and it has the solution
The semiclassical solutions of (50) for and are obtained in Appendices Appendix C and Appendix D, respectively.
The solution for (the vacuum solution) reads
where is the positive constant defined in Appendix C.
The solution for is given by
In Figure 2 and Figure 3, the plots of and are shown for (56) and (57) when , , , and (it corresponds to , , ).
The phase incursion of along the localization curve (Figure 2 and Figure 3) indicates that we have the vortex solution that corresponds to the solutions of the GPE with 10 quantum vortices [48]. However, since we construct asymptotics only in the neighborhood of the circumference, it do not allows us to obtain the complete information about the phase features in the center.
6. Conclusions
We applied the method proposed in [25] to constructing the solutions of the spectral problem for the stationary nonlocal GPE (7). The considered equation has rather general form with the limitation of the smoothness of the kernel of the nonlocal nonlinearity and of the symbol of the linear part. The solutions of the spectral problem were obtained by choosing the special solutions of the non-stationary nonlocal GPE semiclassically concentrated on time-invariant curves. Such solutions are found in terms of periodic solutions of the second order Hamilton–Ehrenfest system (A2) and (A5) and the set of linear independent skew-orthogonal Floquet solutions of the variational system (37).
The generalization of the method [25] for the stationary GPE is based on the following ideas. The time independence of the Equation (8) operator yields the symmetry (33) with respect to the t- and s-shift for the Equation (12) if the system (34) admits periodic solutions. The property (33) provides to the same symmetry for the second order Hamilton–Ehrenfest system, which leads to the Equations (A2) and (A5), and to the periodicity of coefficients in the variational system (23), that allows us to pose the Floquet problem for (37). Then, the periodic solutions of (A2) yield the symmetry for the associated linear GPE (18) and (20) that effectively reduces the number of independent variables in (18) from (n variables , t, and s) to (n variables , and ). This symmetry together with the Floquet solutions of (37) allows us to obtain the particular solutions (40) and (41) of the associated linear GPE (18) and (20) that give the solutions (42) of the stationary GPE (7) subject to the quantization conditions (48).
The feature of our method is that we do not reduce the initial n-dimensional problem to the one-dimensional. It means that, in general, our approach does not require the GPE to have some special symmetries. Nevertheless, various concepts of symmetries may be required for analytical integration of the auxiliary Equations (34), (37) and (A5) (e.g., [49]).
The solution constructed for the simplest form of the external field reminds the giant hole solution that was studied in [50]. It is known that vortices in the BEC are quantized so that each vortex corresponds to the phase incursion of (see, e.g., [48]). The giant hole solutions are characterized by the proximity of vortex centres so that they form one density pit of the BEC [50] like our example. The indistinguishability of the vortex centres is the consequence of the semiclassical approximation. The centres of vortices likely can be identified within our approach by the consideration of a less trivial localization curve that is an open problem.
The solutions obtained here within the semiclassical approach proposed have the limit (the linear case) and decrease rapidly (exponentially) with distance from the localization curve. In the sense of these properties, the solutions of the nonlocal GPE correspond to ones studied in [11] for the one-dimensional local GPE. However, although the solutions localized on a curve cannot be directly reduced to the one-dimensional solutions found in [11], they explicitly reflect topological features of the BEC.
Note that the numerical calculations for the local and nonlocal GPEs show that, for the delta-like interaction kernel, the solutions of these type weakly change with an increase in the nonlocality coefficient , and qualitatively do not differ from the solutions of the local GPE when is small. The solution is sensitive to the variation of when the localization domain of the interaction potential becomes larger than the localization area of the BEC density. It caused by the sharp decrease in the interaction energy. Therefore, the nonlocal GPE can be treated not only as a separate problem but as the deformation of the local GPE for using our approach.
Author Contributions
Conceptualization, A.V.S. and A.Y.T.; methodology, A.E.K. and A.Y.T.; validation, A.E.K., A.V.S. and A.Y.T.; formal analysis, A.E.K.; investigation, A.E.K. and A.Y.T.; writing—original draft preparation, A.E.K. and A.V.S. and A.Y.T.; writing—review and editing, A.E.K. and A.V.S. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by Russian Foundation for Basic Research (RFBR) and Tomsk region grant number 19-41-700004.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
Not applicable.
Acknowledgments
The work supported by Tomsk Polytechnic University under the International Competitiveness Improvement Program; by Tomsk State University under the International Competitiveness Improvement Program.
Conflicts of Interest
The authors declare no conflict of interest.
Appendix A. The Second Order Hamilton–Ehrenfest System
We can reduce the nonlocal GPE (8) to the parametric family of the associated linear GPEs (18) and (20) due to the possibility of obtaining the dispersion matrix and the first order correction with respect to ℏ for the first order moment from the auxiliary system of equations in the semiclassical approximation, which is called the second order Hamilton–Ehrenfest system. This system was derived in the general case in [25]. In this Appendix, we consider the second order Hamilton–Ehrenfest system for the special type of solutions (9) to (8).
The matrix of second order moments accurate to is determined by the leading term of asymptotics since it corresponds to an operator with the estimate . The formula for reads
where , is -dimensional vector of variables that complement the variable to form the coordinate system in , and is the domain of the variables . Here, the operator acts subject to the condition . Note that the function included in the relation (40) affects only the phase of and commutes with . Hence, the matrix does not affect the result of the formula (A1), i.e., this formula is not recursive. The function given by (A1) satisfy the equation
where is the solution of the matrix variational system
and the matrix is given by (A1). In view of (33) and (36), the system of equations for can be written as
where the operator acts on the arguments of the function and , and the function is given by
It was shown in [25] that . Hence, in the system (A5), we can put
In view of (A7), the system (A5) is the system of linear integro-differential equations. Note that, for every , the function and , which are included in (A5), are different. As opposed to , the function is not determined by the leading term of asymptotics. Thus, our formalism relies on the assumption that system (A5) for every admits T-periodic solution under some initial conditions for .
Appendix B. Auxiliary Theorems
TheoremA1.
The matrix given by (26), (24) and (37) is T-periodic, i.e.
where , are scalar coefficients. Let us make the following change of variables corresponding to the rotation of the vector by the angle :
Then, the variational system reads
where . We seek in the form
In view of , this relation can be written as
The characteristic equation
has the roots
Substitution of into this relation yields
It is clear that we have only aperiodic solutions for , which do not satisfy the Floquet condition (24). Hence, we consider the case when . Denote . Then, the function takes the form
Substitution of this relation into the formula for yields the following fundamental system of solutions:
For , we have
The first solution corresponds to the already known solution . The second solutions does not satisfy the Floquet condition. The third and forth solutions are related by . For , , we have:
It is easy to see that
and
Thus, taking into account the normalization condition , we have the solution
with the Floquet exponent . Substitution of and into (26) yields
and
The coefficient N is obtained from the unit normalization condition:
Then, we have
Note that in our example. Hence, it is not necessary to obtain the whole matrix for the system (A5) but it is enough to obtain only , which is given by
Substituting (40) into (A6), the formula for is as follows:
Since does not depend on t, it does not affect the phase of the wave function in the example considered. In view of (47) and (48), the function is as follows:
The relation (A20) for can be rearranged in the form
Periodic solutions of the second order Hamilton–Ehrenfest system for other can be obtained in the same way. For , the function is as follows:
References
Pitaevskii, L.; Stringari, S. Bose–Einstein Condensation and Superfluidity; Oxford University Press: Oxford, UK, 2016. [Google Scholar]
Dalfovo, F.; Giorgini, S.; Pitaevskii, L.P.; Stringari, S. Theory of Bose–Einstein condensation in trapped gases. Rev. Mod. Phys.1999, 71, 463–512. [Google Scholar] [CrossRef] [Green Version]
Kevrekidis, P.G.; Frantzeskakis, D.J.; Carretero-Gonzalez, R. (Eds.) Emergent Nonlinear Phenomena in Bose–Einstein Condensates: Theory and Experiment; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
Kevrekidis, P.G.; Frantzeskakis, D.J.; Carretero-Gonzales, R. The Defocusing Nonlinear Schrödinger Equation: From Dark Solitons to Vortices and Vortex Rings; SIAM: Philadelphia, PA, USA, 2015. [Google Scholar]
Baranov, M.A. Theoretical progress in many-body physics with ultracold dipolar gases. Phys. Rep.2008, 464, 71–111. [Google Scholar] [CrossRef]
Lahaye, T.; Menotti, C.; Santos, L.; Lewenstein, M.; Pfau, T. The physics of dipolar bosonic quantum gases. Rep. Prog. Phys.2009, 72, 126401. [Google Scholar] [CrossRef]
Beaufils, Q.; Chicireanu, R.; Zanon, T.; Laburthe-Tolra, B.; Maréchal, E.; Vernac, L.; Keller, J.; Gorceix, O. All-optical production of chromium Bose–Einstein condensates. Phys. Rev. A At. Mol. Opt. Phys.2008, 77, 061601. [Google Scholar] [CrossRef] [Green Version]
Jackson, A.D.; Kavoulakis, G.M.; Pethick, C.J. Solitary waves in clouds of Bose–Einstein condensed atoms. Phys. Rev. A At. Mol. Opt. Phys.1998, 58, 2417–2422. [Google Scholar] [CrossRef] [Green Version]
Pethick, C.J.; Smith, H. Bose–Einstein Condensation in Dilute Gases; Cambridge University Press: Cambridge, UK, 2008; pp. 1–569. [Google Scholar]
D’Agosta, R.; Malomed, B.A.; Presilla, C. Stationary solutions of the Gross–Pitaevskii equation with linear counterpart. Phys. Lett. Sect. A Gen. At. Solid State Phys.2000, 275, 424–434. [Google Scholar] [CrossRef] [Green Version]
D’Agosta, R.; Presilla, C. States without a linear counterpart in Bose–Einstein condensates. Phys. Rev. A At. Mol. Opt. Phys.2002, 65, 436091–436097. [Google Scholar] [CrossRef] [Green Version]
Antoine, X.; Duboscq, R. GPELab, a Matlab toolbox to solve Gross–Pitaevskii equations I: Computation of stationary solutions. Comput. Phys. Commun.2014, 185, 2969–2991. [Google Scholar] [CrossRef]
Bao, W.; Du, Q. Computing the ground state solution of Bose–Einstein condensates by a normalized gradient flow. SIAM J. Sci. Comput.2004, 25, 1674–1697. [Google Scholar] [CrossRef] [Green Version]
Bao, W.; Tang, W. Ground-state solution of Bose–Einstein condensate by directly minimizing the energy functional. J. Comput. Phys.2003, 187, 230–254. [Google Scholar] [CrossRef] [Green Version]
Carr, L.D.; Clark, C.W.; Reinhardt, W.P. Stationary solutions of the one-dimensional nonlinear Schrodinger equation. I. Case of repulsive nonlinearity. Phys. Rev. A At. Mol. Opt. Phys.2000, 62, 063610. [Google Scholar] [CrossRef] [Green Version]
Carr, L.D.; Clark, C.W.; Reinhardt, W.P. Stationary solutions of the one-dimensional nonlinear Schrodinger equation. II. Case of attractive nonlinearity. Phys. Rev. A At. Mol. Opt. Phys.2000, 62, 063611. [Google Scholar] [CrossRef] [Green Version]
Alfimov, G.L.; Gegel, L.A.; Lebedev, M.E.; Malomed, B.A.; Zezyulin, D.A. Localized modes in the Gross–Pitaevskii equation with a parabolic trapping potential and a nonlinear lattice pseudopotential. Commun. Nonlinear Sci. Numer. Simul.2019, 66, 194–207. [Google Scholar] [CrossRef] [Green Version]
Mallory, K.; Van Gorder, R.A. Stationary solutions for the nonlinear Schrödinger equation modeling three-dimensional spherical Bose–Einstein condensates in general potentials. Phys. Rev. E Stat. Nonlinear Soft Matter Phys.2015, 92, 013201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Sacchetti, A. Solution to the double-well nonlinear Schrödinger equation with Stark-type external field. J. Phys. A Math. Theor.2015, 48, 035303. [Google Scholar] [CrossRef]
Maslov, V. Operational Methods; Mir Publishers: Moscow, Russia, 1976. [Google Scholar]
Maslov, V. The Complex WKB Method for Nonlinear Equations. I. Linear Theory; Birkhauser Verlag: Basel, Switzerland, 1994. [Google Scholar]
Belov, V.V.; Dobrokhotov, S.Y. Semiclassical maslov asymptotics with complex phases. I. General approach. Theor. Math. Phys.1992, 92, 843–868. [Google Scholar] [CrossRef]
Shapovalov, A.V.; Kulagin, A.E.; Trifonov, A.Y. The Gross–Pitaevskii equation with a nonlocal interaction in a semiclassical approximation on a curve. Symmetry2020, 12, 201. [Google Scholar] [CrossRef] [Green Version]
Maslov, V. Complex Markov Chains and the Feynman Path Integral; Nauka: Moscow, Russia, 1976. [Google Scholar]
Maslov, V. Equations of the self-consistent field. J. Sov. Math.1978, 11, 123–195. [Google Scholar] [CrossRef]
Karasev, M.V.; Maslov, V.P. Algebras with general commutation relations and their applications. II. Unitary-nonlinear operator equations. J. Sov. Math.1981, 15, 273–368. [Google Scholar] [CrossRef]
Belov, V.V.; Kondrat’eva, M.F.; Smirnova, E.I. Semiclassical soliton-type solutions of the Hartree equation. Dokl. Math.2007, 76, 775–779. [Google Scholar] [CrossRef]
Anikin, A.Y.; Dobrokhotov, S.Y.; Nazaikinskii, V.E.; Tsvetkova, A.V. Uniform Asymptotic Solution in the Form of an Airy Function for Semiclassical Bound States in One-Dimensional and Radially Symmetric Problems. Theor. Math. Phys. Russ. Fed.2019, 201, 1742–1770. [Google Scholar] [CrossRef]
Shapovalov, A.V.; Rezaev, R.O.; Trifonov, A.Y. Symmetry operators for the Fokker-Plank-Kolmogorov equation with nonlocal quadratic nonlinearity. Symmetry Integr. Geom. Methods Appl. SIGMA2007, 3, 5. [Google Scholar] [CrossRef] [Green Version]
Levchenko, E.A.; Shapovalov, A.V.; Trifonov, A.Y. Asymptotics semiclassically concentrated on curves for the nonlocal Fisher-Kolmogorov-Petrovskii-Piskunov equation. J. Phys. A Math. Theor.2016, 49, 305203. [Google Scholar] [CrossRef]
Levchenko, E.A.; Shapovalov, A.V.; Trifonov, A.Y. Pattern formation in terms of semiclassically limited distribution on lower dimensional manifolds for the nonlocal Fisher-Kolmogorov-Petrovskii-Piskunov equation. J. Phys. A Math. Theor.2014, 47, 025209. [Google Scholar] [CrossRef] [Green Version]
Pereskokov, A.V. Semiclassical asymptotics of the spectrum near the lower boundary of spectral clusters for a Hartree-type operator. Math. Notes2017, 101, 1009–1022. [Google Scholar] [CrossRef]
Pereskokov, A.V. Asymptotics of the Spectrum of a Two-dimensional Hartree Type Operator Near Upper Boundaries of Spectral Clusters. Asymptotic Solutions Located Near a Circle. J. Math. Sci.2017, 226, 517–530. [Google Scholar] [CrossRef]
Lugiato, L.A.; Lefever, R. Spatial dissipative structures in passive optical systems. Phys. Rev. Lett.1987, 58, 2209–2211. [Google Scholar] [CrossRef] [PubMed]
Chembo, Y.K.; Menyuk, C.R. Spatiotemporal Lugiato-Lefever formalism for Kerr-comb generation in whispering-gallery-mode resonators. Phys. Rev. A At. Mol. Opt. Phys.2013, 87, 053852. [Google Scholar] [CrossRef] [Green Version]
Coen, S.; Randle, H.G.; Sylvestre, T.; Erkintalo, M. Modeling of octave-spanning Kerr frequency combs using a generalized mean-field Lugiato-Lefever model. Opt. Lett.2013, 38, 37–39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Godey, C.; Balakireva, I.V.; Coillet, A.; Chembo, Y.K. Stability analysis of the spatiotemporal Lugiato-Lefever model for Kerr optical frequency combs in the anomalous and normal dispersion regimes. Phys. Rev. A At. Mol. Opt. Phys.2014, 89, 063814. [Google Scholar] [CrossRef] [Green Version]
Butov, L.V.; Lai, C.W.; Ivanov, A.L.; Gossard, A.C.; Chemla, D.S. Towards Bose–Einstein condensation of excitons in potential traps. Nature2002, 417, 47–52. [Google Scholar] [CrossRef] [PubMed]
Sob’Yanin, D.N. Bose–Einstein condensation of light: General theory. Phys. Rev. E Stat. Nonlinear Soft Matter Phys.2013, 88, 022132. [Google Scholar] [CrossRef] [Green Version]
Berman, O.L.; Kezerashvili, R.Y. Superfluidity of dipolar excitons in a transition metal dichalcogenide double layer. Phys. Rev. B2017, 96, 094502. [Google Scholar] [CrossRef] [Green Version]
Berman, O.L.; Gumbs, G.; Kezerashvili, R.Y. Bose–Einstein condensation and superfluidity of dipolar excitons in a phosphorene double layer. Phys. Rev. B2017, 96, 014505. [Google Scholar] [CrossRef] [Green Version]
Elizalde, E.; Lidsey, J.E.; Nojiri, S.; Odintsov, S.D. Born-Infeld quantum condensate as dark energy in the universe. Phys. Lett. Sect. B Nucl. Elem. Part. High Energy Phys.2003, 574, 1–7. [Google Scholar] [CrossRef] [Green Version]
Addazi, A.; Capozziello, S.; Odintsov, S. Born–Infeld condensate as a possible origin of neutrino masses and dark energy. Phys. Lett. Sect. B Nucl. Elem. Part. High Energy Phys.2016, 760, 611–616. [Google Scholar] [CrossRef] [Green Version]
Bagrov, V.G.; Belov, V.V.; Trifonov, A.Y. Semiclassical trajectory-coherent approximation in quantum mechanics I. High-order corrections to multidimensional time-dependent equations of Schrödinger type. Ann. Phys.1996, 246, 231–290. [Google Scholar] [CrossRef]
Fetter, A.L.; Svidzinsky, A.A. Vortices in a trapped dilute Bose–Einstein condensate. J. Phys. Condens. Matter2001, 13, R135–R194. [Google Scholar] [CrossRef] [Green Version]
Obukhov, V.V. Algebra of symmetry operators for Klein-Gordon-Fock equation. Symmetry2021, 13, 727. [Google Scholar] [CrossRef]
Kasamatsu, K.; Tsubota, M.; Ueda, M. Giant hole and circular superflow in a fast rotating Bose–Einstein condensate. Phys. Rev. A At. Mol. Opt. Phys.2002, 66, 053606. [Google Scholar] [CrossRef] [Green Version]
Figure 1.
The values of for various quantum numbers l and for , (black), , (red), , (blue).
Figure 1.
The values of for various quantum numbers l and for , (black), , (red), , (blue).
Figure 2.
Density plots of .
Figure 2.
Density plots of .
Figure 3.
Density plots of .
Figure 3.
Density plots of .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Kulagin, A.E.; Shapovalov, A.V.; Trifonov, A.Y.
Semiclassical Spectral Series Localized on a Curve for the Gross–Pitaevskii Equation with a Nonlocal Interaction. Symmetry2021, 13, 1289.
https://doi.org/10.3390/sym13071289
AMA Style
Kulagin AE, Shapovalov AV, Trifonov AY.
Semiclassical Spectral Series Localized on a Curve for the Gross–Pitaevskii Equation with a Nonlocal Interaction. Symmetry. 2021; 13(7):1289.
https://doi.org/10.3390/sym13071289
Chicago/Turabian Style
Kulagin, Anton E., Alexander V. Shapovalov, and Andrey Y. Trifonov.
2021. "Semiclassical Spectral Series Localized on a Curve for the Gross–Pitaevskii Equation with a Nonlocal Interaction" Symmetry 13, no. 7: 1289.
https://doi.org/10.3390/sym13071289
APA Style
Kulagin, A. E., Shapovalov, A. V., & Trifonov, A. Y.
(2021). Semiclassical Spectral Series Localized on a Curve for the Gross–Pitaevskii Equation with a Nonlocal Interaction. Symmetry, 13(7), 1289.
https://doi.org/10.3390/sym13071289
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
Article Metrics
No
No
Article Access Statistics
For more information on the journal statistics, click here.
Multiple requests from the same IP address are counted as one view.
Kulagin, A.E.; Shapovalov, A.V.; Trifonov, A.Y.
Semiclassical Spectral Series Localized on a Curve for the Gross–Pitaevskii Equation with a Nonlocal Interaction. Symmetry2021, 13, 1289.
https://doi.org/10.3390/sym13071289
AMA Style
Kulagin AE, Shapovalov AV, Trifonov AY.
Semiclassical Spectral Series Localized on a Curve for the Gross–Pitaevskii Equation with a Nonlocal Interaction. Symmetry. 2021; 13(7):1289.
https://doi.org/10.3390/sym13071289
Chicago/Turabian Style
Kulagin, Anton E., Alexander V. Shapovalov, and Andrey Y. Trifonov.
2021. "Semiclassical Spectral Series Localized on a Curve for the Gross–Pitaevskii Equation with a Nonlocal Interaction" Symmetry 13, no. 7: 1289.
https://doi.org/10.3390/sym13071289
APA Style
Kulagin, A. E., Shapovalov, A. V., & Trifonov, A. Y.
(2021). Semiclassical Spectral Series Localized on a Curve for the Gross–Pitaevskii Equation with a Nonlocal Interaction. Symmetry, 13(7), 1289.
https://doi.org/10.3390/sym13071289
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.