Publication Cover
Applicable Analysis
An International Journal
Volume 101, 2022 - Issue 15
662
Views
0
CrossRef citations to date
0
Altmetric
Research Article

On global attractors for 2D damped driven nonlinear Schrödinger equations

&
Pages 5490-5503 | Received 26 Oct 2020, Accepted 18 Feb 2021, Published online: 03 Mar 2021

Abstract

The well-posedness and global attraction are proved for nonautonomous 2D nonlinear damped Schrödinger equations with almost periodic external source in a bounded region. The proofs rely on a development of the method of energy equation. The research is inspired by the problems of laser and maser coherent radiation.

2020 Mathematics Subject Classifications:

1. Introduction

We consider the weakly damped driven nonlinear Schrödinger equation in a bounded region ΩR2 with a smooth boundary Ω, (1) iψ˙(x,t)=Δψ(x,t)+f(ψ(x,t))iγψ(x,t)+p(x,t),xΩ,tR,(1) where the ‘friction coefficient’ γ>0, and the ‘pumping’ p(,t) is an almost periodic function. All the derivatives here and below are understood in the sense of distributions. Equations of type (Equation1) arise for instance in plasma physics [Citation1] and in optical fibers models [Citation2].

We impose the Dirichlet boundary condition (2) ψ(x,t)=0,xΩ,tR.(2) Our analysis of the nonautonomous Schrödinger equation is motivated by mathematical problems of laser and maser coherent radiation [Citation3,Citation4]. In particular, the maser action is described by the damped driven nonlinear Maxwell–Schrödinger equations in a bounded cavity, and the Dirichlet boundary condition (Equation2) holds in the case of a solid-state gain medium [Citation3].

We will identify complex numbers zC with two-dimensional real vectors z=(Rez,Imz). We assume that the nonlinear term is potential, i.e. (3) f(ψ)=U(ψ),ψCR2.(3) Here U(ψ)C3(R2) is a real function, and the derivatives mean derivatives with respect to ψ1=Reψ and ψ2=Imψ, i.e. (4) f(ψ):=(1U(ψ),2U(ψ)),ψR2.(4) Accordingly, equation (Equation1) is an identity of R2-valued distributions, and the multiplication by i in (Equation1) means the application of the real 2×2-matrix (0110). Similarly, f(ψ) and f(ψ) are tensors: (5) fij(ψ):=jfi(ψ1,ψ2)=jiU(ψ1,ψ2),fijk(ψ):=jkfi(ψ1,ψ2),i,j,k=1,2.(5) Note that f(ψ) is a real symmetric matrix. The Schrödinger equation (Equation1) under the boundary conditions (Equation2) can be written in ‘Hamiltonian form’ (6) iψ˙(t)=DH(ψ(t))iγψ(t)+p(ωt),tR,(6) with the Hamilton functional (7) H(ψ):=Ω12|ψ(x)|2dx+U(ψ),U(ψ):=ΩU(ψ(x))dx.(7) Our main results. For the damped driven Schrödinger equation (Equation1) with almost-periodic pumping in a bounded region, we establish

(I).

The well-posedness of the Cauchy problem for equation (Equation1) (Theorem 3.1);

(II).

Global attraction in the Sobolev norm H1 of all finite energy solutions to a compact subset in H1 (Theorem 4.2);

Let us comment on previous works in these directions. In the case p(x,t)0 (and γ=0), the well-posedness for nonlinear Schrödinger equations of type (Equation1) was established in [Citation5]. The approach [Citation5] relies essentially on the conservation of charge and energy.

Compact attractors in the energy space H1 for autonomous damped nonlinear Schrödinger equations were constructed (i) on the circle Ω=R/Z or a bounded interval ΩR by Ghidaglia [Citation6], (ii) on a bounded region ΩR2 by Abounouh and Goubet [Citation7,Citation8] and (iii) on the entire space Ω=RN with N3 by Laurençot [Citation9]. Similar results were obtained by Ghidaglia and Heron [Citation10] for the Ginzburg–Landau equations on a bounded region ΩRn with n = 1, 2. The methods of these papers rely on the method of energy equation introduced by Ball [Citation11]. Note that the pumping terms (or ‘external force’) in [Citation7–10] do not depend on time. In [Citation6], bounded absorbing sets are established for the pumping term depending on time, while the attractor is constructed for the autonomous case and for time-periodic pumping. The papers [Citation12–14] deal with smoothness properties of functions from attractors of 1D nonlinear Schrödinger equations.

Tao [Citation15] established the existence of a global attractor for radial solutions of nonlinear autonomous defocusing Schrödinger equation without damping in Rn with n11. The paper [Citation16] concerns the well-posedness and decay of solutions to the 2D damped Schrödinger equation in a bounded region.

Recently Cazenave and Han established the long-time behavior for nonlinear autonomous Schrödinger equation in Rn with a nonlinear subcritical dissipation [Citation17].

