The Existence and Stability of Standing Waves for the Inhomogeneous Nonlinear Schrödinger Equation with Magnetic Potential

Abstract

In this paper, we study the Cauchy problem for the three dimensional inhomogeneous nonlinear Schrödinger equation (INLS) in the presence of a constant magnetic field. We establish the existence and orbital stability of normalized standing waves. Our approach introduces a novel methodology that circumvents the widely utilized concentration-compactness principle.

Share and Cite:

Yan, P. (2026) The Existence and Stability of Standing Waves for the Inhomogeneous Nonlinear Schrödinger Equation with Magnetic Potential. Journal of Applied Mathematics and Physics, 14, 1525-1534. doi: 10.4236/jamp.2026.144072.

1. Introduction

This paper is concerned with the Cauchy problem for inhomogeneous nonlinear Schrödinger equations with a constant magnetic field in three dimensions

{ i t ψ+ ( +iA ) 2 ψ= | x | b | ψ | p ψ, ψ| t=0 = ψ 0 , ( t,x ) + × 3 , (1.1)

where A is a vector-valued potential of the form

A( x )= a 2 ( x 2 , x 1 ,0 ),x=( x 1 , x 2 , x 3 ) 3 (1.2)

modeling the effect of an external magnetic field

curl( A )=( 0,0,a ),a0. (1.3)

As a fundamental equation of quantum mechanics, the time-independent Schrödinger equation Schrödinger [1] successfully describes the quantum evolution of particles without external fields or interparticle interactions. Practical systems are often driven by external sources, and homogeneous equations fail to describe the evolutionary behavior under such external driving. Therefore, it is necessary to investigate the inhomogeneous nonlinear Schrödinger equation with a magnetic potential (see [2]-[5] and references therein). The main objective of this paper is to study the existence and stability of solutions to the inhomogeneous nonlinear Schrödinger equation with a constant magnetic potential (1.1). Moreover, the Schrödinger equation with a uniform magnetic field serves as an effective model for describing the behavior of a single non-relativistic quantum particle subjected to an electromagnetic field (see, for example, [6] [7]). Avron, Herbst, and Simon [8]-[10] conducted a thorough mathematical analysis of the linear Schrödinger operator in the presence of a constant magnetic field.

Furthermore, there are conservation laws of mass and energy, namely

M( ψ( t ) )= ψ( t ) L 2 2 =M( ψ 0 ), (Mass)

and

E( ψ( t ) )= 1 2 ( +iA )ψ( t ) L 2 2 1 p+2 3 | x | b | ψ | p+2 dx =E( ψ 0 ), (Energy)

for all t[ 0, T * ) . Moreover, we have the following useful identity (see e.g., [11], Lemma 2.2)

( +iA )u L 2 2 = u L 2 2 +aR( u )+ a 2 4 ρu L 2 2 , (1.4)

where ρ:= x 1 2 + x 2 2 and

R( u ):=i ( x 2 x 1 u x 1 x 2 u ) u ¯ = L z u u ¯ dx , (1.5)

where L z :=i( x 2 x 1 x 1 x 2 ) . Next, we focus on the existence and stability of standing waves with a specified mass for Equation (1.1). By standing waves, we mean solutions to (1.1) of the form ψ( t,x )= e iωt ϕ( x ) , where ω and ϕ is a solution to

( +iA ) 2 ϕ+ωϕ | x | b | ϕ | p ϕ=0. (1.6)

This study recalls the research on Cauchy solutions for the nonlinear Schrödinger equation. When b>0 , a compact embedding exists (see Lemma 2.3 [12]), which aids in establishing the compactness of Palais-Smale sequences and the existence of solutions. However, the inclusion of the | x | b term for b>0 complicates the analysis of the symmetry and decay properties of solutions, necessitating different approaches and careful treatment. It is important to emphasize that the parameter b>0 plays a crucial and significant role throughout the discussion.

The presence of standing waves for Equation (1.1) can be established by minimizing the energy functional E( u ) subject to the mass constraint

J( m ):={ u H A 1 ( 3 ):M( u )=m },

where m>0 . Specifically, we focus on the minimization problem

I( m ):=inf{ E( u ):uJ( m ) }.

Theorem 1.1 [13] Let 0<p< 42b 3 . For any m>0 , there exists a minimizer for I( m ) .

Theorem 1.2 [14] The set

( m ):={ ϕJ( m ):E( ϕ )=I( m ) }

is orbitally stable with respect to the flow of Equation (1.1), meaning that for any ε>0 , there exists δ>0 such that for any initial condition u 0 H A 1 ( 3 ) satisfying

inf ϕ( c ) ψ 0 ϕ H A 1 δ,

the corresponding solution to (1.1) exists globally in time and satisfies

inf ϕ( c ) inf y 3 e iA( y )x ψ( t,+y )ϕ H A 1 ε,t0.

Remark 1.1 Lee and Seo [15] established the existence of solutions under critical conditions by constructing a complete metric space using the contraction mapping principle and applying the fixed point theorem to obtain local solutions. However, this approach, which circumvents compactness issues by selecting a sufficiently small time interval T, is unable to analyze the long-term energy behavior of solutions or address the lack of compactness. Consequently, it cannot directly resolve problems related to non-compactness.

In contrast, we employ the concentration compactness lemma along with suitable scaling parameters to prove compactness by ruling out vanishing—demonstrated by showing that every minimizing sequence of I(m) has an L p+2 -norm bounded away from zero (see Lemma 3.1)—and dichotomy.

2. Preliminaries

In this section, let us recall some fundamental properties of the magnetic Sobolev space H A 1 ( 3 ) and provide several lemmas to support the subsequent proofs.

Lemma 2.1 [14] Suppose A L loc 2 ( 3 , 3 ) . Then H A 1 ( 3 ) equipped with the inner product

u,v H A 1 := 3 u v ¯ dx + 3 ( +iA )u ( +iA )v ¯ dx

is a Hilbert space.

Lemma 2.2 (Gagliardo-Nirenberg inequality [12]) Let 0<p<42b , and 0<b<2 , then the Gagliardo-Nirenberg inequality

3 | x | b | ψ( x ) | p+2 dx K opt ψ L 2 ( 3 ) 3p 2 +b ψ L 2 ( 3 ) 2 p 2 b (2.1)

holds, and the sharp constant K opt >0 is explicitly given by

K opt = ( 3p+2b 42bp ) 43p2b 4 2p+4 ( 3p+2b ) Q L 2 ( 3 ) p ,

where Q is the unique positive radial solution to

ΔQQ+ | x | b | Q | p Q=0. (2.2)

Moreover, the solution Q satisfies the following relations

Q L 2 ( 3 ) = ( 3p+2b 4p2b ) 1 2 Q L 2 ( 3 ) , (2.3)

and

3 | x | b | Q( x ) | p+2 dx =( 2p+4 4p2b ) Q L 2 ( 3 ) 2 . (2.4)

Lemma 2.3 (Diamagnetic inequality [16]) Let A L loc 2 ( 3 , 3 ) and u H A 1 ( 3 ) . Then | u | H 1 ( 3 ) . In particular, we have

| | u |( x ) || ( +iA )u( x ) | a.e. x 3 .(2.5)

Lemma 2.4 [14] Let A L loc 2 ( 3 , 3 ) . Then the following properties hold:

1) C 0 ( 3 ) is dense in H A 1 ( 3 ) .

