Skip to main content

Blow-up criteria for the inhomogeneous nonlinear Schrödinger equation

Abstract

In this paper, using the variational characteristic of the virial identity and a new estimate of the kinetic energy, we obtain a new sufficient condition for the existence of blow-up solutions.

MSC:35Q55, 35B44.

1 Introduction

In this paper, we study the Cauchy problem of the inhomogeneous nonlinear Schrödinger equation (INSE)

i u t +u+|x | b |u | p 2 u=0,t0,x R N ,
(1.1)
u(0,x)= u 0 ,
(1.2)

where i 2 =1; = j = 1 N 2 x j 2 is the Laplace operator in R N ; u=u(t,x): [0,T)× R N C is the complex valued function and 0<T+; the parameter b0 and 2<p< p ˜ (we use the convention: p ˜ =+ for N=2, p ˜ = 2 N N 1 + 2 b N 1 for N3); N2 is the space dimension. A few years ago, it was suggested that stable high power propagation can be achieved in a plasma by sending a preliminary laser beam that creates a channel with a reduced electron density, and thus reduces the nonlinearity inside the channel (see [1, 2]). In this case, the beam propagation can be modeled by the inhomogeneous nonlinear Schrödinger equation in the following form:

i ϕ t +ϕ+K(x)|ϕ | p 2 ϕ=0,ϕ(0,x)=φ H 1 ( R N ) .
(1.3)

Recently, this type of inhomogeneous nonlinear Schrödinger equations has been widely investigated. When k 1 K(x) k 2 with k 1 , k 2 >0 and p=2+ 4 N , Merle [3] proved the existence and nonexistence of blow-up solutions to the Cauchy problem (1.3). When K(x)=K(ε|x|) C 4 ( R N ) L ( R N ) with small ε and p=2+ 4 N , Fibich, Liu and Wang [2, 4] obtained the stability and instability of standing waves to the Cauchy problem (1.3).

We recall some known results on the blow-up solutions for the classical nonlinear Schrödinger equation

i v t +v+|v | p 2 v=0,v(0,x)= v 0 H 1 ( R N ) .
(1.4)

Ginibre and Velo [5] showed the local well-posedness in H 1 ( R N ). Glassey [6] showed the existence of blow-up solutions when the energy is negative and |x| v 0 L 2 ( R N ). Ogawa and Tsutsumi [7] obtained the existence of blow-up solutions in radial case without the restriction |x| v 0 L 2 ( R N ). Weinstein [8] and Zhang [9] obtained the sharp conditions of global existence for critical and supercritical nonlinearity. Merle and Raphaël [10] showed the existence of blow-up solutions without |x| v 0 L 2 ( R N ) for p=2+ 4 N . Lushnikov [11] and Holmer et al. [12] obtained some sufficient conditions for existence of blow-up for p=4 and N=3 basing on an estimate of the kinetic energy.

In this paper, we study blow-up criteria for the Cauchy problem (1.1)-(1.2) with p> 2 N + 2 b + 4 N , where the nonlinearity |x | b |u | p 2 u includes an unbounded potential |x | b . We note that (1.1) has scaling: u λ (t,x)= λ 2 + b p 2 u( λ 2 t,λx) is a solution if u(t,x) is a solution. The scale-invariant Lebesgue norm for this equation is L p c -norm, where p c = N ( p 2 ) 2 + b . Since p> 2 N + 2 b + 4 N , we have p c = N ( p 2 ) 2 + b >2 and we may call (1.1) a class of Schrödinger equations with L 2 -super critical nonlinearity. Chen and Guo [13] and Chen [14] showed the local well-posedness of the Cauchy problem (1.1)-(1.2) in H r 1 = H r 1 ( R N ), where H r 1 ( R N ) is the set of radially symmetric functions in H 1 ( R N ). Moreover, u(t,x) satisfies the following conservation laws:

R N |u(t,x) | 2 dx= R N | u 0 | 2 dx

and