For nonautonomous PDEs, the theory of attractors was started by Haraux [Citation18,Citation19] in the context of nonlinear parabolic and damped wave equations with almost periodic external force. The method of energy equation [Citation11] was extended to nonautonomous PDEs first in [Citation20] for 1D Schrödinger equations on the circle, and in [Citation21] for nonautonomous KdV and 2D Navier-Stokes equations.

In [Citation22], the theory of attractors of nonautonomous equations was profoundly developed by Chepyzhov and Vishik [Citation22] in the context of Navier–Stokes equations, Ginzburg–Landau equations, damped nonlinear wave equations, and reaction–diffusion equations. A great deal of attention was focused there on the classes of translation-compact external forces. In [Citation23], the theory [Citation22] was extended to more general classes of admissible external forces, such as translation compact, time-regular, space-regular, and normal in the terminology of [Citation23]. Attractors of nonautonomous Schrödinger equations with nonperiodic pumping were not considered previously.

Let us comment on our approach. First, we prove the existence of a weak solution, using the standard Galerkin approximations. Our choice of dimension two is motivated by the Trudinger inequality, which ensures uniqueness of the weak solution by Theorem 3.6.1 of [Citation5]. Further, we develop the method of energy equation [Citation11,Citation20,Citation21] to prove that (i) the weak solution is the strong solution, and (ii) the solution continuously depends on the initial state and the pumping.

To prove the convergence in H1-norm of all finite energy solutions to a compact attractor, we extend Ball's ideas [Citation11] and the general theory [Citation22] to nonlinear Schrödinger equations.

We consider almost-periodic pumping which is translation-compact in the terminology [Citation22,Citation23]. This condition is physically relevant in the problem of laser and maser coherent radiation, where the pumping is a mixture of harmonics with a discrete set of frequencies [Citation3,Citation4]. We suppose that this condition can be relaxed by using methods [Citation23]. We do not touch the continuous dependence of the attractor on the pumping function, that will be considered elsewhere.

2. A priori estimates for smooth solutions

Here we prove a priori estimates for sufficiently smooth solutions to the Schrödinger equation (Equation1). The estimates ensure the existence of bounded absorbing sets in the energy norm. We will get rid of the smoothness assumption in the next section. Our assumptions are as follows.

Nonlinearity. We will assume that the potential U(ψ)C3(C), satisfies the following estimates (8) ϰ1|ψ|4b1U(ψ)C|ψ|4,(8) (9) Ref(ψ)ψ¯ϰ2U(ψ)b2,(9) (10) |f(ψ)|ϰ3(1+|ψ|3),(10) (11) |f(ψ)|ϰ4(1+|ψ|2),(11) (12) f(ψ)b3,(12) (13) |f(ψ)|ϰ5(1+|ψ|),(13) with all ϰl>0 and all blR. Let us recall that f(ψ) in (Equation11) and (Equation12) denotes the matrix (Equation5). For example, the potentials (14) U(ψ)=a2|ψ|4+a1|ψ|2+a0(14) satisfy all these conditions if akR and a2>0. Any smooth potential which differs from (Equation14) only in a bounded domain |ψ|C also satisfies these conditions.

Denote Hs:=Hs(Ω), Lp:=Lp(Ω). Let us introduce the Hilbert phase space (15) E:=H01(Ω),ψE=ψ,(15) where stands for the norm in L2. Then E=H1.

The pumping. We will assume that (16) pCb(R,E) is a uniform almost periodic function oftR,p0:=pCb(R,E)<.(16) By definition (see Appendix of [Citation24] for the case of scalar functions), a function pCb(R,E) is uniformly almost periodic if for any sequence tkR there is a subsequence tk such that the translations pk(t):=p(t+tk) are uniformly converging, (17) pkCb(R,E)p,k.(17)

2.1. Energy estimates

First we prove a priori estimates in the norm of the energy space E.

Lemma 2.1

Let conditions (Equation8)–(Equation12) and (Equation16) hold, and ψ(x,t)C2(Ω¯×R¯) be a solution to (Equation1) with ψ(,0)E. Then (18) ψ(t)E2C0eα±t+D0,±t0,(18) where C0=C0(ψ(0)E), α±>0, and D0 depends on p0 only.

Proof.