2) H A 1 ( 3 ) is continuously embedded in L k ( 3 ) for all 2k62b .

3) Assume that A is linear, i.e., A( x+y )=A( x )+A( y ) for all x,y 3 . Let y 3 , u H A 1 ( 3 ) , and set

u ˜ ( x ):= e iA( y )x u( x+y ),x 3 .

Then ( +iA ) u ˜ ( x )= e iA( y )x ( +iA )u( x+y ) . In particular,

( +iA ) u ˜ L 2 = ( +iA )u L 2 .

4) If A L loc 3 ( 3 , 3 ) , then H A 1 ( 3 ) is continuously embedded in H loc 1 ( 3 ) . In particular, H A 1 ( 3 ) is compactly embedded in L loc k ( 3 ) for all 2k<62b .

Lemma 2.5 Let A W loc 1, ( 3 , 3 ) and j,k{ 1,2,3 } . Then for any u C 0 ( 3 ) , we have

| ( j A k k A j )u u ¯ dx | ( j +i A j )u L 2 2 + ( k +i A k )u L 2 2 .

In particular, if A is as in (1.2), then

| a | u L 2 2 ( +iA )u L 2 2 . (2.6)

Lemma 2.6 Let A L loc 2 ( 3 , 3 ) and ( u n ) n1 be a bounded sequence in H A 1 ( 3 ) . Assume that u n u weakly in H A 1 ( 3 ) . Then we have

