arxiv: v1 [math.ap] 4 Dec 2017

Transcription

1 RESOLVENT ESTIMATES FOR THE 3D SHRÖDINGER OPERATOR WITH INVERSE-SQUARE POTENTIAL A. ADAM AZZAM arxiv: v1 [math.ap] Dec 017 Abstract. We consider the unitary group e ith for the Schrödinger operator with inverse-suare potential H := 1 x. We adapt ombes- Thomas estimates to show that, when restricted to non-radial functions, H enjoys much better estimates that mirror those of the Laplacian. ontents 1. Introduction. Main Results 3 3. Acknowledgements. Preliminaries 5 5. ombes-thomas Estimates 5 6. Kernel Estimates 7 7. Proofs of Main Results 10 References 1 Date: November 9,

2 A. ADAM AZZAM 1. Introduction In this paper, we will consider harmonic analysis uestions concerning the unitary group e ith for the Schrödinger operator with inverse-suare potential (1.1) H := 1 x where is the Laplacian in R 3. The operator H in (1.1) belongs to a family of operators (1.) H a := a a 1 x of interest in uantum mechanics (see [3], [6]). Since the potential term is homogenous of order, the potential and Laplacian scale exactly the same. For this reason, this makes the study of the associated PDEs { (i t H a )u = 0, (t,x) R R 3, (1.3) u(0,x) = u 0 (x) Ḣ1 x(r 3 ) { (i t H a )u = u p u, (t,x) R R 3, (1.) u(0,x) = u 0 (x) Ḣ1 x (R3 ) of particular mathematical interest. Indeed, since the potential and Laplacian are of eual strength at every scale, analysis of (1.3) is immune to standard perturbative methods. Thus, Schrödinger euations with inverse suare potentials have necessitated a more delicate approach. This challenge has been met with a long campaign to understand those properties enjoyed by H a and the local and global behavior of solutions to (1.). The choice of the constant a = 1 in this paper is critical in another sense. This constant corresponds to the sharp constant in the Hardy ineuality. Lemma 1.1 (Hardy Ineuality). If u 0 (R 3 ), then 1 u(x) (1.5) dx u(x) dx. x Moreover, (1.5) fails if the constant 1 is replaced by any constant a > 1. When a 1 and u 0 (R 3 ) we see that (1.6) H a u,u L x (R 3 ) = u(x) dx a R 3 u(x) x dx 0, and so H a is positive semi-definite. When a > 1, all self-adjoint extensions of H a are unbounded below. In this case the spectrum of H a has infinitely many negative eigenvalues that diverge to ([11]), making it impossible to obtain global dispersive estimates. Substantial progress has been made in understanding H a and (1.) in both the critical case a = 1/ and the subcritical case a < 1. The first wave of progress in the linear problem came in the form of studying the full range

3 RESOLVENT ESTIMATES FOR 3D NLS WITH INVERSE-SQUARE POTENTIAL 3 of dispersive and Strichartz estimates for solutions to (1.3). When a 1, [], Bur Planchon Stalker Tahvildar-Zadeh demonstrated local smoothing for solutions to (1.3). The nonlinear problem in this critical regime was advanced in [1], when Suzuki demonstrated global wellposedness for (1.) in R 3 when a 1 and p < 5 for a suitable class of initial data by proving Strichartz estimates for (1.3) avoiding the sharp endpoint (, 6). Recently, in [10], Mizutani recovers the sharp endpoint in a Lorentz space framework. In [8], Killip-Miao-Visan-Zhang-Zheng prove a variety of harmonic analysis tools for studying H a when a 1. In particular, they develop a Littlewood- Paley theory for H a, which allows them to prove in [7] that (1.) with p = 5 is globally well-posed when a < 1 1. This represents the first progress made 5 in treating the energy critical nonlinearity u u in the presence of an inverse suare potential.. Main Results In this paper, we advocate for the idea that any problem involving the critical inverse suare potential should be decomposed into a system of euations: one for the radial part of the solution (which is more susceptible to classical ODE techniues), and another for the remainder. Indeed, what is clear from previous advances in [], [10], and [1] is that data in L L rad enjoy a wealth of useful estimates for H a which in turn yield Strichartz estimates that are not immediately available in the general case. When restricted to higher angular momentum, the operator H a enjoys estimates that mirror those of the Laplacian. To make this discussion more precise, let P : L (R 3 ) L rad (R3 ) be the projection defined by Pf(x) = 1 f( x ω)dσ(ω). π S We define P = I P to be the orthogonal projection onto L := L L rad, the space of angular functions. As P and P commute with and 1, x they commute with H. We will often write H := P H = HP when convenient. The advantage of restricting to functions in L is due, in part, to an improvement to Lemma 1.1 in this setting. For a clear exposition, see []. Lemma.1 (Improved Hardy Ineuality). If f : R 3 is a Schwartz function, then 9 R3 P f (.1) dx P f dx. x R 3 Moreover, (.1) fails if the constant 9 is replaced by any constant a > 9.

4 A. ADAM AZZAM For data in L this improvement allows one to apply techniues available only in the subcritical regime a < 1 at the critical constant a = 1. In this paper, we exploit this improvement to expand upon the harmonic analysis tools developed in [8]. We seek out which embeddings enjoyed by are also enjoyed by H in order to develop a Littlewood-Paley theory for H. Towards this end, we study the resolvent of H and, in particular, how it differs from the resolvent of when projected to angular functions. In this paper we improve the resolvent estimate from [10] to prove the following: Theorem.1. Let (,s) [1, ) [1, ) satisfy (.) If, additionally, (.3) then (.) s 3 > max{ 3 s, 3, 3 5, 3 s 5 }, P ((H +1) 1 ( +1) 1 ) L x (R 3 ) L s x (R3 ) <. This estimate allows us to improve the range of Bernstein estimates enjoyed by H that are found in [10]. Indeed, we obtain that H obeys many of the useful Bernstein estimates enjoyed by the Laplacian. More precisely, we prove the following: orollary. (Sharp Sobolev Embedding). If f : R 3 is Schwartz, then (.5) P (H +1) 1 f L 6 x (R 3 ) f L x (R 3 ). orollary.3. If f : R 3 is Schwartz, then (.6) (.7) P (H +1) 1 f L x (R 3 ) f L x (R 3 ) P (H +1) 1 f L x (R 3 ) f L 1 x (R 3 ). Theorem. (BernsteinEstimates). Theoperator P e H : L 1 (R 3 ) L (R 3 ) is bounded. Remark.1. The estimates in Theorem.1, orollary.3, and Theorem. fail in the absence of the projection P. 3. Acknowledgements It is difficult for me to articulate just how much I appreciate the guidance, support, and friendship of my advisors, Rowan Killip and Monica Visan. I am incredibly grateful to them for introducing me to this problem, for their time in discussing it with me, and for carefully reading this manuscript. This work was supported in part by NSF grant DMS (P.I. Rowan Killip) and NSF grant DMS (P.I. Monica Visan).

5 RESOLVENT ESTIMATES FOR 3D NLS WITH INVERSE-SQUARE POTENTIAL 5. Preliminaries We begin by fixing some notation. We will write X Y if there exists a constant so that X Y. When we wish to stress the dependence of this implicit constant on a parameter ε (say), so that = (ε), we write X ε Y. We write X Y if X Y and Y X. If there exists a small constant c for which X cy we will write X Y. For1 r < werecall thelebesguespace L r x(r 3 ), which isthecompletion of smooth compactly supported functions f : R 3 under the norm ( f r = f L r x (R 3 ) := f(x) r dx R 3 )1 r. When r =, we employ the essential supremum norm. Our convention for the Fourier transform on R 3 is ˆf(ξ) = 1 e ix ξ f(x) dx. (π) 3 R 3 The Fourier transform allows us to define the fractional differentiation operators s f(ξ) = ξ sˆf(ξ) and s f(ξ) = (1+ ξ ) s ˆf(ξ). The fractional differentiation operators give rise to the(in-)homogenous Sobolev spaces. We define Ḣ1,r (R 3 ) and H 1,r (R 3 ) to be the completion of smooth compactly supported functions f : R 3 under the norms f Ḣ1,r (R 3 ) = f L r x and f H 1,r (R 3 ) = f L r x. When r =, we simply write H 1,r = H 1 and Ḣ1,r = Ḣ1. Let us define the operator H = 1 as the Friedrichs extension of the x symmetric and non-negative sesuilinear form (.1) Q(u,v) = u v 1 x u v dx u,v 0 (R 3 ). More explicitly, let Q be the closure of Q with domain D( Q) given by the completion of 0 (R 3 ) with respect to the norm ( u L x (R3 ) + Q(u,u))1. For each element u in D(H) = {u D( Q) : Q(u,v) u v L x (R 3 ) for all v D( Q)}, we define Hu to be the uniue element in L x (R3 ) satisfying Q(u,v) = (Hu,v) for all v D( Q). 5. ombes-thomas Estimates In this section we prove a resolvent estimate for H, of which a special case is anessential stratagem in[10]. To do this, we adopt a techniue due to ombes and Thomas (see [9], [5]) for proving decay estimates for the resolvent.

6 6 A. ADAM AZZAM Lemma 5.1 (ombes-thomas Estimate, [5]). Suppose A is a self-adjoint operator on a Hilbert space H and its spectrum, σ(a), is positive and satisfies d := dist(σ(a),0) > 0. If B is self-adjoint on H, then A+iβB is self adjoint for all β R and (A+iβB) 1 d 1. Proof. If u H, then u (A+iβB)u Re u,(a+iβb)u = Re( u,au iβ u,bu ) d u. To use Lemma 5.1 we will need a useful reformulation of Lemma.1. orollary 5.1. If f : R 3 is Schwartz and P f = f, then (5.1) and, likewise, (5.) (H )f,f 0 (H x )f,f 0. Proof. Suppose f is Schwartz. As P H = HP, by (.1) we see that (H )f,f L x (R3 ) = f,f L x (R 3 ) 1 x f,f L x (R 3 ) f,f L x (R3 ) Similarly, we see that = 1 9 ( f, f L x(r 3 ) 9 x f,f L x (R 3 )) 0. (H x )f,f L x (R 3 ) = f,f L x (R 3 ) 9 x f,f L x (R 3 ) 0. Theorem 5. (ombes-thomas Estimates). For p R, define { } (5.3) A := 1 p (H +1) 1 1 (H +1) 1 x, and { B := ip(h +1) 1 x x (5.) + }(H +1) 1. x x If p <, then A+iB is invertible. Proof. First, note that A and B are self-adjoint and σ(a) = σ(a) R +. By Lemma 5.1, A + ib is invertible provided that d = dist(σ(a),0) > 0. To determine the range of p for which d > 0, it behooves us to compute (H +1) 1 1 x (H +1) 1 L x (R 3 ) L x(r 3 ). Towards this end we note that, by (5.), x (1+H ) 1 (1+H ) 1 x = x (1+H ) x x (1+ x ) x = + x.

7 RESOLVENT ESTIMATES FOR 3D NLS WITH INVERSE-SQUARE POTENTIAL 7 It follows, then, that (5.5) (H +1) 1 x (H +1) 1 L x (R 3 ) L x (R3 ) 1. By a TT argument, it follows that (5.6) x 1 (H +1) 1 L x (R 3 ) L x (R3 ) = (H +1) 1 x 1 L x (R 3 ) L x (R3 ) 1. onseuently, (5.7) (H +1) 1 x x L x (R 3 ) L x (R3 ) 1. Thus, d 1 p so d > 0 provided p <. It follows that A + ib is invertible provided that p <, as desired. orollary 5.3 (Weighted Resolvent Estimate). The operator is bounded for all p (, ). x 1 p (H +1) x 1+p : L x (R3 ) L x (R3 ) Proof. In the notation of Lemma 5., with A as defined in (5.3) and B as defined in (5.), we see that (5.8) x 1 p (H +1) x 1+p = x (H +1) 1 {A+iB}(H +1) 1 x. By (5.6), theoperators(h +1) x and x (H 1 +1) 1 areboundedly invertible. By Lemma 5., A+iB is boundedly invertible provided that p (, ). Thus, x 1 p (H + 1) x 1+p is bounded provided that p (, ), as claimed. 6. Kernel Estimates In this section we prove kernel estimates for the operator P ( +1) 1. Lemma 6.1. Fix y R 3 and suppose x R 3. If x < y, then (6.1) If y < x, then (6.) P ( x 1±p ( +1) 1 δ y ) x ±p y e y /. P ( x 1±p ( +1) 1 δ y ) y x 3±p e x /. Proof. As x 1±p is radial it commutes with P and so (6.3) P ( x 1±p ( +1) 1 δ y ) = x 1±p P (( +1) 1 δ y ). It suffices, then, to bound (6.) P (( +1) 1 δ y )(x) = S e x y x y e ω x y ω x y dσ(ω),

8 8 A. ADAM AZZAM where dσ denotes the uniform probability measure on the unit sphere S. By the mean value theorem, for each ω S there exists some c(ω) between x y and ω x y so that (6.5) (6.6) e x y e ω x y x y ω x y = e c(ω) ( c(ω)+1 )( x y ω x y ) c(ω) e c(ω) ( c(ω)+1 c(ω) ) x ω x y x y + ω x y. Suppose for now that x < y. It follows, then, that c(ω) min{ x y, ω x y } x y > y. Since the function t e t( ) t+1 t is decreasing, we see that ( ) c(ω)+1 ( ) (6.7) e c(ω) e y / y + e y /. c(ω) y y Moreover, we see that (6.8) x ω x y x y x. x y + ω x y y ombining (6.6), (6.7), and (6.8) we get (6.9) e x y x y e y /. e ω x y x y ω x y Integrating over the sphere, we combine (6.9) with (6.3) to deduce (6.1). If instead y < x, then the left hand side of (6.9) is bounded by the right hand side of (6.9) with y replaced with x and vice versa. This gives (6.). Lemma 6.. If p < 3 and y R3, then P ( x 1±p ( +1) 1 δ y ) L x (R ) y ±p min{ y 1, y 1 }. Proof. Fix some y R 3. We split the integral into three regions ( ) (6.10) + P ( x 1±p ( +1) 1 δ y ) dx = I 1 +I +I 3, + A B where A = { x < y }, B = { x > y }, and = { y / x y }. On the region A, we employ (6.1) to see that P ( x 1±p ( +1) 1 δ y ) dx y e y / x ±p dx A provided that p < 3. A e y / y 1±p y ±p ( y y 1 ),

9 RESOLVENT ESTIMATES FOR 3D NLS WITH INVERSE-SQUARE POTENTIAL 9 On the region B, we employ (6.) to see that P ( x 1±p ( +1) 1 δ y ) dx y B B x 6±p e x / dx y e y r ±p dr y y 1±p e y y ±p ( y y 1 ), provided that p < 3. Since P : L x () L x () is bounded, we see that (6.11) P ( x 1±p ( +1) 1 δ y ) dx x ±p ( +1) 1 δ y ) dx (6.1) x ±pe x y dx x y (6.13) y ±p dx. e x y x y We may further split = 1 where 1 = { x y < y } and = \ 1. On 1 we see by a change of coordinates that (6.1) 1 e x y dx x y y 0 e r dr = 1 e y (1 y ). On we know that x y y and so e x y e y (6.15) dx 1 dx y e y (1 y ). x y y It follows, then, that P ( x 1±p ( +1) 1 δ y ) dx y ±p ( y y 1 ). Summing over the regions, we have our desired bound. orollary 6.1. The integral kernel of P ((H +1) 1 ( +1) 1 ), namely, (6.16) obeys the bound (6.17) G(z,y) = δ z,p ((H +1) 1 ( +1) 1 )δ y, G(z,y) for all p (, ). ( z y z + y ) p( z 1 z 1 )( y 1 y 1 ) Proof. Given two invertible operators A, B, we have the resolvent identity (6.18) A 1 B 1 = B 1 (B A)B 1 +B 1 (B A)A 1 (B A)B 1.

10 10 A. ADAM AZZAM As ( +1) (H +1) = 1 x, by (6.18) we see that (6.19) where (6.0) (6.1) P (H +1) 1 P ( +1) 1 = (I)+(II), (I) := 1 P ( +1) 1 x P ( +1) 1 and (II) := 1 16 P ( +1) 1 x P (H +1) 1 x P ( +1) 1. Let G (I) (z,y) and G (II) (z,y) be the integral kernels of (I) and (II), respectively, defined as above. It follows, then, from auchy-schwarz and (6.) that G (I) (z,y) x 1 p P ( +1) 1 δ z, x 1+p P ( +1) 1 δ y L x (R ) (6.) x 1 p P ( +1) 1 δ z L x (R 3 ) x 1+p P ( +1) 1 δ y L x (R 3 ) z p y p ( y 1 y 1 )( z 1 z 1 ). Since (I) is self-adjoint, G (I) is conjugate symmetric and hence obeys the same bounds with the roles of z and y reversed. So (6.3) G (I) (z,y) ( z y z + y ) p( z 1 z 1 )( y 1 y 1 ), provided that p < 3. For p <, we know by orollary 5.3 that p := x 1+p P (H +1) 1 x 1 p L x (R 3 ) L x(r 3 ) p 1. It follows then by (6.), as before, that G (II) (z,y) p x 1 p P ( +1) 1 δ z L x (R 3 ) x 1+p P ( +1) 1 δ y L x (R 3 ) (6.) p z p y p ( y 1 y 1 )( z 1 z 1 ). Again by conjugate symmetry of G (II), we see that G (II) obeys (6.3), as claimed. 7. Proofs of Main Results In this section we begin by proving Theorem.1, and conclude by deducing a number of important Bernstein estimates. Theorem.1. Let (,s) [1, ) [1, ) satisfy (7.1) If, additionally, (7.) then s 3 > max{ 3 s, 3, 3 5, 3 s 5 }, (7.3) P ((H +1) 1 ( +1) 1 ) L x (R 3 ) L s x(r 3 ) <.

11 RESOLVENT ESTIMATES FOR 3D NLS WITH INVERSE-SQUARE POTENTIAL 11 Proof. We will only treat the case,s [1,] [1,]. For if (,s) [, ) [, ) satisfy (7.1) and (7.), then (s, ) belongs to [1,] [1,] and satisfy (7.1) and (7.). If the theorem is known to be true in this case, then P ((H +1) 1 ( +1) 1 ) L s x (R3 ) L x (R 3 ) <, but then by duality it must also follow that P ((H +1) 1 ( +1) 1 ) L x (R 3 ) L s x (R3 ) <. Now, by duality, it suffices to show that sup sup G(x,y)f(x)g(y) dxdy <. g L s f Lx (R 3 ) 1 R 3 R 3 y (R3 ) 1 Suppose, then, that f L x (R 3 ) 1 and g L s y (R3 ) 1. As G depends only on x and y we may, and do, assume without loss of generality that g = g( y ) and f = f( x )are radial. There areprecisely six regions to consider, each case corresponding to an ordering of x, y, and 1 (e.g. x y 1, x 1 y, etc.). By symmetry, it suffices to treat the three cases in which we always have x y. (1) With A = {(x,y) R 3 R 3 : x y 1}, we have by orollary 6.1, for any p (, ), that G(x, y)f(x)g(y) dxdy (7.) A r 1 1 r 1 1 k k=0 j=0 = 0 j k r 1 p ρ 1+p f(r) g(ρ) dρdr r 1 p ρ 1 s +p f(r)r g(ρ)ρ s dρdr k(1 + 1 p) j(1 s + 1 s +p) f L x ( x k ) g L s y ( y j ) ( ) p+ 3 j j( 3 k 3 s ) f L x ( x ) g k L s y ( y j ) We may choose p so that 3 p < 0. Since we see that by s 3 Schur s Test ( )1 ( )1 G(x, y)f(x)g(y) dxdy f L x( x k ) g L s y ( x j ) (7.5) A since,s. k=0 f L x (R 3 ) g L s y (R3 ) 1, k=0

12 1 A. ADAM AZZAM () With B = {(x,y) R 3 R 3 : x 1 y }, we similarly have by orollary 6.1, for any p (, ), that G(x, y)f(x)g(y) dxdy (7.6) B k=0 j=0 r 1 p ρ 3 s +p f(r)r g(ρ)ρ s dρdr k(1 + 1 p) j(3 s + 1 s +p) f L x ( x k ) g L s y ( y j ). hoosing p so that both 5 3 +p > 0 and, simultaneously, 3 p < 0, s by auchy-schwarz, we see that (7.7) G(x,y)f(x)g(y) dxdy f L x (R ) g 3 L s y (R 3 ) 1. B (3) With = {(x,y) R 3 R 3 : 1 x y }, we similarly compute that G(x, y)f(x)g(y) dxdy (7.8) (7.9) 1 r r r 3 p ρ 3 +p f(r) g(ρ) dρdr r 3 p ρ 3 s +p f(r)r g(ρ)ρ s dρdr 0 0 ( ) 5 j 3 k 0 k j s +p k(5 3 s 3 ) f L x ( x k ) g L s y ( y j ) We may choose a suitable p satisfying 5 3 +p > 0. As , by s s Schur s Test we get as desired. G(x,y)f(x)g(y) dxdy 1, We now pause to enjoy a few specific dividends of Theorem.1. orollary. (Sharp Sobolev Embedding). If f : R 3 is Schwartz, then (7.10) Proof. If f : R 3 is Schwartz, then (7.11) (7.1) (7.13) P (H +1) 1 f L 6 x (R 3 ) P (H +1) 1 f L 6 x (R 3 ) f L x (R 3 ). P ( +1) 1 f L 6 x (R 3 ) + P ((H +1) 1 ( +1) 1 )f L 6 x (R 3 ) f L x (R 3 )

13 RESOLVENT ESTIMATES FOR 3D NLS WITH INVERSE-SQUARE POTENTIAL 13 by Theorem.1 since the pair (,6) satisfies = and > max{ 1, 1, 1,0}. orollary.3. If f : R 3 is Schwartz, then (7.1) (7.15) P (H +1) 1 f L x (R 3 ) f L x (R 3 ) P (H +1) 1 f L x (R 3 ) f L 1 x (R 3 ). Proof. To see the first ineuality, note that if f : R 3 is Schwartz, then P (H +1) 1 f L x (R 3 ) P ( +1) 1 f L x (R 3 ) + P ((H +1) 1 ( +1) 1 )f L x (R 3 ) f L x (R 3 ). by Theorem.1 since the pair (, ) satisfies and 3 3 > max{ 1, 1, 1, 1}. The second ineuality follows from duality. Theorem. (Bernstein Estimates). The operator P e H : L 1 L is bounded. Proof. First, write It follows that P e H = (H +1) 1 P e H (H +1) (H +1) 1 P. P e H L 1 x L x = (H +1) 1 P e H (H +1) (H +1) 1 P L 1 x L x (H +1) 1 P L x L x e H (H +1) L x L x (H +1) 1 P L 1 x L x. By (7.1), (H +1) 1 P L x L 1 <. By (7.15), (H x +1) 1 P L 1 x L <. x Lastly, e H (H +1) L x L < by the Spectral Theorem. x

14 1 A. ADAM AZZAM References [1] Jean-Marc Bouclet and Haruya Mizutani. Uniform resolvent and Strichartz estimates for Schrödinger euations with critical singularities [] Nicolas Bur, Fabrice Planchon, John G. Stalker, and A.Shadi Tahvildar-Zadeh. Strichartz estimates for the wave and Schrödinger euations with the inverse-suare potential. Journal of Functional Analysis, 03():519 59, 003. [3] K. M. ase. Singular potentials. Phys. Rev., 80: , Dec [] T. Ekholm and R. L. Frank. On Lieb-Thirring ineualities for Schrödinger operators with virtual level. ommunications in Mathematical Physics, 6(3):75 70, 006. [5] P.D. Hislop. Exponential decay of two-body eigenfunctions: A review. Electronic Journal of Differential Euations, (0):65 88, 000. [6] H. Kalf, U.-W. Schmincke, J. Walter, and R. Wüst. On the spectral theory of Schrödinger and Dirac operators with strongly singular potentials, pages Springer Berlin Heidelberg, Berlin, Heidelberg, [7] R. Killip,. Miao, M. Visan, J. Zhang, and J. Zheng. The energy-critical NLS with inverse-suare potential. Preprint arxiv: , 015. [8] R. Killip,. Miao, M. Visan, J. Zhang, and J. Zheng. Sobolev spaces adapted to the Schrödinger operator with inverse-suare potential. Preprint arxiv: , 015. [9] J.M. ombes L. Thomas. Asymptotic behaviour of eigenfunctions for multiparticle Schrödinger operators. ommun. Math. Phys, 3:51 70, [10] Haruya Mizutani. Remarks on endpoint Strichartz estimates for Schrödinger euations with the critical inverse-suare potential. Preprint arxiv: , 016. [11] Fabrice Planchon, John G. Stalker, and A. Shadi Tahvildar-Zadeh. L p estimates for the wave euation with the inverse-suare potential. Discrete and ontinuous Dynamical Systems, 9():7, 003. [1] Toshiyuki Suzuki. Solvability of nonlinear Schrödinger euations with some critical singular potential via generalized hardy-rellich ineualities. Funkcialaj Ekvacioj, 59(1):1 3, 016.