Differentiating the Hamilton functional H(t):=H(ψ(t)) and using (Equation6), we obtain the energy balance (19) H˙(t)=DH(ψ(t)),ψ˙(t)=DH(ψ(t)),iDH(ψ(t))+Δψ+f(ψ),γψ(t)ip(t),tR,(19) where the brackets stand for duality, which is the inner product in the real Hilbert space L2(Ω)R2. Hence the first term on the right-hand side vanishes. Therefore, (20) H˙(t)=Δψ(t)+f(ψ(t)),γψ(t)ip(t)=γψ(t),ψ(t)γf(ψ(t)),ψ(t)ψ(t),ip(t)f(ψ(t)),ip(t),tR.(20) According to (Equation9), (21) f(ψ(t)),ψ(t)ϰ2U(t)b2|Ω|,U(t):=ΩU(ψ(x,t))dx.(21) Further, (Equation8), (Equation10), and (Equation16) imply that (22) |f(ψ(t)),ip(t)|ϰ3(1+|ψ(x,t)|3)|p(x,t)|dxϰ3p(t)L1(Ω)+(γϰ1ϰ22|ψ(x,t)|4+C|p(x,t)|4)dxγϰ22U(t)+C(p0)(22) since p(τ)EL4(Ω) by the Sobolev embedding theorem in dimension two. Moreover, (Equation16) gives (23) |ψ(t),ip(t)|γ4ψ(t)2+1γp(t)2γ4ψ(t)2+1γp0.(23) Now (Equation20)–(Equation23) imply (24) H˙(t)γψ(t)2γϰ2U(t)+γ4ψ(t)2+γϰ22U(t)+C1(p0)α+H(t)+C1(p0),tR,(24) where α+:=γ2min(3,ϰ2)>0. Hence the Gronwall inequality gives (25) H(t)C(H(0))eα+t+D(p0),t0.(25) Now (Equation8) implies (Equation18) for t0.

For t0, the bound (Equation18) follows by similar arguments applied to the time-reversed equation (Equation1) for the function φ(t)=ψ(t): (26) φ˙(t)=iDH(φ(t))+γφ(t)+ip(t),tR.(26) Now for H(t):=H(φ(t)), similarly to (Equation20), we obtain (27) H˙(t)=γφ(t),φ(t)+γf(φ(t)),φ(t)+φ(t),ip(t)+f(φ(t)),ip(t).(27) In this case, the bound (Equation21) becomes (28) |f(φ(t)),φ(t)|2ϰ3ϰ1U(t)+C|Ω|,(28) which follows from (Equation8) and (Equation10). Now (Equation27) together with (Equation22), (Equation23), and (Equation28) imply (29) H˙(t)αH(t)+C,t>0,(29) with some α>0. Hence the Gronwall inequality gives H(t)C(H(0))eαt+C for t>0, which implies (Equation18) for t<0 by (Equation8).

3. Well-posedness

In this section, we prove the well-posedness of the Cauchy problem for the Schrödinger equation (Equation1) in the energy class. The key role is played by the energy equation (Equation20).

Theorem 3.1

Let conditions (Equation8)–(Equation12), (Equation16) hold. Then

(i)

For any initial state ψ(0)E equation (Equation1) admits a unique solution ψ(t)C(R,E)C1(R,E), and the bounds (Equation18) hold.

(ii)

The energy equation (Equation20) holds, (30) ddtH(ψ(t))=Δψ(t)+f(ψ(t)),γψ(t)ip(t),tR.(30)

(iii)

The map S:(ψ(0),p)ψ() is continuous (31) S:E×C(R,E)C(R,E).(31)

3.1. Weak solutions

First, we construct ‘weak solutions’.

Lemma 3.2

Let the conditions of Theorem 3.1 hold. Then, for any initial state ψ(0)E, equation (Equation1) admits a unique solution ψ(t)C(R,Ew)Lloc(R,E)Wloc1,(R,E), and the bounds (Equation18) hold.

Proof.

The existence of a weak solution ψ(t)C(R,Ew)Lloc(R,E)Wloc1,(R,E) follows by the Galerkin approximations. We recall this construction in Appendix 1, since we will use it in the proof of (Equation30). The uniqueness of this solution in the case p(x,t)0 and γ=0 is deduced in Theorem 3.6.1 of [Citation5] from the Trudinger inequality [Citation25]. The proof of uniqueness for γ>0 and p(x,t) satisfying (Equation16) remains almost unchanged: we give the required modifications in Appendix 2. Uniqueness implies that the convergence (EquationA5) actually holds for the entire sequence of Galerkin approximations: for any s<1 (32) ψm(t)C([T,T],Hs)ψ(t),m,T>0.(32) The uniform bounds (EquationA2) hold for the Galerkin approximations ψm. Hence the bounds (Equation18) hold also for their limit ψ(t).

3.2. Energy equation

Let us prove the energy equation (Equation30).

Proposition 3.3

Let the conditions of Theorem  3.1 hold. Then the energy equation (Equation30) holds for any solution ψ(t)C(R,Ew)Lloc(R,E)Wloc1,(R,E) to (Equation1).

Proof.

The identity (Equation30) for smooth solutions can be rewritten as 12tψ(t)2+γψ(t)2+tU(ψ(t))=ψ(t),ip(t)+f(ψ(t)),γψ(t)ip(t),tR.Equivalently, (33) 12[e2γtψ(t),ψ(t)]|0T=[e2γtU(ψ(t))]|0T+2γ0Te2γtU(ψ(t))dt0Te2γtf(ψ(t)),γψ(t)+ip(t)dt0Te2γtψ(t),ip(t)dt,TR.(33) We will deduce (Equation33) by passage to the the limit in the corresponding equation for the Galerkin approximations. Namely, the Galerkin equation (EquationA1) imply that for any TR (34) 12[e2γtψm(t),ψm(t)]|0T=[e2γtU(ψm(t))]|0T+2γ0Te2γtU(ψm(t))dt0Te2γtf(ψm(t)),γψm(t)+ipm(t)dt0Te2γtψm(t),ipm(t)dt.(34)

