Extended and localized Hopf-Turing mixed-mode in non-instantaneous Kerr cavities

We investigate the spatio-temporal dynamics of a ring cavity filled with a non-instantaneous Kerr medium and driven by a coherent injected beam. We show the existence of a stable mixed-mode solution that can be either extended or localized in space. The mixed-mode solutions are obtained in a regime where Turing instability (often called modulational instability) interacts with self-pulsing phenomenon (Andronov-Hopf bifurcation). We numerically describe the transition from stationary inhomogeneous solutions to a branch of mixed-mode solutions. We characterize this transition by constructing the bifurcation diagram associated with these solutions. Finally, we show stable localized mixed-mode solutions, which consist of time-periodic oscillations that are localized in space.


I. INTRODUCTION
Transition from spatially uniform state to a self-organized or ordered structures is a universal feature in far from equilibrium systems and has been observed in many natural systems (see overviews on this issue [1][2][3][4][5]). This transition is triggered by a few number of modes often called Turing modes that leads to the formation of spatially stationary structures that can be either periodic or localized in space. On the order hand, a transition to a time oscillation or a self-pulsing state through the Andronov-Hopf (termed Hopf in the following) bifurcation has also been observed. These transitions are responsible for the symmetry breaking in space (Turing) and in time (Hopf). The interaction between Turing and Hopf bifurcations concerns almost all fields of nonlinear science such as biology, chemistry, physics, fluid mechanics, and optics. When both instabilities are close one to another, the space-time dynamics of various spatially extended systems may be significantly impacted. In particular, it has been shown that in this regime either a bistable behavior between Turing patterns and self-pulsing states or front waves between Hopf and Turing-type domains may occur [6][7][8][9][10][11]. Moreover, mixed-mode solutions resulting from the interplay between Andronov-Hopf and Turing modes may dominate the spatio-temporal dynamics in optical frequency conversion systems [12,13]. This interaction can also lead to a spatio-temporal chaos [14].
In this paper, we consider a ring cavity filled with a non instantaneous and a nonlocal Kerr medium, and driven by a coherent injected beam. We show that the threshold associated with Turing and Hopf bifurcations can occur arbitrarily close one to another leading to codimension-two point where both bifurcations coincide. In the monostable regime, we show that the Turing pattern becomes unstable and leads to the formation of a Turing-Hopf mixedmode. More importantly, we show that localized mixed-mode solutions can be stabilized when the system exhibits a coexistence between an homogeneous steady state and an extended mixed-mode solution.

II. MODEL EQUATIONS AND A LINEAR STABILITY ANALYSIS
The slowly varying envelope of the electric field E = E(x, t) circulating inside the cavity is described by ∂E/∂z = inE + i∂ 2 E/∂x 2 coupled with an additional equation for the material photoexcitation n = n(x, t): γ −1 ∂n/∂z = |E| 2 + d∂ 2 n/∂x 2 with z is the longitudinal coordinate, and the parameter d is proportional to the ratio between the diffusion length and diffraction coefficients. This parameter describes the degree of nonlocality. Indeed, if d = 0, the above propagation equations have a local nonlinear response. The strong nonlocal response is characterized by a large d. The parameter γ is proportional to the ratio between the characteristics decay times associated with the electric field and the refractive index n, respectively. At each round trip, the light inside the Kerr media is coherently superimposed with the input beam. This can be described by the cavity boundary conditions E(z = 0, x) = rS + τ E(L, x) exp (iφ e ) and E(z = 0, x) = rn(L, x) exp (iφ n ), with φ e,n are linear phase shifts, S is the ampliude of the injected field, and L denotes the cavity length. The letters "r" and "τ " indicate the transmission and the reflexion coefficients. We adopt the well know mean field approach proposed by Lugiato and Lefever [16]. A detailled calculations of the derivation of the mean field model can be found in the appendix of the paper [17]. This approach is valid under the following approximations: the cavity possesses a high Fresnel number and we assume that the cavity length is much shorter than the diffraction, diffusion and the nonlinearity spatial scales. A single longitudinal mode operation is also assumed. The mean field approach applied to our system leads to the adimensional coupled partial differential equations where δ is the frequency detuning between the injected light and the cavity resonance. S is the injected field amplitude considered as being real without loss of generality. Diffraction of the electric field in the transverse x coordinate is modeled by i∂ 2 E/∂x 2 . The density of material photoexcitation is denoted by n that diffuses according to ∂ 2 n/∂x 2 .
In the limit of δ = 0, the model has been derived for a Kerr medium but in a single feedback mirror configuration [15]. For a ring cavities, and in the particular limit where d = 0 and γ = 0, one obtains from Eq. (2) n = |E| 2 . By replacing n by |E| 2 in Eqs.
(1), we recover the well known Lugiato-Lefever (LL) equation [16]. The model Eqs. (1) and 2) may be viewed as an extension of LL equation to include the non instantaneous response of the medium and the nonlocal effects. Homogeneous steady state (HSS) solutions of Eqs. (1) and (2) are found by setting the time and space derivatives equal to zero: n s = I s and The characteristic I s as a function of the input intensity S 2 is monstable when δ < √ 3. The HSS's undergo a bistable behavior when δ > √ 3. The linear stability analysis of the HSS with respect to a finite wavelength perturbation of the form exp (λt + ikx) leads to the thirdorder polynomial characteristic equation: The Turing type of bifurcation occurs when λ = 0 and ∂λ/∂k 2 = 0. The threshold associated with this instability and the critical wavelength at the Turing bifurcation are obtained from the expression a 0 = 0. This instability is characterized by an intrinsic wavelength which is determined by the dynamical parameters and not by the external geometrical constraints such as boundary conditions.
A Hopf bifurcation occurs if a pair of complex-conjugate roots has a vanishing real part and a nonzero imaginary part. In the limit γ = 1 and δ < 2, only one Hopf bifurcation point is possible. The associated critical intensity is explicitly given by I H = 2(2−δd)(d+ √ d 2 + 1). At this bifurcation point the real part vanishes at nonzero wavenumber Note that the HSS's may undergo a homogeneous Hopf bifurcation where the real part of the eigenvalues vanishes for k = 0.
We are interested in the regime where Turing and Hopf instabilities are close one to another. For this purpose, we fix γ = 1 and δ = 1, and we vary the coefficient d and the input field intensity. An important and interesting feature is that the parameter d controls the relative position of the thresholds associated with both Turing and Hopf instabilities. Typical marginal stability curves are shown in Fig. 1. By increasing the value of the parameter d, the first bifurcation is of Turing type, and the Hopf bifurcation appears as secondary instability as shown in Fig. 1(a). There exists a critical value of the parameter d for which both bifurcations coincide, i.e., I H = I T = I c as shown in Fig. 1(b). At this co-dimensional two point, one of the real roots of the characteristic equation vanishes λ = 0 with a non zero finite intrinsic wave length 2π/k c and two other roots are purely imaginary λ = ±iω. Note that the critical mode k c is alsod-dependent. When further increasing d, the first instability is the Hopf bifurcation followed by the Turing instability as shown in Fig. 1(c).

