Generation and evolution of quadratic dark spatial solitons using the wavefront modulation method

We study the evolution of dark spatial solitons by χ:χ cascaded nonlinearity. The generation and evolution of solitons pairs by the optical branching effect is investigated. It is shown that the number of solitons pairs and the transverse velocity of a given pair depend on the parameters of the input beam. We also investigate the hollows generated in a beam with a uniform background when the two sides of the input beam have various phase difference. It is concluded that dark spatial solitons by χ:χ cascaded nonlinearity have analog performance compared with that in χ defocusing Kerr media. ©2005 Optical Society of America OCIS codes: (190.5530) Pulse propagation and solitons; (190.2620) Frequency conversion. References and Links 1. N. N. Akhmediev, “Spatial solitons in Kerr and Kerr-like Media,” Opt. Quantum. Electron. 30, 535-569 (1998). 2. M. Segev, “Optical spatial solitons,” Opt. Quantum. Electron. 30, 503-533 (1998). 3. J. S.Aitchison, “Observation of spatial optical solitons in a nonlinear glass waveguide,” Opt. Lett. 15, 471-473 (1990). 4. P. -Andre Belanger and Pierre Mathieu, “Dark soliton in a Kerr defocusing medium,” Appl. Opt. 26, 111 -113 (1987) 5. G. Bartal, O. Manela, O. Cohen, J. W. Fleischer, and M. Segev, “Observation of 2nd-band vortex solitons in 2D photonic lattices,” Phys. Rev. Lett. 95, 053904 (2005). 6. M. Segev and G. Stegeman, “Self-trapping of optical beams: spatial solitons,” Phys. Today 51, 42-48 (1998). 7. G. I. Stegeman, D. N. Christodoulides, and M. Segev, “Optical spatial solitons: historical perspectives,” Invited Paper, Special Millennium Issue of the IEEE J. Sel Top. Quantum. Electron. 6, 1419-1427 (2000). 8. G. P. Agrawal, Nonlinear Fiber Optics (Academic 1995). 9. A. D. Boardman, L. Pavlov, and S. Tanev, Advanced Photonics with Second-order Optical Nonlinear Processes (Kluwer, 1998). 10. A. D. Boardman, K. Marinov, D. I. Pushkarov, and A Shivarova, “Wave-beam coupling in quadratic nonlinear optical waveguides: effects of nonlinearly induced diffraction,” Phys. Rev. E. 62, 2871-2877 (2000). 11. A. D. Boardman and K. Xie, “Spatial bright-dark soliton steering through waveguide coupling,” Opt. Quantum. Electron. 30, 783-794 (1998). 12. Y. S. Kivshar and B. Luther-Davies, “Dark optical solitons: physics and applications,” Phys. Rep. 298, 81-197 (1998). 13. Y. S. Kivshar and G. P. Agrawal, Optical Solitons (Academic, 2003). 14. V. E. Zakharov and A. B. Shabat, “Exact Theory of two-dimensional self-focusing and one dimensional selfmodulation of waves in nonlinear media,” Sov. Phys. JETP 34, 62-69 (1972). 15. G. I. Stegeman and M. Segev, “Optical spatial solitons and their interactions: universality and diversity,” Science 286, 1518-1523(1999). 16. D. S. Jiang, S. G. Ouyang and W. L. She, “Interaction of photovotaic spatial solitons between dark and dark as well as dark and bright,” Acta. Phys. Sin. 53, 3777-3785 (2004). 17. S. G. Ouyang, X. S. Wang, and W. L. She, “Interaction between photovotaic spatial solitons with difference #8512 $15.00 USD Received 19 August 2005; revised 9 October 2005; accepted 9 October 2005 (C) 2005 OSA 31 October 2005 / Vol. 13, No. 22 / OPTICS EXPRESS 8699 wavelengths,” Acta. Phys. Sin. 53, 761-772(2004). 18. Y. Q. Xie and Q. Guo, “Interaction between optical spatial solitons in nonlocal Kerr media,” Acta. Phys. Sin. 53, 3020-3024 (2004) 19. Y. V. Kartashov, L. Torner, and V. A. Vysloukh, “Parametric amplification of soliton steering in optical lattices,” Opt. Lett. 29, 10-13 (2004). 20. L. A.Ostrovskii, “Self action of light in crystals,” JETP Lett. 5, 272-275 (1967). 21. V. E. Zakharov and S. V. Manakov “The theory of resonance interaction of wave packets in nonlinear media,” Sov. Phys. JETP 42, 842-850 (1976). 22. C. R. Menyuk, R. Schiek, and L. Torner, “Solitary waves due to (2) (2) : χ χ cascading,” Opt. Soc. Am. B 11, 2434-2443 (1994). 23. C. B. Clausen, O. Bang, and Y. S. Kivshar, “Spatial solitons and Induced Kerr effects in quasi-phase-matched quadratic media,” Phys. Rev. Lett. 78, 4749-4752(1997). 24. L. Torner and C. R. Menyuk, W. E. Torruellas, and G. I. Stegeman, “Two-dimensional solitons with second-order nonlinearities,” Opt. Lett. 20, 13-15 (1995). 25. W. E. Torruellas, Z. Wang, D. J. Hagan, E. W. VanStryland and G. I. Stegeman, “Observation of two-dimensional spatial solitary waves in quadratic midium,” Phys. Rev. Lett. 74, 5036-5039 (1995). 26. N. Nikolov, D. Neshev and O. Bang, “Quadratic solitons as nonlocal solitons,” Phys. Rev. E 68, 036614 (2003). 27. G. Assanto and G. I. Stegeman, “Simple physics of quadratic spatial solitons,” Opt. Express 10, 388-396 (2002), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-10-9-388. 28. P. Vidakovic, D. J. Lovering, and J. A. Levenson, “Large nonlinear phase shift owing to cascaded (2) χ in quasi-phase-matched bulk LiNbO3,” Opt. Lett. 22, 277-229 (1997). 29. B. Bourliaguet, V. Couderc, and A. Barthelemy, “Observation of quadratic spatial solitons in periodically poled lithium niobate,” Opt. Lett. 24, 1410-1412 (1999). 30. H. Kim, L. Jankovic, and G. I. Stegeman, “Quadratic spatial solitons in periodically poled KTiOPO4,” Opt. Lett. 28, 640-642 (2003). 31. D. R. Andersen, D. E. Hooton, and G. A. Swartzlander, Jr., “Direct measurement of the transverse velocity of dark spatial solitons,” Opt. Lett. 15, 783-785 (1990). 32. S. A. Gredeskul and Yu. S. Kivshar, “Dark-pulse solitons in nonlinear-optical fibers,” Phys. Rev. A 41, 3994-4008 (1990). 33. B. Luther-Davies and X. Yang, “Steerable optical waveguides formed in self-defocusing media by using dark spatial solitons,” Opt. Lett. 17, 1755-1757 (1992). 34. C. B. Clausen and L. Torner, “Spatial switching of quadratic solitons in engineered quasi-phase-matched structures,” Opt. Lett. 24, 7-10 (1999).

