Complete spatial and temporal locking in phase ... - OSA Publishing

Complete spatial and temporal locking in phasemismatched second-harmonic generation Eugenio Fazio,1,* Federico Pettazzi,1,4 Marco Centini,1 Mathieu Chauvet,2 Alessandro Belardini,1 Massimo Alonzo,1 Concita Sibilia,1 Mario Bertolotti,1 and Michael Scalora3 1

Ultrafast Photonics Laboratory, Dipartimento di Energetica and CNISM, Sapienza Università Roma, via Scarpa 16, I-00161 Roma Italy-UE Département d'Optique, Institut FEMTO-ST, UMR CNRS 6174, Université de Franche-Comté, F-25030 Besançon, France-UE 3 Charles M. Bowden Research Facility, RDECOM, US Army Aviation and Missile Command, Redstone Arsenal, AL35803, USA 4 Currently with the European Space Agency, ESTEC - TEC-MME, Keplerlaan 1, 2200 AG Noordwijk, The Netherlands-EU * Corresponding author: [email protected] 2

Abstract: We experimentally demonstrate simultaneous phase and group velocity locking of fundamental and generated second harmonic pulses in Lithium Niobate, under conditions of material phase mismatch. In phasemismatched, pulsed second harmonic generation in addition to a reflected signal two forward-propagating pulses are also generated at the interface between a linear and a second order nonlinear material: the first pulse results from the solution of the homogeneous wave equation, and propagates at the group velocity expected from material dispersion; the second pulse is the solution of the inhomogeneous wave equation, is phaselocked and trapped by the pump pulse, and follows the pump trajectory. At normal incidence, the normal and phase locked pulses simply trail each other. At oblique incidence, the consequences can be quite dramatic. The homogeneous pulse refracts as predicted by material dispersion and Snell’s law, yielding at least two spatially separate second harmonic spots at the medium’s exit. We thus report the first experimental results showing that, at oblique incidence, fundamental and phase-locked second harmonic pulses travel with the same group velocity and follow the same trajectory. This is direct evidence that, at least up to first order, the effective dispersion of the phase-locked pulse is similar to the dispersion of the pump pulse. 2009 Optical Society of America OCIS codes: (190.0190) Nonlinear optics. (190.2620) Harmonic generation and mixing. (190.4350) Nonlinear optics at surfaces. (190.5530) Pulse propagation. (190.7110) Ultrafast nonlinear optics.

References and links 1. 2. 3. 4. 5. 6. 7.

J. A. Armstrong, N. Bloembergen, J. Ducuing, and P. S. Pershan, “Interactions between Light Waves in a Nonlinear Dielectric”, Phys. Rev. 127, 1918 (1962). N. Bloembergen and P. S. Pershan, “Light Waves at the Boundary of Nonlinear Media”, Phys. Rev. 128, 606 (1962). P. D. Maker, R. W. Terhune, N. Nisenoff, and C. M. Savage, “Effects of Dispersion and Focusing on the Production of Optical Harmonics”, Phys. Rev. Lett. 8, 21 (1962). N. Bloembergen, H. J. Simon, and C. H. Lee, “Total Reflection Phenomena in Second-Harmonic Generation of Light”, Phys. Rev. 181, 1261 (1969). A. Feisst and P. Koidl, “Current induced periodic ferroelectric domain structures in LiNbO3 applied for efficient nonlinear optical frequency mixing”, Appl. Phys. Lett. 47, 1125 (1985). U. Osterberg and W. Margulis, “Dye laser pumped by Nd: YAG laser pulses frequency doubled in a glass optical fiber”, Opt. Lett. 11, 516 (1986). G. A. Magel, M. M. Fejer, and R. L. Byer, “Quasi-phase-matched second-harmonic generation of blue light in periodically poled LiNbO3”, Appl. Phys. Lett. 56, 108 (1990).

#103628 - $15.00 USD

(C) 2009 OSA

