Semidiscrete optical vortex droplets in quasi-phase-matched photonic crystals

Xiaoxi Xu; Feiyan Zhao; Jiayao Huang; Hexiang He; Hexiang He; Li Zhang; Zhaopin Chen; Zhongquan Nie; Boris A. Malomed; Boris A. Malomed; Yongyao Li; Yongyao Li; Yongyao Li

doi:10.1364/OE.506130

1. Introduction

Semidiscrete vortex quantum droplets, a new type of vortices, were initially predicted in binary Bose-Einstein condensates trapped in an array of tunnel-coupled quasi-1D potential wells [1]. Unlike vortex modes in fully continuous or fully discrete systems [2–7], these are stripe-shaped localized states, which are continuous in one direction and discrete in the perpendicular one, and do not exhibit rotational symmetry. It is well known that the stability of self-trapped vortex modes in two-dimensional (2D) and three-dimensional (3D) geometries is a challenging problem because the self-attractive nonlinearity gives rise to strong splitting instability of vortex rings and tori, even if the collapse instability that affects fundamental (zero-vorticity) solitons in the same media may be suppressed [8–13]. Due to the competition between the mean-field (MF) and beyond-MF effects in the bosonic condensate [14–23], semidiscrete vortex quantum droplets may maintain stability in this setting against the azimuthal (splitting) perturbations. The droplets exhibit significant differences from traditional solitons. The droplet state features a flat-top structure, characterized by an almost constant density distribution, bounded by relatively sharp edges. In the field of nonlinear optics, somewhat similar objects in the form of “photon droplets" were experimentally demonstrated in optical media with nonlocal (thermal) nonlinearity [24,25]. Actually, the competition between different nonlinear terms is a common effect [26–32], which occurs in the propagation of high-power laser beams in various media [32,33]. Optical semidiscrete vortex droplets can be maintained by the balance between the competing nonlinearities. In particular, stable self-bound semidiscrete vortex modes in the spatial domain were predicted in coupled planar waveguides with the cubic-quintic nonlinearity [34]. Similarly, self-bound spatiotemporal vortex modes can be predicted in coupled arrays of nonlinear fibers [35].

Recently, patterned quasi-phase matched (QPM) nonlinear photonic crystals in the 3D space have been produced by means of the thermoelectric field polarization [36], laser erasing [37], and femtosecond laser poling technique [38], which provides more possibilities for the creation of vortex states. The QPM technique has developed to a well-known method for achieving accurate phase matching in $\chi ^{(2)}$ crystals for the nonlinear frequency conversion [39–49] and nonlinear beam shaping [50–57] in different dimensions. Very recently, stable vortex solitons were predicted in 3D QPM photonic crystals [58]. The structure of the vortex solitons can be engineered by fixing different phase-matching conditions in different cells of the photonic crystals, thus inducing effective discreteness in this 3D optical medium. This technique offers a possibility of building semidiscrete vortex modes in the QPM-structured bulk photonic crystals. If the input power of the light field increases up to tens of GW/cm$^2$, third-order nonlinear effects are excited [32,33]. However, if one aims to induce fifth-order nonlinear effects, as demonstrated in Ref. [28], a higher power or a careful selection of specific optical materials is necessary.

The objective of this work is the construction of semidiscrete vortex optical droplets. This type of nontrivial self-bound vortices are able to yield more detailed information compared to the vortex self-localization solitons outlined in Ref. [27]. It is relevant to mention that effective 2D (but not 3D) discrete waveguiding structures for optical beams with the extraordinary polarization can be induced by means of a different technique in photorefractive materials illuminated by interfering beams with the ordinary polarization [59]. Such virtual photonic lattices were used for the creation of quasi-discrete 2D vortex solitons [60–62].

In this paper, we propose a scenario for the creation of semidiscrete vortex optical droplets in 3D photonic crystals with quadratic ($\chi ^{(2)}$) and defocusing cubic ($\chi ^{(3)}$) nonlinearities and a striped QPM-induced spatial structure, as shown in Fig. 1(a,b,c). As mentioned above, the $\chi ^{(3)}$ nonlinearity is able to compete with the $\chi ^{(2)}$ interaction if intensities of the light fields reach the level of several GW/cm$^{2}$ [32], making it possible to create self-bound photon droplets. We demonstrate that the phase-matching condition, adjusted to the striped structure, induces effective discreteness between adjacent stripes, which is necessary for the design of semidiscrete states. The balance of the competition between quadratic and cubic nonlinearity allows the self-bound vortex modes to maintain their stability.

Fig. 1. (a,b) The structures corresponding to OC and IC modulations, which are defined as per Eq. (3), $\ell$ and $\mathbf {O}$ representing the stripe’s width and center of the modulation pattern, respectively. The black and gray blocks represent $\sigma =-1$ and $+1$, respectively. (c) The periodic modulation along the $Z$ axis is defined as per Eq. (4), $\Lambda$ being the period of the longitudinal modulation. (d) A schematic of the experimental setup for the creation of the semidiscrete vortex droplets: L1, L2, L3 – lenses, SLM1, SLM2 – spatial light modulators, PPLN – the periodically polarized lithium niobate crystal, FF – the fundamental frequency, SH – the second harmonic. Two bottom plots display the input and output intensity and phase patterns of the FF and SH beams.

Download Full Size | PDF

The subsequent material is arranged as follows. The model is introduced in Section 3. Numerical results and estimates for the experimental setup are presented in Section 4. The paper is concluded by Section 5.

2. Model

The paraxial propagation of light beams through the 3D QPM photonic crystals with the competing $\chi ^{(2)}$ and $\chi ^{(3)}$ nonlinearities is governed by coupled equations for the slowly varying fundamental frequency (FF) and second harmonic (SH) amplitudes, $A_{1}$ and $A_{2}$:

