Next Article in Journal
Marcus Stochastic Differential Equations: Representation of Probability Density
Previous Article in Journal
Reciprocal Hyperbolic Series of Ramanujan Type
Previous Article in Special Issue
Bounds of Different Integral Operators in Tensorial Hilbert and Variable Exponent Function Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Decay in the Energy Space of Solutions to the Damped Magnetic Radial Schrödinger Equation with Non-Local Nonlinearities

1
Institute of Mathematics and Informatics, Bulgarian Academy of Science, Acad. Georgi Bonchev Str., Block 8, 1113 Sofia, Bulgaria
2
Mathematics and Science Department, American University in Bulgaria, 1 Georgi Izmirliev Sq., 2700 Blagoevgrad, Bulgaria
3
Department of Mathematical Analysis and Differential Equations, Faculty of Applied Mathematics and Informatics, Technical University of Sofia, 1756 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2024, 12(19), 2975; https://doi.org/10.3390/math12192975
Submission received: 31 July 2024 / Revised: 29 August 2024 / Accepted: 23 September 2024 / Published: 25 September 2024
(This article belongs to the Special Issue Variational Problems and Applications, 2nd Edition)

Abstract

:
We will explore, in any space dimension d 4 , the decay in the energy space for the damped magnetic Schrödinger equation with non-local nonlinearity and radial initial data in H 1 ( R d ) . We will also display new Morawetz identities and corresponding localized Morawetz estimates.

1. Introduction

We will analyze the following Cauchy problem for the nonlinear defocusing magnetic Schrödinger equation with non-local nonlinearity in the high dimensional frame d 4 :
i t u + Δ x A u + i b ( t ) u k [ | · | ( d γ ) | u | 2 ] u = 0 , ( t , x ) [ 0 , ) × R d , u ( 0 , x ) = f ( x ) H 1 ( R d ) ,
with k R , where u = u ( t , x ) : [ 0 , ) × R d C , x A = i A , A = A 1 , , A d C l o c 1 ( R d { 0 } ; R ) , so that div A = 0 , Δ x A = x A · x A is self-adjoint on L 2 ( R d ) and b : [ 0 , ) C is a measurable function that contains dissipative and oscillatory terms. We shall also assume that
| A | 2 2 i A · L d 2 , ( R d ) , A L d , ( R d ) .
Moreover,
| x | x B L ( R d ) 2 ( d 1 ) ( d 3 ) ,
where the magnetic field B : R d M d × d ( R ) is defined by
B : = D A ( D A ) t ,
with
( D A ) i j = i A j , ( D A ) i j t = ( D A ) j i .
We will impose further the conditions on the nonlinear terms:
d 2 γ < d ,
b ( t ) , b ( t ) C ( [ 0 , ) ) with b ( t ) 0 and
B ( t ) = 0 t b ( s ) d s , inf t > 0 B ( t ) t 0 .
The last of the two conditions above, roughly speaking, means that every global solution of (1) behaves like the solution of the associated free equation (that is, b ( t ) = k = 0 ) as t + . The main goal of this paper is to show the decay of the solutions to (1) in the energy space. More explicitly, we will prove the following theorem.
Theorem 1. 
Let d 4 and k = 1 , and let u C ( [ 0 , ) ; H 1 ( R d ) ) be a global solution to (1) with radial initial data f H 1 ( R d ) such that (2) and the strict inequality in (3)–(5) are satisfied. Then, for 2 < r < 2 d d 2 , one achieves
lim t e B ( t ) u ( t , x ) L r ( R d ) = 0 .
Equation (1) is significant in many mathematical physics models. For instance, it was introduced in quantum mechanics to analyze the behavior of Bose–Einstein condensates by considering the self-interactions of charged particles, as discussed in [1,2,3], and the references therein. This has spurred numerous studies on the Schrödinger–Hartree equation. For example, Ref. [4] demonstrates the asymptotic completeness and the existence of wave operators for both the nonlinear Schrödinger equation with L 2 H 1 intercritical nonlinearity and the Schrödinger–Hartree equation. Subsequent improvements on these results for the Schrödinger–Hartree equation are found in [5]. Additionally, Refs. [6,7] employed the pseudo-conformal transform to study scattering solutions of the Schrödinger–Hartree equation in spaces with higher regularity than H 1 . In the critical case, Ref. [8] established scattering for general data with d 5 . Scattering in the focusing case was achieved in [9,10] for small and radial data. Further references for the NLS in a general setting include [11,12]. A principal tool in studying the dynamics of solutions to (1) is the Morawetz multiplier technique and its associated estimates. In our recent work [13], we developed a method combining Morawetz inequalities, a localization step, and interpolation with a contradiction argument to achieve the decay of solutions for the Schrödinger–Hartree equation. This robust property is crucial in scattering theory, as highlighted in [12,13,14]. Motivated by these developments, we present a generalization of this method for the damped magnetic Schrödinger equation with Hartree-type nonlinearity. The linearly damped nonlinear Schrödinger equation plays a significant role across multiple scientific disciplines, including nonlinear optics, plasma physics, and fluid mechanics. This equation is fundamental for understanding various complex phenomena, such as the propagation of optical pulses in nonlinear media, the behavior of plasma waves in magnetized environments, and the dynamics of fluid flows under certain conditions. We quote here, for example, [15,16]. Our result is novel in the literature, and we make minimal assumptions on the magnetic function A ( x ) . Furthermore, our strategy simplifies and extends the damped magnetic Schrödinger equation to the approach used in [17,18,19,20]. We emphasize also that the approaches previously proposed, for instance, in [21,22,23] (see also references therein), are outperformed since we coped with the complex-valued function b ( t ) in (1).