Received 3 Nov 2008; revised 11 Dec 2008; accepted 20 Jan 2009; published 17 Feb 2009

2 March 2009 / Vol. 17, No. 5 / OPTICS EXPRESS 3141

8. 9. 10. 11.

12. 13.

14.

15. 16. 17. 18.

19.

20.

21.

S. K. Kurtz and T. T. Perry, “A Powder Technique for the Evaluation of Nonlinear Optical Materials”, J. Appl. Phys. 39, 3798 (1968). J. P. van der Ziel, “Phase−matched harmonic generation in a laminar structure with wave propagation in the plane of the layers“, Appl. Phys. Lett. 26, 60 (1976). K. Sakoda and K. Othaka, “Sum-frequency generation in a two-dimensional photonic lattice”, Phys. Rev. B 54, 5742 (1996). R. J. DeSalvo, D. J. Hagan, M. Sheik-Bahae, G. Stegeman, E. W. Van Stryland, and H. Vanherzeele, “Self-focusing and self-defocusing by cascaded second-order effects in KTP”, Opt. Lett. 17, 28 (1992). N. C. Kothari and X. Carlotti, “Transient second-harmonic generation: influence of effective groupvelocity dispersion”, J. Opt. Soc. Am. B 5, 756 (1988). R. Maleck Rassoul, A. Ivanov, E. Freysz, A. Ducasse, and F. Hache, “Second-harmonic generation under phase-velocity and group-velocity mismatch:influence of cascading self-phase and cross-phase modulation”, Opt. Lett. 22, 268 (1997). S. Cussat-Blanc, R. Maleck-Rassoul, A. Ivanov, E. Freysz, and A. Ducasse, “Influence of cascading phenomena on a type I second-harmonic wave generated by an intense femtosecond pulse: application to the measurement of the effective second-order coefficient”, Opt. Lett. 23, 1585 (1998). E. Fazio, M. Zitelli, S. Dominici, C. Sibilia, G. D'Aguanno, and M. Bertolotti, “Phase-driven pulse breaking during perfectly-matched second harmonic generation”, Opt. Commun. 148, 427 (1998). D. Noordam, H.J. Bakker, M.P. de Boer, H.B. van Linden van den Heuvell, “Second-harmonic generation of femtosecond pulses: observation of phase-mismatch effects”, Opt. Lett. 15, 1464 (1990) M. Mlejnek, E. M. Wright, J. V. Moloney, and N. Bloembergen, “Second Harmonic Generation of Femtosecond Pulses at the Boundary of a Nonlinear Dielectric”, Phys. Rev. Lett. 83, 2934 (1999). W. Su, L. Qian, H. Luo, X. Fu, H. Zhu, T. Wang, K. Beckwitt, Y. Chen, and F. Wise, “Induced group-velocity dispersion in phase-mismatched second-harmonic generation”, J. Opt. Soc. Am. B 23, 51 (2006). V. Roppo, M. Centini, C. Sibilia, M. Bertolotti, D. de Ceglia, M. Scalora, N. Akozbek, M. J. Bloemer, J.W. Haus, O. G. Kosareva, V.P. Kandidov, “Role of phase matching in pulsed second-harmonic generation: Walk-off and phase-locked twin pulses in negative-index media”, Phys. Rev. A 76, 033829 (2007) M. Centini, V. Roppo, E. Fazio, F. Pettazzi, C. Sibilia, J.W. Haus, J.V. Foreman, N. Akozbek, M.J. Bloemer, M. Scalora, “Inhibition of Linear Absorption in Opaque Materials Using Phase-Locked Harmonic Generation”, Phys. Rev. Lett. 101, 113905 (2008) P. D. Marker, R. W. Terhune, M. Nisenoff and C. M. Savage, “Effects of Dispersion and Focusing on the Production of Optical Harmonics”, Phys. Rev. Lett. 8, 21 (1962).