(1)$$\begin{aligned} &i\partial _{Z}A_{1}={-}\frac{1}{2k_{1}}\nabla ^{2}A_{1}-\frac{2d(X,Z)\omega _{1}}{cn_{1}}A_{1}^{{\ast} }A_{2}e^{{-}i\Delta k_{0}Z}+\frac{3\chi ^{(3)}\omega _{1}}{2cn_{1}}(|A_{1}|^{2}+2|A_{2}|^{2})A_{1},\\ &i\partial _{Z}A_{2}={-}\frac{1}{2k_{2}}\nabla ^{2}A_{2}-\frac{d(X,Z)\omega _{2}}{cn_{2}}A_{1}^{2}e^{i\Delta k_{0}Z}+\frac{3\chi ^{(3)}\omega _{2}}{ 2cn_{2}}(|A_{2}|^{2}+2|A_{1}|^{2})A_{2}, \end{aligned}$$

where $\nabla ^{2}=\partial _{X}^{2}+\partial _{Y}^{2}$ is the paraxial-diffraction operator, $c$ is the speed of light in vacuum, while $n_{1,2}$, $\omega _{1,2}$ ($\omega _{2}=2\omega _{1}$), and $k_{1,2}$ are, respectively, the refractive indices, carrier frequencies, and wavenumbers of the FF and SH components, and $\Delta k_{0}=2k_{1}-k_{2}$ is the phase-velocity mismatch. $\chi ^{(3)}>0$ is the third-order susceptibility, which accounts for the cubic self-defocusing. The local modulation of the second-order susceptibility $\chi ^{(2)}$ is determined by

(2)$$d(X,Z)=\sigma (X)d(Z),$$

where $\sigma (X)$ is the transverse striped OC (onsite-centered) or IC (intersite-centered) modulation pattern:

(3)$$\sigma (X)= \begin{cases} -\mathrm{sgn}\left[ \cos (\pi X/\ell )\right] & \mathrm{OC}, \\ -\mathrm{sgn}\left[ \sin (\pi X/\ell )\right] & \mathrm{IC}, \end{cases}$$

$\ell$ being the width of a stripe. The OC and IC patterns correspond to the pivot of the vortex beam located, respectively, at the center of a stripe or at the border between two stripes, see Fig. 1(a,b). Further, the factor accounting in Eq. (2) for the modulation in the $Z$ direction, with amplitude $d_{0}$ and period $\Lambda$, is [63–65]

(4)$$d(Z)=d_{0}\mathrm{sgn}\left[ \cos (2\pi Z/\Lambda )\right] \equiv d_{0}\sum_{m\neq 0}\left( \frac{2}{m\pi }\right) \sin \left( \frac{m\pi }{2} \right) \exp \left( i\frac{2\pi m}{\Lambda }Z\right) ,$$

see Fig. 1(c). Actually, only the terms with $m=\pm 1$ are kept in Eq. (4), as they play the dominant role in the QPM effect. Thus, $m=1$ and $-1$ relate to the FF and SH components, respectively.

By means of rescaling [66,67]

(5)$$\begin{aligned} &\zeta =\left( \frac{n_{1}}{\omega _{1}}+\frac{n_{2}}{\omega _{2}}\right) ,\quad z_{d}^{{-}1}=\frac{2}{c\pi }\frac{d_{0}^{2}}{\chi ^{(3)}}\left( \frac{ \omega _{1}^{2}\omega _{2}}{n_{1}^{2}n_{2}}\zeta \right) ^{\frac{1}{2}},\\ &u_{p}=\frac{\chi ^{(3)}}{d_{0}}\sqrt{\frac{n_{p}}{\omega _{p}\zeta }} A_{p}\exp \left[ i(\Delta k_{0}-2\pi /\Lambda )Z\right] ,\quad p=1,2,\\ &z=Z/z_{d},\quad x=\sqrt{k_{1}/z_{d}}X,\quad y=\sqrt{k_{1}/z_{d}}Y,\quad \Omega =\left( \Delta k_{0}-2\pi /\Lambda \right) z_{d},\\ &g_{11}=\frac{3\pi }{4}\sqrt{\frac{\omega _{1}^{2}n_{2}}{n_{1}^{2}\omega _{2}}\zeta },\quad g_{22}=\frac{3\pi }{4}\sqrt{\frac{\omega^{3}_{2}n_{1}^{2}}{ n_{2}^{3}\omega _{1}^{2}}\zeta }\quad g_{12}=\frac{3\pi }{2}\sqrt{\frac{ \omega _{2}}{n_{2}}\zeta }, \end{aligned}$$

Eq. (1), which keep, as said above, the terms with $m=\pm 1$ in Eq. (4), are simplified to the form of

(6)$$\begin{aligned} i\partial _{z}u_{1}={-}\frac{1}{2}\nabla ^{\prime 2}u_{1}-\Omega u_{1}-2\sigma (x)u_{1}^{{\ast} }u_{2}+\left( g_{11}|u_{1}|^{2}+g_{12}|u_{2}|^{2}\right) u_{1}, \end{aligned}$$

(7)$$\begin{aligned} i\partial _{z}u_{2}={-}\frac{1}{2\eta }\nabla ^{\prime 2}u_{2}-\Omega u_{2}-\sigma (x)u_{1}^{2}+\left( g_{22}|u_{2}|^{2}+g_{12}|u_{1}|^{2}\right) u_{2}, \end{aligned}$$