III. EXTENDED AND LOCALIZED MIXED-MODE
In what follows we focus on numerical investigations of the model Eqs. (1) and (2) by using an adaptive step-size Bulirsch-Stoer method [18]. We choose parameters values such that the homogeneous steady state is destabilized first by Turing instability. From this bifurcation point, a branch of stationary spatially periodic solutions appears supercitically with a well defined wavelength. However, when increasing further the injected field intensity a stable mixed-mode solution is spontaneously generated in the system. These solutions correspond to oscillations both in time and in space. A typical example of such a behavior is plotted in Fig. 2. To characterize the transition from stationary periodic patterns to a branch of mixed-mode Turing-Hopf structures, we plot in Fig. 3 the maxima of both Turing structures and mixed-mode solutions together with the homogeneous steady state. As we increase the amplitude of the injected field, the HSS becomes unstable, and a spatially periodic structure emerges from the Turing bifurcation point S = S T . These structures are stable in the parameter range S T < S < S MM . For S > S MM , the Turing structure becomes unstable and bifurcates to a stable mixed-mode solution. The extended solutions either Turing structures or mixed-mode solutions are obtained in the supercritical regime. Besides these extended structures there exists another type of solutions which are aperiodic but localized in space (LS's). The latter are found in the subcritical regime associated with the Turing instability. It is well known that the emergence of these solutions does not necessarily requires a bistable regime [19]. The prerequisite condition for the generation of LS's is the coexistence between a single HSS and the spatially periodic solutions. This solutions have been predicted theoretically [20] and experimentally evidenced in an instantaneous and local optical Kerr medium [21]. Notice that temporal localized solutions known also as temporal solitons have also been experimentally observed in a dispersive Kerr cavity in [22]. Here we show that in the case of a non-instantaneous and a nonlocal Kerr medium, modeled by the Eqs. (1) and (2), support localized structures. An example of such a behavior is shown in Fig. 4. We show in this figure only localized structures with one, two and three peaks. The number of localized peaks and their spatial distributions are determined by the initial conditions [20]. There exists an infinite number of spatially localized solutions if the size of the system is infinite. When increasing the injected field intensity, localized structures exhibit a pulsing phenomenon, i.e., time oscillations in a wide range of parameters leading to the formation of Hopf-Turing mixed-mode solutions. A typical example of such a behavior is shown in Fig. 5(a). The maxima of intensity associated with localized structures and localized mixed-mode solutions are plotted together with the homogeneous steady states in Fig. 5(b). When increasing the input intensity, a single peak stationary localized structure is formed in the range S SL < S < S MM . For S > S MM , stable localized mixed-mode solutions are generated. Their maximum intensity are shown in Fig. 5(b). The width of these time dependent structures varies in the course of time.

IV. CONCLUSIONS
We have shown that extended mixed-mode solutions can be generated in the output of a driven ring cavity filled with a non-instantaneous and a nonlocal Kerr medium. Our investigations are focused on the regime where Turing and Hopf instabilities interact strongly. In the monostable regime, the mixed-mode solutions appear as a transition from spatially periodic structures. We have drawn the bifurcation diagram showing their stability domain. In the bistable regime, we have shown occurrence of stationary localized structures that exhibit multistability in a finite range of parameters. Finally, we have established evidence of localized mixed-mode solutions that emerge from the branch of a single peak localized structures. We have constructed a bifurcation diagram associated with these localized structures.