1. Introduction The second harmonic generation process has been widely investigated beginning with seminal contributions by authors [1-4] that derived the now well-known theory described in every nonlinear optics textbook. Typically, second harmonic generation processes are investigated and used at or near the phase-matching condition, so that the pump and generated frequencies travel inside the nonlinear material experiencing similar propagation velocities, and yielding as large a conversion efficiency as possible. Outside of the phase matching condition conversion efficiencies decrease rapidly [1]. Since the achievement of large conversion efficiency was the most sought-after outcome, the study of propagation phenomena under phase mismatched condition was seldom tackled. Investigations typically focused on the achievement of two main goals: a) compensation of material dispersion by using materials with alternate orientations of poling [5-7], and b) realization of new materials, either organic [8] or structured [9-10], in order to achieve a larger effective nonlinear coefficient. Parametric processes were also investigated for signal processing purposes by observing that a small phase-mismatch between generating and generated beams can induce phasemodulation of the fundamental frequency by means of cascaded parametric generation [11] as well as pulse breaking [12-15]. However, even for large phase-mismatch Noordam et al [16] theoretically and experimentally demonstrated that two second harmonic pulses are generated. Theoretical aspects of such phenomena have also been discussed by M. Mlejnek et al [17] in a birefringent semiconductor and by Su et al [18], both experimentally and theoretically, in BBO crystals. In all cases one finds that the second harmonic signal breaks-up to form two distinct pulses. #103628 - $15.00 USD

(C) 2009 OSA



It was not until recently, however, that the generation of two distinct components was understood as a phase locking mechanism that causes the pump to trap portions of the harmonics [19], an effect later experimentally confirmed even in absorbing materials [20]. In ref.[20] an investigation of harmonic generation was conducted at normal incidence in a GaAs substrate, at wavelengths below the band edge (~890nm). A 100-fs pump pulse tuned to 1300nm generated signals at 650nm (SH) and 433nm (TH). The harmonic signals then became trapped by the pump pulse and subsequently propagated through the 450-micron thick GaAs sample without being absorbed, notwithstanding the fact that typical absorption lengths in GaAs below the band edge are much smaller than 100nm. In contrast, both second and third harmonic homogenous pulses were immediately absorbed. Although at first counter-intuitive, the suppression of absorption that portions of the harmonic signals experience can be understood within the frame-work of the phase locking mechanism. In refs.[19] it was shown theoretically that the real part of the effective index of refraction that the harmonics experience is equal to that of the pump. This may be done by performing a full spectral decomposition of the fields, and by calculating the effective index of refraction neff=c/ as the ratio of expectation values ω =∞

k =∞

< k >=

∫

k | H (ω , k ) | dk and < ω > = 2

k =−∞

∫ ω

ω | H (ω , z ) |2 dω

[19].

The E field

=−∞

may also be used in the calculation. Once the real part of the index neff is mapped onto the pump index of refraction, the imaginary part of the index is obtained by performing a Kramers-Kronig reconstruction, which in turn yields an imaginary part for the second harmonic field that is identical to the index of refraction for the pump field. Therefore, it is sufficient to tune the pump to a region of transparency to guarantee suppression of absorption at the harmonic wavelengths. In the present work we theoretically predict and experimentally demonstrate that under conditions of a phase mismatch and oblique incidence, away from absorption resonances, the phase locked pulses display the same phase and group velocities as the pump field. This determination is made by performing second harmonic generation at 400 nm in a lithium niobate crystal using femtosecond pulses, and by monitoring pulse trajectories and group velocities inside the sample. 2. Numerical simulations To model oblique incidence we use a numerical model similar to that used in ref.[20], with the addition of a transverse coordinate to include diffraction and propagation on a twodimensional plane. We use a Lorentz oscillator model to describe material dispersion, and for simplicity we assume a TM-polarized incident pump field. The fields may be written as a superposition of harmonics as follows: ∞

ℓ =1 ∞

∞