( +iA ) u n L 2 2 = ( +iA )u L 2 2 + ( +iA )( u n u ) L 2 2 + o n ( 1 ),

3 | x | b | u n | k dx = 3 | x | b | u | k dx + 3 | x | b | u n u | k dx + o n ( 1 ),2k62b.

Proof As H A 1 ( 3 ) is continuously embedded in L r ( 3 ) for all 2r62b . The second identity follows directly from the refined Fatou’s lemma established by Brézis and Lieb [17]. Now, we will demonstrate the first identity. Set v n := u n u . We see that v n 0 weakly in H A 1 ( 3 ) . We compute

( +iA ) u n L 2 2 = ( +iA )( u+ v n ) L 2 2 = ( +iA )u L 2 2 + ( +iA ) v n L 2 2 +2Re ( +iA )u ( +iA ) v n ¯ dx .

Let ε>0 . Since C 0 ( 3 ) is dense in H A 1 ( 3 ) , we take ϕ C 0 ( 3 ) so that ( +iA )( uϕ ) L 2 <ε/ ( 2K ) , where K:= sup n1 v n H A 1 < . Since v n 0 weakly in H A 1 ( 3 ) , we see that

| ( +iA )ϕ ( +iA ) v n ¯ dx |0asn.

Thus there exists n 0 such that for n n 0 ,

| ( +iA )u ( +iA ) v n ¯ dx | | ( +iA )( uϕ ) ( +iA ) v n ¯ dx |+| ( +iA )ϕ ( +iA ) v n ¯ dx | ( +iA )( uϕ ) L 2 ( +iA ) v n L 2 +ε/2 <ε.

The proof is complete. □

Lemma 2.7 Let A L loc 3 ( 3 , 3 ) be linear. Let ( u n ) n1 be a bounded sequence in H A 1 ( 3 ) , i.e., sup n1 u n H A 1 < . Assume that there exists ε 0 >0 such that

inf n1 u n L k ε 0 (2.7)

for some 2<k<62b . Then up to a subsequence, there exist u H A 1 ( 3 )\{ 0 } and { y n } n1 3 such that

e iA( y n )x u n ( x+ y n )u weakly in H 1 ( 3 ) .

3. Existence and Stability of Normalized Standing Waves

In this section, we prove the existence and orbital stability of normalized standing waves related to (1.1). To prove it, we present the relevant arguments in this section. A prerequisite for this proof is the following result, which is essential to rule out the possibility of vanishing.

Lemma 3.1 Let A be as in (1.2) and 0<p< 42b 3 . Let m>0 and ( u n ) n1 be a minimizing sequence for I( m ) . Then there exists M>0 such that

liminf n u n L p+2 K>0.

Proof Assume by contradiction that there exists a subsequence still denoted by ( u n ) n1 satisfying lim n u n L p+2 =0 . As a result of Hölder’s inequality, we then know that

3 | x | b | u n | p+2 dx = B R ( 0 ) | x | b | u n | p+2 dx + 3 \ B R ( 0 ) | x | b | u n | p+2 dx ( B R ( 0 ) | x | br r1 dx ) r1 r ( B R ( 0 ) | u n | ( p+2 )r dx ) 1 r + R b 3 \ B R ( 0 ) | u n | p+2 dx ,