where $\nabla ^{^{\prime }2}=\partial _{xx}+\partial _{yy}$ and $\eta =k_{2}/k_{1}$. Neglecting the slight difference in the FF and SH refractive indices, i.e., setting $n_{1}=n_{2}$, results in coefficients $g_{11}=3\pi \sqrt {3}/8$, $g_{22}=g_{12}=4g_{11}$ and $\eta =2$ in Eqs. (6) and (7). In Eqs. (6) and (7), if we transform the fundamental frequency wave and the second harmonic: $u_{1}=\tilde {u}_{1}, u_{2}=\sqrt {2}\tilde {u}_{2}$ and multiply Eq. (7) by 2 in its both sides, the coefficients of third-order nonlinear terms changes to $\tilde {g}_{11}=g_{11}: \tilde {g}_{12}=2g_{12}: \tilde {g}_{22}=4g_{22}$. The ratio of them is $\tilde {g}_{11}:\tilde {g}_{12}:\tilde {g}_{22}=1: 8: 16$, which is consistent with those reported in Ref. [26].

Equations (6) and (7) conserve two dynamical invariants, viz., the Hamiltonian and total power (alias the Manley-Rowe invariant [68–70]):

(8)$$\begin{aligned} H=\int \int \ \left( \mathcal{H}_{k}+\mathcal{H}_{\Omega }+\mathcal{H}_{2}+ \mathcal{H}_{3}\right) dxdy, \end{aligned}$$

(9)$$\begin{aligned} P=\iint {\left( |u_{1}|^{2}+2|u_{2}|^{2}\right) }dxdy\equiv P_{1}+P_{2}, \end{aligned}$$

where $\mathcal {H}_{k}=\frac {1}{2}|\nabla u_{1}|^{2}+\frac {1}{2\eta }|\nabla u_{2}|^{2}$, $\mathcal {H}_{\Omega }=-\Omega \left ( |u_{1}|^{2}+|u_{2}|^{2}\right )$, $\mathcal {H}_{2}=\sigma (x)\left ( u_{1}^{\ast 2}u_{2}+\mathrm {c.c.}\right )$ and $\mathcal {H}_{3}=\frac {1}{2}g_{11}|u_{1}|^{4}+g_{12}|u_{1}|^{2}|u_{2}|^{2}+\frac {1}{2}g_{22}|u_{2}|^{4}$. The power-sharing between the FF and SH components is defined as the ratio $r=P_{1}/P_{2}$. Control parameters for the subsequent analysis are $P$, $\ell$, and $\Omega$ (the total power, the stripe’s width, and the scaled detuning).

3. Results

3.1 Numerical results

Stationary solutions to Eqs. (6) and (7) with propagation constant $\beta$ were looked for as

(10)$$u_{p}\left( x,y,z\right) =\phi _{p}\left( x,y\right) \mathrm{e}^{ip\beta z},\quad p=1,2$$

where $\phi _{1,2}$ are the stationary amplitudes of the FF and SH components. Vortex solutions were generated by means of the imaginary-time (imaginary-$z$) method, applied to Eqs. (6) and (7), with the input taken at $z=0$ as

(11)$$\begin{aligned} \phi _{1}=r^{|S|}\exp \left( -\alpha _{1}r^{2}+iS\theta \right) , \end{aligned}$$

(12)$$\begin{aligned} \phi _{2}=r^{2|S|}\exp \left[ -\alpha _{2}r^{2}+i2S\tilde{\theta}(x)\right] , \end{aligned}$$

where $r$ and $\theta$ are the 2D polar coordinates, and

(13)$$\tilde{\theta}(x)=\theta +\frac{1}{4S}[\sigma (x)-1]\pi ,$$

with $S$ and $2S$ representing the winding numbers of FF and SH components, respectively. In Eq. (13), the matching condition between the phases of the FF and SH components, $\varphi _{1,2}(x,y)\equiv \arg \{\phi _{1,2}\}$, is defined by setting

(14)$$\varphi _{2}(x,y)=2\varphi _{1}(x,y)-\varphi _{d}(x,y),$$

with $\varphi _{d}=0$ and $\pi$ corresponding, respectively, to $\sigma (x)=1$ and $-1$ (i.e., $\varphi _{d}=-[\sigma (x)-1]\pi /2$). According to Eqs. (11) and (12), one has $\varphi _{1}=S\theta$ and $\varphi _{2}=2S\tilde {\theta }(x)$, hence Eq. (13) is derived via Eq. (14).

The stability of the stationary vortex droplets was tested by direct real-$z$ simulations (through the algorithm of split-step Fourier transformation) of the perturbed evolution in the framework of Eqs. (6) and (7) up to $z=10000$. If the intensity pattern of the droplet remains unchanged up to this length (which corresponds to several diffraction lengths), it is categorized as a stable solution. Unstable solutions readily exhibit splitting in the course of the simulations.

As the stripe modulation acts solely in the $x$-direction, the vortex solutions feature similar modulation in the same direction and can be characterized by the number of stripes, $\mathcal {N}$, in the localized solution. According to Eq. (3), the solutions are also categorized into the OC and IC types. For the OC- and IC-type solutions, with the pivot located at the center of a stripe or at the border between adjacent stripes, numbers $\mathcal {N}$ are, respectively, odd or even. Typical examples for the OC-type vortex solution with $\mathcal {N}=3,5,7,9$ and IC-type ones with $\mathcal {N}=2,4,6,8$, all carrying the winding number $S=1$ in their FF component, are shown, respectively, in Figs. 2 and 3. These states feature a flat-top density structure, due to the competition between the quadratic and cubic nonlinear terms, which is a characteristic feature of droplet patterns. All these states are stable, as confirmed by direct simulations up to $z=10000$. Because the modulation is applied only along the $x$-direction, the vortex solutions feature a typical semidiscrete configuration similar to that reported in previous works [1,34]. However, unlike those works, the stable vortex modes are elaborated here in the bulk crystals with the spatially modulated local $\chi ^{(2)}$ susceptibility. In Fig. 1 (a), black and gray stripes represent positive and negative polarizations, respectively. Consequently, the density distribution of optical droplets at these positions on the x-axis induces a phase difference of $\pi /2$ in the second harmonic, leading to a quasi-discrete structure of the field. The vortex phase patterns exhibited in Figs. 2 and 3 of the SH component are striped to obey the matching condition of Eq. (14), hence, the effective angular coordinate of this component is defined by Eq. (13), showing a complicated striped-mixed vorticity pattern of this component.

