Generation of broadband two-dimensional multicolored arrays in a sapphire plate

Broadband 2-D multicolored arrays with more than ten periodic columns and more than ten rows were generated in a sapphire plate using two crossed femtosecond laser beams overlapping in time and space. These multicolored 2-D arrays were sensitive to rotation of the sapphire plate in the plane normal to the incident beams. The broadband spectrum, pulse duration less than 40fs, and less than 1%RMS power stability of the array signals make them well suited for various applications, for example multicolor pump-probe experiments. 2009 Optical Society of America OCIS codes: (190.4380) Nonlinear optics, four-wave-mixing; (320, 2250) Femtosecond phenomena; (190, 7110) Ultrafast nonlinear optics; (190.4420) Nonlinear optics, transverse effects in. References and links 1. A. Saliminia, S. Chin, and R. Vallée, "Ultra-broad and coherent white light generation in silica glass by focused femtosecond pulses at 1.5 um," Opt. Express 13, 5731-5738 (2005). 2. M. Segev and G. Stegeman, “Self-Trapping of Optical Beams: Spatial Solitons,” Phys. Today 51, 43-48 (1998). 3. C. G. Durfee III, S. Backus, M. M. Murnane, and H. C. Kapteyn, "Ultrabroadband phase-matched optical parametric generation in the ultraviolet by use of guided waves," Opt. Lett. 22, 1565-1567 (1997). 4. C. Conti, S. Trillo, P. Di Trapani, G. Valiulis, A. Piskarskas, O. Jedrkiewicz, and J. Trull, “Nolinear electromagnetic X waves,” Phys. Rev. Lett. 90, 170406 (2003). 5. T. Kobayashi, A. Shirakawa, and T. Fuji, “Sub-5-fs transform-limited visible pulse source and its application to real-time spectroscopy,” IEEE J. Sel. Top. Quantum Electron. 7, 525-538 (2001). 6. H. Crespo, J. T. Mendonça, and A. Dos Santos, "Cascaded highly nondegenerate four-wave-mixing phenomenon in transparent isotropic condensed media," Opt. Lett. 25, 829-831 (2000). 7. J. Liu and T. Kobayashi, “Wavelength-tunable multicolored femtosecond laser pulse generation in fused silica glass,” Opt. Lett. 34, 1066-1068 (2009). 8. J. Liu and T. Kobayashi, “Generation of uJ-level multicolored femtosecond laser pulses using cascaded four-wave mixing,” Opt. Express 17, 4984-4990 (2009). 9. J. Liu and T. Kobayashi, “Cascaded four-wave mixing and multicolored arrays generation in a sapphire plate by using two crossing beams of femtosecond laser,” Opt. Express 16, 22119-22125 (2008). 10. M. Zhi and A. V. Sokolov, "Broadband coherent light generation in a Raman-active crystal driven by twocolor femtosecond laser pulses," Opt. Lett. 32, 2251-2253 (2007). 11. E. Matsubara, T. Sekikawa, and M. Yamashita, "Generation of ultrashort optical pulses using multiple coherent anti-Stokes Raman scattering in a crystal at room temperature," Appl. Phys. Lett. 92, 071104 (2008). 12. H. Matsuki, K. Inoue, and E. Hanamura, "Multiple coherent anti-Stokes Raman scattering due to phonon grating in KNbO3 induced by crossed beams of two-color femtosecond pulses," Phys. Rev. B 75, 024102 (2007) 13. K. Inoue, J. Kato, E. Hanamura, H. Matsuki, and E. Matsubara, “Broadband coherent radiation based on peculiar multiple Raman scattering by laser-induced phonon grating in TiO2,” Phys. Rev. B 76, 041101(R) (2007). 14. E. Matsubara, K. Inoue, and E. Hanamura, “Violation of Raman selection rules induced by two femtosecond laser pulses in KTaO3,” Phys. Rev. B 72, 134101 (2005). (C) 2009 OSA 25 May 2009 / Vol. 17, No. 11 / OPTICS EXPRESS 9226 #108884 $15.00 USD Received 17 Mar 2009; revised 21 Apr 2009; accepted 22 Apr 2009; published 18 May 2009 15. J. Liu, J. Zhang, and T. Kobayashi, "Broadband coherent anti-Stokes Raman scattering light generation in BBO crystal by using two crossing femtosecond laser pulses," Opt. Lett. 33, 1494-1496 (2008). 16. D. Faccio, A. Dubieties, G. Tamosauskas, P. Polesana, G. Valiulis, A. Piskarskas, A. Lotti, A. Couairon, and P. Di Trapani, “Phaseand group-matched nonlinear interactions mediated by multiple filamentation in Kerr media,” Phys. Rev. A 76, 055802 (2007). 17. D. Kip, M. Soljacic, M. Segev, E. Eugenieva, and D. N. Christodoulides, “Modulation instability and pattern formation in spatially incoherent light beams,” Science 290, 495-498 (2000). 18. J.W. Fleischer, M. Segev, N. K. Efremidis, and D. N. Christodoulides, “Observation of two-dimensional discrete solutions in optically induced nonlinear photonic lattices,” Nature 422, 147-150 (2003). 19. X. Wang, A. Bezryadina, Z. Chen,K. G. Makris, D. N. Christodoulides, and G. I. Stegeman, “Observation of Two-Dimensional Surface Solitons,” Phys. Rev. Lett. 98, 123903 (2007). 20. H. Zeng, J. Wu, H. Xu, and K. Wu, “Generation and weak beam control of two-dimensional multicolored arrays in a quadratic nonlinear medium,” Phys. Rev. Lett. 96, 083902 (2006). 21. V. Kartazaev, and R. R. Alfano, “ Polarization properties of SC generated in CaF2,” Opt. Commun. 281, 463-468 (2008). 22. P. B. Lundquist, D. R. Andersen, and Y. S. Kivshar, “Multicolor solitons due to four-wave mixing,” Phys. Rev. E 57, 3551–3555 (1998). 23. G. Fanjoux, J. Michaud, M. Delqu ́e, H. Mailotte, and T. Sylvestre, “Generation of multicolor vector Kerr solitons by cross-phase modulation, four-wave mixing, and stimulated Raman scattering,” Opt. Lett. 31, 3480–3482 (2006).


