## Abstract

We demonstrate two alternative techniques for controlling and stabilizing domain walls (DW) in phase-sensitive, nonlinear optical resonators. The first of them uses input pumps with spatially modulated phase and can be applied also to dark-ring cavity solitons. An optical memory based on the latter is demonstrated. Here the physical mechanism of control relies on the advection caused to any feature by the phase gradients. The second technique uses a plane wave input pump with holes of null intensity across its transverse plane, which are able to capture DWs. Here the physical mechanism of control is of topological nature. When distributed as a regular array, these holes delimit spatial optical bits which constitute an optical memory. These techniques are illustrated in a degenerate optical parametric oscillator model, but can be applied to any phase-sensitive nonlinear optical cavity.

©2004 Optical Society of America

## 1. Introduction

Cavity solitons (CS) are robust, self-sustained localized structures that can be excited across the transverse section of the light beam produced by a nonlinear optical resonator with high Fresnel number [1–4]. CSs are interesting objects from the applicative viewpoint as they can serve as bits of an optical memory, while their manipulation and control can allow the parallel processing of information [5]. Especially versatile systems for the generation of CSs are those in which the nature of the nonlinear interaction is phase-sensitive and the generated field can take two possible phase values, differing by π. Thus, such systems display phase bistability. A paradigm of this behavior is represented by the degenerate optical parametric oscillator (DOPO) in which the generated, subharmonic field is phase-locked to the pump field, modulo π. In these phase bistable systems the light field phase can take different values in different regions of the transverse plane, which are separated among them by dark lines. These lines are known as domain walls (DW) or dark cavity-solitons because of their resemblance with the dark solitons of the nonlinear Schrödinger equation. These DWs are essentially 1D structures; related 2D objects are optical droplets (a circular front connecting the two equivalent, oppositely phased homogeneous solutions) [6] and dark-ring cavity solitons (DRCS), where a circular front separates an inner, high intensity peak, from an outer, homogeneous light distribution, of opposite phase [4].

In nonlinear optical cavities DWs were first described in DOPOs [7, 8] and have been experimentally demonstrated in a photorefractive oscillator [9]. In [7] DWs were shown to be stable structures for both one and two transverse dimensions. While for one transverse dimension this seems to be always true irrespectively of the size of the excited region, this is not the case for two transverse dimensions: It has been shown [10] that DWs are only transient solutions of the system that dissapear through the borders of the illuminated region. This would make these structures of reduced practical interest in real systems unless something is done for stabilizing and controlling them. In this sense, the stabilization of DWs in DOPOs pumped by helical waves has been theoretically shown [10].

In this paper two simple methods for controlling and stabilizing DWs are proposed. The methods are illustrated in the case of DOPO but in principle they can be applied to any phase-sensitive nonlinear optical cavity displaying phase bistability. The first method uses an injected field with spatially modulated phase profile [11]. Here the physical mechanism of control relies in the advection exerted by phase gradients on DWs, which tend to be trapped along regions in which the injected field phase is maximum and escape from regions of phase minima, as we show next. The second method uses an injected field with spatially modulated intensity, in the form of a “hole burned” plane wave input. Here the holes of null intensity are able to capture the ends of DWs, and the physical mechanism of control is of topological nature as will be discussed below. When distributed as a regular array these holes delimit spatial optical bits which constitute an optical memory.

It is worth mentioning that although DWs have not yet been observed in DOPOs, their observation in photorefractive oscillators [9] makes possible the implementation of the stabilizing and controlling techniques we present here.

## 2. Model

In the usual paraxial and mean-field approximations, the dynamics of a doubly resonant, phase-matched, type-I DOPO with planar mirrors is governed by the following dimensionless equations for the pump *A*
_{0}(*x*,*y*,*t*) and signal (subharmonic) *A*
_{1}(*x*,*y*,*t*) field envelopes [12]:

Denoting by *γ*
_{1} (*γ*
_{0}) the cavity decay rate for the signal (pump) field, and by *k*
_{1} the signal field wavenumber within the crystal, *t* is time normalized to ${\gamma}_{1}^{-1}$, *γ* = *γ*
_{0}/*γ*
_{1}, (*x*,*y*) are transverse coordinates normalized to $\sqrt{\frac{c}{{\gamma}_{1}}{k}_{1}}$ (*c* is the speed of light in vacuum), ∇^{2} = *∂*
^{2}/*∂X*
^{2}+*∂*
^{2}/*∂y*
^{2}, Δ_{0} and Δ_{1} are the cavity detunings of pump and signal fields normalized to *γ*
_{0} and *γ*
_{1} respectively, and *E*(*x*,*y*) is the complex envelope of the injected pump field.

For a spatially homogeneous input field *E* (in which case *E* can be taken as a positive real without loss of generality), the DOPO oscillation threshold occurs at *E*
^{2} = 1 + ${\mathrm{\Delta}}_{0}^{2}$ for Δ_{1} ≥ 0, giving rise to a spatially uniform signal field, or at *E* = 1 for Δ_{1}< 0, giving rise to a spatially modulated signal field in the form of a roll pattern [12, 13]. Consider the particularly simple, perfectly resonant case Δ_{0} =Δ_{1} = 0. Above threshold the DOPO admits two equivalent nontrivial solutions $\left\{{A}_{0}=1,{A}_{1}=\pm \sqrt{E-1}\right\}$ in which the signal field can take two values (±) differing by π. In this case DWs exist [7,8] and are stable in two transverse dimensions as far as the injected field is transversely unbounded. But a real optical system has always a finite transverse extension which, in most cases, is usually determined by the transverse size of the input beam. Oppo et al. [10] have shown that for a pump with circular cross section DWs disappear after a transient time and a single phase ultimately covers the whole transverse plane. We show next how DWs can be stabilized and controlled.

## 3. Two control methods

#### 3.1. Method based on phase modulated pump

This method was originally introduced by Firth and Scroggie [11] in the context of the “optical bullet holes” (cavity solitons in which a bright light spot is surrounded by a homogeneous background) formed in an optical resonator filled by a saturable absorber. These solitons are of a nature different from that of DWs, as the former do not depend on any phase bistability: They rely on the bistability between a homogeneous solution and an hexagonal pattern. Firth and Scroggie [11] showed that phase gradients are effective in advecting bright structures and used this fact to pin them at desired locations. Oppo et al. [10] also showed that DWs are advected by phase gradients as they found that DWs created in DOPOs pumped by helical waves (with a linearly varying phase in the azymuthal direction) were rotating permanently. Here we demonstrate that the phase gradient method is able to control, select, and pin DWs at desired locations.

In order to get some insight into the effect of a spatially modulated input phase we write *E* = ∣*E*(*x*,*y*)∣exp[*iψ*(*x*,*y*)] and rewrite Eqs. (1) and (2) for the variables ${F}_{0}={A}_{0}\mathrm{exp}[-\frac{1}{2}\mathit{i\psi}\left(x,y\right)]$ [11],

where

is a convective derivative, *κ*
_{0}=1 + ∇^{2}
*ψ*/2*γ*, *κ*
_{1} = 1 + ∇^{2}
*ψ*/2, *δ*
_{0} = Δ_{0} +∣∇ψ∣^{2}/2*γ* and *δ*
_{1} = Δ_{1} +∣∇*ψ*∣^{2}/4. Eqs. (3,4) are formally equivalent to the original DOPO Eqs. (1,2), but now the parameters are, in general, space dependent. On the one hand linear losses (*κ*
_{0,1}) can increase (reduce) for ∇^{2}
*ψ*< 0 (∇^{2}< 0). On the other hand the cavity mistunings are always increased as Δ_{0,1} ≥ Δ_{0,1}. But the most important modification is the appearance of the advection term (**∇**
*ψ*)∙ **∇** accompanying the time derivative in Eq. (5), which indicates that any feature in the transverse plane will move towards phase maxima with velocity **v** =**∇**
*ψ*. This suggests that the use of a convergent (divergent) input field could stop (increase) the natural tendency of DWs to escape the illuminated region.

In order to test these ideas the DOPO Eqs. (1) and (2) were integrated for a top hat shape input pump: Its intensity ∣*E*∣^{2} = *P*
^{2} was uniform inside a circle of radius *r*
_{0}, falling off very quickly to zero outside it. As a consequence the intracavity pump and signal fields were off outside the circle. Its phase *ψ*(*x*,*y*) was allowed to be spatially varying. A split-step algorithm with periodic boundary conditions was used on a square grid of side *L*(>2 *r*
_{0}) with 256×256 points. Since both pump and signal fields were excited only within a finite region of the integration window, the chosen boundary conditions had no effect on the DOPO dynamics. We show next results obtained for *γ*= 1, Δ_{0} =Δ_{1} = 0, *P* = 3, *L* = 100, and *r*
_{0} = 45 (other parameters, whenever Δ_{1} ≥ 0, yield similar results), and for a temporal step δ*t*= 0.01 (smaller δ*t* ’s were also used for asessing the convergence of the results). Also other spatially confined pumps with different (e.g., Gaussian) intensity profiles were used, yielding similar results.

