Vortex Solitons for a Class of Schrödinger Equation with Square Root Nonlinear Term ()

Weikang Chen^{}

College of Sciences, University of Shanghai for Science and Technology, Shanghai, China.

**DOI: **10.4236/apm.2020.104011
PDF HTML XML
210
Downloads
388
Views
Citations

College of Sciences, University of Shanghai for Science and Technology, Shanghai, China.

In this paper, we prove an existence theory for ring-profiled optical vortex solitons via constrained minimization, which are considered in the context of an electromagnetic light wave propagating in a nonlinear media and governed by a nonlinear Schrödinger type equation with square root nonlinear term.

Share and Cite:

Chen, W. (2020) Vortex Solitons for a Class of Schrödinger Equation with Square Root Nonlinear Term. *Advances in Pure Mathematics*, **10**, 174-180. doi: 10.4236/apm.2020.104011.

1. Introduction

In optics research, a fundamental prototype situation is that the light waves are described by a complex-valued wave function governed by nonlinear Schrödinger equations [1] - [8]. These rigorous mathematical treatments of such nonlinear problems provide more possibilities for the existence and properties of optical vortices. Our interest is motivated by the work of Lin, Belić, Petrović, Hajaiej and Chen [9], the mathematical analysis of Lin and Ren [10].

In dimensionless form, consider the following nonlinear Schrödinger equation,

$i{\partial}_{z}E+\frac{1}{2}{\nabla}_{\perp}^{2}E+f\left(I\right)E=\mathrm{0,}$ (1)

where E is the evolution of the slowly varing electric field envelope propagating in the longitudinal z-direction; ${\nabla}_{\perp}^{2}$ is the Laplace operator over the transverse plane of coordinates $\left(x\mathrm{,}y\right)$ which is perpendicular to the z-axis. The function f depends on the total field intensity, I, i.e. $I={\left|E\right|}^{2}$ , and we will concentrate henceforth on the model of the self-focusing square-root nonlinearity.

$f\left({\left|E\right|}^{2}\right)=1-\frac{1}{\sqrt{1+{\left|E\right|}^{2}}}\mathrm{,}$ (2)

which describes narrow-gap semiconductors [11] [12].

We focus on spatial optical solitons. Spatially localized solutions of (1), which do not change their intensity profile during propagation, can be described under the spatial soliton ansatz

$E\left(r\mathrm{,}\theta \mathrm{,}z\right)=u\left(r\right){\text{e}}^{i\left(n\theta +\alpha z\right)}$ (3)

where r and $\theta $ are real polar coordinates over ${\mathbb{R}}^{2}$ , and $r=\sqrt{{x}^{2}+{y}^{2}}$ , $\theta =\mathrm{arctan}\left(y/x\right)$ , $u\left(r\right)$ is the radial profile function which gives rise to the intensity of light waves, $n\in \mathbb{Z}$ is the winding number, and $\alpha \in \mathbb{R}$ is the wave propagation constant. This ansatz describe a vortex wave centered around the z-axis. Inserting (3) into (1) and in a square root nonlinear media, we arrive at the following equation

${\left(r{u}_{r}\right)}_{r}-\frac{{n}^{2}}{r}u+2ru-\frac{2ru}{\sqrt{1+{u}^{2}}}-2\alpha ru=0$ (4)

Due to the presence of the vortex core, in other words, the regularity of u at $r=0$ , we impose the condition $u\left(0\right)=0$ . Besides, such ring-like beams remain localized that allows us to mathematically impose the “boundary” condition $u\left(R\right)=0$ for $R>0$ sufficiently large, where R represents the distance from the vortex core.

Therefore, in view of (4), we can get the n-vortex equation with boundary conditions.