Lemma 3.4

(a)

The right-hand side of (Equation34) converges to the same of (Equation33).

(b)

The left-hand side of (Equation33) is dominated by the limit of the left-hand side of (Equation34).

Proof.

It suffices to consider the case T>0.

(a) The convergence (Equation32) of the Galerkin approximations with any s<1 implies that (35) ψm(t)C([0,T],Lr)ψ(t),m.(35) with any r< by the Sobolev embedding theorem in dimension two. Hence the convergence of the first three terms on the right-hand side of (Equation34) to the same of (Equation33) follows.

The last term of (Equation34) converges to the last term of (Equation33) by the Lebesgue theorem due to (1) the uniform bounds (EquationA2), (2) the strong convergence pm(t)Ep(t) and (3) the weak convergence ψm(t)Eψ(t) for all t0. This weak convergence follows from (Equation32) and (EquationA2).

(b) The second assertion of the lemma follows from the weak convergence ψm(T)Eψ(T), and the strong convergence ψm(0):=Pmψ(0)Eψ(0).

As a result, we obtain in the limit the energy inequality: the left-hand side of (Equation33) is dominated by its right-hand side for every TR.

On the other hand, ψ(t)C(R,Ew) by Lemma 3.2, and hence, ψ(T)E. Therefore, the same lemma implies the existence of a solution ψ~(t)C(R,Ew)Lloc(R,E)Wloc1,(R,E) starting from the initial state ψ~(T)=ψ(T)E. The same arguments of Lemma 3.4 show that the right-hand side of (Equation33), with ψ~ instead of ψ, is dominated by its left-hand side.

Finally, the uniqueness of a solution in Lemma 3.2 implies that ψ~(t)=ψ(t) for tR. Hence the integral energy equations (Equation33) holds for tR. We conclude that (Equation30) holds for a.a. tR and in the sense of distributions. The proposition is proved.

3.3. Strong solutions

Now we can prove that ψ(t)C(R,E)C1(R,E). Indeed, ψ(t)C(R,Ew)Lloc(R,E) by Lemma 3.2. On the other hand, the energy equation (Equation33) holds for all TR, and the right-hand side of (Equation33) is a continuous function of TR by the same arguments as in the proof of Lemma 3.4. Hence (Equation33) implies that the norm ψ(t)E is also a continuous function of tR. Therefore, ψ(t)C(R,E), and the equation (Equation1) together with Lemma A.1, (ii) and condition (Equation16) imply that ψ˙(t)C(R,E).

3.4. Continuous dependence

The continuity of the map (Equation31) follows by similar arguments. Namely, let ψ(0),ψn(0)E, p,pnC(R,E) and (36) supnψn(0)E<,ψn(0)Eψ(0);pnC(R,E)p,n.(36) Consider the corresponding unique solutions ψn(t)C(R,E)C1(R,E) to (Equation1). Similarly to (Equation32), we obtain for any s<1 and any T>0 (37) ψn(t)C([T,T],Hs)ψ(t),n.(37) Hence the Banach theorem on weak compactness implies that ψn(t)Eψ(t),n,tR.On the other hand, the energy equation (Equation33) implies the norm-convergence ψn(t)EC(R)ψ(t)E,nby the same arguments as in the proof of Lemma 3.4, using the uniform bounds (Equation18) and the uniform convergence (Equation37). Hence ψn(t)C(R,E)ψ(t),n.Theorem 3.1 is proved.

3.5. Weak continuity of the process

In conclusion, we prove a technical lemma which will be important in the next section.

Lemma 3.5

Let conditions (Equation8)–(Equation11) and (Equation16) hold, and let 1/3s<1. Then, for any T>0, the map S(T,0):(ψ(0),p)ψ(T) is continuous from Hs×C([0,T],E) to Hs.

Proof.

Let ψ(0),ψn(0)E, p,pnC([0,T],E) and (38) supnψn(0)E<,ψn(0)Hsψ(0);pnC([0,T],E)p,n.(38) Consider the corresponding unique solutions ψn(t)C(R,E)C1(R,E) to (39) iψ˙n(x,t)=Δψn(x,t)iγψn(x,t)+f(ψn(x,t))+pn(x,t),xΩ,tR.(39) These trajectories are equicontinuous on [0,T] with values in Hs. This follows from the interpolation inequality (EquationA4) since supnsupt[0,T]ψn(t)H1< by (Equation18), while supnsupt[0,T]ψ˙n(t)E< by equations (Equation39) together with Lemma A.1, (ii) and (Equation16).