Introduction
When an ultrashort pulse propagates in a nonlinear bulk medium, various nonlinear optical phenomena such as supercontinuum light generation [1], optical solitons [2], optical parametric generation [3], and X-wave generation [4] have been observed.When two ultrashort pulses that overlap in space and time propagate in a nonlinear bulk medium, much more interesting nonlinear optical phenomena have appeared.Optical parametric amplification based on a nonlinear crystal is well established, and has been used for many scientific experiments [5].A 1-D broadband multicolored sideband array was observed when two crossed femtosecond beams were synchronized in BK7 glass [6], fused silica glass [7][8], sapphire [9] and certain crystals, for example PbWO 4 [10], LiNbO 3 [11], KNbO 3 [12], TiO 2 [13], KTaO 3 [14], and BBO [15].However, all these phenomena occurred in one dimension only.There were only few experiments based on 2-D nonlinear optical phenomena.2-D arrays of transverse patterns were experimentally observed in the multifilament interactions [16,17], and 2-D solitons were the subject of considerable research over the last decade.Twodimensional discrete solitons [18] and surface solitons [19] were observed in nonlinear photonic lattices.Regular 2-D multicolored transverse arrays were also observed when two crossed femtosecond laser beams were synchronized in a quadratic nonlinear crystal [20].We recently reported the observation of 2-D multicolored arrays in a sapphire plate [9], while only 5 rows of arrays were generated and only the spectrum property is studied in that paper.The data are not enough to understand this new phenomenon and to make use of the 2-D multicolored sidebands.
In this letter, we reported the generation of stable broadband 2-D multicolored arrays with more than ten periodic columns and more than ten rows in a sapphire plate using two crossed femtosecond laser beams overlapping in time and space.The properties of the spectrum, spatial mode, power stability, polarization, pulse duration and crossing angle between two neighbor sidebands of the 2-D muliticolored arrays were studied in detail.These researches will help people to understand the process and explain the phenomenon in the future.These data also clarify the possibility of the many kinds of applications of the generated sidebands.Moreover, we found that these stable 2-D multicolored arrays were sensitive to the orientation of the optical axis of the sapphire plate with respect to the plane of polarization of the incident beams.This is an interesting phenomenon that makes it possible to control the 2-D multicolor arrays using polarizer optics or rotating the sapphire plate.