where R>0 is a constant. Note that

( B R ( 0 ) | x | br r1 dx ) r1 r <,

then we have lim n 3 | x | b | u n | p+2 dx =0 . Thanks to (2.6), we see that

I( m )= lim n E( u n )= lim n 1 2 ( +iA ) u n L 2 2 lim n | a | 2 u n L 2 2 = | a |m 2 . (3.1)

Denote x=( x , x 3 ) with x =( x 1 , x 2 ) 2 and x 3 , and set v( x ):= | a | 2π e | a | 4 | x | 2 . One can readily check that

v L 2 ( 2 ) =1, v L 2 ( 2 ) 2 + a 2 4 ρv L 2 ( 2 ) 2 =| a |.

Let w C 0 ( ) be such that w L 2 ( ) 2 =m and set

u λ ( x )=v( x ) w λ ( x 3 ), w λ ( x 3 )= λ 1 2 w( λ x 3 ) (3.2)

with λ>0 to be chosen later. We have u λ L 2 2 =c for all λ>0 . Using (1.4), we see that

( +iA ) u λ L 2 2 = u λ L 2 2 +bR( u λ )+ a 2 4 ρ u λ L 2 2 = v L 2 ( 2 ) 2 w λ L 2 ( ) 2 + v L 2 ( 2 ) 2 3 w λ L 2 ( ) 2 +a( 2 L z v v ¯ d x ) w λ L 2 ( ) 2 + a 2 4 ρv L 2 ( 2 ) 2 w λ L 2 ( ) 2 =m( v L 2 ( 2 ) 2 + a 2 4 ρv L 2 ( 2 ) 2 )+ λ 2 3 w L 2 ( ) 2 =m| a |+ λ 2 3 w L 2 ( ) 2 ,

where

2 L z v v ¯ d x =i 2 x 2 v 1 v ¯ d x 1 d x 2 +i 2 x 2 v 2 v ¯ d x 1 d x 2 ,

as v is radially symmetric, we have

1 v ¯ = v ( x ) x 1 x

2 v ¯ = v ( x ) x 2 x ,

then we get L z v v ¯ d x =0 . It follows that

E( u λ )= | a |m 2 + λ 2 2 3 w L 2 ( ) 2 C λ p 2 +b Φ( x ).

where C is a constant, Φ is a function of x . As p<42b , by taking λ>0 sufficiently small, we have E( f λ )< | a |m 2 . In particular, I( m )< | a |m 2 which contradicts (3.1). Then the proof is complete. □

Proof of Theorem 1.1 We first show that I( m ) is well-defined, i.e., I( m )> . Let uJ( m ) and 0<ε< 1 2 , by the Young’s inequality and Gagliardo-Nirenberg inequality (2.1), we have

E( u ) 1 2 ( +iA )u L 2 2 M( ε ) ( +iA )u L 2 2 C( ε ) u L 2 164p2b 43p2b =( 1 2 M( ε ) ) ( +iA )u L 2 2 C>C

for all uJ( m ) , where 0<C< . This shows that I( m )> . Now let ( u n ) n1 be a minimizing sequence for I( m ) . From the above estimate, we have

( 1 2 M( ε ) ) ( +iA ) u n L 2 2 E( u n )I( m )asn.

This shows that ( u n ) n1 is a bounded sequence in H A 1 ( 3 ) .

Moreover, by Lemma 3.1, we see that up to a subsequence,

inf n1 u n L 102b 3 K>0.

By Lemma 2.7, up to a subsequence, there exist u H A 1 ( 3 )\{ 0 } and ( y n ) n1 3 such that

u ˜ n ( x ):= e iA( y n )x u n ( x+ y n )uweaklyin H A 1 ( 3 ).

By the weak convergence in H A 1 ( 3 ) , we have

0< u L 2 2 liminf n u ˜ n L 2 2 = liminf n u n L 2 2 =m

and

( +iA )u L 2 2 liminf n ( +iA ) u ˜ n L 2 2 = liminf n ( +iA ) u n L 2 2 .