Fig. 2. Typical examples of the intensity distribution, $|\phi _{1}(x,y)|^{2}$ (the first row) and $|\phi _{2}(x,y)|^{2}$ (the third row), and phase patterns of $\phi _{1}(x,y)$ (the second row) and $\phi _{2}(x,y)$ (the fourth row) of a semidiscrete vortex optical droplet with $S=1$, which corresponds to point “A-D" in the stable area of the OC-type in Fig. 4(a). The parameters are $(P,\ell )=(80,18)$ in (a), $(80,11)$ in (b), $(80,8)$ in (c), and $(80,6.5)$ in (d). In this case, the effective detuning is fixed as $\Omega =0$.

Download Full Size | PDF

Fig. 3. A typical example of a stable IC-type optical droplet with $(P,\Omega )=(80,0)$ and $S=1$. The first and second rows display the intensity and phase distributions of the FF component, while the third and fourth rows exhibit the same for the SH component. The stripe’s widths in panels (a-d) are $\ell =22,13,9,$ and $7$, respectively, corresponding to points “A-D" in Fig. 4(b).

Download Full Size | PDF

Stability areas for the vortex solutions of the OC and IC types, with different values of $\mathcal {N}$ are shown, in the $(P,\ell )$ plane with $\Omega =0$ and $S=1$, in Fig. 4(a,b). These plots demonstrate that the vortex solutions with larger values of $P$ and smaller values of $\ell$ produce overlaps between different stability areas [see the cyan and yellow areas in Fig. 4(a), and the cyan area in Fig. 4(b)]. The overlaps indicate the presence of multistability in the system. To further illustrate this feature, we select $\ell =10$ and examine solutions with different values of $\mathcal {N}$, up to $P=200$. In this region, we find that four different values of $\mathcal {N}$ stably coexist at $P=200$ for both OC and IC types of solutions. Note also that the solutions with larger values of $\mathcal {N}$ require larger values of $P$ to support the stability. Multistable semidiscrete vortex solutions are characterized by values of $H$, $\beta$ and $r$, which are displayed, as functions of $P$, in Fig. 5, while keeping $\ell$ and $\Omega$ fixed. Notably, curves $H(P)$ and $\beta (P)$ for these solutions overlap almost completely, indicating degenerate solutions. The nearly flat dependences, with $d\beta /dP\approx 0$, indicate the existence of broad states, which may be considered as effectively liquid ones, cf. Reference [14]. A nearly constant value $r(P)\approx 2.7$ in Fig. 4(e,f) implies the domination of the FF component in the system. It is necessary to point out that the number of stripes, $\mathcal {N}$, is not a parameter affecting the droplet´s stabilization and multistability. Rather, this number distinguishes OC and IC (onsite- and intersite-centered) type solutions and to characterizes coexisting states. The key factors affecting the droplet’s stabilization and multistability are the total power ($P$) and period ($\ell$).

Fig. 4. Stability areas of semidiscrete vortex optical droplets of the OC (a) and IC (b) types, in the $(P,\ell )$ plane with $\Omega =0$ and $S=1$. Digits in colored stability areas indicate the number of stripes in the droplets, which are odd in (a) and even in (b). Note the presence of bistability in the cyan and yellow areas. Stability ranges of multistable states are shown in (c) and (d) for the droplets of the OC and IC types, respectively, with $(\ell,\Omega )=(10,0)$ and varying values of the total power, $P$.

Download Full Size | PDF

Fig. 5. (a,b) Hamiltonian $H$, defined as per Eq. (8), (c,d) the FF propagation constant $\beta$, and (e,f) ratio $r=P_{1}/P_{2}$ vs. total power $P$. The lines of blue stars, cyan circles, orange squares, and red balls in (a,c,e) indicate, respectively, semidiscrete vortex states with stripe numbers $\mathcal {N}=5,7,9$, and $11$, while in (b,d,f), they represent $\mathcal {N}=6,8,10$, and $12$. Other parameters are $\ell =10,\Omega =0$, and vorticity of fundamental frequency $S=1$.

Download Full Size | PDF

Finally, the phase-mismatch factor is determined by the modulation period $\Lambda$ (after the normalization, the phase mismatch exhibits by the effective detuning ´$\Omega$´). In Fig. 6 we present the range of the effective detuning for a fixed total power, $P=100$. The results indicate that stable semidiscrete vortex droplets exist for $\Omega \neq 0$.

Fig. 6. (a) and (b): The FF propagation constant $\beta$ for stable vortex optical droplets of the OC and IC types vs. effective detuning $\Omega$. Solid and dashed lines denote the number of stripes $\mathcal {N}=5$ and $\mathcal {N}=7$ for OC, or $\mathcal {N}=6$ and $\mathcal {N}=8$ for IC, respectively. Other parameters are $(P,\ell )=(100,10)$ and $S=1$.

Download Full Size | PDF