Experiment setup
A 1kHz Ti:sapphire regenerative amplifier fs laser system (Micra+Legend-USP, Coherent) with 40 fs pulse duration and 2.5 W average output power was used as a pump source.The laser pulse after the regenerative amplifier system was split into three beams using beam splitters.One of the beams (beam_1) was spectrally broadened in a hollow fiber with a 250 µm inner-diameter and 60 cm in length that was filled with krypton gas.The broadband spectrum after the hollow fiber was dispersion compensated with a pair of chirped mirrors and a pair of glass wedges.The pulse duration after the hollow fiber compressor was about 10 fs.After passing through a bandpass filter at 720 nm, beam_1 was focused into a 2-mm-thick sapphire plate by a concave mirror with 60cm focal length.Another beam (beam_2) passed through a delay stage with less than 3 fs resolution.Beam_2 was first attenuated by a variable neutral density (VND) filter and was then focused into the sapphire plate by a lens with 1m focal length.The third beam (beam_3) was used to generate correlation signals with the input beams and the generated sidebands in a 10-µm-thick BBO crystal, which was used to measure pulse durations.The spectra of array signals on the center row B m, 0 were measured using a multichannel spectrometer (USB4000, Ocean Optics), as shown in Fig. 2(a).For clarity, not all sideband spectra are shown.A broadband spectrum from 400 nm to 1.2 µm with more than 1.5 octaves was generated [9].These generated sidebands were explained to be the result of a cascaded FWM process, which was recently reported in detail [6][7][8][9].The spectra were tunable by changing the crossing angle α between the two input beams, and also by changing the center wavelength of the bandpass filter.The maximum difference between peak wavelengths of the side spots and the center spot on the same column was about 20 nm, as shown in Fig. 2 (b).The spatial profiles of different signals in the arrays were measured using a CCD camera.Figure 3(a) shows spatial profiles of B 0, 0 , B 1, 0 , and B 4, 1 in two dimensions and one dimension.Figure 3(b) shows one-dimensional spatial profiles of B 0, 0 , B 1, 0 , and B 4, 1 with logarithmic scale in the intensity.The spatial profile changed from Gaussian (B 0, 0 ) to a Lorentzian profiles (B 4, 1 ).The cyan dashed line in Fig. 3(b) is the Gaussian fit of B 0, 0 .The magenta dash-dotted line in Fig. 3(b) is the Lorentzian fit of B 4, 1.The pulse duration of the sidebands were measured by cross-correlation with beam_3.Figure 3(c) shows the retrieved XFROG pulse trace and the phase of B 1, 0 with a retrieved error of 0.01022 from the commercial software from Femtosoft Technologies.The retrieved pulse duration was 35 ± 3 fs.The retrieved phase shows that there was some chirp in the pulse due to the positive chirped input pulses and the dispersion of the glass.The pulse duration of B 1, 0 was even shorter than the two incident pulses, the cross-correlation width of which with beam_3 were  We measured the powers of some array signals when the power of the two incident beams, beam_1 and beam_2, were 0.1 mW and 25 mW, respectively, as shown in Fig. 4(a).The signal power decreased rapidly with increasing column order for signals on the center row B m, 0 , as shown by the star symbols in Fig. 4(a).The difference in output power between the side spots and the center spot on the same column decreased continuously as the row order increased.For signals on column 5 (B 5, n ), signal power on the center row B 5, 0 was even smaller than that on rows ±1 and ±2.The dependence of the different sideband output power on the input power of beam_2 is shown in the inset of Fig. 4(a).The increase in signal power with higher column or row order was delayed, but more rapid than that of lower order columns or rows.In this case, the power of beam_1 was very low, and saturation did not take place.In fact, the output power of beam_1 was amplified from 0.1 mW to 0.17 mW in the experiment.The power stabilities of different array signals were monitored by a Si power sensor, as shown in Fig. 4(b).The stability of B 1, 0 was about 1.84% RMS over 200 seconds.Interestingly, the stabilities of the high order sidebands were much better than that of the firstorder, especially for the sidebands beside the beam on the center row.The stabilities of B 4, 0 , B 3, 1 , and B 4, 1 were 1.25%RMS, 0.63% RMS, and 0.97% RMS over 200 seconds, respectively.We observed the multicolored arrays to be sensitive to the rotation of the sapphire plate in the plane normal to the input beams.The brightest and largest number of array signals appeared when the plane of the polarization of the two input beams coincided with one of the crystal axes, as shown in Fig. 5(a).It need to clarify that the white color of the two input beams due to the over exposure of camera.There is no white light generation in the two input beams.It was found that there were four angles 0 o , 90 o , 180 o , and 270 o at which the sidebands were brightest, and four angles 45 o , 135 o , 225 o , and 315 o at which the sidebands were weakest, as shown in Fig. 5(b).The periodic array did not appear continuously but appeared aperiodically i.e. when the angle was rotated 6 o , 14 o , 0 o , -6 o , -11 o , -13 o , -16 o , and -19 o and further.At other angles, these regular arrays were replaced by a noise pattern, as shown in Fig. 5(c).The rotation angle of the sapphire plate also affected the position of the array signals.Figures 5(d) and 5(e) showed photographs when the sapphire plate was rotated by -16 o and 14 o , respectively.It can be seen that the column line was tilted a different direction when the sapphire plate was rotated.The spectrum of the tilted signal was slightly narrower in the shorter wavelength region when the signal was tilted in the high order direction, and vice versa.If the power of beam_2 was increased to 27 mW, bright stable "fish-like" multicolor arrays were observed, as shown in Fig. 5(f).The dependence of the B 1, 0 output power on the rotation angle θ of the sapphire plate was measured, as shown in Fig. 6.Here we only show the evolution in the rotation angle region from 0 o to 210 o .Clearly, the output power changed periodically with rotation of the sapphire plate.This periodic evolution was because of the periodic dependence of χ (3) (θ)∝ cos 2 (2θ) +1of the sapphire plate on the rotation angle θ, as shown by the dashed line in Fig. 6.The peak wavelength of the sidebands shifted continuously to a shorter wavelength by about 20 nm when the sapphire plate was rotated from 0 o to 45 o due to a phase matching condition.The polarization of the sidebands also changed with rotation of the sapphire plate.To detect the polarization rotation of the sidebands, a thin film polarizer was located normal to the sideband beams and placed 50 cm after the sapphire plate.The polarizations of the multicolored arrays were the same as those of the input beams when the plane of polarization of the input beams coincided with one of the crystal axes.Very little orthogonally polarized light was detected.As the sapphire plate was rotated from 0 o to 90 o , the orthogonally polarized light continuously increased from 0 o to 45 o and decreased from 45 o to 90 o , and light with parallel polarization showed the opposite effect.At 45 o , the sidebands were the weakest, and the orthogonally polarized light had power equal to the parallel polarized light.We rotated the thin film polarizer to minimize the intensity of light passing through the film polarizer.The rotation angle of the thin film polarizer β was in accordance with the rotation angle of the sapphire plate θ, and had the same angle when the sapphire plate was rotated from 0 o to 45 o and was 90 oθ when the sapphire plate was rotated from 45 o to 90 o .This rotation angle dependence phenomenon was also recently observed in supercontinuum generation process [21].The physical explanation of the phenomenon is complex because it included combined effects of the cascaded third-order nonlinear processes of FWM, coherent anti-Stokes Raman scattering (CARS), cross-phase modulation (XPM), and self-phase modulation (SPM) [22][23].Much detail physical explanation of this phenomenon will be explored in the future.

