Semiclassical approach to the nonlocal nonlinear Schr¨odinger equation with a non-Hermitian term

The nonlinear Sch¨odinger equation (NLSE) with a non-Hermitian term is the model for various phenomena in nonlinear open quantum systems. We deal with the Cauchy problem for the nonlocal generalization of multidimensional NLSE with a non-Hermitian term. Using the ideas of the Maslov method, we propose the method of constructing asymptotic solutions to this equation within the framework of semiclassically concentrated states. The semiclassical nonlinear evolution operator and symmetry operators for the leading term of asymptotics are derived. Our approach is based on the solutions of the auxiliary dynamical system that effectively linearize the problem under certain algebraic conditions. The formalism proposed is illustrated with the specific example of the NLSE with a non-Hermitian term that is the model of an atom laser. The analytical asymptotic solution to the Cauchy problem is obtained explicitly for this example.


I. INTRODUCTION
A great variety of nonlinear phenomena in matter-wave and optical media is modelled based on generalizations and modifications of the nonlinear Schrödinger equation (NLSE).
A significant share of theoretical research based on NLSE is associated with nonlinear optics [1], the physics of Bose-Einstein condensates (BEC) [2,3], the dynamics of quantum vortices [4,5], and other areas of nonlinear physics.In the BEC theory based on the mean-field approximation, the NLSE, termed the Gross-Pitaevskii equation (GPE) [2], is the base model equation.Mathematically, it is considered in multidimensional space-time with the variable coefficients responsible for the external fields of the traps confining the condensate in some area.In nonlinear optics, the NLSE describes optical solitons [6][7][8], which represent spatially localized perturbations of an electromagnetic field steadily propagating in a nonlinear medium.
These models usually describe conservative systems isolated from the environment, in which dissipative phenomena are not taken into account.However, in real conditions, quantum systems interact with the environment.In many-particle quantum systems, this interaction is dissipative, weakening the coherent effects, and, accordingly, blurring the manifestation of pronounced quantum properties.On the other hand, the combination of dissipative effects with quantum ones gives rise to new scenarios of behavior in nonlinear systems.This encourages the interest in the study of dissipative phenomena in nonlinear quantum systems, primarily BEC and nonlinear optical ones that has been the subject of detailed studies in many publications.Some examples below give an idea of these studies.
In [9] the use of BEC atoms in an atom laser is discussed.This problem is directly related to the interaction of the condensate with noncondensed trapped atoms which is a dissipative process.The theoretical description of this system is constructed in the mean field approximation and GPE with additional dissipative terms.After some simplifications, we obtain a closed description in terms of the complex Ginzburg-Landau equation (GLE) [10], which formally has the NLSE form with a non-Hermitian operator.In the review [11], in particular, the impact on the BEC parameters by Feshbach resonance, which leads to the generation of trains of solitons, was studied.This problem was considered in terms of a dissipative one-dimensional GPE with a time-dependent complex addition to the potential.
The paper [12] considers the dynamics of nonlinear waves in periodic complex PT-potentials, when the NLSE becomes invariant under the parity and time-reversal symmetry.A special kind of non-equilibrium stationary states was introduced and studied in [13] based on the GPE with a stochastic noise term, which is considered by adding terms that violate the hermiticity of the GPE operator.The paper emphasizes that the states under study and the resulting specific phase transitions in BEC are possible only in the presence of dissipation.
The dissipative NLSE also arise in the description of solitons in nonlinear media such as the cavity of the mode-locked lasers (the so-called Haus master equation [14] that is the (1+1)-NLSE with a non-Hermitian term) and related models [15] including multidimensional [16] and nonlocal [17] ones.
Most of the studies devoted to the mentioned nonlinear systems, including the BEC, deal with the local form of the NLSE since it is simpler for the mathematical analysis especially for the numerical one.While such simplification is reasonable for short range interactions, it can not be used for long range interactions.The example of the last ones is dipole-dipole interaction [18,19] that significantly affects the dynamics of the BEC [20,21].In [22], the author studied the influence of the nonlocal interaction of condensate particles and an external periodic field on the BEC dynamics in the framework of the mean field theory and nonlocal generalized GPE without dissipative terms.Therefore, it is of interest to study the nonlocal NLSE as a more general problem since the well studied local form can be treated as the limiting case of the nonlocal one in a manner.
Using the nonlocal form of the model equations, we apply the semiclassical formalism to the problem under consideration.The semiclassical approximation is widely used for linear equations of quantum mechanics.Some modern semiclassical approaches based on ideas of the Maslov method [23] were also applied to some nonlinear problems (see, e.g., [24,25]).In [26][27][28], the formalism of semiclassical asymptotics for a generalized nonlocal GPE in a special class of trajectory concentrated functions is developed that corresponds to closed quantum systems.In [29][30][31], this formalism was applied to kinetic reaction-diffusion equations that correspond to open classical systems.The conception of this work is to combine the ideas of those approaches to solve the nonlinear problem corresponding to open quantum systems.
In this paper, following [26,27,29,30], we extend the method of semiclassical asymptotics as applied to a generalized nonlocal NLSE with a non-Hermitian term that is responsible for the dissipation.In particular cases, the equation considered in the paper transforms into a complex nonlocal GLE [10], as well as into a nonlocal NLSE with a complex potential [12,32].The general formalism is illustrated by an example.
The paper is organized as follows.In Section II, the original nonlinear problem is posed.
In Section III, we introduce the moment of the desired solution and give some additional notations.In Section IV, we explain what we mean by the semiclassically concentrated states, derive the classical equations corresponding to the nonlinear quantum problem, and introduce the class of functions where the asymptotic solutions are sought.In Section V, we deal with the auxiliary dynamical system that allows us to proceed to the linear partial differential equation associated with the original nonlinear problem.In Section VI, we derive this linear partial differential equation and construct the leading term of the asymptotic solution to the Cauchy problem for the original nonlinear equation under some algebraic conditions.The explicit analytical form of the semiclassical nonlinear evolution operator is given.In Section VII, we discuss the semiclassical symmetries for the problem under consideration.Section VIII provides the example for the presented formalism.Here, we apply our method to the specific NLSE that is the model of an atom laser.Analytical asymptotic solutions are obtained for this equation.In Section IX, we conclude with some remarks.
We consider solutions Ψ to the equation (2.1) that belong to the Schwartz space S with respect to ⃗ x and deal with the L 2 (R n x ) scalar product 2) The operators ẑ, ŵ satisfy the following commutation relations: Let us introduce the set S of functions A(z, t, ℏ)=A(⃗ p, ⃗ x, t, ℏ) that satisfy the following conditions for every fixed t ≥ 0: 1) A(⃗ p, ⃗ x, t, ℏ) ∈ C ∞ with respect to ⃗ p and ⃗ x; 2) A(⃗ p, ⃗ x, t, ℏ) and all its derivatives grow not faster that polynomials of |⃗ p| and |⃗ x| as 3) A(⃗ p, ⃗ x, t, ℏ) regularly depends on the parameter ℏ in a neighborhood of ℏ = 0.
We denote by A the set of pseudo-differential operators defined above.
These operators are Hermitian with respect to the scalar product (2.2).
Whenever no confusion arises, we will simplify our notation.We will drop the explicit dependence on ℏ in the functions and indicate it where appropriate, in particular, in solutions Ψ to equation (2.1), Ψ(⃗ x, t, ℏ) = Ψ(⃗ x, t).