As we know, quadratic nonlinearity can effectively mimic cubic properties.The existence of nonlinear phase shift in fundamental waves during SHG was first discussed by Ostrovskii in 1967 [20], and the existence of solitons was predicted in 1974 by V. E. Zakharov and S. V. Manakov [21].In the past 12 years, it was already proved experimentally that bright solitary waves can exist with cascaded second-order nonlinearities ( : χ χ ) in quadratic media [22][23][24][25].An advanced study on the cascading effect and the induced Kerr effect was conducted by Bang.Theoretically, Bang and his colleagues showed that quadratic solitons are equivalent to solitons of a nonlocal Kerr medium [26].Assanto and Stegeman used the concepts of cascading phase shift to interpret the formation and properties of quadratic solitons [27].In quasi-phase-matched (QPM) cascaded SHG processes, both positive and negative values of 2 cascad n are observed and measured [28].Bright spatial solitons have been observed in periodically poled lithium niobate and KTP crystals during second-harmonic generation [29,30].Dark spatial solitons by (2) ( : χ χ nonlinearity in QPM configuration can also be expected.
Since a cascading second-order nonlinearity is much larger than a third-order nonlinearity, lower power is needed to achieve similar effects in cascaded second-order quadratic media.It has been reported that the optical branching effect was experimentally observed in Kerr media [4].The beam propagating in the media splits into several beams.These beams propagate at different angles on each side of the beam.Theoretical calculation and experimental measurement of the transverse velocity of dark solitons in Kerr media are reported by Anderson et al. [31].By the wavefront modulation method, dark spatial solitons and soliton pairs can be easily generated, and they can be steerable by controlling the blocking width or the phase difference at the wavefront of the beam.The purpose of this paper is to investigate the evolution characters of dark solitons in cascaded media and to compare them with those of solitons in Kerr media.We numerically simulate this effect in quadratic media to study its concrete properties and the potential influences on it.In addition, we investigate the generation of a hollow in a beam with a uniform background and the method to control the location of the generated hollow.

Theory
Under the condition of slowly varying envelope approximation, the wave coupling equations in a QPM quadratic nonlinearity media in a reduced normalized form is where 1 2 , a a are the normalized amplitudes of the fundamental and harmonic waves, respectively, and , k k are the wavenumbers at the two frequencies.The Λ is the one-order QMP grating period), and η is the characteristic beam transverse width.ξ is the propagation distance in the unit of Based on the above wave coupling equation, the generation and evolution of dark solitons under various conditions can be numerically simulated.In our calculation, we integrate the wave coupling equation numerically with a split-step approach.The linear part ( 2 ⊥ ∇ ) is integrated in the Fourier space, and the nonlinear part is integrated by a fourth Runge-Kutta algorithm.We divide the propagation process into many steps.In every step, the diffraction effect is first considered exclusively, and then only the nonlinear process is calculated.