If the pump phase is flat (*ψ*= constant) we observe that a straight DW excited in the initial condition leaves the illuminated region after a transient, in agreement with [10]. This behavior is due to local curvature driven dynamics, which develops as the initially straight DW tends to allign itself perpendicularly to the circular boundary. This dynamics is faster for DWs initially excited far from the pump beam center but also occurs for straight DWs initially crossing that center, whenever noise is considered. After our previous discussion one can envisage a way to stabilize this DW: If the pump wavefront is given a cylindrical shape so that it has a maximum along a diameter of the circular input, the DW should move towards that diameter and then remain there forever. Fig. 1 shows this pinning effect for a cylindrical wavefront, convergent along the *x* direction. The initial condition was a straight DW located at *x* = 0.56*r*
_{0}. The DW moves to the left until it finally settles at the vertical line *x* = 0, where *ψ* is maximum.

One can wonder whether this control mechanism will be effective in the presence of noise. We have integrated the DOPO equations with additive noise terms (δ-correlated Gaussian white noise of zero mean) for different values of the noise variance. The results show that the mechanism is robust and works perfectly well for any noise value that allows the formation of DWs. As an example, in a numerical experiment as that of Fig. 1 but with added noise, we arrive at the same final state (left panel in Fig. 2). The intensity trace along a horizontal cut shown in the right panel of Fig. 2 shows that, even for very large noise levels, the DW is effectively stabilized and kept attached to the prescribed position. Let us remark that noise does not introduce any noticeable delay in reaching the final state.

We have also investigated the formation and evolution of DWs starting from a random seed. In this case, if the pump phase is flat, we observe the formation of several phase domains of irregular shapes which in the course of time disappear in a way similar to Fig. 9 in [10]. One could ask whether in this so disordered situation a single DW can be picked and stabilized. The positive answer is given in Fig. 3, where the same pump phase profile as in Fig. 1 was used, but now the simulation started from a random seed. After a short transient a labyrinthic pattern is formed. As time goes on all DWs tend to dissapear except one that tends to align itself with the pump phase maximum. We observe that the dynamics can be very slow especially if DRCSs bond to the DW as it occurs in Fig. 3. But one can remove these defects as well by giving an additional curvature to the pump wavefront in the yet neutral direction *y*, as evidenced in the bottom row of Fig. 3, which forces the defects to exit the illuminated domain, and the final state is a single vertical DW [14]. The previous situation can be generalized to several phase maxima as shown in Fig. 4.

*We have checked also the effect of noise when the initial seed is random. As in the previous case, we have found that the control mechanism is effective in selecting and pinning a particular DW even in the presence of large noise.*

*This method can be used as well for stabilizing DRCSs in much the same way as it was used in [11] for stabilizing “optical bullet holes”. We illustrate this in Fig. 5. The left panel shows the initial condition, where nine DRCSs were written in a somewhat disordered manner. Then an egg-box like phase profile with 3×3 maxima was given to the injected field. The result after a transient is shown in the middle panel where the nine DRCSs are finally fixed at the positions of the input field phase maxima. The fact that this constitutes an optical memory is demonstrated in the right panel where the central “bit” has been removed by injecting a narrow pulse of appropriate phase [15], without disturbing the remaining bits.*

*The calculation was repeated in the presence of additive noise (same variance as in Fig. 2) obtaining the same results. As an example we just show in Fig. 6 the state analogous to the
middle panel in Fig. 5. The right panel in Fig. 6 evidences the high level of noise present in the simulation.*

*3.2. Method based on intensity modulated pump*

*As DWs are topological solitons [7] a different control method can be envisaged that makes use of input pumps whose transverse profile contains zeroes (or regions below DOPO oscillation). In such case there are regions in the transverse plane where the signal field is off, to which a DW can be attached. The point is that a DW cannot be detached from a “hole” as, in that case, a DW with a free end would exist. But this cannot occur in phase-locked systems as a DW separates regions with opposite phases and the only possibility for that would be that a continuous 2π variation of the signal phase would develop around that free end, similarly to a vortex. An example of DW stabilization using this method is given in Fig. 7. In the left panel the initial condition is shown: A pair of DWs (note that an even number of DWs ending in the hole must be excited as the signal phase must be continuous after a circulation around the hole) were written in such a way that they were not perpendicular to the pumped region boundaries. Curvature forces the two DW to move (central panel) until they become perpendicular to the boundaries (right panel) where the DW remain at rest. If the initial condition is such that the DW are perpendicular to the boundaries, they do not move. This method can be combined with the phase gradient method in order to move the DWs at will.*

*The method can be exploited further to constuct an optical memory as shown in Fig. 8. In this case the input pump contains an array of 4×4 holes. In the left panel the signal field phase is spatially uniform, and bits of information can be addressed by applying gaussian pulses whose phases are opposite (with a large tollerance) to the existing one, middle panel. In this way square patches (phase domains) of a given phase can be written. This can be done independently in each of the square regions delimited by holes, as well as use of oppositely phase pulses allows to erase existing domains, right panel.*

*4. Conclusions*