$\{\begin{array}{l}{\left(r{u}_{r}\right)}_{r}-\frac{{n}^{2}}{r}u+2ru-\frac{2ru}{\sqrt{1+{u}^{2}}}-2\alpha ru=0\hfill \\ u\left(0\right)=0,\mathrm{}u\left(R\right)=0\hfill \end{array}$ (5)

In this paper, we treat (5) as a nonlinear eigenvalue problem and prove the existence of positive solution pairs $\left(u\mathrm{,}\alpha \right)$ by a constrained minimization approach, with a prescribed energy flux constrained.

2. Preliminary Setting and Main Theorems

In this section, we give some basic notations and lemmas which will be used in next section. In order to approach the Equation (5), we write down the action functional ${I}_{\alpha}\mathrm{:}H\to \mathbb{R}$ defined as

${I}_{\alpha}\left(u\right)=\frac{1}{2}{\displaystyle {\int}_{0}^{R}}\left\{r{u}_{r}^{2}+\frac{{n}^{2}}{r}{u}^{2}-2\left(1-\alpha \right)r{u}^{2}+4r\sqrt{1+{u}^{2}}\right\}\text{d}r$ (6)

where $\left|n\right|\ge 1$ , H is the completion of

$X=\left\{u\in {C}^{1}\left[0,R\right]|u\left(0\right)=0=u\left(R\right)\right\}$ (7)

equipped with the inner product

$\left(u,v\right)={\displaystyle {\int}_{0}^{R}}\left\{r{u}_{r}{v}_{r}+\frac{1}{r}uv\right\}\text{d}r,\mathrm{}u,v\in H$ (8)

As a Hilbert space, H may be considered as an embedded subspace of ${W}_{0}^{\mathrm{1,2}}\left({B}_{R}\right)$ which is composed of radially symmetric functions such that any element $u\in H$ enjoys the desired property $u\left(0\right)=0$ , where ${B}_{R}\mathrm{:}=\left\{\left(x\mathrm{,}y\right)\in {\mathbb{R}}^{2}\mathrm{:}{x}^{2}+{y}^{2}\le {R}^{2}\right\}$ .

Lemma 2.1. From the inequalities

${\left(\sqrt{1+{u}^{2}}-1\right)}^{2}\le \frac{1}{4}{u}^{4}\mathrm{}\text{for}\text{\hspace{0.17em}}\text{\hspace{0.05em}}u\in \mathbb{R}$ (9)

and

${\int}_{0}^{R}}\text{\hspace{0.05em}}r{u}^{2}\text{d}r\le {R}^{2}{\displaystyle {\int}_{0}^{R}}\frac{{u}^{2}}{r}\text{d}r$ (10)

we get that there exists a constant $C>0$ , such that ${I}_{\alpha}\left(u\right)\le C{\Vert u\Vert}_{H}^{2}$ .

For convenience, we define the “energy” functional as

$\epsilon \left(u\right)=\frac{1}{2}{\displaystyle {\int}_{0}^{R}}\left\{r{u}_{r}^{2}+\frac{{u}^{2}}{r}+4r\sqrt{1+{u}^{2}}\right\}\text{d}r$ (11)

Now, we state our main theorem in this paper.

Theorem 2.2. For any parameters $\left|n\right|\ge 1$ , consider the n-vortex Equation (5) with boundary conditions, describing ring-profile vortex solitons in a square-root nonlinear media, with the prescribed energy flux $\Phi \left(u\right)={\Phi}_{0}>0$ , and $R>0$ .

1) There exists a solution pair $\left(u\mathrm{,}\alpha \right)$ with $u\left(r\right)>0$ , $r\in \left(\mathrm{0,}R\right)$ and $\alpha \in \mathbb{R}$ .

2) For $r\in \left[\mathrm{0,}R\right]$ , the energy flux $\Phi \left(u\right)={\Phi}_{0}\le \frac{1}{4}$ , and there exists no nontrivial solution, if ${n}^{2}+2{r}^{2}\alpha >0$ .

3. Existence of Vortices via Constrained Minimization

In this section, we consider the wave propagation constant $\alpha $ as a Lagrange multiplier, we prove the existence of solution of the Equation (5) with constrained minimization approach.

We rewrite the n-vortex Equation (5) as

