Amplitude and phase reconstruction of ... - OSA Publishing

0 downloads 0 Views 392KB Size Report
presence of a photogalvanic nonlinearity or an externally applied electric field,” Phys. ... been observed in LiNbO3, due to the defocusing nature of the photovoltaic effect [8]. .... Fig.4(c) it is noticeable a central circular area in the phase map.
Amplitude and phase reconstruction of photorefractive spatial bright-soliton in LiNbO3 during its dynamic formation by digital holography M. Paturzo, L. Miccio, S. De Nicola, P. De Natale and P. Ferraro Istituto Nazionale di Ottica Applicata - CNR and LENS-European Laboratory for Nonlinear Spectroscopy, Via Campi Flegrei 34, 80078 Pozzuoli (Na), Italy [email protected]

Abstract: The time behaviour of bright spatial solitons in congruent undoped lithium niobate crystal is experimentally investigated. Full field characterization of the optical wavefront emerging from the crystal during the soliton formation process is performed by digital holographic method. Experimental results of the amplitude and phase maps of the field distribution at the exit face of the crystal allow the real-time monitoring of the evolution of the soliton beam from the application of the external field to the end of the process when the generation of the channel waveguide appears to be stable. The features of the dynamics of the soliton formation are visualized, analyzed and compared to a time-dependent numerical model. ©2007 Optical Society of America OCIS codes: (190.5330 ) Photorefractive nonlinear optics; interferometry;

(090.2880) Holographic

References and links 1.

M Segev, B Crossignani, A Yariv, and B Fisher, “Spatial solitons in photorefractive media,” Phys. Rev. Lett. 68, 923-926 (1992). 2. G. C. Duree, J. L. Shultz, G. J. S. Alamo, M. Segev, A. Yariv, B. Crossignani, P. Di Porto, E. Sharp and R. R. Neurgaonkar, “Observation of self-trapping of an optical beam due to the photorefractive effect,” Phys. Rev.Lett. 71, 533-536 (1993). 3. M. F. Shih, M. Segev, G. C. Valley, G. Salamo, B. Crossignani, and P. Di Porto, “Observation of 2dimensional steady-state photorefractive screening solitons,” Electron Lett. 31, 826-827 (1995) 4. J. Petter and C. Denz, “Guiding and dividing waves with photorefractive solitons,” Opt. Commun, 188, 5561 (2001). 5. J. W. Fleiscer, T. Carmon, M. Segev, N. K. Efremedis, and D. N. Christodoulides, “Observation of discrete solitons in optically induced real time waveguide arrays,” Phys Rev. Lett 90, 023901 (2003). 6. E. Fazio, F. Renzi, R. Rinaldi, M. Bertolotti, M. Chauvet, W.Ramadan, A. Petris, and V. I. Vlad, “Screeningphotovoltaic bright solitons in lithium niobate and associated single-mode waveguides,” Appl. Phys. Lett. 85, 2193–2195 (2004). 7. G. Couton, H. Maillotte, R. Giust, and M. Chauvet, “Formation of reconfigurable singlemode channel waveguides in LiNbO3 using spatial solitons,” Electron. Lett. 39, 286-287 (2003). 8. M. Taya, M. C. Bashaw, M. M. Fejer, M. Segev, and G. C. Valley, “Observation of dark photovoltaic spatial solitons,” Phys. Rev. A 52, 3095-3100 (1995). 9. S. Mailis, C. Riziotis, I. T. Wellington, P. G. R. Smith, C. B. E. Gawith, and R. W. Eason, "Direct ultraviolet writing of channel waveguides in congruent lithium niobate single crystals ," Opt. Lett. 28, 1433-1435 (2003). 10. I. E. Barry, G. W. Ross, P. G. R. Smith, and R. W. Eason, “Ridge waveguides in lithium niobate fabricated by differential etching following spatially selective domain inversion” Appl. Phys. Lett. 74, 1487-1488 (1999). 11. J. L. Jackel, C. E. Rice, and J. J. Veselka, “Proton exchange for high-index waveguides in LiNbO3” Appl. Phys. Lett. 41, 607-608 (1982). 12. S. Grilli, P. Ferraro, S. De Nicola, A. Finizio, G. Pierattini and R. Meucci, “Whole optical wavefields reconstruction by digital holography,” Opt. Express 9, 294-302 (2001).