Exploring vortex states with $S>1$ is a challenging problem, as they may be stable only for sufficiently large values of $P$. Two typical examples of stable vortex solutions of the OC and IC types with $S=2$, and $\mathcal {N}=5$ and $6$, are shown in Fig. 7. The stable solution with $S=2$ exists at $P>280$. This threshold is much higher than its counterpart for $S=1$, in which case the stable solutions are found at $P>40$.

Fig. 7. (a1)-(a4) A typical example of a stable semidiscrete vortex droplet of the OC type, with vorticity $S=2$ of the FF component. Other parameters are $(P,\ell,\Omega )=(300,15,0)$. (b1)-(b4) An example of a stable semidiscrete vortex droplet of the IC type, with vorticity $S=2$ of the FF component, for the same parameters as in (a1)-(a4). The four columns of panels from left to right display the intensity and phase patterns of the FF and SH components. These droplets with $S=2$ remain stable for the propagation distance $z=10000$, up to which the simulations were running.

Download Full Size | PDF

3.2 Outline of the experimental setup

The fabrication of 3D photonic crystal by means of femtosecond laser pulses is a mature technology [71,72]. To estimate parameters of the setting under the consideration, we consider the nonlinear photonic crystal implemented in LiNbO$_{3}$, which has the second-order nonlinearity coefficient $d_{0}=d_{22}=2.1$ pm/V [73], and the third-order one $\chi ^{(3)}=6.6\times 10^{-22}$ m$^{2}$/V$^{2}$ [74]. The wavelengths of the FF and SH components are selected as $1064$ nm and $532$ nm, respectively. The relations between the scaled units, in which Eqs. (6) and (7) are written, and their physical counterparts can be established by means of Eq. (5), as summarizes in Table 1.

Table 1. Relations between scaled and physical units of the coordinates, stripe’s width, total power, and intensity

View Table

According to the results of the simulations (see Figs. 2 and 3), the peak droplet’s intensity in the FF and SH components are $I_{\mathrm {FF}}\approx 4.32$ GW/cm$^{2}$ and $I_{\mathrm {SH}}\approx 0.432$ GW/cm$^{2}$, respectively. If we select the pulse width as $200$ ps, the energy densities for the FF and SH components are $0.86$ J/cm$^{2}$ and $0.086$ J/cm$^{2}$, which are lower than the damage threshold ($\sim$1.5 J/cm $^{2}$ [75,76]) of the PPLN crystals. The characteristic propagation distance, which is $z=10000$, amounts to $2.8$ m, which is several times the underlying diffraction length. Therefore, the stability of the droplets, predicted by the simulations, is a reliable prediction.

The sketch for the experimental observation of these droplets is shown in Fig. 1(d): A high-power sub-nanosecond laser (e.g., one with the power exceeding $0.6$ mJ per pulse, $200$ ps pulse width, and $400$ Hz repetition rate) may be a suggested light source in the proposed experiment. The input patterns on the front surface of the PPLN can be generated by a spatial light modulator (SLM) and a Lens. Input FF and SF components can be selected with energies $0.37$ mJ and $0.16$ mJ per pulse, which are close to the FF/SF power ratio $r\approx 2.7$ for stationary semidiscrete vortex droplets in Figs. 6(e) and (f). If the input’s power ratio is essentially different from this value, it naturally gives rise to strong oscillations between the FF and SH components. The power and phase patterns of the FF component of the input, which are shown in the left inset of Fig. 1(d), can be produced by a properly designed SLM1 and L1, and the input SH patterns with $0.14$ mJ per pulse can be produced by SLM2 and L2, respectively. The FF and SH beams are coupled by the dichroic mirror, which sends them onto the PPLN coaxially. The necessary PPLN crystal may be $4$ cm long along the axial direction. Finally, the beams transmitted through the PPLN are imaged on to a camera by Lens L3. The simulated output patterns at the back surface of the crystals are shown in the right inset in Fig. 1(d).

4. Conclusion

We have proposed 3D photonic crystals with the striped structure and the combination of the $\chi ^{(2)}$ and defocusing $\chi ^{(3)}$ nonlinearities. The results of the analysis predict the creation of two types of stable semidiscrete vortex solutions, OC and IC (onsite- and intersite-centered ones), which exhibit an odd or even number of stripes in their structure, respectively. The smallest number of the stripes is $\mathcal {N}_{\text {OC}}=3$ or $\mathcal {N}_{\text {IC}}=2$. The setting admits with different numbers of stripes for the same parameters. Unlike the multistability in systems with the competing cubic-quintic nonlinear interactions, the Hamiltonian and propagation constant of these coexisting states are equal, thus featuring degeneracy. The range of the multistability has been found. The stable solutions for semidiscrete vortices exist within a certain range of positive or negative values of the phase mismatch of the $\chi ^{(2)}$ interaction in the medium.

The scheme proposed in this paper can be developed in other settings, such as photonic crystals with ring- or fan-shaped structures. Those settings can be used, in particular, for the implementation of various scenarios of beam shaping. It is worth noting that it is feasible to generalize the model introduced in this paper to include the sum-frequency generation and three-wave mixing by quadratic terms. In particular, it is relevant to analyze how soliton states change into droplets under the action of the three-wave mixing. Moreover, Furthermore, as detailed in Ref. [28], the ability to generate distinctive vortex self-localized states, such as ones with hidden vorticity, can be considered in the process of three-wave mixing. These findings are expected to extend into the investigation of vortex states in the realm of structured light fields [77,78].

Funding

National Natural Science Foundation of China (11874112, 12274077); Israel Science Foundation (1695/22); Ciuangdong-Hong Kong-Macao Joint l.aboratory for lntelligent Micro-Nano Optoelectroniclechnology (2020B1212030010).

Acknowledgment