Now the Arzelà–Ascoli theorem implies that there exists a subsequence ψn converging in C([0,T],Hs): (40) ψnC([0,T],Hs)ψ~,n.(40) The limit function ψ~(t)C([0,T],Hs) satisfies the Schrödinger equation (Equation1). This follows by the passage to the limit in the equation (Equation39), taking into account Lemma A.1, i) and the uniform convergence of pn in  (Equation38). Moreover, ψ~(0)=ψ(0) by the convergence of ψn(0) in (Equation38). Hence the uniqueness arguments as above imply that ψ~(t)ψ(t) for all t[0,T]. Therefore, the continuity of S(T,0) follows from (Equation40) since the limit does not depend on the subsequence n.

4. Global attractor

In this section, we prove the convergence in the H1-norm of all finite energy solutions to a compact attractor. Proposition 3.1,(i) allows us to define the continuous process Sp(t,τ) in E acting by (41) Sp(t,τ):ψ(τ)ψ(t),(41) where ψ(t)C(R,E)C1(R,E) is a solution to (Equation1). The bounds (Equation18) imply that the process admits a bounded uniform absorbing set in H1. Indeed, let BR denote the ball {ψE:ψER}. Denote B^:=BD+1, where D is the constant from (Equation18). Then (Equation18) implies that Sp(t,τ)BRB^,tτ>logA(R)α.Now we prove the existence of a uniform compact attractor in H1. Let us recall the definition of uniform attractor for nonautonomous equations, see, e.g. [Citation19], [Citation26, Definition A2.3], [Citation27, Definition (7)].

Definition 4.1

A set AE is called a uniform attractor of the process Sp(t,τ) if it is a uniformly attracting set and if any other closed attracting set contains A.

The term ‘uniformly attracting’ means that, for any bounded subset BE, (42) distE(Sp(t,τ)B,A)0,tτ.(42)

Theorem 4.2

Let the assumptions of Theorem 3.1 hold , γ>0, and let the pumping p(t) be an almost periodic function with values in E. Then the process Sp(t,τ) admits a compact uniform attractor in E.

For the proof, we develop the strategy [Citation11] for a nonautonomous equation. Namely, we define the functionals (cf. (2.3) and (2.4) of [Citation9] and (2.12), (2.25) of [Citation20]) (43) Φ(ψ,t):=H(ψ)+p(t),ψ,Ψ(ψ,t):=U(ψ)12f(ψ),ψ+12p(t),ψ,ψE,t0.(43)

Lemma 4.3

For any s1/2 and any t0, the functionals U() and Ψ(,t) are continuous on Hs and bounded on bounded subsets.

Proof.

Using (Equation10) and (Equation11), we obtain |U(ψ1)U(ψ2)|CΩ(1+|ψ1(x)|3+|ψ2(x)|3)|ψ1(x)ψ2(x)|dx,|f(ψ1),ψ1f(ψ2),ψ2||f(ψ1)f(ψ2),ψ1|+|f(ψ2),ψ1ψ2|C(Ω(1+|ψ1(x)|2+|ψ2(x)|2)|ψ1(x)ψ2(x)||ψ1|dx+Ω(1+|ψ2(x)|3)|ψ1(x)ψ2(x)|dx)C1Ω(1+|ψ1(x)|3+|ψ2(x)|3)|ψ1(x)ψ2(x)|dxApplying the Hölder inequality and the Sobolev embedding theorem, we obtain for any s1/2 |U(ψ1)U(ψ2)|+|f(ψ1),ψ1f(ψ2),ψ2|C(1+ψ1L43+ψ2L43)ψ1ψ2L4C(1+ψ1Hs3+ψ2Hs3)ψ1ψ2Hs.Similarly, |p(t),ψ1ψ2|p(t)ψ1ψ2Cp0ψ1ψ2Hs,for any s0.

Let ψ(t)C(R,E)C1(R,E) be a solution to (Equation1). Now the energy equation (Equation30) implies (44) ddtΦ(ψ(t),t)=2γΦ(ψ(t),t)+2γΨ(ψ(t),t)+p˙(t),ψ(t),t>0,(44) where the last term is well defined in the sense of distributions: for any test function φC(R+¯) (45) 0tp˙(s),ψ(s)φ(s)ds:=[p(s),ψ(s)φ(t)]|0t0t[p(s),ψ˙(s)φ(s)+p(s),ψ(s)φ˙(s)]ds.(45) To prove (Equation44), we differentiate using (Equation1) and (Equation30), ddtΦ(ψ(t),t)p˙(t),ψ(t)=H˙(ψ(t))+ψ˙(t),p(t)=γψ(t),ψ(t)γf(ψ(t)),ψ(t)+iψ(t),p(t)+if(ψ(t)),p(t)+iΔψ(t)γψ(t)if(ψ(t))ip(t),p(t)=2γH(ψ(t))+2γU(ψ(t))γf(ψ(t)),ψ(t)γψ(t),p(t)=2γΦ(ψ(t),t)+2γΨ(ψ(t),t).Integrating (Equation44), we obtain (46) Φ(ψ(t),t)=e2γtΦ(ψ(0),0)+0te2γ(tτ)[2γΨ(ψ(τ),τ)+p˙(τ),ψ(τ)]dτ,t0.(46)