(

) = xˆ ∑ (H

(

)

H = xˆ ∑ H x lω ( z , y, t ) + c.c

lω x

( z , y, t ) eiℓ ( kz −ωt ) + c.c

ℓ =1 ∞

(

E = yˆ ∑ E y lω ( z, y, t ) + c.c + zˆ ∑ E z lω ( z , y, t ) + c.c ℓ =1

ℓ =1

∞

)

) = yˆ ∑ ( E ℓ =1

(1) ℓω y

)

∞

(

( z , y, t ) + c.c + zˆ ∑ Ezℓω ( z , y, t ) + c.c ℓ =1

)

where Eyℓω ( z , y , t ), H yℓω ( z , y , t ) are generic, spatially- and temporally-dependent, complex envelope functions; k and ω are carrier wave vector and frequency, respectively, and l is an integer. Eqs.(1) are a convenient representation of the fields, and no a priori assumptions are made about the field envelopes. We have also assumed that a TM-polarized incident field generates similarly polarized harmonics. The linear response of the medium is described by a ω 2p Lorentz oscillator model: ε (ω ) = 1 − , and µ (ω ) = 1 , where γ , ω p , and ωr are the ω 2 + iγω − ωr2

damping coefficient, the plasma and resonance frequencies, respectively. The second order nonlinear polarization is: PNL = χ (2) E 2 . Expanding the field into its components yields nonlinear polarization terms at the fundamental and second harmonic frequencies: #103628 - $15.00 USD

(C) 2009 OSA



(2) 2 Pω ( z , t ) = 2 χ ω(2) Eω* E2ω and P2ω ( z , t ) = χ 2ω Eω .

Assuming that polarization and currents may be decomposed as in Eqs.(1), we obtain the following Maxwell-Lorentz system of equations for the l th field components, in the scaled two-dimensional space ( ~y , ξ ) plus time ( τ ) coordinate system: ℓω ∂H xℓω ∂E ℓω ∂E y = iβ ℓ (H xℓω + E zℓω sin θ i + E yℓω cos θi ) − z + ∂τ ∂yɶ ∂ξ ∂E yℓω ∂τ

(

)

ω = i β ℓ E yℓω + H xℓω cos θ i − 4π ( J yℓω − iβ ℓ Pyℓω ) + i 4πβℓ P yℓ, NL − 4π

ω ∂P yℓ, NL

∂τ

+

∂H xℓω ∂ξ

∂P ℓω ∂H xℓω . ∂E ω = i β ℓ E zℓω + H xℓω sin θi − 4π ( J zℓω − i βℓ Pz ℓω ) + i 4πβ ℓ Pz ,ℓNL − 4π z , NL − ∂τ ∂τ ∂yɶ ℓω z

∂J yℓω

(

)

(

)

πω p2 ℓω E ω02 y

(

)

πω p2 ℓω E ω02 z

= 2i βℓ − γ

(

)J

ℓω y

+ βℓ2 + iγβ ℓ − β r2 P yℓω +

∂J zℓω = 2i βℓ − γ ∂τ

)J

ℓω z

+ βℓ2 + iγβ ℓ − β r2 Pz ℓω +

∂τ