$\{\begin{array}{l}{\left(r{u}_{r}\right)}_{r}-\frac{{n}^{2}}{r}u+2ru-\frac{2ru}{\sqrt{1+{u}^{2}}}=2\alpha ru\hfill \\ u\left(0\right)=0,\mathrm{}u\left(R\right)=0\hfill \end{array}$ (12)

Define the function I and the soliton energy flux as

$I\left(u\right)=\frac{1}{2}{\displaystyle {\int}_{0}^{R}}\left\{r{u}_{r}^{2}+\frac{{n}^{2}}{r}{u}^{2}-2r{u}^{2}+4r\sqrt{1+{u}^{2}}\right\}\text{d}r$

$\Phi \left(u\right)={\displaystyle {\int}_{0}^{2\pi}}\text{\hspace{0.05em}}\text{d}\theta {\displaystyle {\int}_{0}^{R}}\text{\hspace{0.05em}}r{u}^{2}\text{d}r=2\pi {\displaystyle {\int}_{0}^{R}}\text{\hspace{0.05em}}r{u}^{2}\text{d}r$

Thus, to get a solution of (12), it suffices to show that a solution to the following exists:

$\mathrm{min}\left\{I\left(u\right)|u\in \Lambda ,\Phi \left(u\right)={\Phi}_{0}\right\},{\Phi}_{0}>0$ (13)

where the nonempty admissible class $\Lambda $ is defined by

$\Lambda =\left\{u\left(r\right)\mathrm{}\text{isabsolutelycontinuousover}\mathrm{}\left[0,R\right],u\left(0\right)=u\left(R\right)=0,\epsilon \left(u\right)<\infty \right\}$ (14)

with $\epsilon \left(u\right)$ being defined by (11).

The proof of Theorem 2.2. 1) Using the energy flux ${\Phi}_{0}$ , we have

$I\left(u\right)\ge \frac{1}{2}\left({\displaystyle {\int}_{0}^{R}}\text{\hspace{0.05em}}r{u}_{r}^{2}\text{d}r+{n}^{2}{\displaystyle {\int}_{0}^{R}}\frac{{u}^{2}}{r}\text{d}r\right)-\frac{{\Phi}_{0}}{2\pi}$ (15)

Let be a minimizing sequence of (13). Then (15) gives the bound

(16)

where is a constant independent of m. We know the fact that the distributional derivative of u must satisfy, and the functionals I and are even. Thus, we may assume that the sequence consists of non-negative valued functions. Therefore, it is clear that we may view these functions as radially symmetric over the disk and vanishing on its boundary. Moreover, with (16) and (10), it can be seen that belongs in under the reduced norm,

(17)

Therefore, is bounded in. Without loss of generality, we get the weak convergence of to an element. Using the compact embedding for, strongly in as. Hence, u is radially symmetric as well with.

In view of (16) and Fatou’s lemma, Let be a measure space and a sequence of nonnegative measurable functions. Then the function is measurable and

(18)

we have

(19)

(20)

(21)

Therefore, from (10) and (19)-(21), we get

Following as in [13]. Let be a sequence in where. It is clear that for any, is bounded in. We may get that uniformly over as applying the compact embedding. Thus, we have for any pair, with (16),

(22)

Taking, we get

(23)

Since, the right-hand side of (23) tends to zero as. Hence,

As a consequence, the boundary condition is achieved. With (13), u is a solution to (13), and there is a real number such that satisfies (12).

Moreover, we may suppose that there is a point such that, then since is a minimum point for the function. By the uniqueness theorem of the initial value problem of ordinary differential equations, we have for all, thus contradicting the fact. Hence, for all.

2) We establish

(24)

Suppose otherwise that (24) is not valid, equivalently, , then there is a and so that for all. However,

(25)

which contradicts with. So, (24) is valid. From (24), we can find a sequence such that as and

(26)

Multiplying (5) by u, integrating over, letting. Appealing to (26), we obtain

(27)

Using, we have that

(28)

We may treat u as a radially symmetric function defined over with its support contained in the disk. Hence, from the classical inequality over, we deduce

with, we have

Therefore, when for,. as claimed.

4. Conclusion

Through the prove of the theorem 2.2, we get that the existence of positive solution pairs by a constrained minimization approach. In other words, we get the existence of ring-profiled optical vortex solitons propagating in a square-root nonlinear media. Moreover, we obtain that there is no nontrivial small-energy-flux solution satisfying, if for.