4.1. Strong convergence

The following proposition is the key step in the proof of Theorem 4.2.

Proposition 4.4

For any sequences tk,τk with tkτk and ϕkB^, there is a subsequence k and ϕE such that (47) Sp(tk,τk)ϕkEϕ,k.(47)

Proof.

The sequence Sp(tk,τk)ϕk is bounded in E by (Equation18). Let us fix an arbitrary s[1/2,1). The inclusion EHs is compact by the Sobolev embedding theorem. Hence there is a subsequence strongly converging in Hs: (48) Sp(tk,τk)ϕkHsϕ,k.(48) Now to prove (Equation47) it suffices to check that for a subsequence {k}{k}, we have (49) lim supkSp(tk,τk)ϕkEϕE.(49) Recall that tkτk. Hence similarly to (Equation48), for any T>0, there exists a subsequence {k}{k} and an element ϕ(T)E such that (50) Sp(tkT,τk)ϕkHsϕ(T),k.(50) We set pk(t):=p(tkT+t) and denote (51) ψk(t):=Sp(tkT+t,τk)ϕk=Spk(t,0)Sp(tkT,τk)ϕk,t[0,T].(51) In particular, for t = T the convergence (Equation48) implies that (52) ψk(T):=Sp(tk,τk)ϕkHsϕ,k.(52)

Now we are going to apply Lemma 3.5(iv). At this step, we need the almost-periodicity (Equation17) of the pumping p(t). It implies that for a subsequence {k}{k}, we have (53) pkCb(R,E)p,k.(53) Hence (Equation50), (Equation51) and the continuity of the map (Equation31) imply the convergence for t[0,T], (54) ψk(t)=Spk(t,0)Sp(tkT,τk)ϕkHsψ(t):=Sp(t,0)ϕ(T),k.(54) In particular, for t = T the convergence (Equation52) gives that (55) ψ(T):=Sp(T,0)ϕ(T)=ϕ.(55) Now we apply the integral identity (Equation46) to solutions ψk(t) and ψ(t) of Equation (Equation1) with the pumping pk(t) and p(t) respectively. We obtain (56) Φ(ψk(T),T)=e2γTΦ(ψk(0),0)+0Te2γ(Tτ)[2γΨ(ψk(τ),τ)+p˙k(τ),ψk(τ)]dτ(56) (57) Φ(ψ(T),T)=e2γTΦ(ψ(0),0)+0Te2γ(Tτ)[2γΨ(ψ(τ),τ)+p˙(τ),ψ(τ)]dτ,(57) where the integrals with p˙k(τ) and p˙(τ) are definite according to (Equation45). Letting k in (Equation56), we obtain (58) lim supkΦ(ψk(T),T)Ce2γT+0Te2γ(Tτ)[2γΨ(ψ(τ),τ)+p˙(τ),ψ(τ)]dτ(58) by the following arguments:

  1. supkψk(0)E<, while the functionals Φ(,0) are bouded on bounded subsets of E,

  2. ψk(τ)Hsψ(τ) by (Equation54), while the functionals Ψ(,τ) are continuous on Hs by Lemma 4.3, since s[1/2,1), and

  3.  (Equation53) implies the convergence of the integrals with p˙k(τ) by their definition (Equation45).

Now eliminating the integral term from (Equation57) and (Equation58), we obtain (59) lim supkΦ(ψk(T),T)Ce2γT+Φ(ψ(T),T)e2γTΦ(ψ(0),0).(59) By definitions (Equation7) and (Equation43), Φ(ψk(T),T)=12ψk(T)E+U(ψk(T))+pk(T),ψk(T),Φ(ψ(T),T)=12ψ(T)E+U(ψ(T))+p(T),ψ(T),Now Lemma 4.3 together with (Equation48), (Equation52), and (Equation55) imply that limkU(ψk(T))=U(ϕ)=U(ψ(T)),limkpk(T),ψk(T)=p(T),ϕ=p(T),ψ(T).Hence (Equation59) becomes 12lim supkψk(T)ECe2γT+12ϕEe2γTΦ(ψ(0),0).According to the definition (Equation52), we can replace ψk(T) by Sp(tk,τk)ϕk. Now letting T, we obtain (Equation49) since γ>0.

4.2. Compactness

Let us denote by A the set of all points ϕ from (Equation47). This set is obviously closed. Now Theorem 4.2 will follow from the next lemma.

Lemma 4.5

  1. The set A is compact in E, and

  2. the set A is uniformly attracting in E.

Proof.

(i) Let us consider a sequence ϕnA. Then, for each n, (60) Sp(tnk,τnk)ϕnkEϕn,tnkτnk(60) as k, where ϕnkB^. Hence there exists a sequence k(n) such that (61) Sp(tnk(n),τnk(n))ϕnk(n)ϕnE0,tnk(n)τnk(n)(61) as n. However, Proposition 4.4 implies that for a subsequence n, we have (62) Sp(tnk(n),τnk(n))ϕnk(n)Eϕ¯,n.(62) Hence ϕ¯A, and (Equation61) implies that ϕnEϕ¯. Therefore, the first assertion of the lemma is proved.