Optical branching effect in cascaded quadratic media
In Kerr media, it has been observed experimentally that beam splitting occurs after a certain value of the input beam intensity has been reached.The dark solitons are evolved and generated in pairs [4].The nth pair from the center of the beam is called the nth-order pair.These phenomena result from the inherent repulsive character of their interaction: it is energetically favorable to break a soliton into two pieces traveling away from each other.Naturally, a problem is raised: does the optical branching effect occur in the cascaded quadratic media?
A box-like dark spatial pulse as an initial input condition is suitable for analyzing the branching effect.An ideal box-like pulse is given by 0 1 If above ideal box-like pulse is employed as the initial condition, the simulation results are singular due to abrupt discontinuity at the points x a = ± .To avoid this issue, an input beam described by the function below is introduced to simulate the box-like dark solitons: In the calculation, we found that the soliton pairs do not appear simultaneously when the beam enters the media.Along with the evolution of dark pulses in the media, soliton pairs appear one after another symmetrically beside the center of the hollow and evolve to the edges.For fixed intensity and as time increases, more branches appear until a steady state is reached where the number of branches is fixed.The above behaviors are the same as those in Kerr media [4].To observe the steady state, we studied the evolution of box-like dark solitons after propagating various diffraction lengths, and we find that after a certain number of diffraction lengths, the soliton pair number does not increase.
Figure 1 shows the generation and evolution process of dark soliton pairs at different diffraction lengths.At 15 diffraction lengths [Fig.1(b)], only the first-order pair is completely produced; at 20 diffraction lengths [Fig.1(c)], the second-order pair is found; at 35 diffraction lengths [Fig.1(d)], the third-order pair is produced.After 35 diffraction lengths, no higher-order pair appears and the interval between the existing three pairs become large.In addition, it is shown that the low-order pairs are darker than the high-order ones: the intensity of bottoms of high-order pairs is no longer zero, and it increases with the order of the pairs.The steady soliton pair number is governed by the parameters of the input beam and the cascading effect.As the input intensity increases, the generation of multiple solitons is favored.We investigated the steady branch numbers for input beams with different intensities.In this process, since the required diffraction length for stable branch number varies for different intensity, and since much longer evolution length will arouse reflections from the beam edges, repeated calculations are conducted to ensure that the branch number for each input beam is what it can finally engender.The results are shown in Fig. 2. When the input amplitude is 2, the steady state is reached at about 100 diffraction lengths and only the first-order pair is produced.When the amplitude is 4, 80 diffraction lengths are needed and three pairs are produced.While the amplitude is 7, only 60 diffraction lengths are needed and four pairs are produced.These results confirm the prediction that the higher the intensity of the input beam, the more soliton pairs are generated.The number of branches also differs when the initial width of the input hollow is changed.
With the increase of initial width, the branch number increases notably.We changed ω [the width of the input hollow described by function (2)] to ½ω , 2ω , 4ω and simulated the evolution process.The results are shown in Fig. 3.It can be found that, when the initial width is large, more pairs are produced and they are closer to each other than the pairs produced with a narrow input hollow.In addition, due to the conservation of energy, when a number of dark soliton pairs are produced, the intensity of the background increases obviously: in Fig. 3(d), the normalized intensity of the background reaches 20, while its initial value is 16.From the inverse scattering theory, in Kerr media the eigenvalue of a box-like dark soliton pair is determined by , where v is the transverse velocity measured experimentally and E 0 represents the amplitude of the input beam.This shows that for a given pair, the transverse velocity is directly proportional to the amplitude of the input beam.As in cascaded media, it is relatively difficult to obtain the analytic solution from the wave coupling equation.Through our numerical calculation, we proved that the soliton pairs in cascaded media have similar performance to those in Kerr media.With different fundamental wave amplitudes, we recorded the transverse displacements of the first-order pair at same diffraction length to analyze the relation between the amplitude and its transverse velocity.The results shown in Fig. 4 demonstrate exactly the linear relation between the displacement (indicating the transverse velocity) and the input beam amplitude, which is just as predicted in Kerr media.From the above discussions we can conclude that the optical branching effect in cascaded quadratic media proves to have analog performance compared with that in Kerr media: the soliton pairs are evolved and generated gradually through the propagation in the media; there is a steady state where the number of soliton pairs does not increase; the steady number of soliton pairs increases with the intensity of the input beam and the initial width of the hollow; the transverse velocity of a given pair is directly proportional to the input beam amplitude.