2. Preliminaries

Before outlining our main achievements, we will unveil some necessary notations and several useful results. For any two positive real numbers a , b , we write a b (resp. a b ) to denote a C b (resp. C a b ), with C > 0 , and we unravel the constant only when it is necessary. We introduce the Banach space L r ( R d ) = L x r for 1 r . In addition, we introduce
H 1 , r ( R d ) = 1 Δ x 1 2 L r ( R d ) , H 1 , r ( R d ) = H x 1 , r ,
and denote it with H 1 , 2 ( R d ) = H 1 ( R d ) = H x 1 . Given any Banach space X, we define
f L t X = ess sup t R f ( x ) X .
We adopt the notation L T X when one restricts t [ 0 , T ) , for T > 0 . The following results are also useful (see [9,17,19], respectively).
Lemma 1. 
Let f be a radial function in H x 1 . Then,
| x | d 1 2 f L x 2 f L x 2 x f L x 2 .
Proposition 1. 
Let A be as in (2) and (3). For any 1 < r < d , one obtains
( Δ x A ) 1 2 f L x r ( Δ x ) 1 2 f L x r
and
( Δ x ) 1 2 f L x 2 ( Δ x A ) 1 2 f L x 2 .
We also have the following maximal estimate (see, for example, [24]), as a straightforward consequence of the Hardy inequality.
Proposition 2. 
Let 0 < γ < d . We have
[ | · | d γ | u | 2 ] L x C ( d , γ ) u H ˙ x d γ 2 2 .
We recall also that the solutions to (1) satisfy the conservation laws. We summarize them in the following.
Proposition 3. 
Let d 1 . Then, a sufficiently smooth solution to (1) satisfies the following identities:
u ( t ) L x 2 = e B ( t ) f L x 2 , H ( u ( t ) ) = H ( f ) ,
where
H ( u ( t ) ) = e 2 B ( t ) R d | x A u ( t , x ) | 2 d x + k e 2 B ( t ) R d R d | u ( t , x ) | 2 | u ( t , y ) | 2 | x y | d γ d x d y + 2 k 0 t R d R d b ( s ) e 2 B ( t ) | u ( t , x ) | 2 | u ( t , y ) | 2 | x y | d γ d x d y .
Proof. 
We utilize the change in variable
v ( t , x ) : = e B ( t ) u ( t , x )
and see that u satisfies (1) if v solves
i t v + Δ x A v = k e 2 B ( t ) [ | · | d + γ | v | 2 ] v , ( t , x ) [ 0 , ) × R d , v ( 0 , x ) = u 0 ( x ) .
We multiply the above equation by u ¯ ( t , x ) , integrate with respect to the x-variable, and take the imaginary part, obtaining the following, since div A = 0 :
i 1 2 d d t R d | v ( t , x ) | 2 d x + R d v ¯ ( t , x ) Δ u ( t , x ) + | A | 2 u ( t , x ) 2 i A · v ( t , x ) d x + R d k e 2 B ( t ) [ | x | d + γ | v ( t , x ) | 2 ] | v ( t , x ) | 2 d x = i 1 2 d d t R d | v ( t , x ) | 2 d x 2 i R 2 A · | v ( t , x ) | 2 d x = 0 .
Thus, solutions local in time satisfy the conservation of mass
v ( t ) L 2 2 = f L 2 2 .
that is, the first identity in (11). We multiply now Equation (14) by u ¯ ( t , x ) , integrate with respect to the x-variable, and take the imaginary equation part. We have
R d A v ( t , x ) A t v ¯ ( t , x ) + k e 2 B ( t ) [ | x | d + γ | v ( t , x ) | 2 ] v ( t , x ) t v ¯ ( t , x ) d x = 0 .
The previous identity is enhanced to
R d 1 2 t | A v ( t , x ) | 2 + 1 2 k e 2 B ( t ) [ | x | d + γ | v ( t , x ) | 2 ] t | v ( t , x ) | 2 d x = 0
and then to
t R d 1 2 | A v ( t , x ) | 2 + 1 2 k e 2 B ( t ) [ | x | d + γ | v ( t , x ) | 2 ] | v ( t , x ) | 2 d x = k R d b ( t ) e 2 B ( t ) [ | x | d + γ | v ( t , x ) | 2 ] | v ( t , x ) | 2 d x .
Integrating with respect to the t-variable identity (15), we obtain
R d | A v ( t , x ) | 2 + k e 2 B ( t ) [ | x | d + γ | v ( t , x ) | 2 ] | v ( t , x ) | 2 d x + 2 k 0 t R d b ( s ) e 2 B ( s ) [ | x | d + γ | v ( s , x ) | 2 ] | v ( s , x ) | 2 d x d s = R d | A v ( 0 , x ) | 2 + k [ | x | d + γ | v ( 0 , x ) | 2 ] | v ( 0 , x ) | 2 d x .
The above relation suggests that the quantity
H ˜ ( v ( t ) ) = R d | A v ( t , x ) | 2 + k e 2 B ( t ) [ | x | d + γ | v ( t , x ) | 2 ] | v ( t , x ) | 2 d x + 2 k 0 t R d b ( s ) e 2 B ( s ) [ | x | d + γ | v ( s , x ) | 2 ] | v ( s , x ) | 2 d x d s
is conserved. Hence, this implies the local conservation of the Hamiltonian in (11) with H ( u ( t ) ) as in (12). □