Experimental results and discussion
A half-wave plate or a quartz-wave plate can also be used in the path of one of the input beams to safely control the multicolored arrays.The polarization, intensity, and position of the multicolored arrays can be controlled by rotating the sapphire plate.Note that this phenomenon did not appear when a fused silica glass was used as the medium, due to its symmetric structure.This phenomenon was very easily repeated in the experiment, and the sapphire plate was not damaged over the course of the experiment.

Conclusion
In conclusion, an interesting 2-D nonlinear optical phenomenon was observed.Broadband femtosecond 2-D multicolored arrays with more than ten periodic columns and more than ten rows were generated in a sapphire plate.These stable 2-D multicolored arrays can be controlled by rotating the sapphire plate, a half-wave plate, or a quarter wave plate.The properties of the spectrum, spatial modal, power stability, and pulse duration of the 2-D multicolored arrays show they could be used in various applications, for example 2-D alloptical switching devices or multicolor pump-probe experiments.

Fig. 1 .
Fig. 1.(a) Schematics of the experimental setup for multicolored array generation.α is the crossing angle between the two input beams, beam_1 and beam_2 in the air.θ is the rotation angle of the sapphire plate.(b) A photograph of the 2-D multicolored array on a UV light sensitive plate.(c) Definition of 2-D multicolored arrays, where B0, 0 and B-1, 0 refer to two incident beams, beam_1 and beam_2, respectively.Schematics of the experimental setup for multicolored array generation are shown in Fig.1(a).The sapphire plate was a-cut and there were two orthogonal optical axes in the plane of the sapphire plate.The plane formed by the two orthogonal crystal axes was set normal to the two input beams.The two input beams had perpendicular polarizations and coincided with one of the crystal axes.The diameters of both incident beams on the surface of sapphire plate were 300 µm, as initially measured by a CCD camera (BeamStar FX 33, Ophir Optronics).The crossing angle α between the two input beams was 1.8 o ± 0.02 o in the air.When beam_1 and beam_2 were synchronously focused on the sapphire plate in both time and space, stable separate 2-D transverse multicolored array signals at different wavelengths were generated.The polarizations of the multicolored arrays were the same as those of the input beams as