Xiaoxi Xu and Feiyan Zhao, appreciate the Graduate Innovative Talents Training Program of Foshan University.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. X. Zhang, X. Xu, Y. Zheng, et al., “Semidiscrete quantum droplets and vortices,” Phys. Rev. Lett. 123(13), 133901 (2019). [CrossRef]

2. B. Malomed and P. Kevrekidis, “Discrete vortex solitons,” Phys. Rev. E 64(2), 026601 (2001). [CrossRef]

3. F. Lederer, G. I. Stegeman, D. N. Christodoulides, et al., “Discrete solitons in optics,” Phys. Rep. 463(1-3), 1–126 (2008). [CrossRef]

4. B. A. Malomed, “Vortex solitons: Old results and new perspectives,” Phys. D 399, 108–137 (2019). [CrossRef]

5. B. A. Malomed, “Nonlinearity and discreteness: Solitons in lattices,” Emerging Frontiers in Nonlinear Science 32, 81–110 (2020). [CrossRef]

6. Z. Zhao, G. Chen, B. Liu, et al., “Discrete vortex quantum droplets,” Chaos, Solitons Fractals 162, 112481 (2022). [CrossRef]

7. X. Xu, F. Zhao, Y. Zhou, et al., “Vortex gap solitons in spin–orbit-coupled bose–einstein condensates with competing nonlinearities,” Communications in Nonlinear Science and Numerical Simulation 117, 106930 (2023). [CrossRef]

8. L. Bergé, “Wave collapse in physics: principles and applications to light and plasma waves,” Phys. Rep. 303(5-6), 259–370 (1998). [CrossRef]

9. P. Di Trapani, W. Chinaglia, S. Minardi, et al., “Observation of quadratic optical vortex solitons,” Phys. Rev. Lett. 84(17), 3843–3846 (2000). [CrossRef]

10. D. Mihalache, D. Mazilu, B. Malomed, et al., “Stable vortex solitons in a vectorial cubic–quintic model,” J. Opt. B: Quantum Semiclass. Opt. 6(5), S341–S350 (2004). [CrossRef]

11. C. Sulem and P.-L. Sulem, “The nonlinear Schrödinger equation: self-focusing and wave collapse139, (Springer Science & Business Media, 2007).

12. G. Fibich, “The nonlinear Schrödinger equation192, (Springer, 2015).

13. B. A. Malomed, “Multidimensional solitons, (AIP Publishing New York, 2022).

14. D. Petrov, “Quantum mechanical stabilization of a collapsing bose-bose mixture,” Phys. Rev. Lett. 115(15), 155302 (2015). [CrossRef]

15. D. Petrov and G. Astrakharchik, “Ultradilute low-dimensional liquids,” Phys. Rev. Lett. 117(10), 100401 (2016). [CrossRef]

16. M. Schmitt, M. Wenzel, F. Böttcher, et al., “Self-bound droplets of a dilute magnetic quantum liquid,” Nature 539(7628), 259–262 (2016). [CrossRef]

17. C. Cabrera, L. Tanzi, J. Sanz, et al., “Quantum liquid droplets in a mixture of bose-einstein condensates,” Science 359(6373), 301–304 (2018). [CrossRef]

18. Y. V. Kartashov, B. A. Malomed, L. Tarruell, et al., “Three-dimensional droplets of swirling superfluids,” Phys. Rev. A 98(1), 013612 (2018). [CrossRef]

19. Y. Li, Z. Chen, Z. Luo, et al., “Two-dimensional vortex quantum droplets,” Phys. Rev. A 98(6), 063602 (2018). [CrossRef]

20. Z.-H. Luo, W. Pang, B. Liu, et al., “A new form of liquid matter: Quantum droplets,” Front. Phys. 16(3), 32201 (2021). [CrossRef]

21. Z. Lin, X. Xu, Z. Chen, et al., “Two-dimensional vortex quantum droplets get thick,” Communications in Nonlinear Science and Numerical Simulation 93, 105536 (2021). [CrossRef]

22. Y. Zheng, S. Chen, Z. Huang, et al., “Quantum droplets in two-dimensional optical lattices,” Front. Phys. 16(2), 22501 (2021). [CrossRef]

23. G. Li, X. Jiang, B. Liu, et al., “Two-dimensional anisotropic vortex quantum droplets in dipolar Bose-Einstein condensates,” Front. Phys. 19, 22202 (2024).

24. K. E. Wilson, N. Westerberg, M. Valiente, et al., “Observation of photon droplets and their dynamics,” Phys. Rev. Lett. 121(13), 133903 (2018). [CrossRef]

25. N. Westerberg, K. E. Wilson, C. W. Duncan, et al., “Self-bound droplets of light with orbital angular momentum,” Phys. Rev. A 98(5), 053835 (2018). [CrossRef]

26. A. V. Buryak, Y. S. Kivshar, S. Trillo, et al., “Optical solitons supported by competing nonlinearities,” Opt. Lett. 20(19), 1961–1963 (1995). [CrossRef]

27. D. Mihalache, D. Mazilu, B. A. Malomed, et al., “Stable vortex solitons supported by competing quadratic and cubic nonlinearities,” Phys. Rev. E 69(6), 066614 (2004). [CrossRef]

28. A. S. Desyatnikov, D. Mihalache, D. Mazilu, et al., “Two-dimensional solitons with hidden and explicit vorticity in bimodal cubic-quintic media,” Phys. Rev. E 71(2), 026615 (2005). [CrossRef]

29. Z. Wu, Y. Zhang, C. Yuan, et al., “Cubic-quintic condensate solitons in four-wave mixing,” Phys. Rev. A 88(6), 063828 (2013). [CrossRef]