(ii) Let us assume the contrary. Then there exists a sequence S(tk,τk)ϕk such that (63) distE(S(tk,τk)ϕk,A)ε>0,tkτk.(63) On the other hand, Proposition 4.4 implies that for a subsequence k, we have (64) S(tk,τk)ϕkEϕA,k.(64) This contradiction proves the second assertion of the lemma.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Additional information

Funding

The research is supported by the Austrian Science Fund (FWF) [grant numbers P28152–N35 and P34177–N].

References

  • Nozaki K, Bekki N. Low-dimensional chaos in a driven damped nonlinear Schrödinger equation. Phys D. 1986;21:381–393.
  • Blow KJ, Doran NJ. Global and local chaos in the pumped nonlinear Schrödinger equation. Phys Rev Lett. 1984;52:526–529.
  • Haken H. Laser theory. Berlin: Springer; 1984.
  • Nussenzveig H. Introduction to quantum optics. London: Gordon and Breach; 1973.
  • Cazenave T. Semilinear Schrödinger equations. New York: AMS; 2003.
  • Ghidaglia JM. Finite-dimensional behaviour for weakly damped driven Schrödinger equations. Ann Inst Henri Poincaré. 1988;5:365–405.
  • Abounouh M. Asymptotic behavior for a weakly damped Schrödinger equation in dimension two. Appl Math Lett. 1993;6:29–32.
  • Abounouh M, Goubet O. Attractor for a damped cubic Schrödinger equation on a two-dimensional thin domain. Differ Integral Equ. 2000;13(1–3):311–340.
  • Laurençot P. Long-time behaviour for weakly damped driven nonlinear Schrödinger equations in RN, N≤3. Nonlinear Differ Equ Appl. 1995;2(3):357–369.
  • Ghidaglia JM, Heron B. Dimension of the attractors associated to the Ginzburg–Landau partial differential equation. Physica D. 1987;28:282–304.
  • Ball J. Global attractors for damped semilinear wave equations. Discrete Cont Dyn Syst. 2004;10(1/2):31–52.
  • Goubet O. Regularity of the attractor for the weakly damped nonlinear Schrödinger equation. Appl Anal. 1996;60:99–119.
  • Goubet O. Regularity of the attractor for a weakly damped nonlinear Schrödinger equation in R2. Adv Differ Equ. 1998;3:337–360.
  • Goubet O, Molinet L. Global attractor for weakly damped nonlinear Schrödinger equations in L2(R). Nonlinear Anal Theory Methods Appl. 2009;71:317–320. hal-00421278.
  • Tao T. A global compact attractor for high-dimensional defocusing non-linear Schrödinger equations with potential. Dyn PDE. 2008;5:101–116.
  • Saanouni T. Nonlinear damped Schrödinger equation in two space dimensions. Electron J Differ Equ. 2015;2015(121):1–9.
  • Cazenave T, Han Z. Asymptotic behavior for a Schrödinger equation with nonlinear subcritical dissipation. Discrete Cont Dyn Syst A. 2020;40(8):4801–4819.
  • Haraux A. Two remarks on dissipative hyperbolic problems. In: Brezis H, Lions JL, editors. Nonlinear partial differential equations and their applications. College de France Seminar. Vol. VII. Research Notes in Mathematics. Vol. 122. London: Pitman; 1985. p. 161–179.
  • Haraux A. Systèmes dynamiques dissipatifs et applications, R.M.A. 17. Collection dirigée par Ph. Ciarlet et Lions JL. Paris: Masson; 1991.
  • Wang X. An energy equation for the weakly damped driven nonlinear Schrödinger equations and its application to their attractors. Phys D. 1995;88:167–175.
  • Moise I, Rosa R, Wang X. Attractors for noncompact nonautonomous systems via energy equations. Discrete Cont Dyn Syst. 2004;10(1):473–496.
  • Chepyzhov VV, Vishik MI. Attractors for equations of mathematical physics. Providence (RI): AMS; 2002.
  • Zelik SV. Strong uniform attractors for non-autonomous dissipative PDEs with non translation-compact external forces. Discrete Cont Dyn Syst Ser B. 2015;20(3):781–810.
  • Hale J. Ordinary differential equations. Malabar (FL): Krieger; 1980.
  • Adams RA. Sobolev spaces. New York: Academic; 1975.
  • Chepyzhov VV, Vishik MI. Non-autonomous dynamical systems and their attractors. In: Vishik M, editor. Asymptotic behaviour of solutions of evolutionary equations. Cambridge: Cambridge University Press; 1992. p. 111–150.
  • Chepyzhov VV, Vishik MI. A Hausdorff dimension estimate for kernel sections of non-autonomous evolution equations. Indiana Univ Math J. 1993;42(2):1057–1076.
  • Lions J-L. Quelques méthodes de résolution des problèmes aux limites non linéaires. Gauthier-Villars, Paris: Dunod; 1969.