#81427 - $15.00 USD

(C) 2007 OSA

Received 26 Mar 2007; revised 3 May 2007; accepted 8 May 2007; published 18 Jun 2007

25 June 2007 / Vol. 15, No. 13 / OPTICS EXPRESS 8243

13. R. Osellame, N. Chiodo, V. Maselli, A. Yin, M. Zavelani-Rossi, G. Cerullo, P. Laporta, L. Aiello, S. De Nicola, P. Ferraro, A. Finizio, and G. Pierattini, "Optical properties of waveguides written by a 26 MHz stretched cavity Ti:sapphire femtosecond oscillator," Opt. Express 13, 612-620 (2005) 14. S. Grilli, P. Ferraro, M. Paturzo, D. Alfieri, P. De Natale, M. de Angelis, S. De Nicola, A. Finizio, and G. Pierattini, "In-situ visualization, monitoring and analysis of electric field domain reversal process in ferroelectric crystals by digital holography," Opt. Express 12, 1832-1842 (2004). 15. M. Paturzo, P. Ferraro, S. Grilli, D. Alfieri, P. De Natale, M. de Angelis, A. Finizio, S. De Nicola, G. Pierattini, F. Caccavale, D. Callejo, and A. Morbiato, "On the origin of internal field in Lithium Niobate crystals directly observed by digital holography," Opt. Express 13, 5416-5423 (2005). 16. N. Fressengeas, J. Maufoy and G. Kugel, “Temporal behaviour of bidimensional photorefractive bright spatial solitons,” Phys. Rev E. 54, 6866-6875 (1996). 17. A. A. Zozulya and D. Z. Anderson, “Propagation of an optical beam in a photorefractive medium in the presence of a photogalvanic nonlinearity or an externally applied electric field,” Phys. Rev. A 51, 1520-1531 (1995). 18. A. A. Zozulya, D. Z. Anderson, A. V. Mamaev, and M. Saffman, “Solitary attractors and low-order filamentation in anisotropic self-focusing media,” Phys. Rev. A 57, 522-534 (1998). 19. M Chauvet, V Coda, M. Maillotte, E. Fazio, and G. Salamo “Large self-deflection of soliton beams in LiNbO3,” Opt. Lett. 30, 1977-1979 (2005).

1. Introduction Photorefractive nonlinear medium provides an intensity-dependent, saturable nonlinear response which can be conveniently exploited for the purpose of spatial-soliton generation. The existence of photorefractive (PR) spatial soliton, first predicted in Ref. [1] has stimulated substantial research interest because of an attractive property of these spatial solitons in photorefractive materials, that is the very low laser power, in the range of mW and lower, necessary for their generation and manipulation. The investigation of PR spatial solitons is thus accessible with conventional continuouswave laser sources and allows to foresee a number of realist applications in waveguiding, switching and optical information processing systems, owing to the compactness and relatively low cost of the experimental devices. PR spatial solitons are based on an optical nonlinearity which results from charge excitation, transport and trapping and gives rise to a space-charge filed that in turn modulates the refractive index by means of the electro-optic (i.e. Pockels) effect. Although the first observations [2, 3] and advanced experiments [4, 5] were realized in strontium barium niobate (SBN) crystal, it was recently demonstrated that lithium niobate (LiNbO3 ) can be efficiently employed for forming narrow bright screening spatial solitons [6]. At first planar and channel waveguides in LiNbO3 has been produced by dark solitons, that leave the material modified for a long time [7]. Dark photovoltaic (PV) solitons have been observed in LiNbO3, due to the defocusing nature of the photovoltaic effect [8]. However has been observed that bright solitons are much more stable than dark solitons and offer the simplest way to photo-induce slab and channel waveguides. Lithium niobate is available in good optical quality and it is a widely used material for optoelectronic applications, with large electro-optic, nonlinear and acusto-optic coefficients. The possibility of generating and employing spatial solitons in LiNbO3 for writing optimum single-mode waveguide memorized for long time, represents a convenient low-power light technique [6], advantageous compared to other light-induced techniques such as direct UV writing [9], micromachining [10] or ion and proton exchange [11]. An important issue for complete solitonic waveguide characterization in photorefractive media, such as lithium niobate, and a prerequisite for understanding the physical mechanisms underlying the writing process is an accurate investigation of the dynamics of the soliton formation leading to the refractive index change across the waveguide. Complete characterization of the dynamics of the soliton development requires measuring not only the output light intensity distribution but also the wavefield phase at the exit plane of the photorefractive crystal. To the best of our knowledge, a technique for the simultaneous measurement of the transient dynamics of both the intensity and phase of the output beam propagating in the photorefractive medium has not been presented before in the literature. We

