Next Article in Journal
A New Feature Selection Method Based on a Self-Variant Genetic Algorithm Applied to Android Malware Detection
Previous Article in Journal
General Norm Inequalities of Trapezoid Type for Fréchet Differentiable Functions in Banach Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Semiclassical Spectral Series Localized on a Curve for the Gross–Pitaevskii Equation with a Nonlocal Interaction

by
Anton E. Kulagin
1,2,*,†,
Alexander V. Shapovalov
3,*,† and
Andrey Y. Trifonov
1,†
1
Department of Mathematics and Computer Sciences, Tomsk Polytechnic University, 30 Lenin A., 634050 Tomsk, Russia
2
V.E. Zuev Institute of Atmospheric Optics, Siberian Branch of the Russian Academy of Sciences, 1 Academician Zuev Sq., 634055 Tomsk, Russia
3
Department of Theoretical Physics, Tomsk State University, 1 Novosobornaya Sq., 634050 Tomsk, Russia
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2021, 13(7), 1289; https://doi.org/10.3390/sym13071289
Submission received: 17 June 2021 / Revised: 9 July 2021 / Accepted: 11 July 2021 / Published: 17 July 2021
(This article belongs to the Section Physics)

Abstract

:
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.

1. Introduction

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]):
i t 2 2 m Δ + V t r ( x ) + ϰ W ^ [ Ψ ] ( x , t ) Ψ ( x , t ) = 0 .
Here, t = t , Δ is the Laplace operator, m is the particle mass, ϰ is the nonlinearity parameter, V t r ( x ) is the external field (trap) potential, the function Ψ ( x , t ) is the macroscopic wavefunction of the condensate (the order parameter) such that | Ψ ( x , t ) | 2 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 Ψ ( x , t , ) in formulas where it does not cause the confusion.The nonlinear interaction operator W ^ [ Ψ ] ( x , t ) in (1) is usually chosen as
W ^ [ Ψ ] ( x , t ) = W ^ l o c a l [ Ψ ] ( x , t ) = g | Ψ ( x , t ) | 2
and it describes the local interaction with g = 4 π 2 a / m where a is the s-wave scattering length. However, to take into account long-range interaction, the nonlocal operators are considered here
W ^ [ Ψ ] ( x , t ) = W ^ n o n l o c a l [ Ψ ] ( x , t ) = R n W ( x , y ) | Ψ ( y , t ) | 2 d y .
The local operator W ^ l o c a l of (2) can be treated as the limit of the nonlocal version (3) with a delta-like function W ( x , y ) 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 W ^ l o c a l the local GPE, and with W ^ n o n l o c a l the nonlocal GPE.
Substitution of Ψ ( x , t ) = exp i μ t ψ ( x ) into Equation (1) yields the time- independent GPE
2 2 m Δ + V t r ( x ) + ϰ W ^ [ ψ ] ( x ) ψ ( x ) = μ ψ ( x ) ,
which describes the stationary states of the BEC. If solutions ψ ( x ) of (4) are restricted to the unit norm
| | ψ ( x ) | | 2 = R n | ψ ( x ) | 2 d x = 1 ,
then eigenvalues μ of the nonlinear eigenvalue problem (4) and (5) correspond to the chemical potential and eigenfunctions ψ ( x ) minimize the energy functional
E [ ψ ] = R n 2 2 m | ψ ( x ) | 2 + V t r ( x ) | ψ ( x ) | 2 + ϰ 2 W ^ [ ψ ] ( x ) | ψ ( x ) | 2 d x ,
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 ϰ = 0 . 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 ( ϰ = 0 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 W ( x , y ) = W ( x y ) , 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 W ( x , y ) 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 V t r .
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
H ( z ^ ) [ ψ ] μ ( ) ψ ( x ) = 0 , H ( z ^ ) [ ψ ] = V ( z ^ ) + ϰ R n W ( x , y ) | ψ ( y ) | 2 d y .
Here, z ^ = ( p ^ , x ) , x R n , p ^ = i x , and V ( z ^ ) is a linear partial differential operator depending on z ^ with a symbol V ( z ) = V ( p , x ) smoothly depending on p and x ( p R n ). 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:
i t + H ( z ^ ) [ Ψ ] Ψ ( x , t ) = 0 , Ψ ( x , t ) | t = 0 = ψ ( x ) .
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 Ψ ( x , t ) to the Cauchy problem (8), found for various initial functions ψ ( x ) , those Ψ ( x , t ) that have the form
Ψ ( x , t ) = exp i μ t · ψ ( x ) .
In the process of finding the functions (9), we come to ψ ( x ) and μ which are the eigenfunctions and eigenvalues of the problem (7).
Consider solutions of (8) concentrated in a vicinity of the one-dimensional manifold Λ t 1 (the curve) in the 2 n -dimensional phase space M 2 n = R n × R n , p R n , x R n :
Λ t 1 = { z = ( p , x ) = Z ( t , s , ) = P ( t , s , ) , X ( t , s , ) | s [ 0 , ω T ] , t [ 0 , T ] } , T , ω > 0 .
Here, s is the curve Λ t 1 parameter and:
Z ( t , s , ) = Z ( t , s ) + Z ( 1 ) ( t , s ) + O ( 3 / 2 ) .
For the functions Ψ semiclassically concentrated on the curve Λ t 1 , the dynamics of the curve in the limits 0 is determined by the ( 1 , 1 ) -type Hamilton–Ehrenfest system (the term ( 1 , 1 ) implies that the parameter s is a one-dimensional variable and the vector Z ( s , t ) corresponds to the first order moments of the function Ψ ( x , t ) , see [25] for details):
X ˙ ( t , s ) = V p Z ( t , s ) , P ˙ ( t , s ) = V x Z ( t , s ) ϰ | | Ψ | | 2 ω T 0 ω T W x X ( t , s ) , X ( t , r ) d r .
Here, | | Ψ | | is the L 2 -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 | | Ψ | | = 1 .
The class J τ of functions semiclassically concentrated on curves Λ 1 is given by
J τ = Φ : Φ ( x , t , ) = χ ( x , t , s , ) | s = ω τ ( x , t ) , χ ( x , t , s , ) P t ( s ) ,
where P t is a class of trajectory-concentrated functions [47]:
P t ( s ) = χ : χ ( x , t , s , ) = φ Δ x , t , s , exp i S ( t , s , ) + P ( t , s ) , Δ x .
In this definition, the function φ ( ξ , t , s , ) belongs to a Schwarz space S in the variable ξ R n , smoothly depends on t and s, and regularly depends on as 0 ; the function S ( t , s , ) together with functions Z ( t , s ) = ( P ( t , s ) , X ( t , s ) ) define the class and are to be determined; the family of hypersurfaces s = ω τ ( x , t ) are given by relation [23,24]
X s ( t , s ) , Δ x = 0 , Δ x = x X t , s .
Hereinafter, we assume that
rank X s ( t , s ) = 1 .
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
Ψ ( x , t ) = χ ( x , t , s ) | s = ω τ ( x , t )
to the non-stationary GPE (8) in the class of functions J τ in the semiclassical approximation with the accuracy of O ( 3 / 2 ) can be reduced to the solution of the associated linear GPE:
i t + H ( t , s , ) [ g ( t , s ) ] + H z ( t , s ) , Δ z ^ + 1 2 Δ z ^ , H z z ( t , s ) Δ z ^ χ ( x , t , s ) = 0 ,
with the additional algebraic constraint
a ^ 0 ( t , s ) χ ( x , t , s ) = 0 .
Here, we denoted
H t , s , [ g ( t , s ) ] = V Z ( t , s ) + ϰ ω T 0 ω T d r { W X ( t , s ) , X ( t , r ) + + W y X ( t , s ) , X ( t , r ) , X ( 1 ) ( t , r ) + 1 2 Sp W y y X ( t , s ) , X ( t , r ) · σ x x ( t , r ) } , H z ( t , s ) = V z Z ( t , s ) + ϰ ω T 0 ω T d r W z X ( t , s ) , X ( t , r ) , W z = 0 , W x , H z z ( t , s ) = V z z Z ( t , s ) + ϰ ω T 0 ω T d r W z z X ( t , s ) , X ( t , r ) , W z z = 0 0 0 W x x , a ^ 0 ( t , s ) = P s ( t , s ) , Δ x X s ( t , s ) , Δ p ^ , Δ z ^ = ( Δ p ^ , Δ x ) = z ^ Z ( t , s ) ,
and functions g ( t , s ) = ( Z ( 1 ) ( t , s , ) , Δ 2 ( t , s ) ) , where Δ 2 = σ p p σ p x σ x p σ x x , are the solutions of the second order Hamilton–Ehrenfest system with integration constants determined from the initial condition for the function χ ( x , t , s ) . In general case, these functions are given by:
Z ( 1 ) ( t , s ) = 1 Δ z ^ , Δ 2 k m ( t , ) = 1 2 Δ z ^ k Δ z ^ m + Δ z ^ m Δ z ^ k , k , m = 1 , 2 n ¯ .
Here, A ^ is the mean value of the operator with a Weyl symbol A ( z , t ) in the class J τ :
A ^ ( t , s ) = 1 ω T D χ * ( x , t , s ) A ( z ^ , t ) χ ( x , t , s ) d ξ ,
where ξ is ( n 1 ) -dimensional vector of variables that complement the variable s = ω τ ( x , t ) to form the coordinate system in R n , D is the domain of the variables ξ , and the operator z ^ acts subject to the condition s = const .
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:
a ˙ ( t , s ) = J H z z ( t , s ) a ( t , s ) ,
where a ( t , s ) C 2 n , J = 0 I n × n I n × n 0 is unit symplectic matrix, ( I n × n is the unit n × n -matrix). The function a 0 ( t , s ) = Z ˙ ( t , s ) is the particular solutions of (23). Let ( n 1 ) linear independent solutions of (23) be known and satisfy the conditions
{ a 0 , a k } = 0 , { a 0 , a k * } = 0 , { a l , a k } = 0 , { a l * , a k } = 2 i δ l k , l , k = 1 , ( n 1 ) ¯ .
Here, δ l k is the Kronecker symbol, { · , · } is the skew-orthogonal product.
Denote a k ( t , s ) = W k ( t , s ) , Z k ( t , s ) . Then, the particular solutions of (18) reads [47]
χ 0 ( x , t , s ) [ g ( t , s ) ] = N det C ( t , s ) exp { i [ 0 t ( P ( τ , s ) , X ˙ ( τ , s ) H ( τ , s , ) [ g ( t , s ) ] ) d τ + P ( t , s ) , Δ x + 1 2 Δ x , Q ( t , s ) Δ x ] } ,
where N is the normalization coefficient, and
Q ( t , s ) = B ( t , s ) C 1 ( t , s ) , B ( t , s ) = P s ( t , s ) , W 1 ( t , s ) , . . . , W n 1 ( t , s ) C n × n , C ( t , s ) = X s ( t , s ) , Z 1 ( t , s ) , . . . , Z n 1 ( t , s ) C n × n .
The conditions (16) and (24) ensures the non-degeneracy of matrix C ( t , s ) . In order to make sure that function (40) satisfies the auxiliary constraint (19) it is enough to check the validity of the identity
P s ( t , s ) Q ( t , s ) X s ( t , s ) = P s ( t , s ) B ( t , s ) C 1 ( t , s ) X s ( t , s ) = 0 ,
which results from the property
C 1 ( t , s ) X s ( t , s ) = ( 1 , 0 , . . . , 0 ) ,
following from the definition of the matrix C ( t , s ) .
Let us introduce first order operators of the form
a ^ k ( t , s ) = 1 2 a k ( t , s ) , J Δ z ^ , a ^ k + ( t , s ) = 1 2 a k * ( t , s ) , J Δ z ^ , k = 1 , n 1 ¯ .
Due to (24), the following commutation relations hold:
[ a ^ 0 , a ^ k ] = [ a ^ 0 , a ^ k + ] = [ a ^ k , a ^ m ] = [ a ^ k + , a ^ m + ] = 0 , a ^ k , a ^ m + ] = δ k m , k , m = 1 , n 1 ¯ .
These operators are the symmetry operators of (18) that preserve the property (19) of the solutions. Thus, the operators a ^ k + , k = 1 , n 1 ¯ are creation operators of (18), and the operators a ^ k , k = 0 , n 1 ¯ are annihilation operators. The function | 0 , t (40) is the “vacuum” solution. Let us construct the “occupation number representation” for (18) by determination of the following countable set of solutions:
χ ν ( x , t , s ) = k = 1 n 1 1 ν k ! a ^ k + ( t , s ) ν k χ 0 ( x , t , s ) .
Here, ν Z + n is a ( n 1 ) -dimensional multiindex, ν = ( ν 1 , ν 2 , . . . , ν n 1 ) , ν k = 0 , ¯ , k = 1 , n 1 ¯ :
ν ! = k = 1 n 1 ν k ! , a ν = k = 1 n 1 a k ν k .
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 Z ( t , s ) :
Z ( t , s ) = Z ( t + s / ω ) , Z ( t + T ) = Z ( t ) ,
Here, a parameter ω is the ratio between the periods by s and by t. The curve Λ t 1 determined by such the functions is mapped into itself with change of t, i.e., Λ t 1 = Λ 1 does not depend on t. By the change of variables ( t + s / ω ) t and ( t + r / ω ) r , the system (12) can be written in the following form:
X ˙ ( t ) = V p Z ( t ) , P ˙ ( t ) = V x Z ( t ) ϰ T 0 T W x X ( t ) , X ( r ) d r .
In this case, the equation of hypersurfaces s = ω τ ( x , t ) has the form s = ω τ ( x ) + t or
t + s / ω = τ ( x ) .
Let the solutions of the second order Hamilton–Ehrenfest satisfy the relations g ( t , s ) = g ( t + s / ω ) , g ( t + T ) = g ( t ) , i.e.
Z ( 1 ) ( t , s ) = ( P ( 1 ) ( t , s ) , X ( 1 ) ( t , s ) ) = Z ( 1 ) ( t + s / ω ) , Z ( 1 ) ( t + T ) = Z ( 1 ) ( t ) , Δ 2 ( t , s ) = Δ 2 ( t + s / ω ) , Δ 2 ( t + T ) = Δ 2 ( t ) .
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 ( 1 , 1 ) -type Hamilton–Ehrenfest system satisfying (33), the variational system (23) reads
a ˙ ( t ) = J H z z ( t ) a ( t ) ,
where we put
a ( t , s ) = a ( t + s / ω ) ,
and made the change of the variable ( t + s / ω ) t . Since H z z ( t ) 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 a 0 ( t ) = Z ˙ ( t ) (the T-periodic solution). Let there be ( n 1 ) more linear independent solutions of the Floquet problem for (37) that satisfy the following Floquet conditions:
a k ( t + T ) = exp [ i Ω k T ] a k ( t ) , Im Ω k = 0 , k = 1 , ( n 1 ) ¯
and the skew-orthogonality condition (24). Since only linearly stable Λ 1 are of interest for us, we consider the solutions of (37) with the real Floquet exponent Ω k .
If conditions (33), (36) and (38) are met, the particular solutions (25), (31) of (18) read χ ν ( x , t , s ) [ g ν ( t + s / ω ) ] = exp i μ ν t χ ν ( x , t + s / ω , μ ν ) [ g ν ( t + s / ω ) ] , where
χ 0 ( x , t , μ 0 ) [ g 0 ( t ) ] = N det C ( t ) exp { i [ 0 t ( P ( r ) , X ˙ ( r ) H ( r , ) [ g 0 ( t ) ] + μ 0 ( ) ) d r + P ( t ) , Δ x + 1 2 Δ x , Q ( t ) Δ x ] } ,
χ ν ( x , t , μ ν ) [ g ν ( t ) ] = k = 1 n 1 1 ν k ! a ^ k + ( t ) ν k χ 0 ( x , t , μ ν ) [ g ν ( t ) ] .
Then, considering (9) and (35), the relation (17) is as follows:
ψ ν ( x ) = χ ν ( x , t , μ ν ) [ g ν ( t ) ] | t = τ ( x ) .
The transition (42) from the extended ( n + 1 ) -dimensional space to the n-dimensional configuration space can give the multi-valued function ψ ( x ) . However, the function ψ ( x ) is the single-valued one if the following periodicity condition holds:
χ ν ( x , t + T , μ ν ) [ g ν ( t ) ] = χ ν ( x , t , μ ν ) [ g ν ( t ) ] .
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]
0 T P ( t ) , X ˙ ( t ) H ( t , ) [ g ν ( t ) ] + μ ν ( ) d t = 2 π l + T k = 0 n 1 Ω k ν k + 1 2 , l Z .
Let us denote H ( t , ) [ g ( t ) ] = H ( 0 ) ( Z ( t ) ) + H ( 1 ) ( t ) [ g ( t ) ] , where
H ( 0 ) Z ( t ) = V Z ( t ) + ϰ T 0 T W X ( t ) , X ( r ) d r , H ( 1 ) ( t ) [ g ( t ) ] = ϰ T 0 T d r ( W y X ( t ) , X ( r ) , X ( 1 ) ( r ) + + 1 2 Sp W y y X ( t ) , X ( r ) · σ x x ( r ) ) ,
and assume
μ ( ) = μ ( 0 ) + μ ( 1 ) ( ) ,
where μ ( 1 ) ( ) regularly depends on .
We consider the periodic trajectories of the ( 1 , 1 ) -type Hamilton–Ehrenfest system with the period T = T ( μ ( 0 ) ) determined by the following integral of motion of (34):
H ( 0 ) Z t , T ( μ ( 0 ) ) = μ ( 0 ) .
Then, the quantization condition (44) is satisfied if
0 T ( μ ( 0 ) ) P t , T ( μ ( 0 ) ) , X ˙ t , T ( μ ( 0 ) ) d t = 2 π l , μ ν ( 1 ) = 1 T 0 T H ( 1 ) ( t ) [ g ν ( t ) ] d t + k = 0 n 1 Ω k ν k + 1 2 .
The first equation in (48) implicitly determines μ ( 0 ) . The energy per particle for constructed semiclassical spectral series is given by
E ν = E ( 0 ) + E ν ( 1 ) + O ( 3 / 2 ) , E ( 0 ) = 1 T 0 T d t V Z ( t ) + 1 2 ϰ T 0 T W X ( t ) , X ( r ) d r , E ν ( 1 ) = 1 T 0 T d t [ V z Z ( t ) , Z 1 ( t ) + 1 Sp V z z Z ( t ) · Δ 2 ( t ) + + 1 2 ϰ T 0 T d r ( W x X ( t ) , X ( r ) , Z ( 1 ) ( t ) + 1 Sp W x x X ( t ) , X ( r ) · σ x x ( t ) + + W y X ( t ) , X ( r ) , X ( 1 ) ( r ) + 1 2 Sp W y y X ( t ) , X ( r ) · σ x x ( r ) ) ] .
Note that the quantization condition (48) describes the spectrum well in the linear theory when l 1 / .

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:
H ( z ^ ) [ ψ ] μ ψ ( x ) = V ( z ^ ) + ϰ R 2 W ( x , y ) | ψ ( y ) | 2 d y μ ψ ( x ) = 0 , V ( z ^ ) = 1 2 p ^ 2 + x 2 , W ( x , y ) = exp ( x y ) 2 2 γ 2 .
Here, the units of length and energy are given by such the scales that the term V ( z ^ ) corresponds to 2 2 m Δ + 1 2 m ω h a r m 2 x 2 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 ϰ > 0 , the nonlinear term in (50) corresponds to the repulsive interaction. For the unit normalization of the wave function ψ ( x ) , the coefficient ϰ is proportional to the number of condensate particles. Since the equation has the axial symmetry, the ( 1 , 1 ) -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
X ( t + s / ω ) = R cos ( ω t + s ) , sin ( ω t + s ) ,
where the parameter s is the angle in polar coordinates. The ( 1 , 1 ) -type Hamilton–Ehrenfest system (34) for the GPE (50) reads
X ˙ ( t ) = P ( t ) , P ˙ ( t ) = X ( t ) ϰ T 0 T W x X ( t ) , X ( r ) d r , W x ( x , y ) = y x γ 2 exp ( x y ) 2 2 γ 2 .
From the first equation of the system (52) and (51), we have
P ( t ) = X ˙ ( t ) = ω R · ( sin ω t , cos ω t ) .
For the function (51), the integrals in the Equations (52) are as follows:
1 T 0 T W x 1 ( X ( t ) , X ( r ) ) d r = Γ R e Γ · I 0 ( Γ ) I 1 ( Γ ) cos ω t ,
1 T 0 T W x 2 ( X ( t ) , X ( r ) ) d r = Γ R e Γ · I 0 ( Γ ) I 1 ( Γ ) sin ω t ,
where I 0 , I 1 are modified Bessel functions of the first kind, Γ = R 2 γ 2 . Here, we have taken into consideration that ω = 2 π T and have used the property
I n ( z ) = 1 2 π π π e z cos φ cos n φ d φ , n = 0 , ¯ .
As a result, the ( 1 , 1 ) -type Hamilton–Ehrenfest system yields the following equation for R:
ω 2 = 1 ϰ 1 γ 2 exp [ Γ ] · I 0 ( Γ ) I 1 ( Γ ) .
Consider the approximation R γ . It meets the strongly local interaction. The following asymptotic estimate holds for the modified Bessel functions:
I 0 ( x ) = 1 2 π x e x ( 1 + 1 8 x + O ( 1 / x 2 ) ) , I 1 ( x ) = 1 2 π x e x ( 1 3 8 x + O ( 1 / x 2 ) ) .
Then, we have
R 3 = ϰ γ 2 2 π ( 1 ω 2 ) ( 1 + O ( γ 4 / R 4 ) ) .
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 ω = 1 for ϰ = 0 , i.e., the period T depends on the chemical potential μ ( 0 ) only in the presence of the nonlinearity in (50). The first quantization condition in (48) for the solutions (53) reads
ω l R l 2 = l ( ) , μ l ( 0 ) = 1 2 R l 2 ( 1 + ω l 2 ) + ϰ exp R l 2 γ 2 I 0 R l 2 γ 2 , l Z .
The quantized values of μ ( 0 ) for = 0.1 and for various values of ϰ , γ are shown in Figure 1. Note that the eigenvalue μ ( 0 ) linearly depends on l in the linear case ( ϰ = 0 ). 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 μ ( 0 ) for the given range of quantum numbers l. For the large values of ϰ , there are not closed trajectories z = Z ( t ) .
The energy per particle the zeroth order approximation with respect to is as follows:
E l ( 0 ) = 1 2 R l 2 ( 1 + ω l 2 ) + ϰ 2 exp R l 2 γ 2 I 0 R l 2 γ 2 .
The Equation (15) for the function t + s / ω = τ ( x ) for the curve (51) reads
x sin ( ω t + s ) + y cos ( ω t + s ) = 0 ,
and it has the solution
t + s / ω = τ ( x ) = 1 ω Arg ( x 1 + i x 2 ) = φ ω .
The semiclassical solutions of (50) for ν = 0 and ν = 1 are obtained in Appendices Appendix C and Appendix D, respectively.
The solution for ν = 0 (the vacuum solution) reads
ψ 0 ( x ) = N ( 2 θ ) 1 4 ω R exp θ ( ρ R ) 2 + i · l ( ) φ , x = ( ρ cos φ , ρ sin φ ) ,
where θ is the positive constant defined in Appendix C.
The solution for ν = 1 is given by
ψ 1 ( x ) = N 1 2 ( θ 3 R 2 + 1 ) θ 5 2 R 2 ( 2 θ ) 1 4 ω R ( ρ R ) exp θ ( ρ R ) + i · l ( ) φ .
In Figure 2 and Figure 3, the plots of | ψ ν ( x ) | 2 and Arg [ ψ ν ( x ) ] are shown for ν = 0 (56) and ν = 1 (57) when γ = 1 , ϰ = 1 , = 0.1 , and l = 10 (it corresponds to ω 0.883 , R 1.06 , θ 0.528 ).
The phase incursion of 20 π 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 + 2 ) (n variables x , t, and s) to ( n + 1 ) (n variables x , and ( t + s / ω ) ). 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 2 π (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 ϰ = 0 (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 σ x x and the first order correction with respect to for the first order moment X ( 1 ) 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 Δ 2 ( t ) accurate to O ( 3 / 2 ) is determined by the leading term of asymptotics since it corresponds to an operator with the estimate O ( ) . The formula for Δ 2 ( t ) reads
Δ 2 l k ( t ) [ χ ν ] = 1 2 Δ z ^ l Δ z ^ k + Δ z ^ k Δ z ^ l = = 1 T D χ ν * ( x , t ) 1 2 ( Δ z ^ l Δ z ^ k + Δ z ^ k Δ z ^ l ) χ ν ( x , t ) d ξ + O ( 3 / 2 ) ,
where l , k = 1 , 2 n ¯ , ξ is ( n 1 ) -dimensional vector of variables that complement the variable t = τ ( x ) to form the coordinate system in R n , and D is the domain of the variables ξ . Here, the operator Δ z ^ acts subject to the condition t = const . Note that the function H ( t , ) [ g ν ( t ) ] included in the relation (40) affects only the phase of χ ν ( x , t ) and commutes with Δ z ^ . Hence, the matrix σ x x ( t ) does not affect the result of the formula (A1), i.e., this formula is not recursive. The function given by (A1) satisfy the equation
Δ 2 ˙ ( t ) = J H z z ( t ) Δ 2 ( t ) Δ 2 ( t ) H z z ( t ) J .
In general case, the solution of (A2) reads [25]
Δ 2 ( t ) = A ( t ) D 0 A ( t ) ,
where A ( t ) is the solution of the matrix variational system
A ( t ) = J H z z ( t ) A ( t ) , A ( 0 ) = I 2 n × 2 n ,
and the matrix D 0 = Δ 2 ( t ) | t = 0 is given by (A1). In view of (33) and (36), the system of equations for Z ( 1 ) ( t ) = P ( 1 ) ( t ) , X ( 1 ) ( t ) can be written as
d d t Z ( 1 ) ( t ) Π ( t ) = J z ( V z Z ( t ) , Z ( 1 ) ( t ) + + 1 2 Sp [ V z z Z ( t ) · Δ 2 ( t ) ] + ϰ T 0 T d r ( W y X ( t ) , X ( r ) , X ( 1 ) ( r ) + + 1 2 Sp W y y X ( t ) , X ( r ) · σ x x ( r ) + W x x X ( t ) , X ( r ) · σ x x ( t ) + + W x X ( t ) , X ( r ) , X ( 1 ) ( t ) ) ) + O ( 3 / 2 ) , Π ( t ) = π 0 ( t ) X ˙ ( t ) 2 X ˙ ( t ) , 0 ,
where the operator z = z = p ; x acts on the arguments of the function V ( z ) and W ( x , · ) , and the function π 0 ( t ) is given by
π 0 ( t ) [ χ ν ] = T Re i D χ ν * ( x , t ) t χ ν ( x , t ) d ξ .
It was shown in [25] that π 0 ˙ ( t , ) = O ( 3 / 2 ) . Hence, in the system (A5), we can put
π 0 ( t ) = const = π 0 ( t 0 ) = π 0 .
In view of (A7), the system (A5) is the system of linear integro-differential equations. Note that, for every ν , the function Δ 2 ( t ) [ χ ν ] and π 0 [ χ ν ] , which are included in (A5), are different. As opposed to Δ 2 ( t ) , the function Z ( 1 ) ( t ) 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 Z ( 1 ) ( t ) .

Appendix B. Auxiliary Theorems

Theorem A1.
The matrix Q ( t ) given by (26), (24) and (37) is T-periodic, i.e.
Q ( t + T ) = Q ( t ) .
Proof. 
In view of (26), (24) can be written as
C ( t + T ) = C ( t ) · Λ ( T ) , Λ ( T ) = diag 1 , exp { i Ω 1 T } , exp { i Ω 2 T } , . . . , exp { i Ω n 1 T } .
This relation is also valid for B ( t )
B ( t + T ) = B ( t ) · Λ ( T ) .
Then,
Q ( t + T ) = B ( t + T ) C 1 ( t + T ) = B ( t ) Λ ( T ) Λ 1 ( T ) C 1 ( t ) = B ( t ) C 1 ( t ) = Q ( t ) .
 □
Theorem A2.
For the matrix C ( t ) given by (26), (24) and (37), the following relation holds:
det C 1 ( t + T ) = exp i T 2 k = 0 n 1 Ω k det C 1 ( t )
Proof. 
The statement of this theorem directly follows from (A9). □

Appendix C. The Solution of (50) for ν = 0

The variational system for the function a ( t ) = W ( t ) , Z ( t ) in view of (53) takes the form
Z ˙ ( t ) = W ( t ) , W ˙ ( t ) = ω 2 Z ( t ) ϰ Γ R 2 e Γ · D ( t ) · Z ( t ) , D ˜ ( t ) = 1 2 Γ I 0 ( Γ ) 2 I 1 ( Γ ) + I 2 ( Γ ) · I 2 × 2 + D ( t ) , D ( t ) = cos 2 ω t sin 2 ω t sin 2 ω t cos 2 ω t .
For (A13), we seek the Floquet solutions satisfying (24).
The system (A13) can be written as
W ( t ) = Z ˙ ( t ) , Z ¨ ( t ) = g Z ( t ) + κ · D ( t ) · Z ( t ) .
where g = κ ω 2 , κ = ϰ Γ 2 γ 2 e Γ I 0 ( Γ ) 2 I 1 ( Γ ) + I 2 ( Γ ) are scalar coefficients. Let us make the following change of variables corresponding to the rotation of the vector Z ( t ) by the angle ω t :
Z ( t ) = M ( ω t ) · U ( t ) , M ( t ) = cos t sin t sin t cos t .
Then, the variational system reads
W ( t ) = ω · M ω t + π 2 · U ( t ) + M ( ω t ) · U ˙ ( t ) , D ^ 2 ( g + κ + ω 2 ) 2 ω D ^ 2 ω D ^ D ^ 2 ( g κ + ω 2 ) U ( t ) = 0 ,
where D ^ = t . We seek U ( t ) in the form
U ( t ) = 2 ω D ^ D ^ 2 ( g + κ + ω 2 ) u ( t ) .
In view of g = κ ω 2 , this relation can be written as
U ( t ) = 2 ω D ^ D ^ 2 2 κ u ( t ) .
The characteristic equation
( D 2 ( g + κ + ω 2 ) ) ( D 2 ( g κ + ω 2 ) ) + 4 ω 2 D 2 = = D 4 + 2 ( ω 2 g ) D 2 + g 2 + 2 g ω 2 + ω 4 κ 2 = 0
has the roots
D = ± g ω 2 ± κ 2 4 ω 2 g .
Substitution of g = κ ω 2 into this relation yields
D 1 , 2 = 0 , D 3 , 4 = ± 2 κ 4 ω 2 .
It is clear that we have only aperiodic solutions for κ 2 ω 2 , which do not satisfy the Floquet condition (24). Hence, we consider the case when κ < 2 ω 2 . Denote θ = ω 2 κ 2 . Then, the function u ( t ) takes the form
u ( t ) = C 1 + C 2 t + C 3 e 2 i θ t + C 4 e 2 i θ t .
Substitution of this relation into the formula for U ( t ) yields the following fundamental system of solutions:
U 1 ( t ) = 0 1 , U 2 ( t ) = ω κ t , U 3 ( t ) = i θ e 2 i θ t ω e 2 i θ t , U 4 ( t ) = i θ e 2 i θ t ω e 2 i θ t
For Z ( t ) , we have
Z 1 ( t ) = sin ω t cos ω t , Z 2 ( t ) = ω cos ω t + κ t sin ω t ω sin ω t κ t cos ω t ,
Z 3 ( t ) = e 2 i θ t i θ cos ω t ω sin ω t i θ sin ω t + ω cos ω t , Z 4 ( t ) = e 2 i θ t i θ cos ω t + ω sin ω t i θ sin ω t ω cos ω t .
The first solution corresponds to the already known solution a 0 ( t ) = Z ˙ ( t ) . The second solutions does not satisfy the Floquet condition. The third and forth solutions are related by Z 3 ( t ) * = Z 4 ( t ) . For W 3 ( t ) , W 4 ( t ) , we have:
W 4 ( t ) = Z ˙ 4 ( t ) = e 2 i θ t ( ω 2 2 θ 2 ) cos ω t + i θ ω sin ω t ( ω 2 2 θ 2 ) sin ω t i θ ω cos ω t , W 3 ( t ) = W 4 ( t ) * .
It is easy to see that
Z ˙ ( t ) , a ˜ 4 ( t ) = 0 , a ˜ 4 ( t ) = W 4 ( t ) , Z 4 ( t ) ,
and
a ˜ 4 * ( t ) , a ˜ 4 ( t ) = 4 i θ 3 .
Thus, taking into account the normalization condition { a 1 ( t ) , a 1 * ( t ) } = 2 i , we have the solution
a 1 ( t ) = 1 2 θ 3 e 2 i θ t ( ω 2 2 θ 2 ) cos ω t + i θ ω sin ω t ( ω 2 2 θ 2 ) sin ω t i θ ω cos ω t i θ cos ω t + ω sin ω t i θ sin ω t ω cos ω t
with the Floquet exponent Ω = 2 θ = 2 ω 2 κ / 2 . Substitution of a 0 ( t ) and a 1 ( t ) into (26) yields
det C ( t ) = i e 2 i θ t ω R 2 θ , Re Q ( t ) = ω sin 2 ω t cos 2 ω t cos 2 ω t sin 2 ω t , Im Q ( t ) = θ 2 cos 2 ω t sin 2 ω t sin 2 ω t 2 sin 2 ω t .
and
| χ 0 ( x , t ) | 2 = N 2 2 θ ω R exp 1 Δ x , Im Q ( t ) Δ x .
The coefficient N is obtained from the unit normalization condition:
R 2 | χ 0 ( x , τ ( x ) | 2 d x = 1 , τ ( x ) = 1 ω Arg ( x + i y ) .
N 2 = 2 π π ω 2 + O ( ) .
Then, we have
| χ 0 ( x , t ) | 2 = ω θ R π 2 π exp 1 Δ x , Im Q ( t ) Δ x .
Note that V z z z ( z , t ) = 0 in our example. Hence, it is not necessary to obtain the whole matrix Δ 2 ( t ) for the system (A5) but it is enough to obtain only σ x x ( t ) , which is given by
σ x x ( t ) = T 0 r ( r R ) 2 I 2 × 2 + D ( t ) ω θ 2 R π 2 π · · exp θ ( r R ) 2 R 2 X ( t ) , I 2 × 2 + D ( t ) X ( t ) d r = = T I 2 × 2 + D ( t ) 2 ω θ R π 2 π 0 r ( r R ) 2 exp 2 θ ( r R ) 2 d r = = 4 2 θ ω R I 2 × 2 + D ( t ) .
Substituting (40) into (A6), the formula for π 0 is as follows:
π 0 = H ( 1 ) ( t ) Re Q ( t ) X ˙ ( t ) , Δ x + 1 2 Δ x , Re Q ˙ ( t ) , Δ X .
In view of (27), we have
π 0 = H ( 1 ) ( t ) P ˙ ( t ) , X ( 1 ) ( t ) + 1 2 Sp Re Q ˙ ( t ) σ x x ( t ) .
Substitution of (45) into this formula yields
π 0 = ϰ T 0 T d r W y X ( t ) , X ( r ) , X ( 1 ) ( r ) + 1 2 Sp W y y X ( t ) , X ( r ) · σ x x ( r ) P ˙ ( t ) , X ( 1 ) ( t ) + 1 2 Sp Re Q ˙ ( t ) σ x x ( t ) .
The identity (A14) is the integral of motion of the system (A5) for ν = 0 . The system (A5) for the coordinate part of Z ( 1 ) ( t ) yields
X ˙ ( 1 ) ( t ) = P ( 1 ) ( t ) ,
For the pulse part Z ( 1 ) ( t ) , we have
P ˙ ( 1 ) ( t ) + π 0 R 2 X ( t ) = X ( 1 ) ( t ) + + ϰ T γ 2 exp [ Γ ] 0 T d r exp [ Γ cos ( ω ( r t ) ) ] I n × n 1 γ 2 X ( t ) X ( t ) X ( 1 ) ( r ) + + ϰ 8 2 θ ω γ 2 R 3 exp [ Γ ] ( 2 Γ + 4 Γ 2 ) I 0 ( Γ ) ( 1 + 4 Γ 2 ) I 1 ( Γ ) X ( t ) ϰ γ 2 exp [ Γ ] Γ I 0 ( Γ ) ( 1 + Γ ) I 1 ( Γ ) · I 2 × 2 + D ( t ) · X ( 1 ) ( t ) + + ϰ 2 γ 2 exp [ Γ ] I 0 ( Γ ) I 1 ( Γ ) X ( 1 ) ( t ) .
The integral of the system (A14) can be written as
π 0 = ϰ T γ 2 exp [ Γ ] 0 T d r exp [ Γ cos ( ω ( r t ) ) ] X ( r ) X ( t ) , X ( 1 ) ( r ) + + ω 2 X ( t ) , X ( 1 ) ( t ) + ω 8 2 θ R ϰ 16 2 θ ω R γ 2 exp [ Γ ] ( 1 + 2 Γ ) I 1 ( Γ ) ( 2 Γ 1 ) I 0 ( Γ ) .
We seek the solution of (A15) in the form
X ( 1 ) ( t ) = M ( ω t ) U ( 1 ) ( t ) , U ( 1 ) = U 1 ( 1 ) U 2 ( 1 ) .
Then, for (A16), we get
π 0 = ϰ Γ T γ exp [ Γ ] 0 T d r exp [ Γ cos ( ω ( r t ) ) ] Y ( t , r ) , U ( 1 ) ( r ) + + ω 2 R U 1 ( 1 ) ( t ) + ω 8 2 θ R ϰ 16 2 θ ω R γ 2 exp [ Γ ] ( 1 + 2 Γ ) I 1 ( Γ ) ( 2 Γ 1 ) I 0 ( Γ ) , Y ( t , r ) = cos ( ω ( r t ) ) 1 sin ( ω ( r t ) ) .
The constant vector U ( 1 ) ( t ) = U ( 1 ) satisfies (A17). Let us show that it also satisfies the Equation (A15). The integral (A17) with U ( 1 ) ( t ) = U ( 1 ) is as follows:
π 0 = ϰ Γ γ exp [ Γ ] ( I 1 ( Γ ) I 0 ( Γ ) ) U 1 ( 1 ) + ω 2 R U 1 ( 1 ) + + ω 8 2 θ R ϰ 16 2 θ ω R γ 2 exp [ Γ ] ( 1 + 2 Γ ) I 1 ( Γ ) ( 2 Γ 1 ) I 0 ( Γ ) .
Substituting (A18) into (A15), taking into consideration P ˙ ( 1 ) ( t ) = X ¨ ( 1 ) ( t ) = M ( ω t ) U ( 1 ) and multiplying the equation by the matrix M ( ω t ) on the left, we get
π 0 R 1 0 = = ϰ 2 γ 2 exp [ Γ ] I 2 ( Γ ) I 0 ( Γ ) 4 I 1 ( Γ ) 0 0 I 0 ( Γ ) + I 2 ( Γ ) U ( 1 ) + + ϰ γ 3 Γ 1 8 2 θ ω R exp [ Γ ] ( 2 Γ + 4 Γ 2 ) I 0 ( Γ ) ( 1 + 4 Γ 2 ) I 1 ( Γ ) 1 0 ϰ γ 2 exp [ Γ ] ( Γ I 0 ( Γ ) ( 1 + Γ ) I 1 ( Γ ) ) 1 0 0 0 U ( 1 ) + + ϰ 2 γ 2 exp [ Γ ] ( I 0 ( Γ ) I 1 ( Γ ) ) U ( 1 ) .
The Equation (A19) yields U 2 ( 1 ) = 0 and U 1 ( 1 ) given by
U 1 ( 1 ) = ω 8 2 θ R 2 + ϰ γ 2 exp [ Γ ] 16 2 θ ω R 2 ( 1 6 Γ + 8 Γ 2 ) I 0 ( Γ ) + ( 1 + 2 Γ 8 Γ 2 ) I 1 ( Γ ) ω 2 + ϰ 2 γ 2 exp [ Γ ] ( 2 + 2 Γ ) I 0 ( Γ ) + ( 5 2 Γ ) I 1 ( Γ ) I 2 ( Γ ) .
Finally, we obtain
X ( 1 ) ( t ) = U ( 1 ) R · X ( t ) , P ( 1 ) ( t ) = U ( 1 ) R · P ( t ) .
For H ( 1 ) ( t ) , we have
H ( 1 ) ( t ) = ϰ exp [ Γ ] γ 2 [ R U ( 1 ) I 0 ( Γ ) I 1 ( Γ ) + + 1 4 2 θ ω R ( 1 2 Γ ) I 0 ( Γ ) + ( 1 + 2 Γ ) I 1 ( Γ ) ] .
Since H ( 1 ) ( t ) = H ( 1 ) 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 χ 0 is as follows:
χ 0 ( x , t ) = N 0 det C ( t ) exp i ω 2 R 2 t + θ t + P ( t ) , Δ x + 1 2 Δ x , Q ( t ) Δ x ,
Then, ψ 0 ( x ) = χ 0 ( x , t ) | t = φ / ω yields
ψ 0 ( x ) = N ( 2 θ ) 1 4 ω R exp θ ( ρ R ) 2 + i · l ( ) φ , x = ( ρ cos φ , ρ sin φ ) .

Appendix D. The Solution of (50) for ν = 1

For ν 0 , the function χ ν ( x , t ) of (41) can be written as
χ ν ( x , t ) = N 1 · H μ ( ν ) ξ ( t ) , t det C ( t ) × × exp i ω 2 R 2 t + 2 θ ( 1 2 + ν ) t + P ( t ) , Δ x + 1 2 Δ x , Q ( t ) Δ x , ξ ( t ) = i 2 C * ( t ) 1 Δ x ,
where μ ( ν ) = ( 0 , ν ) is n-dimensional multiindex, H μ ( ν ) are multidimensional Hermitian polynomials that are generated by the matrix W ( t ) = C + ( t ) C 1 ( t ) and are given by
H μ ( ν ) ( ξ , t ) = ( 1 ) ν ν ! ξ W ( t ) ξ μ ( ν ) · 1 .
Here,
W ( t ) = 1 0 2 R θ 3 e 2 i θ t e 4 i θ t .
For ν = 1 , we have
H ( 0 , 1 ) ( ξ , t ) = W ( t ) ξ ( 0 , 1 ) = 2 e 2 i θ t θ 5 / 2 R 3 ω ω ( θ 3 R 2 + 1 ) X ( t ) i θ ω P ( t ) , Δ x .
Substitution of this relation into (A1) yields
σ x x ( t ) [ χ 1 ] = 3 2 2 θ ω R cos 2 ω t cos ω t sin ω t cos ω t sin ω t sin 2 ω t .
The integral π 0 [ χ ν ] is given by
π 0 [ χ 1 ] = ϰ T 0 T d r ( W y X ( t ) , X ( r ) , X ( 1 ) ( r ) [ χ 1 ] + + 1 2 Sp W y y X ( t ) , X ( r ) · σ x x ( r ) [ χ 1 ] ) P ˙ ( t ) , X ( 1 ) ( t ) [ χ 1 ] + + 1 2 Sp Re Q ˙ ( t ) σ x x ( t ) [ χ 1 ] θ 2 θ 3 R 3 + R + 1 θ 3 R 3 + R
or
π 0 [ χ 1 ] = ϰ T γ 2 exp [ Γ ] 0 T d r exp [ Γ cos ( ω ( r t ) ) ] X ( r ) X ( t ) , X ( 1 ) ( r ) [ χ 1 ] θ 2 θ 3 R 3 + R + 1 θ 3 R 3 + R + ω 2 X ( t ) , X ( 1 ) ( t ) [ χ 1 ] + 3 ω 8 2 θ R 3 ϰ 16 2 θ ω R γ 2 exp [ Γ ] ( 1 + 2 Γ ) I 1 ( Γ ) ( 2 Γ 1 ) I 0 ( Γ ) .
Then, X ( 1 ) ( t ) [ χ 1 ] = U 1 ( 1 ) [ χ 1 ] R · X ( t ) where
U 1 ( 1 ) [ χ 1 ] = 3 U 1 ( 1 ) [ χ 0 ] + + θ ( 2 θ 3 R 3 + R + 1 ) θ 3 R 4 + R 2 ω 2 + ϰ 2 γ 2 exp [ Γ ] ( 2 + 2 Γ ) I 0 ( Γ ) + ( 5 2 Γ ) I 1 ( Γ ) I 2 ( Γ ) .
The relation (A20) for ν = 1 can be rearranged in the form
H ( 1 ) ( t ) [ χ ν ] = ϰ exp [ Γ ] γ 2 [ R U ( 1 ) [ χ ν ] I 0 ( Γ ) I 1 ( Γ ) + + 3 4 2 θ ω R ( 1 2 Γ ) I 0 ( Γ ) + ( 1 + 2 Γ ) I 1 ( Γ ) ] .
Periodic solutions of the second order Hamilton–Ehrenfest system for other ν can be obtained in the same way. For ν = 1 , the function ψ 1 ( x ) = χ 1 ( x , t ) | t = φ / ω is as follows:
ψ 1 ( x ) = N 1 2 ( θ 3 R 2 + 1 ) θ 5 2 R 2 ( 2 θ ) 1 4 ω R ( ρ R ) exp θ ( ρ R ) + i · l ( ) φ , x = ( ρ cos φ , ρ sin φ ) .

References

  1. Pitaevskii, L.; Stringari, S. Bose–Einstein Condensation and Superfluidity; Oxford University Press: Oxford, UK, 2016. [Google Scholar]
  2. 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]
  3. 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]
  4. 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]
  5. Baranov, M.A. Theoretical progress in many-body physics with ultracold dipolar gases. Phys. Rep. 2008, 464, 71–111. [Google Scholar] [CrossRef]
  6. 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]
  7. Griesmaier, A.; Werner, J.; Hensler, S.; Stuhler, J.; Pfau, T. Bose–Einstein condensation of chromium. Phys. Rev. Lett. 2005, 94, 160401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. 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]
  9. 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]
  10. Pethick, C.J.; Smith, H. Bose–Einstein Condensation in Dilute Gases; Cambridge University Press: Cambridge, UK, 2008; pp. 1–569. [Google Scholar]
  11. 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]
  12. 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]
  13. 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]
  14. 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]
  15. 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]
  16. 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]
  17. 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]
  18. 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]
  19. Zezyulin, D.A.; Lebedev, M.E.; Alfimov, G.L.; Malomed, B.A. Symmetry breaking in competing single-well linear-nonlinear potentials. Phys. Rev. E 2018, 98, 042209. [Google Scholar] [CrossRef] [Green Version]
  20. 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]
  21. 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]
  22. Maslov, V. Operational Methods; Mir Publishers: Moscow, Russia, 1976. [Google Scholar]
  23. Maslov, V. The Complex WKB Method for Nonlinear Equations. I. Linear Theory; Birkhauser Verlag: Basel, Switzerland, 1994. [Google Scholar]
  24. 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]
  25. Shapovalov, A.V.; Kulagin, A.E.; Trifonov, A.Y. The Gross–Pitaevskii equation with a nonlocal interaction in a semiclassical approximation on a curve. Symmetry 2020, 12, 201. [Google Scholar] [CrossRef] [Green Version]
  26. Maslov, V. Complex Markov Chains and the Feynman Path Integral; Nauka: Moscow, Russia, 1976. [Google Scholar]
  27. Maslov, V. Equations of the self-consistent field. J. Sov. Math. 1978, 11, 123–195. [Google Scholar] [CrossRef]
  28. 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]
  29. 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]
  30. 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]
  31. 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. SIGMA 2007, 3, 5. [Google Scholar] [CrossRef] [Green Version]
  32. 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]
  33. 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]
  34. Pereskokov, A.V. Semiclassical asymptotics of the spectrum near the lower boundary of spectral clusters for a Hartree-type operator. Math. Notes 2017, 101, 1009–1022. [Google Scholar] [CrossRef]
  35. 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]
  36. Lugiato, L.A.; Lefever, R. Spatial dissipative structures in passive optical systems. Phys. Rev. Lett. 1987, 58, 2209–2211. [Google Scholar] [CrossRef] [PubMed]
  37. 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]
  38. 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]
  39. 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]
  40. Butov, L.V.; Lai, C.W.; Ivanov, A.L.; Gossard, A.C.; Chemla, D.S. Towards Bose–Einstein condensation of excitons in potential traps. Nature 2002, 417, 47–52. [Google Scholar] [CrossRef] [PubMed]
  41. Deng, H.; Haug, H.; Yamamoto, Y. Exciton-polariton Bose–Einstein condensation. Rev. Mod. Phys. 2010, 82, 1489–1537. [Google Scholar] [CrossRef] [Green Version]
  42. 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]
  43. Berman, O.L.; Kezerashvili, R.Y. Superfluidity of dipolar excitons in a transition metal dichalcogenide double layer. Phys. Rev. B 2017, 96, 094502. [Google Scholar] [CrossRef] [Green Version]
  44. Berman, O.L.; Gumbs, G.; Kezerashvili, R.Y. Bose–Einstein condensation and superfluidity of dipolar excitons in a phosphorene double layer. Phys. Rev. B 2017, 96, 014505. [Google Scholar] [CrossRef] [Green Version]
  45. 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]
  46. 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]
  47. 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]
  48. Fetter, A.L.; Svidzinsky, A.A. Vortices in a trapped dilute Bose–Einstein condensate. J. Phys. Condens. Matter 2001, 13, R135–R194. [Google Scholar] [CrossRef] [Green Version]
  49. Obukhov, V.V. Algebra of symmetry operators for Klein-Gordon-Fock equation. Symmetry 2021, 13, 727. [Google Scholar] [CrossRef]
  50. 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 μ ( 0 ) for various quantum numbers l and for ϰ = 1 , γ = 1 (black), ϰ = 1 2 , γ = 1 (red), ϰ = 1 2 , γ = 2 (blue).
Figure 1. The values of μ ( 0 ) for various quantum numbers l and for ϰ = 1 , γ = 1 (black), ϰ = 1 2 , γ = 1 (red), ϰ = 1 2 , γ = 2 (blue).
Symmetry 13 01289 g001
Figure 2. Density plots of | ψ ν ( x ) | 2 .
Figure 2. Density plots of | ψ ν ( x ) | 2 .
Symmetry 13 01289 g002
Figure 3. Density plots of Arg [ ψ ν ( x ) ] .
Figure 3. Density plots of Arg [ ψ ν ( x ) ] .
Symmetry 13 01289 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. Symmetry 2021, 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

Back to TopTop