Next we claim that

u L 2 2 =m. (3.3)

Let’s delay verifying (3.2) for now and complete the proof of Theorem 1.1. By the weak convergence in H A 1 ( 3 ) and (3.2), we infer that u ˜ n f strongly in L 2 ( 3 ) . Using this strong convergence and the magnetic Gagliardo-Nirenberg inequality

3 | x | b | u | p+2 dx K opt ( +iA )u L 2 ( 3 ) 3p 2 +b u L 2 ( 3 ) 2 p 2 b ,

we see that lim n 3 | x | b | u ˜ n | p+2 dx = 3 | x | b | u | p+2 dx . Thus we get

I( m )E( u ) liminf n E( u ˜ n )= liminf n E( u n )=I( m ),

hence E( u )=I( m ) or u is a minimizer for I( m ) . This also implies that u ˜ n f strongly in H A 1 ( 3 ) .

It remains to prove (3.2). Assume by contradiction that it is not true, i.e., 0< u L 2 2 <m . We have for any λ>0 ,

E( λu )= λ 2 E( u )+ λ 2 ( 1 λ p ) p+2 3 | x | b | u | p+2 dx

or

E( f )= 1 λ 2 E( λf )+ λ p 1 p+2 3 | x | b | u | p+2 dx .

Set λ 0 = m u L 2 >1 . We have λ 0 u L 2 2 =m and

E( u )= u L 2 2 m E( λ 0 u )+ λ 0 p 1 p+2 3 | x | b | u | p+2 dx > u L 2 2 m I( m )

as u0 and λ 0 >1 . Similarly, set λ n := m u ˜ n u L 2 . By Lemma 2.6, we have u ˜ n u L 2 2 m u L 2 2 as n , hence λ n m m u L 2 2 >1 as n . In particular, we have

lim n E( u ˜ n u )= lim n ( 1 λ n 2 E( λ n ( u ˜ n u ) )+ λ n α 1 α+2 u ˜ n u L α+2 α+2 ) m u L 2 2 m I( m ).

Using the refined Fatou’s lemma (see Lemma 2.6), we get

I( m )= lim n E( u n )= lim n E( u ˜ n )=E( u )+ lim n E( u ˜ n u ) > u L 2 2 m I( m )+ m u L 2 2 m I( m )=I( m ),

which is a contradiction. The proof is completed. □

Proof of Theorem 1.2 Let us now demonstrate that the set of minimizers ( m ) is orbitally stable as described in Theorem 1.1. We will use an argument from [13]. Assume by contradiction that it is not true. Then there exist ε 0 >0 , ϕ 0 ( m ) , and a sequence of initial data ( ψ 0,n ) n1 H A 1 ( 3 ) such that

lim n ψ 0,n ϕ 0 H A 1 =0 (3.4)

and a sequence of time ( t n ) n1 [ 0, ) such that

inf ϕ( m ) inf y 3 e iA( y ) ψ n ( t n ,+y )ϕ H A 1 ε 0 , (3.5)

where ψ n is the solution to (1.1) with initial data ψ n | t=0 = ψ 0,n . Since ϕ 0 ( m ) , we have E( ϕ 0 )=I( m ) . From (3.3) and the Sobolev embedding, we infer that

u 0,n L 2 2 ϕ 0 L 2 2 =m,E( u 0,n )E( ϕ 0 )=I( m )asn.

By the conservation laws of mass and energy, we have

ψ n ( t n ) L 2 2 m,E( ψ n ( t n ) )I( m )asn.

In particular, ( ψ n ( t n ) ) n1 is a minimizing sequence for I( m ) .As argued in Step 1, we see that up to a subsequence, there exist ϕ( m ) and ( y n ) n1 3 such that

e iA( y n ) ψ n ( t n ,+ y n )ϕ H A 1 0asn.