Appendices

Appendix 1. Galerkin approximations

Denote by Em, m=1,2,, the linear span of the first m eigenfunctions of the Laplacian Δ on the region Ω with the Dirichlet boundary conditions (Equation2). We define the Galerkin approximations ψm(t) as solutions to the following finite-dimensional dissipative equations in Em (A1) iψ˙m(t)=Δψm(t))iγψm(t)+Pmf(ψm(t))+pm(t),tR;ψm(0)=Pmψ(0).(A1) Here Pm stands for the orthogonal projection of L2 onto Em, and pm(t):=Pmp(t). Applying the calculations from the proof of Lemma 2.1 to Equation (EquationA1), we obtain the uniform estimates of type (Equation18) for ψm with the same constants: (A2) ψm(t)E2C(ψ(0)E)eα±t+D,tR.(A2) Hence, the Galerkin approximations ψm(t) exist globally in time.

Let us show that these uniform estimates imply the existence of a limiting function ψ(t) of the Galerkin approximations. For this purpose, we will use the known continuity property of the nonlinear term:

Lemma A.1

Let nonlinear function f(ψ) satisfy condition (Equation11). Then

  1. For any r1 and s123r, the nonlinearity N:ψ()f(ψ()) is a continuous map HsLr.

  2. N:ψ()f(ψ()) is a continuous map EE.

  3. N is a continuous map E1L2.

Proof.

(i) Condition (Equation11) implies that |f(ψ1(x))f(ψ2(x))|C(1+|ψ1(x)|2+|ψ2(x)|2)|ψ1(x)ψ2(x)|.Hence, the Hölder inequality and the Sobolev embedding theorem give, for r1, (A3) f(ψ1)f(ψ2)LrC(1+ψ1L3r2+ψ1L3r2)ψ1ψ2L3rC(s)(1+ψ1Hs2+ψ2Hs2)ψ1ψ2Hs.(A3) (ii) In particular, (EquationA3) implies that N is continuous H2/3L2. It remains to note that EH2/3, while L2E. (iii) The map N:E1C(Ω¯) is continuous by the Sobolev embedding theorem.

Further, the Galerkin approximations are equicontinuous in Hs with any s<1 by the Dubinsky ‘Theorem on Three Spaces’ ([Citation28, Theorem 5.1]). Namely, this equicontinuity follows from the interpolation inequality: for any δ>0 (A4) ψm(t1)ψm(t2)Hsδψm(t1)ψm(t2)H1+Cδψm(t1)ψm(t2)E,t1,t2R.(A4) Here the first term on the right is small for sufficiently small δ>0, since supmsupt0ψm(t)E< by (EquationA2), while the second term is small for |t1t2|1, since supmsupt0ψ˙m(t)E< by the Galerkin equation (EquationA1) together with Lemma A.1, (ii), and (Equation16).

Now the Arzelà–Ascoli theorem implies that there exists a subsequence ψm converging in C(R,Hs): (A5) ψm(t)C([T,T],Hs)ψ(t),m.(A5) The limit function ψ(t)C(R,Hs) satisfies the estimates (Equation18) by the uniform bounds (EquationA2), and hence, ψ(t)C(R,Ew)Lloc(R,E).

The limit function ψ(t) satisfies the Schrödinger equation (Equation1) in the sense of distributions. This follows from the convergence (EquationA5) with s<1 close to 1 by passage to the limit in the Galerkin equation (EquationA1) taking into account Lemma A.1. Finally, the convergence ψm(0)ψ(0) as m implies that ψ(0)=ψ(0).

Finally, the equation (Equation1) implies that ψ(t)Wloc1,(R,E) by Lemma A.1, (ii).

Appendix 2. Uniqueness

Here we reduce the proof of the uniqueness of the solution to (Equation1) with γ>0 and fixed ψ(0)E to the case γ=0 considered in Section 3.6 of [Citation5]. Namely, let ψ1,ψ2Lloc(R,E)Wloc1,(R,E) be two solutions to (Equation1) with the same initial state ψ1(0)=ψ2(0). Then the difference z=ψ1ψ2 satisfies the equation iz˙(x,t)=Δz(x,t)iγz(x,t)+f(ψ1(x,t))f(ψ2(x,t)),xΩ,tR.This equation implies (A6) ddtz(t)L22+2γz(t)L22=2f(ψ1(t))f(ψ2(t)),z(t),tR.(A6) Equivalently, we have the equation (A7) ddt[e2γtz(t)L22]=2e2γtf(ψ1(t))f(ψ2(t)),z(t),tR.(A7) Finally, denoting h(t)=|ψ1(t)|+|ψ2(t)|, w(t)=eγtz(t) and using (Equation11), we obtain the inequality (A8) |ddtw(t)L22|CΩ(1+h2(t))|w(x,t)|2dx,tR,(A8) which coincides with the corresponding inequality in Section 3.6 of [Citation5]. Further, using the Trudinger inequality as in [Citation5], we obtain w(t)0.