E( u 0 )=E(u):= 1 2 R N |u | 2 dx 1 p R N |x | b |u | p dx.

Chen and Guo [13] also showed the sharp conditions of blow-up and global existence of solutions to the Cauchy problem (1.1)-(1.2) by the cross-constrained variational arguments. On the other hand, letting J(t):= R N |x | 2 |u(t,x) | 2 dx, this can be interpreted as the average width of the initial distribution |u|. It follows from Chen and Guo’s results in [13] that we have the following proposition.

Proposition 1.1 Assume that u 0 H r 1 , |x| u 0 L 2 and the corresponding solution u(t,x) of the Cauchy problem (1.1)-(1.2) on the interval [0,T). Then, for all t[0,T), one has J(t):=|x | 2 |u(t,x) | 2 dx<+,

J (t)=4 R N x u ¯ udx
(1.5)

and

J (t)=4(Np2N2b)E( u 0 )2(Np2N2b4) R N |u | 2 dx.
(1.6)

In the case Np2N2b40, it follows from the last relation that

J (t)4(Np2N2b)E( u 0 ),

and for E( u 0 )<0 the positive-definite quantity J(t) becomes negative over a finite time by virtue of the above inequality. This means that a singularity appears in the solution of the given INSE. Indeed, applying Weinstein’s arguments [8] and the classical analysis identity

J(t)=J(0)+ J (0)t+ 0 t J (ts)sds,
(1.7)

one has the following theorem (see also Chen and Guo [13]).

Theorem 1.2 Let N2, 0<b<N2 and 2 N + 2 b + 4 N <p< p ˜ (where p ˜ =+ for N=2, p ˜ = 2 N N 2 + 2 b N 1 for N3). Assume u 0 H 1 ( R N ) and |x| u 0 L 2 ( R N ) is radially symmetric. If the initial data satisfies either

  1. (i)
    E( u 0 )<0,
    (1.8)
  2. (ii)
    E( u 0 )=0and J (0)<0,
    (1.9)
  3. (iii)

    E( u 0 )>0 and

    J (0)<2 2 ( N p 2 N 2 b ) E ( u 0 ) J(0),
    (1.10)

then there exists 0<T<+ such that the corresponding solution u(t,x) blows up in finite time T.

We remark that in the case E( u 0 )>0, Np2N2b40 both collapse and spreading of the initial disturbance are possible. Although the INSE is no longer applicable near the formation point of a singularity and dissipative or some other limiting mechanism come to play. It is very important to be able to predict the presence or absence of collapse for different classes of initial conditions. The sufficient conditions for existence of blow-up solutions are given in [13] if either E( u 0 )<0 or J (0)<0. A natural question arises whether there is a sufficient condition for existence of blow-up solutions with E( u 0 )>0 and J (0)>0.

In the present paper, motivated by the studies of the classical nonlinear Schrödinger equation (see [8, 11, 12]), we use variational characteristic of second-order derivatives of the virial identity to catch up with the information of u L 2 , and we obtain a new sufficient condition for the existence of blow-up solutions to the inhomogeneous nonlinear Schrödinger equation (1.1). More precisely, let