Generation of dark solitons in uniform background
In nonlinear-optical fibers, the generation of dark solitons in uniform background with unequal phase at t = ±∞ has been theoretically studied [32].Dark solitons can be generated by an input pulse with a phase difference at its edges that equals π.We demonstrate that, in cascaded quadratic media, spatial dark solitons can also be generated in a uniform background.In our calculation, the uniform background is given by where A represents the amplitude of input beam.A factor of ( ) i e δ is introduced in the calculation to introduce different phases.As we know, spatial dark solitons can be used to form steerable optical waveguides to be applied in all-optical apparatuses [33][34].In cascaded media, the steerablity of dark solitons by changing the phase difference at the center of the uniform background is investigated.The simulation results are shown in Fig. 5; dark solitons in a uniform background are generated with different transverse velocity, depending on the phase difference δ.In Fig. 5(b), when the two sides of the uniform background have a phase difference ofπ , a dark hollow in the center of the uniform background is generated.When the left side of the hollow has a phase of +π/2 ahead of the right side, the soliton hollow is generated at the left side of the center as shown in Fig. 5(c).In this case, the soliton hollow sinks less than the hollow generated in the center of background, which is called the gray soliton.In Fig. 5(d), the left side of the hollow has a phase of +3π/2 ahead of the right side and the hollow is generated right to the center.

Conclusion
In summary, the results of the numerical simulation show that optical branching phenomena also occur in cascaded media and the characteristics of the effect are similar to those in Kerr media.In addition, we demonstrate that hollows in a uniform background can be generated in cascaded media and the position of the hollow are steerable by means of changing the phase difference of two sides of the input beam.The dark solitons caused by the cascading effect may have potential application in all-optical apparatuses because lower power is required.

2 1 k
η .δ accounts for the Poynting vector walk-off when propagation is not along the crystal optical axes.We can set 0 δ = because Poynting vector walk-off is absent in typical QPM geometries.The function ( ) d ξ stands for the effective nonlinear coefficients involved in QPM. of one dimension; here x is the normalized transverse coordinate in units of η .

Fig. 4 .
Fig. 4. Relation between the y and the displacements of the first-order pairs at 100 diffraction lengths.

Fig. 5 .
Fig. 5. Generation of dark solitons in a uniform background.(a) Input uniform background, (b) hollow generated in the center of background with a phase difference of π at two edges of the uniform background, (c) hollow generated at left side of the background when the left side of the hollow has a phase of +π/2 ahead of the right side, (d) hollow generated at the right side of the background when the left side of the hollow has a phase of +3π/2 ahead of the right side.