Nonlinear Differ. Equ. Appl. (3), 5 5 c Springer Basel AG -97/3/35- published online April 3, DOI.7/s3--3- Nonlinear Differential Equations and Applications NoDEA Nonlinear Schrödinger problems: symmetries of some variational solutions Christopher Grumiau Abstract. In this paper, we are interested in the nonlinear Schrödinger problem Δu + Vu= u p u submitted to the Dirichlet boundary conditions. We consider p> and we are working with an open bounded domain R N (N ). Potential V satisfies max(v,) L N/ () and min(v,) L + (). Moreover, Δ+V is positive definite and has one and only one principal eigenvalue. When p, we prove the uniqueness of the solution once we fix the projection on an eigenspace of Δ+V.It implies partial symmetries (or symmetry breaking) for ground state and least energy nodal solutions. In the literature, the case V has already been studied. Here, we generalize the technique at our case by pointing out and explaining differences. To finish, as illustration, we implement the (modified) mountain pass algorithm to work with V negative, piecewise constant or not bounded. It permits us to exhibit direct examples where the solutions break down the symmetries of V. Mathematics Subject Classification. Primary 35J; Secondary 35A3. Keywords. Nonlinear Schrödinger problems, Ground state solutions, Least energy nodal solutions, (nodal) Nehari set, Mountain pass algorithm.. Introduction Let N,p>,λ> and an open bounded domain R N. We study the nonlinear Schrödinger problem Δu(x)+V (x)u(x) =λ u(x) p u(x) () submitted to the Dirichlet boundary conditions (DBC). We are interested in the symmetry of solutions. C. Grumiau is partially supported by a grant from the National Bank of Belgium and by the program Qualitative study of solutions of variational elliptic partial differential equations. Symmetries, bifurcations, singularities, multiplicity and numerics of the FNRS, project..55..f of the Fonds de la Recherche Fondamentale Collective.
5 C. Grumiau NoDEA When V belongs to L N/ (), the solutions can be defined as the critical points of the energy functional E p : H () R : u u + Vu λ u p. p Clearly, is solution. Concerning other solutions, if we assume that Δ+ V is positive definite and V := min(v,) L + (), then the norm u = u +Vu defined on H () is equivalent to the traditional norm u = H u (see Proposition.). By working in the same way as in [8], it directly implies the existence of ground state solutions (g.s.) and least energy nodal solutions (l.e.n.s.); i.e. one-signed (resp. sign-changing) solutions with minimal energy. These solutions are characterized as minima of E p respectively on the (resp. nodal) Nehari set { N p := u H () \{} } u + Vu = λ u p (resp. M p := {u : u ± N p }). The Morse index is (resp. ). In this paper, we study the structure of these two types of solutions. We verify whether they are odd or even with respect to the hyperplanes leaving V invariant (i.e. V respects an orthogonal symmetry with respect to the hyperplane). When it is the case, we say that the solution respects the symmetries of V. When V, this type of questions has already been studied. First, on the square in dimension, we can mention a result of Arioli and Koch (see []). They proved the existence of a positive symmetric C -function w such that Δu = wu 3 possesses a non-symmetric positive solution (with as Morse index). The same kind of result has also been obtained for a solution with as Morse index. The proof is partially computer-assisted. Second, in collaboration with Bonheure, Bouchez, Troestler and Van Schaftingen (see [, 3, 7]), we proved for p close to that the symmetries are related to the symmetries of eigenfunctions of Δ. We generalize here the technique at some non-zero potentials V and we make numerical experiments to illustrate it. In Sect., by denoting λ i (resp. E i ) the distinct eigenvalues (resp. eigenspaces) of Δ+V with DBC in H (), we prove the following Theorem.. For this, we assume that λ is the unique principal eigenvalue, i.e. an eigenvalue with a related eigenspace of dimension possessing an one-signed eigenfunction. We also require that eigenfunctions in E have a nodal line of measure. By using a maximum principle, these assumptions are satisfied at least when V L + () (see []). Theorem.. When V L N/ (),V L + () and Δ +V is positive definite such that λ is the unique principal eigenvalue, for p close to, the ground state (resp. least energy nodal) solutions respect the symmetries of their orthogonal projections in H () on E (resp. E ). In particular, when the eigenspace has a dimension, the solutions respect the symmetries of V. As we assumed that λ is the unique principal eigenvalue, ground state solutions respect the symmetries of V.
Vol. (3) Nonlinear Schrödinger problems 53 Depending on the structure of E, some symmetry breaking exist for l.e.n.s. (see Sect..3). In fact, by a traditional bootstrap, a family of ground state (resp. least energy nodal) solutions (u p ) p> converges for C -norm to functions in E (resp. E ). So, for l.e.n.s., u p does not respect the symmetries of V for p small when the projection is not symmetric in E (see Sect. 3 for an example). For larger p, it is depending on the case. In Sect..3, we exhibit rectangles and V (such that the eigenfunctions in E are symmetric) where l.e.n.s. do not respect symmetries of V for p large enough. So, the result. cannot be extended to all p. In Sect. 3, as illustration, we implement the (modified) mountain pass algorithm (see [, 9, ]) to study the cases of V negative constant, piecewise constant or singular. We exhibit direct examples such that the solutions break down the symmetries of V.. Main results The proofs are related to the technique defined in [3]. This is why we just point out and explain the differences and we do not make all the details. The first result implies that the traditional Poincaré s and Sobolev s inequalities are available for := + V ( ). Proposition.. If Δ +V is positive definite, V + L N/ () and V L + (), the norm u := u + Vu and the traditional norm u := H u are equivalent. Proof. Using the Sobolev s inequalities on Vu and as V L N/ (), C > such that u u + C u. Using the Poincaré s inequalities and as V L + (), C > and a real K such that u = ε u +( ε) u + ε Vu ε u + (( ε)c + εk) u ε u, where the last inequality is obtained for ε small enough. Then, the proof of Theorem. is based on two main results. The first one shows that, for p, a priori bounded solutions can be distinguished by their projections on E i... Abstract symmetry Lemma.. There exists ε> such that if a(x) λ i L N/ <εand u solves Δu + Vu= a(x)u with DBC then u =or P Ei u. Proof. Similar as in Lemma 3. in[3], Poincaré s and Sobolev s inequalities are adapted using Proposition..
5 C. Grumiau NoDEA Then, we directly obtain our abstract symmetry result as in the proof of Proposition 3. in[3]. Let us remark that the result holds for any i and not just for i = or as stated in [3]. We denote by B(,M) the ball in H () centered at and radius M. Proposition.3. Let M>. Fori N, p > such that, for p (, p), if u p,v p {u B(,M):P Ei u / B(, M )} solve the boundary value problem with DBC Δu + Vu= λ i u p u then P Ei u p = P Ei v p implies u p = v p. These two results permit us to conclude as in Theorem 3. in[3]. Theorem.. Let (G α ) α E with E = E i be a group acting on H () such that, for g G α and u H (), g(e) =E, g(e )=E, gα = α and E p (gu) =E p (u). For any M >, p > such that, for any family of solutions (u p ) p>p> {u B(,M):P Ei u / B(, M )} of the boundary value problem with DBC Δu + Vu= λ i u p u, u p belongs to the invariant set of G αp where α p is the orthogonal projection P E u p. Theorem. can be used for any bounded family of solutions staying away from. To apply Theorem. at a family (u p ) p> of ground state (resp. least energy nodal) solutions for the problem (), we study the asymptotic behavior when p. We prove that the expected upper and lower bounds are fine if and only if λ = λ (resp. λ ). In some sense, λ (resp. λ )isthe natural rescaling to work with ground state (resp. least energy nodal) solutions of problem (). Let us remark that this condition is not a restriction. Indeed, by homogeneity of λ u p u in Eq. (), the symmetries of ground state (resp. least energy nodal) solutions are independent of λ... Asymptotic behavior Let us denote (u p ) p> a family of ground state (resp. least energy nodal) solutions for the problem (). We consider λ = λ for g.s. (resp. λ n the first not principal eigenvalue for l.e.n.s.) and E = E (resp. E n ). Lemma.5. Concerning the upper bound, lim sup p u p = lim sup ( ) p Ep(u p) / /p u where u E minimizes the limit functional E : E R : u u log u. Proof. The proof is inspired by Lemma. in[3]. First, we define v p := u + (p )w where w H () solves the problem Δw+Vw λ w = λ u log u with P E w =. Then, we prove that the projection of v p on N p (resp. M p ) converges when p. Concerning l.e.n.s., in [3], the result has been stated for n =. Here, let us remark that it works with E n which is not specially E. We just need to ensure that v p is sign-changing for p close to. Nevertheless, we need to assume n = to obtain the lower bound. It is explained in the next Lemma.
Vol. (3) Nonlinear Schrödinger problems 55 Lemma.. Concerning the lower bound, if n =then lim inf p u p >. Proof. The proof is inspired by Lemma. in[3]. Concerning l.e.n.s. (the argument is easier for g.s.), let e be a first eigenfunction in E. By considering s p := u+ p e and s+ up e p := s p, we show the existence of t p > such that v p = t p (s + p u + p + s p u p ) belongs to M p E. Then, we prove that v p stays away from zero using Poincaré s and Sobolev s embeddings, which concludes the proof. For this part, we need to require that λ is the unique principal eigenvalue, i.e. n =. Otherwise, we should prove that v p (E E n ), which cannot be assumed. The two previous results imply Theorem.7. Theorem.7. Assume that Δ+V is positive definite and possesses one and only one principal eigenvalue (n =), V + L N/ () and V L + (). If (u p ) p> is a family of ground state (resp. least energy nodal) solutions for Eq. () then C > such that u p H C ( ) λ i p λ, where i =(resp. ). If p n and ( ) λ i pn λ u pn u in H (), then ( ) λ i pn λ u pn u in H (),u satisfies Δu + Vu = λ i u and E (u ) = inf{e (u) : u E i \{}, de (u),u =}, where E : E i R : u λi u u log u. Remark.8. [(i)] () By a traditional bootstrap, ( ) λ i pn λ u pn u implies ( ) λ i pn λ u pn u for C -norm (see [7]). () If λ<λ (resp. λ ), a family of g.s. (resp. l.e.n.s.) blows up in H (). If λ>λ (resp. λ ), it goes to. So, a family of g.s. (resp. l.e.n.s.) is bounded and stays away from if and only if λ = λ (resp. λ ). (3) By homogeneity of λ u p u, the study of symmetries for only one value of λ is enough to conclude symmetries for any λ>. () By combining Theorems. and.7, we obtain that ground state solutions for p close to respect the symmetries of their projection on E.As first eigenfunctions are unique up to a constant, they keep symmetries of V for p close to. (5) By combining Theorems. and.7, we obtain that l.e.n.s. for p close to respect the symmetries of their projection on E..3. Symmetry breaking for least energy nodal solutions For p, previous results showed that the structure of l.e.n.s. are related to the symmetries of u verifying E (u )=inf{e (u) :u E i \{}, de (u),u =}. In [3] (see Section ), on the square and for V, it is proved that if u does not respect the symmetries of the rectangle, i.e. u is not odd or even with respect to the medians (which is numerically observed), then there exists a symmetry breaking on rectangles sufficiently close to the square and p sufficiently large. In our case, this property can be stated as follows. Theorem.9. Let us work on a square. If V is odd or even with respect to a median but u does not respect this symmetry, then there exist some rectangles and p such that l.e.n.s. u p does not respect the symmetries of V.
5 C. Grumiau NoDEA Moreover, as u p converges for C -norm, we are able to directly construct V such that the l.e.n.s. break down the symmetry of V. It will happen once u is not symmetric. In the next section numerical experiments illustrate this interesting case. 3. Numerical illustrations: non-zero potentials V In this section, we compute the (resp. modified) mountain pass algorithm to approach one-signed (resp. sign-changing) solutions (see [, 5, 9, ]). While it is not sure that approximate solutions have least energy, all the other solutions that we have found numerically have a larger energy. So, we will assume that the approximations are ground state (resp. least energy nodal) solutions. We also give some level curves: and for g.s., ± and± for l.e.n.s. Numerically, we study p =. Let us remark that we always obtain the same kind of symmetry for smaller values of p. We work with p = to illustrate that the result of Theorem. seems to hold at least for a non-negligible interval. 3.. Negative constant potential on a square As first example, we consider a constant V such that λ >, i.e. V > λ where λ is the first eigenvalue of Δ. So, the required assumptions on V are clearly satisfied. Theorem. holds. In particular, concerning symmetries, we obtain () for p close to, on convex domains, ground state solutions are even with respect to each hyperplane leaving invariant; () for p close to, l.e.n.s. on a rectangle are even and odd with respect to a median; (3) for p close to, l.e.n.s. on radial domains are even with respect to N orthogonal directions and odd with respect to the orthogonal one; () for p close to, l.e.n.s. on a square are odd with respect to the barycenter. Numerically, we consider Δu π u = u3 defined on the square = (, ) in R. First and second eigenvalues of Δ are given by π and 5π. On the following graph, one-signed (resp. nodal) numerical solutions have the expected symmetries. Ground state solutions respect the symmetries of the square and l.e.n.s. are odd with respect to the center. Moreover, the nodal line of the sign-changing solutions seems to be a diagonal, as for V (see [3]).
Vol. (3) Nonlinear Schrödinger problems 57 y x For g.s.: max(u) =.8, E (u) =.5 For l.e.n.s.: min(u) =., max(u) =., E (u) =33. Starting function for g.s.: (x )(y )(x + )(y +) Starting function for l.e.n.s.: sin(π(x + )) sin(π(y + )) 3.. Piecewise constant potential on a rectangle As second example, V is piecewise constant on (, ) (, ). In [], it is proved that there exists just one principal eigenvalue. So, our assumptions are satisfied and Theorem. is available. For λ =,p =andv (x, y) =V := when x<(resp.v + := otherwise), following graphs indicate that approximations are just even with respect to a direction. Ground state solutions are more 8 5 3..8..8.....8 8.5..8....8. y x
58 C. Grumiau NoDEA or less located in x< (the side minimizing energy) and respect symmetries of V. l.e.n.s. seem to be formed by g.s. on each nodal domains and are even with respect to a direction. For g.s.: max(u) =5.98, E (u) =3.98 For l.e.n.s.: min(u) = 8.7, max(u) =.53, E (u) =7.3 Starting function for g.s.: (x )(y )xy Starting function for l.e.n.s.: sin(π(x + )) sin(π(y + )) If V =andv + = 35, we get the same symmetry for g.s. but l.e.n.s. are not symmetric. The mass is more or less located in the square defined by x<. So, we obtain a direct symmetry breaking. To minimize the energy, the difference in the potential is so large that it is better to locate the mass in one side of the rectangle. On a square and for V =, it is conjectured that l.e.n.s. is odd with respect to a diagonal. It explains the structure of the approximation. 8 7 5 3..8..8...8.. 8.5..8....8. y x For g.s.: max(u) =.9, E (u) =33. For l.e.n.s.: min(u) = 9.8, max(u) =9.7, E (u) = 89 Starting function for g.s.: (x )(y )xy Starting function for l.e.n.s.: sin(π(x + )) sin(π(y + )) 3.3. A singular potential on a ball As last example, we study singular potentials. First, λ =,p = and V (x, y) = on the ball B(, ). Approximations show as expected that x +y the ground state solutions are radial and the l.e.n.s. are odd and even with respect to a diagonal. We obtain the same symmetry as for the potential V = (see [3]).
Vol. (3) Nonlinear Schrödinger problems 59 8.5. 3.5 3..5..5.5 8 y x For g.s.: max(u) =.5, E (u) =9.9 For l.e.n.s.: min(u) =.3, max(u) =.3, E (u) =7. Starting function for g.s.: cos(π(x + y ).5 /) Starting function for l.e.n.s.: cos(π(x + y ).5 /) cos(π(x + y ).5 cos(π(x + y ).5 Second, V (x, y) = (x.5) +y on the ball B(, ). Ground state solutions seem to be even with respect to a direction but are not radial. One can remark the work of the singularity on the level curve. l.e.n.s. are just odd with respect to a direction. The mass is a little bit attracted by the side x< (the side minimizing the energy). 8.5. 3.5 3..5..5.5 8 y (.5,) x
5 C. Grumiau NoDEA For g.s.: max(u) =., E (u) =8.7 For l.e.n.s.: min(u) =.5, max(u) =.5, E (u) =7.3 Starting function for g.s.: cos(π(x + y ).5 /) Starting function for l.e.n.s.: cos(π(x + y ).5 /) cos(π(x + y ).5 cos(π (x + y ).5 References [] Arioli, G., Koch, A.: Non-symmetric low-index solutions for a symmetric boundary value problem. Preprint () [] Bonheure, D., Bouchez, V., Grumiau, C.: Asymptotics and symmetries of ground-state and least energy nodal solutions for boundary-value problems with slowly growing superlinearities. Differ. Integral Equ. (9 ), 7 7 (9) [3] Bonheure, D., Bouchez, V., Grumiau, C., Van Schaftingen, J.: Asymptotics and symmetries of least energy nodal solutions of Lane-Emden problems with slow growth. Commun. Contemp. Math. (), 9 3 (8) [] Choi,.S., McKenna, P.J.: A mountain pass method for the numerical solution of semilinear elliptic problems. Nonlinear Anal. (), 7 37 (993) [5] Costa, D.G., Ding,., Neuberger, J.M.: A numerical investigation of sign-changing solutions to superlinear elliptic equations on symmetric domains. J. Comput. Appl. Math. 3( ), 99 39 () [] Godoy, T., Gossez, J.-P., Paczka, S.R.: A minimax formula for the principal eigenvalues of Dirichlet problems and its applications. In: Proceedings of the International Conference in honor of Jacqueline Fleckinger. Electronic Journal of Differential Equations, Conference, vol., pp. 37 5. Texas State University-San Marcos, Department of Mathematics, San Marcos, T (7) [7] Grumiau, C., Troestler, C.: Oddness of least energy nodal solutions on radial domains. Electron. J. Differ. Equ. Conf. 8, 3 3 () [8] Neuberger, J.M.: A numerical method for finding sign-changing solutions of superlinear Dirichlet problems. Nonlinear World (), 73 83 (997) [9] hou, J., Li,.: A minimax method for finding multiple critical points and its applications to semilinear elliptic PDE s. SIAM Sci. Comput. 3, 8 85 () [] hou, J., Li,.: Convergence results of a local minimax method for finding multiple critical points. SIAM Sci. Comput., 85 885 ()
Vol. (3) Nonlinear Schrödinger problems 5 Christopher Grumiau Institut de Mathématique Université de Mons Le Pentagone, Place du Parc 7 Mons Belgium e-mail: Christopher.grumiau@umons.ac.be Received: October. Accepted: March.