3. Well-Posedness

Here, we present the following existence and uniqueness result, which is crucial for the proof of (6). Specifically, we prove the following proposition.
Proposition 4. 
Let d 3 . Assume that (2)–(5) are satisfied. Then, for all f H x 1 , there exists T > 0 such that problem (1) has a unique local solution u C ( [ 0 , T ) ; H x 1 ) with
e B ( t ) u ( t , x ) u L T H x 1 f H x 1 .
Moreover, the solution can be extended globally in time if k > 0 .
Proof. 
We shall accomplish a fixed-point argument. Namely, consider the integral operator associated with (14) to be defined for all f H x 1 as
T f ( e B ( t ) u ) = e i t Δ x A + B ( t ) f + k 0 t e 2 B ( t ) e i ( t τ ) Δ x A [ | · | ( d γ ) | e B ( t ) u | 2 ] e B ( t ) u ( τ ) d τ .
We need to show that it is possible to find a T = T f H x 1 > 0 and a unique
e B ( t ) u ( t , x ) L T H x 1
satisfying the property
T f ( e B ( t ) u ( t ) ) = e B ( t ) u ( t ) ,
for t [ 0 , T ) . For the sake of simplicity, we will divide the proof into different steps.
  • Step One: For any e B ( t ) u H x 1 , there exist T = T f H x 1 > 0 and R = R f H x 1 > 0 such that
    T f ( B L T H x 1 ( 0 , R ) ) B L T H x 1 ( 0 , R ) ,
    for any T < T .
By the classical Hardy–Littlewood–Sobolev inequality combined with (8) and (9), we have
T f ( e B ( t ) u ) L T L x 2 + x T f ( e B ( t ) u ) L T L x 2 T f ( e B ( t ) u ) L T L x 2 + x A T f ( e B ( t ) u ) L T L x 2 f H x 1 + 0 T e 2 B ( · ) [ | · | d + γ | e B ( · ) u | 2 ] e B ( · ) u L T H x 1 d τ .
At this point, by condition (5), the last term in the above chain of inequalities can be controlled as follows:
f H x 1 + 0 T e 2 B ( · ) [ | · | d γ | e B ( · ) u | 2 ] e B ( · ) u L T H x 1 d τ f H x 1 + T e 2 B ( · ) L t [ | · | d + γ | e B ( · ) u | 2 ] e B ( · ) u L T H x 1 f H x 1 + T [ | · | d + γ | e B ( · ) u | 2 ] L T L x e B ( · ) u L T H x 1 + T [ | · | d + γ | e B ( · ) u | 2 ] L T H x 1 , 2 d d γ e B ( · ) u L T L x 2 d γ f H x 1 + T e B ( · ) u L T H ˙ x d γ 2 2 e B ( · ) u L T H x 1 + T | e B ( · ) u | 2 L T H x 1 , 2 d d + γ e B ( · ) u L T L x 2 d γ f H x 1 + T e B ( · ) u L T H ˙ x d γ 2 2 e B ( · ) u L T H x 1 + T e B ( · ) u L T H x 1 e B ( · ) u L T L x 2 d γ 2 f H x 1 + T e B ( · ) u L T H ˙ x d γ 2 2 e B ( · ) u L T H x 1 f H x 1 + T R 3 .
By selecting R and T so that
2 f H x 1 = R , 2 C T R 2 1 ,
we finish the proof of this step.
  • Step Two: Let T , R > 0 be as in the above step. Then, there exists T ¯ = T ¯ f H x 1 < T such that T f is a contraction on B L T ¯ H x 1 ( 0 , R ) , equipped with the norm . L T ¯ L x 2 .
Consider e B ( · ) v 1 , e B ( · ) v 2 B L T H x 1 ( 0 , R ) . We obtain, by arguing as in the previous lines, the following chain of inequalities:
T f e B ( · ) v 1 T f e B ( · ) v 2 L T L x 2 T e 2 B ( · ) L t [ | · | d + γ | e B ( · ) v 1 | 2 ] e B ( · ) v 1 [ | · | d + γ | e B ( · ) v 2 | 2 ] e B ( · ) v 2 L T L x 2 T [ | · | d + γ | e B ( · ) v 1 | 2 ] ( e B ( · ) v 1 e B ( · ) v 2 ) L T L x 2 + T [ | · | d + γ ( | e B ( · ) v 1 | 2 | e B ( · ) v 2 | 2 ) ] e B ( · ) v 2 L T L x 2 T e B ( · ) v 1 L T H x d + γ 2 2 e B ( · ) v 1 e B ( · ) v 2 L T L x 2 + T [ | · | d + γ ( | e B ( · ) v 1 | 2 | e B ( · ) v 2 | 2 ) ] L T L x 2 d d γ e B ( · ) v 2 L T L x 2 d γ T R 2 e B ( · ) v 1 e B ( · ) v 2 L T L x 2 + R v 1 + v 2 L T L x 2 d γ e B ( · ) v 1 e B ( · ) v 2 L T L x 2 T R 2 v 1 v 2 L T L x 2 .
Then, we arrive at
T f e B ( · ) v 1 T f e B ( · ) v 2 L T L x 2 T R 2 e B ( · ) v 1 e B ( · ) v 2 L T L x 2 .
This inequality allows us to say that T f is a contraction on B L T H x 1 ( 0 , R ) if T is chosen in a suitable manner.
  • Step Three: The solution exists and is unique in L T ¯ H x 1 , where T ¯ is as in the above step.