#81427 - $15.00 USD

(C) 2007 OSA

Received 26 Mar 2007; revised 3 May 2007; accepted 8 May 2007; published 18 Jun 2007

25 June 2007 / Vol. 15, No. 13 / OPTICS EXPRESS 8244

show here how reconstruction of amplitude and phase of the wavefront is possible as it emerges from the crystal exit face while the writing process of the soliton formation takes place. We were able to follow the dynamics of the process of the wavefront at the output of an externally biased photorefractive LN crystal by means of an interferometric method based on digital holography (DH) technique [12]. Recently, DH has been employed as a full field, high spatial resolution technique for measuring stationary refractive index profiles of femtosecond laser written waveguides [13] or for in-situ characterization of electric poling of LN z-cut crystals or for measuring internal fields [14,15]. For studying the temporal behaviour of photorefractive spatial solitons sequences of digital holograms have been recorded and numerically reconstructed to obtain quantitatively amplitude and phase maps evolution during the process of formation. Comparison of the experimental results with numerical simulations within the framework of a time-dependent one dimensional band transport model is also presented and discussed. 2. Experimental setup The experimental set-up to generate and to study the dynamics of the soliton formation is shown in Fig. 1. It is based on a classical Mach-Zehnder configuration in which a laser source emitting at a wavelength of λ = 514 nm from a Ar+ laser is split into two beams, the reference and the object beam respectively. The object beam is directly focused with beam-waist of 29 μm onto the input face of a congruent undoped LiNbO3 crystal with dimensions 10mm × 7mm × 0.5(c) mm (Altechna Co. Ltd.). The object beam is linearly polarized along the crystallographic Z-axis and its intensity on the entrance face of the crystal is Imax = 2.13 W/cm2. The crystal is externally biased along the Z direction by a static, high voltage field of 3.5 kV/mm and illuminated by a background blue laser emitting at a wavelength of 473nm, along the Y crystallographic direction. The background beam provides a uniform illumination of the sample with an intensity Ib =1.2 mW/cm2. The object beam propagates trough the crystal length along the X-axis. The wave front emerging at the exit face of the crystal, orthogonal to the X axis, is imaged by a 10 × microscope objective (MO) in a plane at a distance d in front of the CCD array plane which records the interference pattern between the plane wave reference beam and the object beam. PBS

PH

M

Laser

λ/2

λ/2

MO Reference beam

collimating lens

background beam

Object beam

Image Plane

CCD Plane

S M

MO HV

BS

d

Fig. 1. Scheme of set-up used for soliton formation

We made use of DH in microscope configuration to reconstruct the complex wave field UL(Y,Z)= |UL(Y,Z,t)|exp[i φL(Y,Z,t)] at the exit plane X=L=7mm of the crystal as function of time t. Whole wave field reconstruction during the formation process of the photorefractive bright soliton means that both the time dependent intensity IL(Y,Z,t)=|UL(Y,Z,t)|2 and phase φL(Y,Z,t) distributions are numerically reconstructed from a sequence of digital holograms recorded at different instant of time, after the externally applied voltage is suddenly switched