Fig. 2 .
Fig. 2. (a) The spectra of array signals on the center row Bm, 0, where beam_1 and beam_2 are two incident beams.(b) The spectra of array signals on the second column B2, n.

Fig. 3 .
Fig. 3. (a)The spatial profiles of B0, 0, B1, 0, and B4, 1 in one dimension.The inset patterns are spatial modals of B0, 0, B1, 0, and B4, 1 from bottom to top, measured by a CCD camera.(b) The spatial profiles of B0, 0, B1, 0, and B4, 1 in one dimension with logarithmic scale in the intensity.The cyan dashed line is the Gaussian fit of B0, 0. The magenta dash-dotted line is the Lorentzian fit of B4, 1. (c) The retrieved XFROG pulse trace and phase of the B1, 0 with a retrieved error of 0.01022.The retrieved pulse duration is 35 ± 3fs.The inset pattern is the measured XFROG trace.

Fig. 4 .
Fig. 4. (a) The output power of array signals on the 0, ±1, and ±2 rows when the power of the two incident beams, beam_1 and beam_2, were 0.1mW and 25mW, respectively.Only the signals in the center row are marked with star symbols and dashed line, as shown on the left.The inset figure shows the dependence of the output power of different sidebands on the input power of beam_2.(b) The power stabilities of B4, 0, B3, 1, and B1, 0 monitored for 200 seconds, which were 1.25%RMS, 0.63%RMS, and 1.84% RMS, respectively.

Fig. 5 .
Fig. 5. Photographs of 2-D multicolored arrays on a sheet of white paper when (a) the plane of polarization of the two input beams coincided with one of the crystal axes; (b) the sapphire plate rotated for 45 o ; (c) noise pattern; (d) and (e) show photographs with the sapphire plate rotated by -16 o and 14 o , respectively.(f) A photograph of the 2-D multicolored arrays on a UV light sensitive plate when the input power of beam_2 was increased to 27mW.