g(y)={ y + 4 N p 2 b 2 N 4 y N p 2 b 2 N 4 4 N p 2 b 2 N N p 2 b 2 N 4 if  0 < y < 1 , y + 4 N p 2 b 2 N 4 y N p 2 b 2 N 4 4 N p 2 b 2 N N p 2 b 2 N 4 if  y 1 .
(1.11)

Then we have the following theorem.

Theorem 1.3 Let N2, 0<b<N2 and 2 N + 2 b + 4 N <p<min{ 2 b + 2 N + 12 N , p ˜ }. Assume that u 0 H 1 ( R N ) and |x| u 0 L 2 ( R N ) is radially symmetric. If

J (0)<2 N p 2 b 2 N 4 N p 2 b 2 N N u 0 L 2 2 g ( 8 ( N p 2 b 2 N ) E ( u 0 ) N 2 ( N p 2 b 2 N 4 ) u 0 L 2 4 J ( 0 ) ) ,
(1.12)

where g(x) is defined by (1.11), then there exists 0<T<+ such that the corresponding solution u(t,x) blows up in finite time T.

2 Notations and preliminaries

In this paper, we denote L q ( R N ), L q ( R N ) , H s ( R N ) and R N dx by L q , L q , H s and dx, respectively. z and z are the real part and imaginary part of the complex number z, respectively. z ¯ is denoted the complex conjugate of the complex number z. The various positive constants will be simply denoted by C.

For the Cauchy problem (1.1)-(1.2), the space we work in is

H r 1 := { u H 1 | u ( x ) = u ( r ) } ,where r= x 1 2 + x 2 2 + + x N 2 ,

which is a Hilbert space. Moreover, we define the energy functional E(u(t)) in H r 1 by

E ( u ( t ) ) := 1 2 |u(t,x) | 2 dx 1 p |x | b |u(t,x) | p dx.

The functional E(u) is well-defined according to the Sobolev embedding theorem (see [15]). Chen and Guo [13] and Chen [14] showed the local well-posedness for the Cauchy problem (1.1)-(1.2) in H r 1 , as follows.

Proposition 2.1 Let N2, b0 and 2+2b/(N1)<p< p ˜ (where p ˜ =+ for N=2, p ˜ = 2 N N 2 + 2 b N 1 for N3). For any u 0 H r 1 , there exists a unique solution u(t,x) of the Cauchy problem (1.1)-(1.2) on the maximal time interval [0,T) such that u(t,x)C([0,T); H r 1 ) and either T=+ (global existence), or T<+ and lim t T u ( t , x ) H r 1 =+ (blow-up). Furthermore, for all t[0,T), u(t,x) satisfies the following conservation laws:

  1. (i)

    Conservation of mass: u ( t ) 2 = u 0 2 .

  2. (ii)

    Conservation of energy: E(u(t,x))=E( u 0 ).

In addition, by some basic calculations, we have the following lemma, which gives further insight in the dynamic criterion for collapse proposed by Lushnikov in [11].

Lemma 2.2 If V(t)>0 is the positive solution of the following differential equation:

d 2 d t 2 V=f(V)
(2.1)

and there exists 0< T 1 <+ such that lim t T 1 V(t)=0, then for the solution V ˜ >0 of the following differential equation:

V ˜ t t =f( V ˜ ) h 2 (t)f( V ˜ ),
(2.2)

there exists 0< T 2 T 1 <+ such that lim t T 2 V ˜ (t)=0.

Proof Since the function h 2 (t) is non-positive, which pulls V ˜ (t) to zero more quickly than V(t) (see also [11]), one sees that the conclusion in Lemma 2.2 is true by the classical analysis identity (1.7). □

3 Proof of Theorem 1.3

Since |x| u 0 L 2 , we have |x|u(t,x) L 2 by the local well-posedness. Taking J:=J(t)=|x | 2 |u | 2 dx, by Proposition 1.1, we get J (t)=4xu u ¯ dx and

J (t)=4(Np2N2b)E( u 0 )2(Np2N2b4)|u | 2 dx.
(3.1)

It follows from some calculations that

|u | 2 dx= 1 N x|u | 2 dx= 2 N xu u ¯ dx.

Then we get N 2 u L 2 2 x u L 2 u L 2 , and by the fact that |z | 2 =|z | 2 +|z | 2 we deduce that

u L 2 2 x u L 2 2 | x u u ¯ d x | 2 = N 2 4 u L 2 4 + | x u u ¯ d x | 2 = N 2 4 u L 2 4 + | J ( t ) | 2 16 .
(3.2)

Injecting (3.2) into (3.1), by the conservation laws, we deduce that

J ( t ) 4 ( N p 2 N 2 b ) E ( u 0 ) N 2 [ N p 2 N 2 b 4 ] 2 J u 0 L 2 4 N p 2 N 2 b 4 8 J J t 2 .
(3.3)

Letting J= A 8 N p 2 N 2 b + 4 and rewriting (3.3) to remove the last term with J t 2 , we have

A t t ( N p 2 b 2 N ) ( N p 2 b 2 N + 4 ) 2 E ( u 0 ) A N p 2 b 2 N 4 N p 2 b 2 N + 4 N 2 ( N p 2 b 2 N 4 ) ( N p 2 b 2 N + 4 ) 16 u 0 L 2 4 A 12 + 2 b + 2 N N p N p 2 b 2 N + 4 ,
(3.4)

which has a simple mechanism analogy. Multiplying A t in (3.4) and integrating with the time variable t, we get the corresponding mechanical energy

ε(t)= 1 2 A t 2 +U ( A ( t ) ) ,
(3.5)

where

U ( A ) = ( N p 2 b 2 N + 4 ) 2 4 E ( u 0 ) A 2 ( N p 2 b 2 N ) N p 2 b 2 N + 4 + N 2 ( N p 2 b 2 N + 4 ) 2 32 u 0 L 2 4 A 2 ( N p 2 b 2 N 4 ) N p 2 b 2 N + 4 .

Restricting ourselves to the case E>0, and according to the assumptions on p, b and N we see that U(A) achieves its maximum U max at A max with

A max = ( N 2 ( N p 2 b 2 N 4 ) u 0 L 2 4 8 ( N p 2 b 2 N ) E ( u 0 ) ) N p 2 b 2 N + 4 8

and

U max = ( N p 2 b 2 N + 4 ) 2 N p 2 b 2 N 4 E( u 0 ) ( N 2 ( N p 2 b 2 N 4 ) u 0 L 2 4 8 ( N p 2 b 2 N ) E ( u 0 ) ) N p 2 b 2 N 4 .

To facilitate the rest of the analysis, we introduce a rescaling. Define J ˜ (s) and ε ˜ (s) by the relations

A ( t ) = A max A ˜ ( α t ) , ε ( t ) = U max ε ˜ ( α t ) , α = N p 2 b 2 N N p 2 b 2 N 4 16 E ( u 0 ) N u 0 L 2 2 , s = α t .

Thus, by (3.4), A ˜ satisfies the following differential inequality:

A ˜ s s N p 2 b 2 N + 4 64 ( A ˜ N p 2 b 2 N 4 N p 2 b 2 N + 4 A ˜ 12 + 2 b + 2 N N p N p 2 b 2 N + 4 ) .
(3.6)

Applying Lemma 2.2, if we show that for the positive solution of the following differential equation:

A ˜ s s = N p 2 b 2 N + 4 64 ( A ˜ N p 2 b 2 N 4 N p 2 b 2 N + 4 A ˜ 12 + 2 b + 2 N N p N p 2 b 2 N + 4 ) ,
(3.7)

there exists a time 0< s 1 <+ and lim s s 1 A ˜ (s)=0, then for the positive solution of the differential inequality (3.6), there exists a time 0< s 2 s 1 <+ and lim s s 2 A ˜ (s)=0. Indeed, setting

U ˜ ( A ˜ )= N p 2 b 2 N 4 4 A ˜ 2 ( N p 2 b 2 N ) N p 2 b 2 N + 4 + N p 2 b 2 N 4 A ˜ 2 ( N p 2 b 2 N 4 ) N p 2 b 2 N + 4 ,
(3.8)

we see that (3.5) converts to

ε ˜ (s)= 16 ( N p 2 b 2 N 4 ) ( N p 2 b 2 N + 4 ) 2 A ˜ s 2 (s)+ U ˜ ( A ˜ ( s ) ) .
(3.9)

It is obvious that the maximum of U ˜ ( A ˜ (s)) is 1, which is attained by the maximum at A ˜ (s)=1. By the variational characteristic of U ˜ ( A ˜ (s)), we claim that under one of the following conditions, A ˜ (s) vanishes in finite time, so does J(t) (which implies the solution u(t,x) blows up in finite time):

  1. (a)

    ε ˜ (0)<1 and A ˜ (0)<1,

  2. (b)

    ε ˜ (0)1 and A ˜ s (0)<0.

Indeed, it follows from (3.7) and (3.9) that ε ˜ (s)= ε ˜ (0). (a) If ε ˜ (0)<1, then ε ˜ (t)<1 and U ˜ ( A ˜ (s))<1 for tI (maximal existence interval). By the assumption A ˜ (0)<1, we deduce that A ˜ (s)<1 for tI. Moreover, we have

A ˜ N p 2 b 2 N 4 N p 2 b 2 N + 4 A ˜ 12 + 2 b + 2 N N p N p 2 b 2 N + 4 <0.

Using (3.7) and the classical analysis identity

A ˜ (s)= A ˜ (0)+ A ˜ (0)s+ 0 s A ˜ (sτ)τdτ,
(3.10)

we see that A ˜ (s) vanishes in finite time. On the other hand, for the second case (b), if ε ˜ (0)1, then ε ˜ (t)1. If follows from (3.9) that U ˜ ( A ˜ (s))=1, which implies that A ˜ (s)=1 and

A ˜ N p 2 b 2 N 4 N p 2 b 2 N + 4 A ˜ 12 + 2 b + 2 N N p N p 2 b 2 N + 4 =0.

Using (3.7) and the classical analysis identity (3.10), we deduce that if A ˜ s (0)<0, then A ˜ (s) vanishes in finite time. This completes the proof of claim (a) and (b). Now, we return to the proof of Theorem 1.3. If we define J ˜ = A ˜ 8 N p 2 b 2 N + 4 , then (3.9) is equal to

ε ˜ (s)= N p 2 b 2 N 4 4 J ˜ N p 2 b 2 N 4 4 [ J ˜ s 2 J ˜ + N p 2 b 2 N N p 2 b 2 N 4 ] .

Taking

h(y)= y + 4 N p 2 b 2 N 4 y N p 2 b 2 N 4 4 N p 2 b 2 N N p 2 b 2 N 4 ,

we see that

ε ˜ <1| J ˜ s |<h( J ˜ )and ε ˜ 1| J ˜ s |h( J ˜ ).

Thus condition (a) is true if and only if

J ˜ (0)<1andh ( J ˜ ( 0 ) ) < J ˜ s (0)<h ( J ˜ ( 0 ) )
(3.11)

and condition (b) is true if and only if J ˜ s (0)<h( J ˜ (0)). Collecting the above two conditions, we deduce that the solution u(t,x) blows up in finite time 0<T<+ provided

J ˜ s (0)<{ h ( J ˜ ( 0 ) ) if  0 < J ˜ ( 0 ) < 1 , h ( J ˜ ( 0 ) ) if  J ˜ ( 0 ) 1 .
(3.12)

Finally, substituting back J(t), we see that

J ˜ (0)= 8 ( N p 2 b 2 N ) E ( u 0 ) N 2 ( N p 2 b 2 N 4 ) u 0 L 2 4 J(0), J ˜ s (0)= N p 2 b 2 N N p 2 b 2 N 4 J t ( 0 ) 2 N u 0 L 2 2 ,

which implies that (3.12) is equivalent to

J t (0)<{ 2 N p 2 b 2 N 4 N p 2 b 2 N N u 0 L 2 2 h ( 8 ( N p 2 b 2 N ) E ( u 0 ) N 2 ( N p 2 b 2 N 4 ) u 0 L 2 4 J ( 0 ) ) if  0 < J ˜ ( 0 ) < 1 , 2 N p 2 b 2 N 4 N p 2 b 2 N N u 0 L 2 2 h ( 8 ( N p 2 b 2 N ) E ( u 0 ) N 2 ( N p 2 b 2 N 4 ) u 0 L 2 4 J ( 0 ) ) if  J ˜ ( 0 ) 1 .
(3.13)

References

  1. Gill TS: Optical guiding of laser beam in nonuniform plasma. Pramana J. Phys. 2000, 55: 842–845.

    Google Scholar 

  2. Liu CS, Tripathi VK: Laser guiding in an axially nonuniform plasma channel. Phys. Plasmas 1994, 1: 3100–3103. 10.1063/1.870501

    Article  Google Scholar 

  3. Merle F:Nonexistence of minimal blow up solutions of equations i u t =uK(x) | u | 4 / N u in R N . Ann. IHP, Phys. Théor. 1996, 64: 33–85.

    MathSciNet  MATH  Google Scholar 

  4. Liu Y, Wang XP, Wang K: Instability of standing waves of the Schrödinger equations with inhomogeneous nonlinearities. Trans. Am. Math. Soc. 2006, 358: 2105–2122. 10.1090/S0002-9947-05-03763-3

    Article  MATH  Google Scholar 

  5. Ginibre J, Velo G: On a class of nonlinear Schrödinger equations. I. The Cauchy problem, general case. J. Funct. Anal. 1979, 32: 1–32. 10.1016/0022-1236(79)90076-4

    Article  MathSciNet  MATH  Google Scholar 

  6. Glassey RT: On the blowing up of solutions to the Cauchy problem for nonlinear Schrödinger equations. J. Math. Phys. 1977, 18: 1794–1797. 10.1063/1.523491

    Article  MathSciNet  MATH  Google Scholar 

  7. Ogawa T, Tsutsumi Y:Blow-up of H 1 solution for the nonlinear Schrödinger equation. J. Differ. Equ. 1991, 92: 317–330. 10.1016/0022-0396(91)90052-B

    Article  MathSciNet  MATH  Google Scholar 

  8. Weinstein MI: Nonlinear Schrödinger equations and sharp interpolation estimates. Commun. Math. Phys. 1983, 87: 567–576. 10.1007/BF01208265

    Article  MATH  Google Scholar 

  9. Zhang J: Sharp conditions of global existence for nonlinear Schrödinger and Klein-Gordon equations. Nonlinear Anal. 2002, 48: 191–207. 10.1016/S0362-546X(00)00180-2

    Article  MathSciNet  MATH  Google Scholar 

  10. Merle F, Raphaël P: Blow-up dynamic and upper bound on the blow-up rate for critical nonlinear Schrödinger equation. Ann. Math. 2005, 16: 157–222.

    Article  MATH  Google Scholar 

  11. Lushnikov PM: Dynamic criterion for collapse. JETP Lett. 1995, 62: 461–467.

    Google Scholar 

  12. Holmer J, Platte R, Roudenko S: Blow-up criteria for the 3D cubic nonlinear Schrödinger equation. Nonlinearity 2010, 23: 977–1030. 10.1088/0951-7715/23/4/011

    Article  MathSciNet  MATH  Google Scholar 

  13. Chen JQ, Guo BL: Sharp global existence and blowing up results for inhomogeneous Schrödinger equations. Discrete Contin. Dyn. Syst., Ser. B 2007, 8: 357–367.

    Article  MathSciNet  MATH  Google Scholar 

  14. Chen JQ: On a class of nonlinear inhomogeneous Schrödinger equation. J. Appl. Math. Comput. 2010, 32: 237–253. 10.1007/s12190-009-0246-5

    Article  MathSciNet  MATH  Google Scholar 

  15. Cazenave T Courant Lecture Notes in Mathematics 10. In Semilinear Schrödinger Equations. Am. Math. Soc., Providence; 2003.

    Google Scholar 

Download references

Acknowledgements

This work is supported partly by the National Natural Science Foundation of P.R. China grants 11226162 and 11371267.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shihui Zhu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

HY participated in the design of the study, SZ studied the virial identity, participated in the sequence alignment and drafted the manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Yang, H., Zhu, S. Blow-up criteria for the inhomogeneous nonlinear Schrödinger equation. J Inequal Appl 2014, 55 (2014). https://doi.org/10.1186/1029-242X-2014-55

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1029-242X-2014-55

Keywords