#81427 - $15.00 USD

(C) 2007 OSA

Received 26 Mar 2007; revised 3 May 2007; accepted 8 May 2007; published 18 Jun 2007

25 June 2007 / Vol. 15, No. 13 / OPTICS EXPRESS 8245

on. The numerical reconstruction process of the digital holograms is based on the Fresnel transformation method [12]. DH is an imaging method in which the hologram resulting from the interference between the reference and the object complex fields, w1(x, y) and w2(x, y) respectively, is recorded with a CCD camera and reconstructed numerically. The hologram is multiplied by the reference wavefield in the hologram plane, namely the CCD plane, to calculate the diffraction pattern in the image plane (see Fig.1). The reconstructed field Γ(ν , μ ) in the image plane, at a distance d from the CCD plane along the beam path, is obtained by using the Fresnel approximation of the Rayleigh-Sommerfield diffraction formula ⎡ iπ

Γ(ν , μ ) ∝ ∫∫ h(ξ ,η )r (ξ ,η )exp⎢ ⎣

(ξ λd

2

)⎤

+ η 2 ⎥ exp[− 2iπ (ξν + ημ )]dξ dη

(1)



where r (ξ , η ) ≡ w1 (ξ , η ) is the reference wave which in the case of a plane wave is simply given by a constant value, h(ξ , η ) = w1 (ξ , η ) + w 2 (ξ , η ) is the hologram function, λ is the laser source wavelength and d is the reconstruction distance, namely the distance measured between the object and the CCD plane along the beam path. The coordinates (ν,μ) are related to the image plane coordinates (x’, y’) by ν = x' λ d and μ = y' λ d . It is clear that the object field Γ(ν,μ) is the Fourier transform of the hologram h(ξ , η ) multiplied by the reference wave 2

r(ξ,η) and the chirp function exp[(iπ λ d )(ξ 2 + η 2 )] . The pixel size (Δx' , Δy ') in the image plane is related to that (Δξ , Δη ) of the CCD array through the equations Δx1 =λd / N∆ξ and Δy1 =λd / N∆η where N is the pixel number of the CCD array. The great advantage of this technique is that it allows to reconstruct numerically the complex field of the object beam. The two-dimensional amplitude A(x’, y’) and phase φ(x’, y’) distributions of the object wavefield can be re-imaged by using one hologram acquisition and performing simple calculations on the object wavefield Γ(ν,μ) reconstructed from the numerical solution of the diffraction problem in the computer: A(x' , y ') = abs[Γ(x' , y ')]; φ (x' , y ') = arctan

Im[Γ(x' , y ')] Re[Γ(x' , y ')]

(2)

3. Experimental results Figures 2-4 show the temporal evolution of the field intensity and phase. The spatial resolution of the reconstructed wave front is determined by the size λd / N∆ξ of the reconstruction pixel, where the ∆ξ = 6.7 µm is the pixel size of the CCD detector and N × N = 512 × 512 is the number of pixels composing the digital holograms. For a reconstruction distance d = 100 mm, the corresponding size of the reconstruction pixel is 18.50 µm, that scaled for the actual magnification M = 24.5 gives a final value of 0.75 µm for the spatial resolution of the reconstructed maps. It important to note that by means of DH it is possible to extract information both in amplitude and phase of the reconstructed optical field under investigation. In Fig. 2(a) is displayed the Gaussian beam intensity map after propagation inside the crystal before the external voltage is applied. Figures 2(b)-2(c) shows intermediate steps at t=84 s and t=172s of the soliton formation during the application of the bias field. The voltage was switched off at a time t=194s. In order to show that the written waveguide lasted longer than the writing time, we characterized the modal profile of the solitonic waveguide at about 20 minutes after the external filed was turned off. Figure 2(d) shows the intensity distribution form retrieved from numerically reconstructed hologram that was recorded at t =1482 s.