*We have demonstrated two robust methods for stabilizing and controlling DWs. The first one is based on the advection introduced by pump phase gradients, so that DWs tend to be trapped along the lines of pump phase local maxima. The use of divergent or tilted wavefronts in order to get rid of defects has been demonstrated as well. The method is also useful for pinning DRCSs at desired locations. The robustness of the method has been demonstrated in the presence of (high levels of) noise. The second method is based on the topological nature of DWs and uses input pumps with zeroes in its transverse spatial distribution which trap DWs. Both methods have been shown to permit the development of optical memories. The techniques have been applied to DOPOs, but they should allow the stabilization and control of DWs and DRCSs in any nonlinear optical resonator exhibiting phase bistability like, e.g., degenerate four-wave mixing oscillators [9] and vectorial Kerr cavities [16,17].*

*Acknowledgments*

*Financial support from the Spanish Ministerio de Ciencia y Tecnología and the European Union FEDER (Project No. BFM2002-04369-C04-01), and from the ESF network PHASE is acknowledged.*

*References and links*

**1. **P. Mandel, *Theoretical Problems in Cavity Nonlinear Optics* (Cambridge University Press, Cambridge, 1997). [CrossRef]

**2. **L. A. Lugiato, M. Brambilla, and A. Gatti, “Optical Pattern Formation,” Adv. At. Mol. Opt. Phys. **40**, 229 (1999). [CrossRef]

**3. **S. Trillo and W. Torruellas, eds., *Spatial Solitons* (Springer-Verlag, Berlin, 2001).

**4. **K. Staliunas and V. J. Sánchez-Morcillo, *Transverse Patterns in Nonlinear Optical Resonators* (Springer-Verlag, Berlin, 2003).

**5. **W. J. Firth and C. O. Weiss, “Cavity and Feedback Solitons,” Opt. Photon. News **13**, 55 (2002). [CrossRef]

**6. **D. Gomila, P. Colet, G. L. Oppo, and M. San Miguel, “Stable Droplets and Growth Laws Close to the Modulational Instability of a Domain Wall,” Phys. Rev. Lett. **87**, 194101 (2001). [CrossRef] [PubMed]

**7. **S. Trillo, M. Haelterman, and A. Sheppard, “Stable Topological Spatial Solitons in Optical Parametric Oscillators,” Opt. Lett. **22**, 970 (1997). [CrossRef] [PubMed]

**8. **K. Staliunas and V. J. Sánchez-Morcillo, “Spatial-Localized Structures in Degenerate Optical Parametric Oscillators,” Phys. Rev. A **57**, 1454 (1998). [CrossRef]

**9. **V. B. Taranenko, K. Staliunas, and C. O. Weiss, “Pattern Formation and Localized Structures in Degenerate Optical Parametric Mixing,” Phys. Rev. Lett. **81**, 2236 (1998). [CrossRef]

**10. **G. L. Oppo, A. J. Scroggie, and W. J. Firth, “Characterization, Dynamics and Stabilization of Diffractive Domain Walls and Dark Ring Cavity Solitons in Parametric Oscillators,” Phys. Rev. E **63**, 066209 (2001). [CrossRef]

**11. **W. J. Firth and A. J. Scroggie, “Optical Bullet Holes: Robust Controlable Localized States of a Nonlinear Cavity,” Phys. Rev. Lett. **76**, 1623 (1996). [CrossRef] [PubMed]

**12. **G. L. Oppo, M. Brambilla, and L. A. Lugiato, “Formation and Evolution of Roll Patterns in Optical Parametric Oscillators,” Phys. Rev. A **49**, 2028 (1994). [CrossRef] [PubMed]

**13. **G. J. de Valcárcel, K. Staliunas, E. Roldán, and V. J. Sánchez-Morcillo, “Transverse Patterns in Degenerate Optical Parametric Oscillation and Degenerate Four-Wave Mixing,” Phys. Rev. A **54**, 1609 (1996). [CrossRef] [PubMed]

**14. **We checked that a similar effect is achieved by introducing a linearly varying phase in the *y* direction that physically corresponds to tilting vertically the input field.

**15. **M. Brambilla, L. A. Lugiato, and M. Stefani, “Interaction and Control of Optical Localized Structures,” Europhys. Lett. **34**, 109 (1996). [CrossRef]

**16. **M. Hoyuelos, P. Colet, M. San Miguel, and D. Walgraef, “Self-Similar Domain Growth, Localized Structures, and Labyrinthine Patterns in Vectorial Kerr Resonators,” Phys. Rev. E **61**, 2241 (2000) [CrossRef]

**17. **V. J. Sánchez-Morcillo, I. Pérez-Arjona, F. Silva, G. J. de Valcárcel, and E. Roldán, “Vectorial Kerr Cavity Solitons,” Opt. Lett. **25**, 957 (2000). [CrossRef]