30. R. DeSalvo, D. J. Hagan, M. Sheik-Bahae, et al., “Self-focusing and self-defocusing by cascaded second-order effects in ktp,” Opt. Lett. 17(1), 28–30 (1992). [CrossRef]

31. M. Chen, D. Kaup, B. A. Malomed, et al., “Three-wave solitons and continuous waves in media with competing quadratic and cubic nonlinearities,” Phys. Rev. E 69(5), 056605 (2004). [CrossRef]

32. C. Bosshard, R. Spreiter, M. Zgonik, et al., “Kerr nonlinearity via cascaded optical rectification and the linear electro-optic effect,” Phys. Rev. Lett. 74(14), 2816–2819 (1995). [CrossRef]

33. A. S. Reyna and C. B. de Araújo, “High-order optical nonlinearities in plasmonic nanocomposites-a review,” Adv. Opt. Photon. 9(4), 720–774 (2017). [CrossRef]

34. X. Xu, G. Ou, Z. Chen, et al., “Semidiscrete vortex solitons,” Adv. Photonics Res. 2(3), 2000082 (2021). [CrossRef]

35. H. Leblond, B. A. Malomed, D. Mihalache, et al., “Spatiotemporal vortices in optical fiber bundles,” Phys. Rev. A 77(6), 063804 (2008). [CrossRef]

36. T. Xu, K. Switkowski, X. Chen, et al., “Three-dimensional nonlinear photonic crystal in ferroelectric barium calcium titanate,” Nat. Photonics 12(10), 591–595 (2018). [CrossRef]

37. D. Wei, C. Wang, H. Wang, et al., “Experimental demonstration of a three-dimensional lithium niobate nonlinear photonic crystal,” Nat. Photonics 12(10), 596–600 (2018). [CrossRef]

38. S. Liu, L. Wang, L. Mazur, et al., “Highly efficient 3d nonlinear photonic crystals in ferroelectrics,” Adv. Opt. Mater. 11(14), 2300021 (2023). [CrossRef]

39. S. Zhu, Y.-y. Zhu, N. Ming, et al., “Quasi-phase-matched third-harmonic generation in a quasi-periodic optical superlattice,” Science 278(5339), 843–846 (1997). [CrossRef]

40. C. B. Clausen, O. Bang, Y. S. Kivshar, et al., “Spatial solitons and induced kerr effects in quasi-phase-matched quadratic media,” Phys. Rev. Lett. 78(25), 4749–4752 (1997). [CrossRef]

41. O. Bang, C. B. Clausen, P. L. Christiansen, et al., “Engineering competing nonlinearities,” Opt. Lett. 24(20), 1413–1415 (1999). [CrossRef]

42. J. F. Corney and O. Bang, “Solitons in quadratic nonlinear photonic crystals,” Phys. Rev. E 64(4), 047601 (2001). [CrossRef]

43. A. Arie and N. Voloch, “Periodic, quasi-periodic, and random quadratic nonlinear photonic crystals,” Laser & Photon. Rev. 4(3), 355–373 (2010). [CrossRef]

44. J. Thomas, V. Hilbert, R. Geiss, et al., “Quasi phase matching in femtosecond pulse volume structured x-cut lithium niobate,” Laser & Photonics Reviews 7(3), L17–L20 (2013). [CrossRef]

45. X. Hu, P. Xu, S. Zhu, et al., “Engineered quasi-phase-matching for laser techniques,” Photonics Res. 1(4), 171–185 (2013). [CrossRef]

46. A. Karnieli and A. Arie, “Frequency domain stern–gerlach effect for photonic qubits and qutrits,” Optica 5(10), 1297–1303 (2018). [CrossRef]

47. A. Karnieli, S. Trajtenberg-Mills, G. Di Domenico, et al., “Experimental observation of the geometric phase in nonlinear frequency conversion,” Optica 6(11), 1401–1405 (2019). [CrossRef]

48. J. Lin, N. Yao, Z. Hao, et al., “Broadband quasi-phase-matched harmonic generation in an on-chip monocrystalline lithium niobate microdisk resonator,” Phys. Rev. Lett. 122(17), 173903 (2019). [CrossRef]

49. H. Li and B. Ma, “Research development on fabrication and optical properties of nonlinear photonic crystals,” Front. Optoelectron. 13(1), 35–49 (2020). [CrossRef]

50. N. V. Bloch, K. Shemer, A. Shapira, et al., “Twisting light by nonlinear photonic crystals,” Phys. Rev. Lett. 108(23), 233902 (2012). [CrossRef]

51. A. Shapira, R. Shiloh, I. Juwiler, et al., “Two-dimensional nonlinear beam shaping,” Opt. Lett. 37(11), 2136–2138 (2012). [CrossRef]

52. A. Shapira, I. Juwiler, A. Arie, et al., “Tunable nonlinear beam shaping by non-collinear interactions,” Laser & Photonics Reviews 7(4), L25–L29 (2013). [CrossRef]

53. S. Trajtenberg-Mills, I. Juwiler, A. Arie, et al., “On-axis shaping of second-harmonic beams,” Laser & Photonics Reviews 9(6), L40–L44 (2015). [CrossRef]

54. Z. Xu, Z. Lin, Z. Ye, et al., “Control the orbital angular momentum in third-harmonic generation using quasi-phase-matching,” Opt. Express 26(13), 17563–17570 (2018). [CrossRef]

55. D. Wei, C. Wang, X. Xu, et al., “Efficient nonlinear beam shaping in three-dimensional lithium niobate nonlinear photonic crystals,” Nat. Commun. 10(1), 4193 (2019). [CrossRef]