We can exhibit the existence and uniqueness of the solution using the contraction principle for the map T f defined on the complete metric space B L T ¯ H x 1 ( 0 , R ) , endowed with the topology induced by . L T ¯ L x 2 .
  • Step Four: The solution can be extended globally.
We obtain, by conservation laws (11), (8), and (9), that
e B ( t ) u H x 1 e B ( t ) u L x 2 + x A e B ( t ) u L x 2 H ( u ( 0 ) ) + f L x 2 .
The previous bound leads to the global well-posedness for (1). □

4. Morawetz Identities and Inequalities

Our first contribution is the Morawetz equalities associated to (1). They are presented in the following.
Lemma 2. 
Let d 1 and u C ( [ 0 , ) ; H x 1 ) denote a global solution to (1) with radial initial data f H x 1 such that (2) and the strict inequality in (3)–(5) are satisfied. Moreover, let ψ = ψ ( x ) : R d R be a sufficiently regular and decaying function, denoted by
V ( t ) : = R d ψ ( x ) | e B ( t ) u ( t , x ) | 2 d x .
Then, the following identities hold:
V ˙ ( t ) = 2 R d e 2 B ( t ) u ¯ ( t , x ) x ψ ( x ) · x A u ( t , x ) d x
and
V ¨ ( t ) = R d Δ x 2 ψ ( x ) | e B ( t ) u ( t , x ) | 2 d x + 4 R d e 2 B ( t ) x A u ( t , x ) D x 2 ψ ( x ) · x A u ( t , x ) ¯ d x 4 R d e 2 B ( t ) u ( t , x ) x ψ ( x ) · B ( x ) x A u ( t , x ) ¯ d x 2 k R d e 2 B ( t ) x ψ ( x ) · x | x | ( d γ ) | e B ( t ) u ( t , x ) | 2 | e B ( t ) u ( t , x ) | 2 d x ,
where D x 2 ψ M d × d ( R ) is the Hessian matrix of ψ, and Δ x 2 ψ = Δ x ( Δ x ψ ) , the bi-Laplacian operator.
Proof. 
We will prove the identities for a smooth, rapidly decreasing solution u = u ( t , x ) , recovering the general case e B ( t ) u C ( R ; H x 1 ) via a density argument. The proof of (18) is similar to the one given in [20], since we can use transformation (13) and then Equation (14). We present details for obtaining (19). We have the following identity for the linear terms, using Theorem 1.2 in [20] and the v ( t , x ) defined in (13) and (14) again. We obtain
2 R d Δ x A v ( t , x ) Δ x ψ ( x ) v ¯ ( t , x ) + 2 x ψ ( x ) · x A v ( t , x ) ¯ d x = R d Δ x 2 ψ ( x ) | v ( t , x ) | 2 d x 4 R d v ( t , x ) x ψ ( x ) · B ( x ) x A v ( t , x ) ¯ d x + 4 R d x A v ( t , x ) D x 2 ψ ( x ) x A v ( t , x ) ¯ d x .
In addition, for the nonlinear terms, one has
2 R d e 2 B ( t ) [ | x | ( d γ ) | | 2 ] | v ( t , x ) | 2 Δ x ψ ( x ) d x + 4 R d e 2 B ( t ) [ | x | ( d γ ) | v ( t , x ) | 2 ] v ( t , x ) x ψ ( x ) · x A v ( t , x ) ¯ d x = 2 R d e 2 B ( t ) [ | x | ( d γ ) | v ( t , x ) | 2 ] | v ( t , x ) | 2 Δ x ψ ( x ) d x + 4 R d e 2 B ( t ) [ | x | ( d γ ) | v ( t , x ) | 2 ] v ( t , x ) x ψ ( x ) · x v ¯ ( t , x ) d x .
The last term of the above identity is equal to
2 R d e 2 B ( t ) [ | x | ( d γ ) | v ( t , x ) | 2 ] | v ( t , x ) | 2 Δ x ψ ( x ) d x + 2 R d e 2 B ( t ) [ | x | ( d γ ) | v ( t , x ) | 2 ] x ψ ( x ) · x | v ( t , x ) | 2 d x .
Then, through integration by parts of the second term in the last line above, one arrives at
2 R d e 2 B ( t ) [ | x | ( d γ ) | v ( t , x ) | 2 ] | v ( t , x ) | 2 Δ x ψ ( x ) d x + 4 R d e 2 B ( t ) [ | x | ( d γ ) | v ( t , x ) | 2 ] v ( t , x ) x ψ ( x ) · x A v ( t , x ) ¯ d x = 2 R d e 2 B ( t ) x ψ ( x ) · x | x | ( d γ ) | v ( t , x ) | 2 | v ( t , x ) ) | 2 d x .
Combining now identities (20) and (21) and turning back to e B ( t ) u ( t , x ) , we obtain (19). □