(

∂P yℓω

(2)

∂Pz ℓω ∂τ ∂τ In Eqs.(2), the functions J,

J yℓω =

J zℓω =

;

P, PNL refer to linear electric currents, polarization, and nonlinear polarization, respectively. The coordinates are scaled so that ξ = z / λ0 , yɶ = y / λ0 , τ = ct / λ0 , ω0 =

2π c , where λ0 = 1µ m is just reference wavelength; γ , βℓω = 2π ℓω / ω0 ,

λ0

β r = 2πωr / ω0 , ω p , are the scaled damping coefficient, wave-vector, resonance and electric plasma frequencies for the l th harmonic, respectively. θ i is the angle of incidence of the pump field. The equations are solved using a split-step, fast Fourier transform-based pulse propagation algorithm [20] that advances the fields in time. The basis for understanding the generation of two distinct signals can be found in reference [2], and as mentioned earlier it is based on the mathematical solution of the homogenous and inhomogeneous wave equations at the second harmonic frequency. Continuity of the tangential components of all the fields at the boundary, a condition automatically and dynamically included in Eqs.(2), guarantees the generation of the two forward-propagating components that interfere in the vicinity of the entry surface and give rise to Maker fringes [21]. In Fig.1 we show the results of a simulation performed in a generic Lorentz medium 80microns thick, where a pump pulse approximately 100fs in duration generates two distinct second harmonic pulses that follow independent trajectories: one follows the pump (yellow arrow), the other refracts according to Snell’s law (red arrow), so that at the medium’s exit two spatially separate spots are clearly visible. From the figure it is easy to see that the SH pulse described by the yellow arrow spatially and temporally overlaps the pump pulse at all times. In contrast, the homogenous pulse, described by the red arrow, propagates much more slowly and exits the medium long after the phase locked pulse finds its way out.

#103628 - $15.00 USD

(C) 2009 OSA



Transverse Coordinate (µ µ m)

110 90 70 50 30 10

(a)

1060

110

160

110 90 70 50 30

(b)

10 1060

110

160

Longitudinal Coordinate (µ µm) Fig. 1. Numerical simulation of pulsed second-harmonic in a generic Lorentz medium, having γ=10-8, ωp=4, ω0=4, under phase-mismatched conditions. The small γ keeps absorption at negligible levels. The yellow box delineates the medium. (a) pump pulse. (b)SH pulses. Two SH signals are discernable, one that tracks the pumps pulse (yellow arrow), the other that refracts according to material dispersion (red arrow). These two components travel with different group velocities and thus tend to separate as distance is gained inside the sample. The leading pulse is phase-locked and trapped by the pump; the second pulse propagates freely and at a group velocity approximately three times smaller compared to the pump.

3. Experiment

The pulse locking was experimentally monitored by performing a SHG experiment using the set-up shown in fig.2. 110-fs light pulses from an amplified Ti-Sapphire laser system with energies of few ten’s of micro-joule are focused onto the XZ face of a congruent lithium niobate crystal (from a commercial z-cut wafer) and propagated along the Y direction. At the entrance face the beam waist was about 600 µm, with intensities of the order of a few GW/cm2. A special prismatic shape of the sample ensured the output face to be 20° tilted from the propagation direction. Thus, this sample acts simultaneously as a generator and a spectrum analyser, as each wavelength refracts at different angle and according to Snell’s law.

Fig. 2. The experimental set-up. A prismatic lithium niobate crystal was used to generate the second harmonic 400 nm signal from a 800 nm pump pulse. The titled output face of the crystal forced the pulses to exit at different angles, according to their refractive indices. The locked arm was then refocused onto a thin BBO crystal at perfect phase-matching for ω+2ω=3ω interaction.

#103628 - $15.00 USD

(C) 2009 OSA



4. Phase-velocity locking

The phase-velocity locking was initially monitored by recording the output angular distribution with a CCD camera in the reference position. Three cases were recorded according to the injected input beam.

Fig. 3. Experimental images of the fundamental (800 nm) and second harmonic (400 nm) signals exiting the tilted output surface. Each position corresponds to a different output angle according to the Snell’s law.

In Figure 3(I) we describe the case of injecting a pump e-polarised beam (A), and generating two e-polarised second-harmonic pulses, thus achieving type-0 polarisation coupling ee-e (∆kee-e=2.45·104 cm-1). One of the beams (B) exits at the angle that corresponds to ne(400nm), while the other (C) exits exactly overlapped with the pump beam, i.e. at the angle that corresponds to ne(800nm). This result demonstrates that the generated second harmonic signal separates in two different parts, one freely propagating without further interaction with the pump pulse, the other locked with the pump, propagating exactly at the same phase-velocity (A-C), i.e. having exactly the same refractive index and following exactly the same trajectory as the pump. It should be mentioned here that an unexpected high conversion efficiency was observed for the locked second-harmonic beam, whose average power was of the order of 10-3-10-2 of the pump one. In Figure 3(II) we show the case of an o-polarised fundamental frequency injected in the crystal (D), resulting in type I coupling oo-e (∆koo-e=1.20·104 cm-1). The generated SH signal is once again broken in two fragments, one freely propagating (E), and the other (F) once again phase-locked to the pump (D-F) that exits the lithium niobate crystal at the same angle as the fundamental frequency. Also in this case there is a transfer of refractive index from the generating to the generated wavelength: in fact the e-polarised locked pulse at 400nm experiences exactly the ordinary refractive index of 800nm pump pulses. Finally in Figure 3(III) we display the case when both e- and o-polarisations are injected inside the sample (G, H), such that type-0, type-I and type-II (oe-e, ∆koe-e=1.8·104 cm-1 , oe-o, ∆koe-o=3.5·104 cm-1) couplings are obtained simultaneously. The fundamental e-polarisation (G) generates two e-polarised second harmonic pulses (M, O) according to the type-0 coupling. The fundamental o-polarisation generates two e-polarised second harmonic pulses according to type-I coupling, (M, N). The temporal overlapping of the fundamental o- and epolarisations generates a second harmonic signal (L) by type-II coupling, which is the only second harmonic o-polarisation. Indeed, due to group velocity mismatch of the two pumps, #103628 - $15.00 USD

(C) 2009 OSA



the phase-locked component cannot be sustained because the two cross-polarised pump pulses rapidly separate and the inhomogeneous (forced) term in the second harmonic field equation vanishes. Thus, only the pulse associated with the homogenous term is recorded at the output. 5. Group-velocity locking

The group velocity locking was monitored in the case of type 0 ee-e second-harmonic generation, the most efficient one. The CCD camera was removed from the reference position and, by using a spherical mirror, the arm containing the locked pulses was refocused inside a 100 µm thick BBO crystal. The BBO was cut in order to achieve perfect phase-matching for the interaction ω+2ω=3ω.This interaction is possible only if the ω and 2ω pulses arrive simultaneously onto the BBO crystal, and its detection would certainly demonstrate that group-velocity locking was taking place inside the LNB crystal. Indeed, using our setup a generated ultraviolet (3ω) pulse train was monitored by the PMT on the scope screen (Figure 4(I)).

Fig. 4. The generated 3ω signal as recorded by the photomultiplier tube (I), with a thin (II) or a thick (III) glass plate at the reference position of the experimental set-up in fig.2.

To be sure that the 3ω signal was not generated inside the LNB crystal, a 100 µm thick glass plate was positioned at the reference position in order to remove any UV light present there. The PMT was still able to register a pulse train (figure 4(II)). In order to be sure of the pulse temporal overlapping as well as of the absence of any 3ω generated inside the BBO crystal, an 8 mm thick glass plate was positioned at the reference position. The glass plate separates the ω and 2ω pulses temporally, so that the two signals arrive at the BBO trailing each other. In that case no 3ω light was observed (figure 4(III)). 6. Conclusions

In conclusion we have experimentally demonstrated that pulsed second harmonic generation under phase and group velocity mismatch conditions gives rise to phase-locked fundamental and second harmonic pulses that propagate with similar effective dispersive properties and locked group-velocities. The locking occurs with type-0 and type-I polarisation couplings, while it is not present with type-II coupling. The absence of a phase locked signal in the latter case can be attributed to the lack of any sustained pump fields overlap. Thus, in the locked cases we have examined a complete imprinting of refractive index and group velocity occurs between the pump and the second harmonic pulse, with the locked pulse at 400 nm experiencing exactly the refractive index of its pump at 800 nm. Acknowledgments

This work has been partially supported by the National PRIN-05 project “sorgenti e amplificatori laser in guide solitoniche in niobato di litio”; Mathieu Chauvet was partially supported by La Sapienza University with the project “professori visitatori”. A.M.D.G.

#103628 - $15.00 USD

(C) 2009 OSA