#81427 - $15.00 USD

(C) 2007 OSA

Received 26 Mar 2007; revised 3 May 2007; accepted 8 May 2007; published 18 Jun 2007

25 June 2007 / Vol. 15, No. 13 / OPTICS EXPRESS 8246

(a)

(b)

(c)

(d)

Fig. 2. Maps of beam intensity during the soliton formation

Figure 3 shows a comparison between the initial and final amplitude profile, along z axis. In fact, the beam profile on the crystal exit face for t=0 s (blue line) and t =1482 s (black line) is shown. Experimental results are fitted with expected theoretical curve (red lines). The experimental curve for t=0 is fitted with a Gaussian curve shape giving the expected information that the actual profile is Gaussian before the soliton formation process takes place. On contrary at t =1482 the fitting curve is in good approximation a squared-hyperbolic secant confirming that the soliton regime has been achieved [6].

Fig. 3 Intensity profile of light beam after and before soliton formation

Figures 4(a)-4(d) display the corresponding time evolution for the phase distribution at the exit face of the crystal. The retrieved phase maps are shown at the same time as in Fig. (3). In Fig.4(c) it is noticeable a central circular area in the phase map. The optical field is clearly confined where the soliton waveguide is written in respect to the surround area. Moreover, it

#81427 - $15.00 USD

(C) 2007 OSA

Received 26 Mar 2007; revised 3 May 2007; accepted 8 May 2007; published 18 Jun 2007

25 June 2007 / Vol. 15, No. 13 / OPTICS EXPRESS 8247

is important to note that the phase map remains stable even after the voltage is switched off [see Fig 4(d)].

(a)

(b)

(c)

(d)

Fig. 4. Phase maps of laser beam during soliton formation

However we have not already investigated completely how much time it remains stable and if changes occurs during the expected decay process. Nevertheless that topic will be the subject of further studies. Figures 5(a) and 5(b) report the corresponding movies displaying the evolution of the soliton formation respectively in amplitude and phase.

Fig. 5. (780 KB- 225KB ) Movie of intensity (a) and phase (b) maps during soliton formation

4. Numerical simulations and comparison with the experimental data Numerical simulations of the dynamics of the soliton formation, with the beam propagation method, were based on the non-stationary 1+1D ( one propagation direction and one diffraction direction ) band transport model developed in ref. [16]. Compared to other approaches, such as the general treatment of Zozulya and Anderson [17, 18] that requires solving globally an elliptic-type equation for an anisotropic electrostatic potential in the

#81427 - $15.00 USD

(C) 2007 OSA

Received 26 Mar 2007; revised 3 May 2007; accepted 8 May 2007; published 18 Jun 2007

25 June 2007 / Vol. 15, No. 13 / OPTICS EXPRESS 8248

photorefractive medium, this simplified model leads to an explicit expression for the change of the induced refractive index Δnb(t) = -(1/2) nb3r33Esc(t) in terms of the time dependent local space-charge field Esc(t),namely



Esc (t ) = ⎜⎜

(E0 + E ph )I + In



1 ⎛ kbT ∂ I ⎞ ⎞ ⎜ ⎟ ⎟ × exp[− (I nt I d Td )] I n ⎝ e ∂Z ⎠ ⎟⎠

(3)

I + Ib 1 ⎛ k T ∂I ⎞ I + E0 d − E ph − ⎜ b ⎟ In I n I n ⎝ e ∂Z ⎠

where nb =2.2 is the extraordinary refractive index, r33 =32pm/V the effective electro-optic coefficient of LiNbO3 crystal for the geometry under consideration, Eph is the photovoltaic field, assumed to be directed along the crystal Z axis, Eo is the amplitude of the external electric field, kb the Bolzmann constant, T the crystal temperature, e the electron charge, Td the dielectric response time of the crystal medium in the dark, Id the equivalent dark irradiance, with Id