A Localized Morawetz Inequality

We start this section with a result that is a consequence of Lemma 2. More precisely, we have the following lemma
Lemma 3. 
Assume d 4 and let u C ( [ 0 , ) ; H x 1 ) be a global solution to (1) with radial initial data f H x 1 such that (2) and the strict inequality in (3)–(5) are satisfied. Then, it holds that
R d 1 | x | 3 | e B ( t ) u ( t , x ) | 2 d x R d e 2 B ( t ) u ¯ ( t , x ) x ψ ( x ) · x A u ( t , x ) d x .
Proof. 
We pick ψ = ψ ( x ) = | x | . This gives
x ψ = x | x | , Δ x ψ = d 1 | x | , Δ x 2 ψ = ( d 1 ) ( d 3 ) | x | 3 ,
if d 4 . A change in variable (13), Equation (14), and an application of identity (19) allow us to write the following:
2 t R d v ¯ ( t , x ) x ψ ( x ) · x A v ( t , x ) d x = 2 R d R d Δ x 2 ψ ( x ) | v ( t , x ) | 2 d x 4 R d v ( t , x ) x ψ ( x ) · B ( x ) x A v ( t , x ) ¯ d x + 4 R d x A v ( t , x ) D x 2 ψ ( x ) x A v ( t , x ) ¯ d x + k ( d γ ) e 2 B ( t ) R d R d 1 | x z | d γ + 2 | v ( t , x ) | 2 | v ( t , z ) | 2 K ( x , z ) d x d z ,
with
K ( x , z ) = ( x z ) · x | x | z | z | .
By the elementary inequality
( x z ) · x | x | z | z | = | x | | z | ( x ) · ( z ) | x | + | z | | x | | z | 0 ,
we have that K ( x , z ) 0 . Therefore, one can drop the last term on the right-hand side of (24). We shall focus now on the linear terms in (24), following the method utilized in [20]. Observe that the relations (23) relate to
x A v ( t , x ) D 2 ψ ( x ) x A v ( t , x ) ¯ = A τ v ( t , x ) 2 | x | ,
(see identity (3.9) in [20]) where the operator A τ is defined as
A τ v ( t , x ) = x A v ( t , x ) x A v ( t , x ) · x | x | x | x | .
Therefore, utilizing (23), we have the following identity:
2 R d R d Δ x 2 ψ ( x ) | v ( t , x ) | 2 d x 4 R d v ( t , x ) x ψ ( x ) · B ( x ) x A v ( t , x ) ¯ d x + 4 R d x A v ( t , x ) D x 2 ψ ( x ) x A v ( t , x ) ¯ d x = 4 R d A τ v ( t , x ) 2 | x | d x + ( d 1 ) ( d 3 ) R d | v ( t , x ) | 2 | x | 3 d x + 4 R d v ( t , x ) x | x | B ( x ) · x A v ( t , x ) ¯ d x .
The last term of the identity above can be estimated as
R d v ( t , x ) x | x | B ( x ) · x A v ( t , x ) ¯ d x R d | v ( t , x ) | 2 | x | 3 d x 1 2 R d | x | 2 x B ( x ) 2 A τ v ( t , x ) 2 d x 1 2 C * R d | v ( t , x ) | 2 | x | 3 d x 1 2 R d A τ v ( t , x ) 2 | x | d x 1 2 ,
where
C * 2 = ( d 1 ) ( d 3 ) .
As a result, the right-hand side of (26) can be bounded as
4 R d A τ v ( t , x ) 2 | x | d x + ( d 1 ) ( d 3 ) R d | v ( t , x ) | 2 | x | 3 d x + 4 R d v ( t , x ) x | x | B ( x ) · x A v ( t , x ) ¯ d x 4 R d A τ v ( t , x ) 2 | x | d x + ( d 1 ) ( d 3 ) R d | v ( t , x ) | 2 | x | 3 d x 4 C ˜ R d | v ( t , x ) | 2 | x | 3 d x 1 2 R d A τ v ( t , x ) 2 | x | d x 1 2 > 0 .
Notice also that the previous inequality and a continuity argument guarantee that
4 R d A τ v ( t , x ) 2 | x | d x + ( d 1 ) ( d 3 ) R d | v ( t , x ) | 2 | x | 3 d x + 4 R d v ( t , x ) x | x | B ( x ) · x A v ( t , x ) ¯ d x > η ˜ ( d 1 ) ( d 3 ) R d | v ( t , x ) | 2 | x | 3 d x ,
for a η ˜ > 0 . The above bound in combination with (24), the fact that K ( x , z ) 0 , and (26) give the proof of (22). □
We have the following corollary, which is a consequence of (22).
Corollary 1. 
Let u C ( [ 0 , ) ; H x 1 ) be a global solution to (1) with radial initial data f H x 1 such that (2) and the strict inequality in (3)–(5) are satisfied. Moreover, let Q x ˜ d ( r ) = x ˜ + [ r , r ] d , with r > 0 and x ˜ R d . Hence, one obtains
0 Q x ˜ d ( r ) 1 | x | 3 | e B ( t ) u ( t , x ) | 2 d x d t < .
Proof. 
By integrating (22) with ψ ( x ) as in (23) with respect to the time variable on the interval J = [ t 1 , t 2 ] , with t 1 , t 2 [ 0 , ) , one arrives at
R d e 2 B ( t ) u ¯ ( t , x ) x ψ ( x ) · x A u ( t , x ) d x t = t 1 t = t 2 t 1 t 2 R d 1 | x | 3 | e B ( t ) u ( t , x ) | 2 d x d t t 1 t 2 Q x ˜ d ( r ) 1 | x | 3 | e B ( t ) u ( t , x ) | 2 d x d y d t .
Applying the Cauchy–Schwartz inequality and Proposition 1, we also infer that
2 R d e 2 B ( t ) u ¯ ( t , x ) x ψ ( x ) · x A u ( t , x ) d x t = t 1 t = t 2 f H x 1 2 < ,
since the H x 1 -norm is a quantity conserved by (17). Finally, we obtain (30) when t 1 = 0 , t 2 . □