Conflicts of Interest

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

[1] |
Adhikari, S.K. (2010) Localization of a Bose-Einstein Condensate Vortex in a Bichromatic Optical Lattice. Physical Review A, 81, Article ID: 043636. https://doi.org/10.1103/PhysRevA.81.043636 |

[2] |
Davydova, T.A. and Yakimenko, A.I. (2004) Stable Multi-Charged Localized Optical Vortices in Cubicquintic Nonlinear Media. Journal of Optics A: Pure and Applied Optics, 97, S197-S201. https://doi.org/10.1088/1464-4258/6/5/010 |

[3] |
Kartashov, Y.V., Vysloukh, V.A. and Torner, L. (2004) Rotary Solitons in Bessel Optical Lattices. Physical Review Letters, 93, Article ID: 093904. https://doi.org/10.1103/PhysRevLett.93.093904 |

[4] |
Kartashov, Y.V., Vysloukh, V.A. and Torner, L. (2005) Stable Ring Vortex Soltions in Bessel Optical Lattices. Physical Review Letters, 94, Article ID: 043902. https://doi.org/10.1103/PhysRevLett.94.043902 |

[5] |
Mamaev, A.V., Saffman, M. and Zozulya, A.A. (1996) Propagation of Stripe Beams in Nonlinear Media: Snake Instability and Creation of Optical Vortices. Physical Review Letters, 76, 2262-2265. https://doi.org/10.1103/PhysRevLett.76.2262 |

[6] |
Neshev, D., Alexander, T.J., Ostrovskaya, E.A., Kivshar, Y.S., Martin, H., Makasyuk, I. and Chen, Z. (2004) Observation of Discrete Vortex Solitons in Optically-Induced Photonic Lattices. Physical Review Letters, 92, Article ID: 123903. https://doi.org/10.1103/PhysRevLett.92.123903 |

[7] |
Rozas, D., Law, C.T. and Swartzlander Jr., G.A. (1997) Propagation Dynamics of Optical Vortices. Journal of the Optical Society of America B, 14, 3054-3065. https://doi.org/10.1364/JOSAB.14.003054 |

[8] |
Skryabin, D.V. and Firth, W.J. (1998) Dynamics of Self-Trapped Beams with Phase Dislocation in Saturable Kerr and Quadratic Nonlinear Media. Physical Review E, 58, 3916-3930. https://doi.org/10.1103/PhysRevE.58.3916 |

[9] |
Lin, T.C., Belić, M.R., Petrović, M.S., Hajaiej, H. and Chen, G. (2017) The Virial Theorem and Ground State Energy Estimates of Nonlinear Schrödinger Equations in with Square Root and Saturable Nonlinearities in Nonlinear Optics. Calculus of Variations and Partial Differential Equations, 56, Article No. 147. https://doi.org/10.1007/s00526-017-1251-4 |

[10] |
Liu, C. and Ren, Q. (2015) On the Steady-State Solutions of a Nonlinear Photonic Lattice Model. Journal of Mathematical Physics, 53, Article ID: 031501. https://doi.org/10.1063/1.4914333 |

[11] |
Petroskia, M.M., Petrovica, M.S. and Belica, M.R. (2007) Quasi-Stable Propagation of Vortices and Soliton Clusters in Previous Saturable Kerr Media with Square-Root Nonlinearity. Optics Communications, 279, Article ID: 196202. https://doi.org/10.1016/j.optcom.2007.07.006 |

[12] |
Skarka, V., Berezhiani, V.I. and Boucher, V. (2000) Self-Trapping of Electromagnetic Pulses in Narrow-Gap Semicon-Ductors. Physical Review B, 61, Article ID: 10201. https://doi.org/10.1103/PhysRevB.61.10201 |

[13] |
Yang, Y. and Zhang, R. (2014) Existence of Optical Vortices. SIAM Journal on Mathematical Analysis, 46, 484-498. https://doi.org/10.1137/120894105 |

Journals Menu

Copyright © 2020 by authors and Scientific Research Publishing Inc.

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