However, this contradicts (3.4). Thus, the proof is finished. □

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1] Tsutsumi, Y. (1987) Schrödinger Equation. Funkcialaj Ekvacioj, 30, 115-125.
http://fe.math.kobe-u.ac.jp/FE/FE_pdf_with_bookmark/FE21-30-en_KML/fe30-115-125/fe30-115-125.pdf
[2] Hefter, E.F. (1985) Application of the Nonlinear Schrödinger Equation with a Logarithmic Inhomogeneous Term to Nuclear Physics. Physical Review A, 32, 1201-1204.[CrossRef] [PubMed]
[3] Yang, R., Li, L., Hao, R., Li, Z. and Zhou, G. (2005) Combined Solitary Wave Solutions for the Inhomogeneous Higher-Order Nonlinear Schrödinger Equation. Physical Review E, 71, Article ID: 036616.[CrossRef] [PubMed]
[4] Tang, X. and Shukla, P.K. (2007) Solution of the One-Dimensional Spatially Inhomogeneous Cubic-Quintic Nonlinear Schrödinger Equation with an External Potential. Physical Review A, 76, Article ID: 013612.[CrossRef]
[5] Schrödinger, E. (1926) An Undulatory Theory of the Mechanics of Atoms and Molecules. Physical Review, 28, 1049-1070.[CrossRef]
[6] Lions, P.L. (1984) The Concentration-Compactness Principle in the Calculus of Variation. The Locally Compact Case, Part I. Annales de lInstitut Henri PoincaréAnalyse Non Linéaire, 1, 109-145, 223-283.
[7] Lions, P.L. (1985) The Concentration-Compactness Principle in the Calculus of Variation. The Locally Compact Case, Part II. Annales de lInstitut Henri PoincaréAnalyse Non Linéaire, 1, 45-121.
[8] Avron, J.E., Herbst, I.W. and Simon, B. (1981) Schrödinger Operators with Magnetic Fields. III. Atoms in Homogeneous Magnetic Fields. Communications in Mathematical Physics, 79, 529-572.[CrossRef]
[9] Avron, J., Herbst, I. and Simon, B. (1978) Schrödinger Operators with Magnetic Fields. I. General Interactions. Duke Mathematical Journal, 45, 847-883.[CrossRef]
[10] Avron, J.E., Herbst, I.W. and Simon, B. (1978) Separation of Center of Mass in Homogeneous Magnetic Fields. Annals of Physics, 114, 431-451.[CrossRef]
[11] Dinh, V.D. (2022) Existence and Stability of Standing Waves for Nonlinear Schrödinger Equations with a Critical Rotational Speed. Letters in Mathematical Physics, 112, Article No. 53.[CrossRef]
[12] Gou, T. (2023) Standing Waves with Prescribed L2-Norm to Nonlinear Schrödinger Equations with Combined Inhomogeneous Nonlinearities. Letters in Mathematical Physics, 114, Article No. 7.[CrossRef]
[13] Cazenave, T. and Esteban, M.J. (1988) On the Stability of Stationary States for Nonlinear Schrödinger Equations with an External Magnetic Field. Matematica Applicata e Computazione, 7, 155-168.
https://mathscinet.ams.org/mathscinet/article?mr=994761
[14] Esteban, M.J. and Lions, P. (1989) Stationary Solutions of Nonlinear Schrödinger Equations with an External Magnetic Field. In: Colombini, F., Marino, A., Modica, L. and Spagnolo, S., Eds., Partial Differential Equations and the Calculus of Variations, Birkhäuser Boston, 401-449.[CrossRef]
[15] Lee, Y. and Seo, I. (2021) The Cauchy Problem for the Energy-Critical Inhomogeneous Nonlinear Schrödinger Equation. Archiv der Mathematik, 117, 441-453.[CrossRef]
[16] Lieb, E.H. and Loss, M. (2001) Analysis (Graduate Studies in Mathematics, Vol. 14). 2nd Edition, American Mathematical Society.
[17] Brézis, H. and Lieb, E. (1983) A Relation between Pointwise Convergence of Functions and Convergence of Functionals. Proceedings of the American Mathematical Society, 88, 486-490.[CrossRef]

Copyright © 2026 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.