III. EXPECTATION OF AN OPERATOR OVER FUNCTIONS FROM S
A non-Hermitian term in the operator of the equation ( where σΨ = dσ/dt, and H(ẑ, t)[Ψ] is given by (2.1).We will consider the solutions with σ Ψ (t, ℏ) = O(1) as ℏ → 0. Actually, this condition depends on the definition of the nonlinearity coefficient κ.We will limit our consideration to the case κ = O(1) since it is of the greatest interest from the physical point of view.It corresponds to the situation when the linear and nonlinear parts of the equation (2.1) are comparable.In this case, regularly perturbation theory in nonlinearity parameter does not yield qualitative result.
For the operator Â(t) ∈ A determined by its Weyl symbol Â(t) = A(ẑ, t), we define the expectation for the solution Ψ(⃗ x, t, ℏ) to (2.1) by the following relation: The evolution equation for the expectation ⟨ Â(t)⟩ Ψ on solutions to (2.1) reads The equations (3.4) can be significantly simplified in the semiclassical approximation and such simplified equations will determine the semiclassical evolution of the solutions to (2.1).In the next Section, we will clarify how we interpret the semiclassical limit within the framework of our approach.

IV. CLASS OF SEMICLASSICALLY CONCENTRATED FUNCTIONS
Definition IV.1.The function Ψ(⃗ x, t, ℏ) belongs to a class T t ℏ (Z(t), σ(t)) of functions semiclassically concentrated in a neighborhood of a trajectory z = Z(t) with a weight of σ(t) if for any operator Â = A(ẑ, t, ℏ) ∈ A with the Weyl symbol A(z, t, ℏ) the following relations hold Here, the functions Z(t) and σ(t) are functional parameters of the class T t ℏ (Z(t), σ(t)).The relation (4.1) is similar to the definition of the Dirac δ-function in terms of δ-sequence.
As is known, δ(x − x 0 ) is determined as a linear functional that maps any function f (x) from the specific space (e.g.smooth functions with compact support) to its value at the point x 0 .In (4.1), the analog of the function f (x) is an operator Â from the set A. The relation (4.1) maps this operator to its Weyl symbol A(z, t, ℏ) on the trajectory z = Z(t) as ℏ → 0.
Thus, the curve z = Z(t) acts as x 0 , the function Ψ acts as δ-sequence, and the operator Â along with its Weyl symbol A(z, t, 0) acts as a test function.Summing up, for every given moment t, the point z = Z(t) is determined in R 2n , and a function Ψ ∈ T t ℏ (Z(t), σ(t)) is semiclassically concentrated in a neighborhood of these points that constitute a trajectory In order to obtain the semiclassical evolution of the expectation ⟨ Â(t)⟩ Ψ for a solution Ψ(⃗ x, t, ℏ) to (2.1) from the class T t ℏ (Z(t), σ(t)), we go to the limit ℏ → 0 in (3.4).In view of (4.1) and properties of Weyl symbols, we obtain Here, A(z), B(z) , z = (⃗ p, ⃗ x), stands for Poisson bracket: In (4.3), we have used the property of the Weyl ordered pseudo-differential operators A(ẑ), , respectively, where that is as follows [34]: Let us consider two important corollaries of (4.3).
By analogy with [27], the system (4.7),(4.8) is termed the Hamilton-Ehrenfest system with dissipation (HESD) of the first order for the equation (2.1) in the class T t ℏ (Z(t), σ(t)).Note that the system (4.7),(4.8) contains the partial information about localization properties of the solutions to (2.1) in the space of the dynamical system of the first and zeroth moments of the function Ψ.From the derived equations (in particular, (4.7), (4.8)) it is clear that it makes sense to consider only the case Λ = O(1) within the framework of the semiclassical approximation.If Λ → ∞ as ℏ → 0, the solution would rapidly damp so that its dynamics would not be observed at damping times.If Λ → 0 as ℏ → 0, we could ignore the non-Hermitian part in the semiclassical approximation that corresponds to the case Λ = 0 in our approach.Note that the last case will be naturally included in our general formalism.
In order to derive explicit analytical expressions for approximate solutions to (2.1) in the form of semiclassical asymptotics with respect to the parameter ℏ, ℏ → 0, let us consider the class P t ℏ Z(t), S(t), σ(t) of functions that singularly depend on asymptotics parameter ℏ [26]: Here, Φ(⃗ x, t, ℏ) is a general element of the class P t ℏ (Z(t), S(t), σ(t)); the real functions Z(t) = ( ⃗ P (t), ⃗ X(t)), S(t), and positively defined function σ(t) are functional parameters of the class ; the function φ( ⃗ ξ, t, ℏ) belongs to the Schwartz space S with respect to the variables ⃗ ξ ∈ R n and regularly depends on ℏ in a neighborhood of ℏ = 0; the functions Z(t), S(t), φ( ⃗ ξ, t, ℏ), and σ(t) smoothly depend on t.Note that in [26] this class was introduced without normalization factor σ(t) since the solutions with conserved norm were considered there.Here, we have added this factor for convenience without loss of generality of the family of classes P t ℏ Z(t), S(t), σ(t) .
For ease of notation, we will abbreviate P t ℏ (Z(t), S(t), σ(t)) by P t ℏ when no confusion will arise.

Functions of the class P t
ℏ are concentrated at a point moving over the trajectory z = Z(t) in the sense (4.1), (4.2) [26].They have the special form determined by the ansatz (4.9).
We will term it as the semiclassical ansatz (4.9).
It was proved in [26,27] that, for functions from the class on a finite time interval t ∈ [0; T ], the following asymptotic estimates hold: Here, the following notations are used.The estimate Â = Ô(ℏ s ), s ≥ 0, in (4.10) stands for an operator Â such that 2n .The functions ∆ In particular, we have We use the contracted notation j = 1, n for j = 1, 2, . . ., n in (4.13) and below.Hereinafter, all calculations and commentaries are given for t ∈ [0; T ] where T < ∞.

V. MOMENTS OF FUNCTIONS FROM THE CLASS P t ℏ
In order to construct semiclassical asymptotics for solutions to (2.1) in the class (4.9), let us introduce the moments of a function Ψ(⃗ x, t, ℏ) ∈ P t ℏ as it follows Here, σ Ψ (t, ℏ) is the zeroth order moment of Ψ(⃗ x, t, ℏ), Z Ψ (t, ℏ) is the first order moment, and the central moments ∆ (α) Ψ (t, ℏ) are given by (4.11).Let us present the operators V (ẑ, t), V (ẑ, t), W (ẑ, ŵ, t), and W (ẑ, ŵ, t) from (2.1) in the form of formal power series in ∆ẑ and ∆ ŵ in a neighborhood of the trajectory z = Z(t).
Then, we have Hereinafter, the following notations are used: and so on by analogy.The expansions for V (ẑ, t) and W (ẑ, ŵ, t) are similar to (5.2).In view of (4.10), the part of the operators dropped (5.3) is estimated as Ô(ℏ (m+1)/2 ) where m is the greatest power of a polynomial of ∆ẑ, ∆ ŵ that was taken into consideration in the asymptotic expansion.In particular, if we limit ourselves to terms that were explicitly written in (5.2), then m = 2.
Let us denote the aggregate vector of moments of the function Ψ(⃗ x, t, ℏ) by It is known that g Ψ (t, ℏ) determines the function Ψ(⃗ x, t, ℏ) (see, e.g., [35]).However, the system of equations for (5.4) is an infinite set of ODEs.
Let us also introduce the aggregate vector g In (5.5), we suppose that all moments included in g Ψ (t, ℏ) are determined within accuracy of O(ℏ (M +1)/2 ) where M indicates both the greatest order of moments included in g Ψ (t, ℏ), t, ℏ . (5.6) Finally, let us introduce the aggregate vector g (M ) (t, ℏ, C) that is a particular solution of with integration constant C. The system (5.7) is termed the HESD of the M -th order.
Note that, in view of properties of the class P t ℏ , the following relations hold [26]: Ψ (t) + O(ℏ 3/2 ). (5.8) From (4.1), (4.2), it also follows that The asymptotic solutions can be obtained in those classes P t ℏ Z(t), S(t), σ(t) that are defined by functions determined by equations (4.7), (4.8) and a function S(t) given by (5.10) Due to the uniqueness of the solution to the Cauchy problem for (5.6) and (5.7), we have In view of (5.9), the algebraic condition (5.11) is degenerate for the integration constants of the HESD of the first order (4.7), (4.8) in the class P t ℏ with the given Z(t) and σ(t).Hence, the functions Z(t) = Z (0) (t) and σ(t) = σ (0) (t) do not depend on C. Thus, equations (4.7), (4.8) form the system (5.7) for M = 1.Next, we will show that we have to solve at least the HESD of the second order to construct the leading term of asymptotics for the solution Ψ(⃗ x, t) to (2.1).In order to construct higher order asymptotics, we must solve the HESD of a higher order.The equations of the HESD of the second order are given explicitly in Appendix A.
Thus, (6.5) can be written as follows: Grouping terms of various orders with respect to ℏ in (6.7), we derive the system of equations for functions Ψ (k) , k = 0, M − 2: In particular, for M = 2, one readily gets that (6.8) yields the following single equation for the leading term of asymptotics: Following [26], we term Ψ (0) (⃗ x, t, ℏ) by the solution with right-hand side accuracy of O(ℏ 3/2 ) in the sense that it is generated by (6.7) with accuracy of Ō(ℏ 3/2 ).Hereinafter, we suppose that L2 is the operator from the expansion (6.4) satisfying (6.6).
Let us introduce the auxiliary linear equation with a solution Φ(⃗ x, t, ℏ, C) ∈ P t ℏ (Z(t), S(t), σ(t)).Then, the following proposition holds [26] Proposition VI.1.Let Ψ(⃗ x, t, ℏ) ∈ P t ℏ (Z(t), S(t), σ(t)) be a solution to the Cauchy problem for (2.1) with the initial condition Then, we have This proposition follows from the uniqueness of the solution to the Cauchy problems for (5.7) and (6.10).
To put it differently, the leading term of the asymptotic solution to the Cauchy problem for the original nonlinear equation (2.1) in the class P t ℏ can be found among solutions to C-parametric family of linear equations (6.10).The search of the appropriate solution is reduced to the algebraic conditions The parameters C (integration constants for (5.7)) that meet the condition (6.14) can be determined by initial conditions for the moments of the function Ψ, g Ψ (t, ℏ) t=0 .Following [27], we term (6.10) as the associated linear Schrödinger equation with dissipation (ALSED).It can be written as follows: where H(t, C, ℏ) = H (0) (t) + ℏH (1) (t, C), ⟨∆ẑ i ∆ẑ j + ∆ẑ j ∆ẑ i ⟩ is the dispersion matrix (matrix of the second order central moments).

VII. SEMICLASSICAL SYMMETRY OPERATORS
The solutions to the equation (6.15), (6.16) that determine the asymptotic solutions to (2.1) can be generated using symmetry operators.
Let a(t) ∈ C 2n be a solution to ȧ = JH zz (t)a.(7.1) Then, the operator is a symmetry operator of the first order for the equation (6.15), (6.16).
In view of the explicit form for the evolution operator (6.17), (6.19), an asymptotic solution to (2.1) in the class P t ℏ (Z(t), S(t), σ(t)) can be written as follows: Then, the function Ψ(0) (⃗ x, t) given by is also an asymptotic solution to (2.1).In terms of the evolution operator Û (t) (6.19), the relation (7.4) reads The relation (7.5) indicates clearly that the functions Ψ(0) and Ψ (0) are related nonlinearly due to the nonlinearity of Û (t).
Let a k (t) ∈ C 2n , k ∈ N, be linearly independent solutions to (7.3) satisfying the skew- where {a 1 , a 2 } = ⟨a 1 , Ja 2 ⟩, and δ km the Kronecker delta.Then, the respective symmetry operators for the ALSE form Heisenberg's Lie algebra: Thus, the linear symmetry operators (7.7) for ALSE that form the Lie algebra (7.8) generate the nonlinear approximate symmetry operators for the original nonlinear equation (2.1).Using the set of operators (7.7), one can construct analogs of the well-known Fock states [37] for the nonlinear equation (2.1).

VIII. EXAMPLE
Hereunder, we illustrate the formalism proposed with the simple but nontrivial example.
Let us consider the model equation It was derived in [9] for description of the field of the BEC in an atom laser that is a fundamentally open system.This equation is the reduction of the system of two related equations.The first one, the GPE, describes the field Ψ of condensed atoms.The second one, the reaction-diffusion equation, describes the density of uncondensed atoms.As we noted earlier, we operate with a nonlocal form of the nonlinearity within the framework of our formalism.The cases of c 1 > 0, ϵ > 0, Λ > 0, and c 2 ≷ 0 were considered in [9].
These coefficients depend on the coupling constant between condensed and uncondensed atoms [38], parameters of the laser pumping, parameters of a trap, and properties of atoms themselves such as the effective mass and self interaction strength.
In our notations, we have Let us pose the initial condition of the form that implies that some amount of condensate is trapped at the initial moment of time.
In view of (8.4) and symmetries in coefficients (8.2), one readily gets The Cauchy problem for (4.7) is given by σ and its solutions reads Note that, in our formalism, the function σ(t) must be positively defined.Hence, for Λ > 0, ϵ < 0, and N κ < ϵ, the asymptotics can be constructed only for 0 The consideration of a finite time interval when constructing asymptotics is reasonable since the exact solutions to the original nonlinear equation do not necessarily exist for infinite time interval.
The condition ϵ > 0 means that the threshold condition (gain=losses) is met [9].The formula (8.7) shows that, for great t, in zeroth approximation by ℏ, the evolution of the condensate in a trap is affected by the effective pump (effective implies that ϵ takes account of both the gain and losses) and nonlinearity factor.The positive nonlinearity factor, κ > 0, corresponds to the respective interatomic interaction.The increase in the effecting pump leads to the increase in the amount of condensate, while the increase in the respectively interaction strength leads to the decrease in the amount of condensate.
Next, we calculate the matrix coefficient In view of the symmetry of the matrix ∆ (1) 2 (t), we denote Then, equation (A3) reads while equation (A5) reads (0) = 0. ( The solution to (8.11) is given by The system(6.18)can be written as Note that the solutions to equation (A3) (that takes the form of (8.9), (8.10) in this particular example) can be expressed via the solutions to (6.18) (that is given by (8.13) in this example) as follows: 2 (0)M (t).(8.14)The exact solutions to (8.13)The substitution of (8.15) and ( 8.3) into (6.19)yields We also compare the analytical asymptotic solution Ψ (0) (x, t) with the numerical solution to (8.1) in Fig. 1 and 2. The numerical solution was obtained using the difference scheme that is common for the NLSE [39].We made the second-order spatial discretization using the method of lines with 1000 points along x and the Dirichlet conditions at x = ±2.The time integration on the spatial mesh was based on the Strang-Marchuk splitting method [40] according to the following scheme: where Î is the identity operator corresponding to the identity matrix on the spatial mesh and τ = t n+1 − t n is the time step.The operator Â = 1 iℏ c 1 p2 − iℏΛp 2 is the differential operator from (8.1) on the spatial mesh, and the rest terms of (8.1) are included in B[Ψ(t)].
The initial condition was equal for both the asymptotic and numerical solutions.
In our method, κ and Λ can be exactly equal to zero.In such way, we can construct asymptotic solutions to the linear Schödinger equation with a non-Hermitian term and to the NLSE without non-Hermitian term, respectively.However, we do not assume these parameters to be small in general.Hence, the solution behaviour can drastically change depending on κ and Λ. Fig. 1 and 2 show that the dynamics of solutions significantly changes when we put one of parameters κ or Λ as well as both of them to be not equal to zero.Judging by the numerical solutions, our analytical asymptotic solutions are reasonably accurate for small ℏ regardless of the presence of a non-Hermitian part.In this paper, we apply the semiclassical asymptotic approach to the NLSE with nonlocal cubic nonlinearity and a non-Hermite operator in the L 2 space (2.1).The non-Hermitian part of the equation accounts the environment impact on the system.We introduce the class of semiclassically concentrated functions T t ℏ in which we deduce evolution equations for the squared norm of the solution σ Ψ (t) and for the first moments Z(t) describing a localization of the semiclassical solutions.These equations already provide partial but important information about the solutions to the equation.They form a dynamical system that can be considered as a nonlinear analogue of the equations of classical mechanics for a linear quantum mechanical equation.However, unlike the linear case, the dynamic moment equations depend both on the symbol of the equation operator and on the class of functions in which solutions of the equation are constructed.
To obtain the leading term of the semiclassical asymptotics of the Cauchy problem for (2.1) explicitly within accuracy of O(ℏ 3/2 ), we follow the approach developed earlier for the GPE with a Hermitian equation operator (see [27] and references therein) in the class P t ℏ of trajectory concentrated functions (4.9).The semiclassical approach required significant modification for the non-Hermitian NLSE compared to the Hermitian one due to non-conservation of the square modulus ||Ψ|| 2 of the solution to (2.1).Although the general scheme of the method for constructing semiclassical asymptotics remains the same as for the nonlocal NLSE with a Hermitian operator, its implementation is subject to change.Following this scheme, we obtain a higher-order dynamical moment system (Hamilton-Ehrenfest system with dissipation) with accuracy of O(ℏ 3/2 ) and an associated linear Schrödinger equation.Together with the algebraic conditions (6.14), these equations allow us to construct a solution to the Cauchy problem in terms of the leading term of the semiclassical asymptotics with accuracy of O(ℏ 3/2 ).The construction of higher corrections also does not cause principal difficulties since they can be obtained using the evolution operator (6.17), (6.19), which is given explicitly in semiclassical approximation.Also, the semiclassical symmetry operators are constructed.The general results are illustrated by the particular example of the non-Hermitian NLSE that admits explicit analytical solutions in the semiclassical approximation.The example is based on the model equation of an atom laser [9] that is a fundamentally open system.
The approach proposed is a new tool for the analytical study of open quantum systems.
The semiclassical approximation was well studied for closed quantum systems using various approaches.One of their common drawbacks is that the error of semiclassical solutions (compared to the exact one) usually grows over time.However, for dissipative systems, the error can not grow indefinitely within the framework of our semiclassical approach since the dynamics of both the exact and asymptotic solutions damps over time due to the dissipation.
Hence, the error of semiclassical approximation should be bounded function with respect to time for such systems.It means that the time-limited semiclassical approaches like ours can be even more natural and useful for study of open quantum systems compared to closed ones.This encourages us to further develop our method for systems with a more complex geometries of the localization domain in the future.For c 1 > 0, c 2 > 0, ϵ > 0, the solution to (8.13) is as follows: (B1) Elements of the matrix Q(τ ) read Here, 2 F 1 (a, b; c; z) is the hypergeometric function [41], and G m,n p,q   z a 1 , ..., a p b 1 , ..., b q   is the Meijer G-function [42].


is the 2n × 2n symplectic identity matrix, I n is the n × n identity matrix.The commutators and anticommutators of operators Â and B are denoted by [ Â, B] = Â B − B Â and [ Â, B] + = Â B + B Â, respectively.The scalar products of vectors from R n and R 2n are denoted by ⟨a, b⟩ = j a j b j .

Figure 1 . 2 Figure 2 .
Figure 1.Dependence of |Ψ (0) (x, t)| 2 on x for various t and ℏ = 0.2 (solid lines).Dashed lines are for the respective numerical solution Appendix B: Solution to the system(8.13) can be expressed via the Meijer G-function.It is quite cumbersome and, for that reason, is given in Appendix A1.