56. J. Imbrock, L. Wesemann, S. Kroesen, et al., “Waveguide-integrated three-dimensional quasi-phase-matching structures,” Optica 7(1), 28–34 (2020). [CrossRef]

57. P. Chen, C. Wang, D. Wei, et al., “Quasi-phase-matching-division multiplexing holography in a three-dimensional nonlinear photonic crystal,” Light: Sci. Appl. 10(1), 146 (2021). [CrossRef]

58. F. Zhao, X. Xu, H. He, et al., “Vortex solitons in quasi-phase-matched photonic crystals,” Phys. Rev. Lett. 130(15), 157203 (2023). [CrossRef]

59. N. K. Efremidis, S. Sears, D. N. Christodoulides, et al., “Discrete solitons in photorefractive optically induced photonic lattices,” Phys. Rev. E 66(4), 046602 (2002). [CrossRef]

60. J. W. Fleischer, G. Bartal, O. Cohen, et al., “Observation of vortex-ring discrete solitons in 2d photonic lattices,” Phys. Rev. Lett. 92(12), 123904 (2004). [CrossRef]

61. D. N. Neshev, T. J. Alexander, E. A. Ostrovskaya, et al., “Observation of discrete vortex solitons in optically induced photonic lattices,” Phys. Rev. Lett. 92(12), 123903 (2004). [CrossRef]

62. B. Terhalle, T. Richter, A. S. Desyatnikov, et al., “Observation of multivortex solitons in photonic lattices,” Phys. Rev. Lett. 101(1), 013903 (2008). [CrossRef]

63. H. Suchowski, G. Porat, A. Arie, et al., “Adiabatic processes in frequency conversion,” Laser & Photonics Reviews 8(3), 333–367 (2014). [CrossRef]

64. A. Karnieli and A. Arie, “Fully controllable adiabatic geometric phase in nonlinear optics,” Opt. Express 26(4), 4920–4932 (2018). [CrossRef]

65. A. Karnieli, Y. Li, A. Arie, et al., “The geometric phase in nonlinear frequency conversion,” Front. Phys. 17(1), 12301 (2022). [CrossRef]

66. Y. Li, O. Yesharim, I. Hurvitz, et al., “Adiabatic geometric phase in fully nonlinear three-wave mixing,” Phys. Rev. A 101(3), 033807 (2020). [CrossRef]

67. F. Zhao, J. Lü, H. He, et al., “Geometric phase with full-wedge and half-wedge rotation in nonlinear frequency conversion,” Opt. Express 29(14), 21820–21832 (2021). [CrossRef]

68. G. Luther, M. Alber, J. Marsden, et al., “Geometric analysis of optical frequency conversion and its control in quadratic nonlinear media,” J. Opt. Soc. Am. B 17(6), 932–941 (2000). [CrossRef]

69. G. Porat and A. Arie, “Efficient, broadband, and robust frequency conversion by fully nonlinear adiabatic three-wave mixing,” J. Opt. Soc. Am. B 30(5), 1342–1351 (2013). [CrossRef]

70. C. Phillips, C. Langrock, D. Chang, et al., “Apodization of chirped quasi-phasematching devices,” J. Opt. Soc. Am. B 30(6), 1551–1568 (2013). [CrossRef]

71. S. Keren-Zur and T. Ellenbogen, “A new dimension for nonlinear photonic crystals,” Nat. Photonics 12(10), 575–577 (2018). [CrossRef]

72. A. Arie, “Storing and retrieving multiple images in 3d nonlinear photonic crystals,” Light: Sci. Appl. 10(1), 202 (2021). [CrossRef]

73. V. G. Dmitriev, G. G. Gurzadyan, D. N. Nikogosyan, et al., “Nonlinear optical properties of crystals,” Handbook of Nonlinear Optical Crystals 64, 53–127 (1991). [CrossRef]

74. R. Ganeev, I. Kulagin, A. Ryasnyanskii, et al., “The nonlinear refractive indices and nonlinear third-order susceptibilities of quadratic crystals,” Opt. Spectrosc. 94(4), 561–568 (2003). [CrossRef]

75. D. H. Titterton and J. A. Terry, “Study of laser-induced damage in ppln,” in Laser-Induced Damage in Optical Materials: 2001, Vol. 4679 (SPIE, 2002), pp. 312–320.

76. D. H. Titterton, J. A. Terry, D. H. Thorne, et al., “Observation of damage in ppln,” in Growth, Fabrication, Devices, and Applications of Laser and Nonlinear Materials, Vol. 4268 (SPIE, 2001), pp. 5–13.

77. S. Fu, C. Guo, G. Liu, et al., “Spin-orbit optical hall effect,” Phys. Rev. Lett. 123(24), 243904 (2019). [CrossRef]

78. W. Zhen, X. Wang, J. Ding, et al., “Controlling the symmetry of the photonic spin hall effect by an optical vortex pair,” Phys. Rev. A 108(2), 023514 (2023). [CrossRef]

$x = 1 & y = 1$	4.64 $μ$ m
$ℓ = 10 \sim 15$	$46.4 \sim 69.6$ $μ$ m
$z = 1$	280 $μ$ m
$P = 1$	31 kW
$\| u_{1} \|^{2} = 0.01 & \| u_{2} \|^{2} = 0.01$	1.44 GW/cm $^{2}$ & 0.72 GW/cm $^{2}$

Semidiscrete optical vortex droplets in quasi-phase-matched photonic crystals

Abstract

1. Introduction

2. Model

3. Results

3.1 Numerical results

3.2 Outline of the experimental setup

4. Conclusion

Funding

Acknowledgment

Disclosures

Data availability

References

Data availability

Cited By

Figures (7)

Tables (1)

Equations (14)

Optics Express