5. The Decay of Solutions

This section is devoted to demonstrating the main Theorem 1.
Proof. 
It is sufficient to prove property (6) for a suitable 2 < q < 2 d d 2 because the thesis for the general case follows conservation law (11) and interpolation. More precisely, it is enough to show that
lim t ± e B ( t ) u ( t , x ) L x 2 + d = 0 .
Then, property (6) follows for all 2 < q < 2 d d 2 by combining (32) with
sup t R e B ( t ) u ( t , x ) H x 1 < .
We recall the following localized Gagliardo–Nirenberg inequality (see [13]):
ζ L x 2 d + 4 d 2 d + 4 d C sup x R d ζ L 2 ( Q x ( 1 ) ) 4 d ζ H x 1 2 ,
where Q x d ( r ) = x + [ 1 , 1 ] d . Next, assume by contradiction that (32) is not true; then, by (33) and by (34), we deduce the existence of a sequence ( t n , x n ) R × R d with t n and ϵ 0 > 0 such that
inf n e B ( t n ) u ( t n , x ) L 2 ( Q x n ( 1 ) ) 2 = ϵ 0 2 .
Notice that by (18) in conjunction with (33), we obtain
sup n , t d d t ϕ ( x x n ) | e B ( t ) u ( t , x ) | 2 d x < ,
where χ ( x ) is a smooth and non-negative cut-off function, such that ϕ ( x ) = 1 for x Q 0 ( 1 ) = [ 1 , 1 ] d and ϕ ( x ) = 0 for x Q 0 ( 2 ) = [ 2 , 2 ] d . Consequently, by the Fundamental Theorem of calculus, we deduce the inequality
R d ϕ ( x x n ) | e B ( σ ) u ( σ , x ) | 2 d x R d ϕ ( x x n ) | e B ( t ) u ( t , x ) | 2 d x C ˜ | t σ | ,
for a C ˜ > 0 that does not depend on n. By choosing t = t n , we have
R d ϕ ( x x n ) | e B ( σ ) u ( σ , x ) | 2 d x R d ϕ ( x x n ) | e B ( t n ) u ( t n , x ) | 2 d x C ˜ | t n σ | ,
which implies the following, considering the support property of function ϕ :
Q x n d ( 2 ) | e B ( σ ) u ( σ , x ) | 2 d x Q x n d ( 1 ) | e B ( t n ) u ( t n , x ) | 2 d x C ˜ | t n σ | ,
for a C ˜ > 0 independent form n. By combining this fact with (35), we have the existence of T > 0 such that
inf n inf t ( t n , t n + T ) e B ( t ) u ( t , x ) L 2 ( Q x n ( 2 ) ) 2 ϵ 1 2 ,
for some ϵ 1 > 0 . Notice that the previous estimate (39) provides, in combination with the Strauss radial inequality (7), that the sequence of centers ( x n ) n N is uniformly bounded. Observe also that since t n , we can assume that the intervals (eventually passing to a subsequence) ( t n , t n + T ) are disjoint. In particular, we acquire the following for d 4 :
n T ϵ 1 2 n t n t n + T Q x n d ( 2 ) | e B ( t ) u ( t , x ) | 2 d x d t
0 sup x ˜ R d Q x ˜ d ( 2 ) 1 | x | 3 | e B ( t ) u ( t , x ) | 2 d x d t ,
So, we obtain a contradiction because the right-hand side of the above (40) is bounded by (30). □
Remark 1. 
Note that (10) and (12) introduce
e B ( t ) u ( t ) H 1 f H 1 + 0 t e B ( τ u ( τ ) H ˙ γ 2 2 e B ( τ ) u ( τ ) H 1 d τ e B ( t ) u ( t ) H 1 + C 0 t e B ( τ u ( τ ) H 1 d τ .
Then, by Gronwall’s inequality, we have
e B ( t ) u ( t ) H x 1 f H x 1 e K t ,
with K > 0 depending on f L x 2 and H ( u ( 0 ) ) . Then, the Sobolev embedding and interpolation with the conservation of mass in (11) lead to
e B ( t ) u ( t ) L x r e ε K t ,
with 0 < ε < 1 and 2 < r 2 d d 2 , which is not sufficient to guarantee a behavior like the one disclosed by (6) in Theorem 1 upon letting t + .
Remark 2. 
We also highlight that one can achieve an exponential decay just by interpolating the conservation of mass in (11) and the estimate arising from the Sobolev embedding and (17). Namely, one has
e B ( t ) u ( t ) L x r 1 .
However, it is not enough to ensure a behavior such as in (6) in Theorem 1, which is a stronger property of the solutions to (1). Moreover, the case B ( t ) = 0 , that is, when i b ( t ) is a real function, cannot be included in the previous analysis.
Remark 3. 
It is important to notice that our results can be used to deal with a class of damped nonlinearities fulfilling (5), particularly when
b ( t ) a ˜ 1 + t , for t > 1 , a ˜ > 0 ,
as considered in [25]. Also, we can address more general damping terms, leading to equations of the form
i t u + Δ x A u + i a ˜ ( 1 + t ) α u + b ˜ ( t ) k [ | · | ( d γ ) | u | 2 ] u = 0 .
with α 0 and where b ˜ ( t ) C ( [ 0 , ) ) is a real-valued function (see [22]).

6. Conclusions

We broaden the outcomes achieved in [17,18,19,20] to the damped scenario. The assumptions formulated for the time-depending function i b ( t ) are more general than the ones found for the example in [21,22,23]. This is because we include an oscillatory part in the perturbed propagator e B ( t ) + i t Δ x A , which can not be treated if one uses the techniques developed in the aforementioned works. We underline that assumptions (2) and (3) related to the operator x A and the function A ( x ) are less restrictive than those imposed in [17,19,20]. This is because our well-posedness analysis relies solely on the energy estimate for (1). We are not using any Strichartz estimates here, forcing the constraint d 2 γ < d . The obstacle here is that the multipliers utilized in the papers previously mentioned are not well suited to handle a non-local nonlinearity because their application cannot guarantee the non-negativity of the last term in (24). This aspect also determines the radial assumption we made on the initial data. A second problem is the equivalence of the norm result of Proposition 1, valid only in the L 2 framework. We are confident that we will overcome all these issues and shed light on the lower regularity frame d 4 γ < d 2 in a future paper. We are also confident that our decay result can greatly simplify the scattering theory associated with (1), as well as shed light on the case B ( t ) < 0 .

7. Open Problems and Further Developments

The theory established in this paper is general and allows us to obtain, in a straightforward manner, the decay in the energy spaces of the solutions to the damped magnetic Schrödinger equation with non-local nonlinearity. We believe that it can be used for the following open problems:
  • The analysis of the scattering in the energy space for the solution to (1);
  • The investigation of the decay properties (and eventually scattering) for the solutions to the generalized Schrödinger–Hartree equation, that is
    i t u + Δ x A u + i b ( t ) u + u [ | · | ( d γ ) | u | p ] | u | p 2 u = 0 ,
    where the positive nonlinear parameter is either p = 2 or satisfies d + γ + 2 d < p < d + γ d 2 ;
  • The exploration of the decaying and scattering properties of the solutions on other nonlinear dispersive equations such as the nonlinear Beam Equation
    t t u + ( Δ x A ) 2 u + i b ( t ) u + u [ | · | ( d γ ) | u | p ] | u | p 2 u = 0 ,
    where ( Δ x A ) 2 = Δ x A ( Δ x A ) is the magnetic bi-Laplacian operator, or the nonlinear Klein–Gordon equation
    t t u Δ x A u + i b ( t ) u + u + [ | · | ( d γ ) | u | p ] | u | p 2 u = 0 ,
    with the nonlinear parameter p defined as above, including the special case when A ( x ) = 0 .

Author Contributions

Conceptualization, T.S., M.T., and G.V.; methodology, T.S., M.T., and G.V.; formal analysis, T.S., M.T., and G.V.; investigation, T.S., M.T., and G.V.; writing—original draft preparation, T.S., M.T., and G.V.; writing—review and editing, T.S., M.T., and G.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Elgart, A.; Schlein, B. Mean field dynamics of boson stars. Comm. Pure Appl. Math. 2007, 60, 500–545. [Google Scholar] [CrossRef]
  2. Lenzmann, E. Well-posedness for semi-relativistic Hartree equations of critical type. Math. Phys. Anal. Geom. 2007, 10, 43–64. [Google Scholar] [CrossRef]
  3. Lewin, M.; Rougerie, N. Derivation of Pekar’s polarons from a microscopic model of quantum crystal. SIAM J. Math. Anal. 2013, 45, 1267–1301. [Google Scholar] [CrossRef]
  4. Ginibre, J.; Velo, G. Quadratic Morawetz inequalities and asymptotic completeness in the energy space for nonlinear Schrödinger and Hartree equations. Quart. Appl. Math. 2010, 68, 113–134. [Google Scholar] [CrossRef]
  5. Nakanishi, K. Energy scattering for Hartree equations. Math. Res. Lett. 1999, 6, 107–118. [Google Scholar] [CrossRef]
  6. Ginibre, J.; Ozawa, T. Long range scattering for nonlinear Schrödinger and Hartree equations in space dimension n ≥ 2. Comm. Math. Phys. 1993, 151, 619–645. [Google Scholar] [CrossRef]
  7. Nawa, H.; Ozawa, T. Nonlinear scattering with nonlocal interactions. Comm. Math. Phys. 1992, 146, 259–275. [Google Scholar] [CrossRef]
  8. Miao, C.; Xu, G.; Zhao, L. Global well-posedness and scattering for the energy-critical, defocusing Hartree equation for radial data. J. Funct. Anal. 2007, 252, 605–627. [Google Scholar] [CrossRef]
  9. Arora, A. Scattering of radial data in the focusing NLS and generalized Hartree equation. Discrete Contin. Dyn. Syst. 2019, 39, 6643–6668. [Google Scholar] [CrossRef]
  10. Arora, A.; Roudenko, S. Global behavior of solutions to the focusing generalized Hartree equation. Michigan Math. J. 2021, 71, 619–672. [Google Scholar] [CrossRef]
  11. Tarulli, M. H2-scattering for Systems of Weakly Coupled Fourth-order NLS Equations in Low Space Dimensions. Potential Anal. 2019, 51, 291–313. [Google Scholar] [CrossRef]
  12. Tarulli, M.; Venkov, G. Decay in energy space for the solution of fourth-order Hartree-Fock equations with general non-local interactions. J. Math. Anal. Appl. 2022, 516, 126533. [Google Scholar] [CrossRef]
  13. Tarulli, M.; Venkov, G. Decay and scattering in energy space for the solution of weakly coupled Schrödinger-Choquard and Hartree-Fock equations. J. Evol. Equ. 2021, 21, 1149–1178. [Google Scholar] [CrossRef]
  14. Cazenave, T. Semilinear Schrödinger Equations; Courant Lecture Notes in Mathematics, 10; New York University Courant Institute of Mathematical Sciences: New York, NY, USA, 2003. [Google Scholar]
  15. Chen, G.; Zhang, J.; Wei, Y. A small initial data criterion of global existence for the damped nonlinear Schrödinger equation. J. Phys. A Math. Theor. 2009, 42, 055205. [Google Scholar] [CrossRef]
  16. Goldman, M.V.; Rypdal, K.; Hafizi, B. Dimensionality and dissipation in Langmuir collapse. Phys. Fluids 1980, 23, 945–955. [Google Scholar] [CrossRef]
  17. Colliander, J.; Czubak, M.; Lee, J.J. Interaction Morawetz estimate for the magnetic Schrödinger equation and applications. Adv. Differ. Equ. 2014, 1, 805–832. [Google Scholar] [CrossRef]
  18. Nikolova, E.; Tarulli, M.; Venkov, G. On the magnetic radial Schrödinger-Hartree equation. Int. J. Appl. Math. 2022, 35, 795–809. [Google Scholar] [CrossRef]
  19. D’Ancona, P.; Fanelli, L.; Vega, L.; Visciglia, N. Endpoint Strichartz estimates for the magnetic Schrödinger equation. J. Funct. Anal. 2020, 258, 3227–3240. [Google Scholar] [CrossRef]
  20. Fanelli, L.; Vega, L. Magnetic virial identities, weak dispersion and Strichartz inequalities. Math. Ann. 2009, 344, 249–278. [Google Scholar] [CrossRef]
  21. Dinh, V.D. Blow-up criteria for linearly damped nonlinear Schrödinger equations. Evol. Equ. Control Theory 2021, 10, 599–617. [Google Scholar] [CrossRef]
  22. Hamouda, M.; Majdoub, M. Energy scattering for the unsteady damped nonlinear Schrödinger equation. arXiv 2024, arXiv:2311.14980. [Google Scholar]
  23. Inui, T. Asymptotic behavior of the nonlinear damped Schrödinger equation. Proc. Amer. Math. Soc. 2019, 147, 763–773. [Google Scholar] [CrossRef]
  24. Kato, T. Pertubation Theory for Linear Operators, 2; Springer: Berlin, Germany, 1980. [Google Scholar]
  25. Bamri, C.; Tayachi, S. Global existence and scattering for nonlinear Schrödinger equations with time-dependent damping. Commun. Pure Appl. Anal. 2023, 22, 2365–2399. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Saker, T.; Tarulli, M.; Venkov, G. On the Decay in the Energy Space of Solutions to the Damped Magnetic Radial Schrödinger Equation with Non-Local Nonlinearities. Mathematics 2024, 12, 2975. https://doi.org/10.3390/math12192975

AMA Style

Saker T, Tarulli M, Venkov G. On the Decay in the Energy Space of Solutions to the Damped Magnetic Radial Schrödinger Equation with Non-Local Nonlinearities. Mathematics. 2024; 12(19):2975. https://doi.org/10.3390/math12192975

Chicago/Turabian Style

Saker, Taim, Mirko Tarulli, and George Venkov. 2024. "On the Decay in the Energy Space of Solutions to the Damped Magnetic Radial Schrödinger Equation with Non-Local Nonlinearities" Mathematics 12, no. 19: 2975. https://doi.org/10.3390/math12192975

APA Style

Saker, T., Tarulli, M., & Venkov, G. (2024). On the Decay in the Energy Space of Solutions to the Damped Magnetic Radial Schrödinger Equation with Non-Local Nonlinearities. Mathematics, 12(19), 2975. https://doi.org/10.3390/math12192975

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