Copyright 2001 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in Appl. Phys. Lett. 78(13), 1861 (2001), and may be found at http://link.aip.org/link/?apl/78/1861. To properly read some greek symbols on this page, you will need to have a Symbol font installed. Piezoelectric Contributions to Pulsed Degenerate FourWave MixingSwiss Federal Institute of Technology, ETHHönggerberg, CH8093 Zürich, Switzerland
(Received 6 November 2000)
In solid state materials lacking a center of symmetry, cascaded second order nonlinear optical effects [1,2,3], especially the combination of optical rectification and the linear electrooptic effect, are an important contribution to Degenerate FourWave Mixing (DFWM) [4,5], an experimental tool widely used to measure the thirdorder nonlinear optical susceptibility. This quantity determines the magnitude of several interesting phenomena [6], among which the optical Kerr effect, selfphase modulation, soliton formation, and the interaction of different optical waves, with applications like optical limiting and alloptical switching. In order to extract an accurate value of the thirdorder susceptibility of noncentrosymmetric materials from a DFWM experiment, it is necessary to take into account the peculiar geometry dependences introduced by second order effects [4,5,7]. They can both lead to misleading results, if ignored, or to a alternative possibility for calibrating DFWM experimental setups, if properly taken into account [5]. Some predicted effects of the electrooptic and optical rectification contributions to DFWM have been observed in Refs. 4,5,7. This letter now reports the experimental observation of another, associated effect: the nanosecond fast piezoelectric relaxation whose contributions to DFWM was calculated in Ref. 5. This is important because the possibility of piezoelectric relaxation introduces additional dependencies of the DFWM signal from the angle between the interacting beams and from the pulselength. This must be taken into account in order to be able to compare experimental results obtained with different DFWM setups. At the origin of the piezoelectric contributions discussed in this work are two spacemodulated rectified polarizations induced by two pairs of interacting waves in DFWM. For the experimental configuration shown schematically in Fig. [1], the rectified polarizations can be described by two static (on the time scale of the pulse length used in DFWM) plane waves with wavevectors k_{a} and k_{b}. [4,5] Approximating the interacting beams and the generated polarizations by plane waves with real amplitude (1/2)A(w, k) exp[i ( t  k x)] +c.c. and complex amplitude A(w, k), one can write the complex amplitude of the rectified polarization with wavevector k_{a} as
where e_{0} is the permittivity of vacuum, E_{n}( w, k) is the component of the electric field amplitude along n, E_{n}( w,  k) is its complex conjugate, and c_{pjl}^{(2)}(0,  w, w) is the second order susceptibility describing optical rectification and the linear electrooptic effect. It is related to the standard electro optic coefficients r_{ijk} by c_{kij}^{(2)}(0,  w, w) =  (1/2) n_{i}^{2} n_{j}^{2} r_{ijk}, where n_{i} is the refractive index for light polarized along the iaxis. An analogous expression gives the rectified polarization with wavevector k_{b}. [5] The electrooptic interaction of one of the input beams with these optical rectification gratings contributes to DFWM, and can be represented by two effective nonlinear optical susceptibility tensors, c^{ ka}_{ijkl} and c^{ kb}_{ijkl}, which must be added to the genuine thirdorder susceptibility c^{(3)}_{ijkl} to get an effective susceptibility tensor c^{(3), EFF}_{ijkl} that describes DFWM in noncentrosymmetric materials. More details about the nomenclature used here for these tensors can be found in Ref. 5. The results of Ref. 5 can be cast in the form:
where r_{ijk} and e_{ij} are effective electrooptic and dielectric tensors to be discussed below, d_{qp}=0, q<>p;d_{qp}=1,q = p, and z_{qp}(k) is
In Eqs. (2)(3), z_{qp} can be substituted with z_{qp}=0 for a transverse polarization (perpendicular to k) and with z_{qp}=d_{qp}/e_{pp} for a longitudinal polarization (parallel to k) oriented along the paxis.
Piezoelectric effects become important when the propagation time of an acoustic elastic wave over a distance equal to the modulation period of the rectified polarization (L = 2 p/k_{a} for Eq. (1) is shorter than the laser pulse length used. In this case a strain pattern can be established in the crystal while the interacting pulses are still present, and it leads to additional contributions to the DFWM signal. These additional contributions can be included in the limit of long wavevectors by using, in Eqs. (2)(4), effective electrooptic and dielectric tensors that take into account that the crystal textemdash\ because of the sinusoidal spatial modulation of the rectified polarizations textemdash\ is neither in a completely strainfree case, nor in a completely stressfree case. These effective tensors can be calculated by [5]
where k is the wavevector of the corresponding optical rectification grating (see Eq. (1)), r^{S}_{ijk} is the electrooptic tensor at constant strain, e^{S}_{ij} is the dielectric tensor at constant strain, k_{i}=k_{i}/k, p_{ijkl} is the elastooptic tensor, e_{ijk} is the piezoelectric tensor, A_{ik} = C^{E}_{ijkl} k_{j} k_{l}, and C^{E}_{ijkl} is the elastic stiffness tensor. For short wavevectors the crystal has no time to elastically relax, and r_{ijk}(k)=r^{S}_{ijk}, e_{ij}(k)=e^{S}_{ij}. This leads to a dependence of the thirdorder susceptibility determined in DFWM on the length of the wavevector of the rectified polarization grating, a fact that influences DFWM measurements in various ways. As an example, in the arrangement of Fig. 1 the wavevector k_{b} is generally quite long, and piezoelectric relaxation can already become possible in some picoseconds. The k_{b} contribution will thus be different for experiments performed, e.g., with 100 ps pulses and for experiments performed with 1 ps pulses. However, to directly observe such an effect, one would need to tune the pulse length, which is impractical. To experimentally demonstrate the existence of the piezoelectric contribution, and to relate its magnitude to the predicted theoretical value, it is better to focus on the rectified polarization grating with the shorter wavevector k_{a}. When the duration of the laser pulse is similar to the acoustic wave propagation time discussed above, the effect of piezoelectric relaxation can be observed by using a nanosecond laser in the experimental configuration of Fig. 2 and measuring the DFWM signal as a function of the angle b between the interacting beams (which is related to the spatial period L = 2 p/k_{a} of the optical rectification grating by L = l/[2 sin (b/2)], with l the laser wavelength in vacuum). Barium Titanate (BaTiO_{3}) is a good candidate to observe these effects because its electrooptic, piezoelectric, and elastic properties are well characterized, and because it has large electrooptic coefficients. All the relevant material parameters have been tabulated in Ref. 8 and the dispersion of the electrooptic coefficients is given in Ref. 9. BaTiO_{3} has tetragonal symmetry with point group mm4, with the caxis parallel to the ferroelectric polarization. The piezoelectric contribution is expected to be particularly strong for c^{ ka}  c_{1133}, i.e. for the contribution of the optical rectification grating with wavevector k_{a} when the caxis of the crystal is parallel to k_{a}. Numerical values are given below. They were all calculated using the material parameters in Ref. 8 and in Ref. 9 at a wavelength of 1.06 mum.
The DFWM experiments were realized with a Qswitched laser producing 7 ns long pulses at a wavelength of 1.064 mum and the experimental setup shown in Fig. 2. This setup has been especially constructed to allow a relatively easy adjustment of the angle between the beams while keeping the path lengths of all beams, and therefore the relative coherence of the interacting pulses, constant. To measure the effective susceptibility c^{(3), EFF}_{1133}, beams 1 and 2 were polarized in the plane of incidence, while beams 3 and 4 were polarized perpendicular to it (by the halfwave plate marked ``HW'' in Fig. 2. For the BaTiO_{3} crystal oriented as in Fig. 2, with its c axis in the beam incidence plane, c^{ k}_{a}  c_{1133} contributes to the signal, and it is proportional to (r_{131})^{2}/(e_{11}+2) (see Eq. (1)). It can be calculated to be c^{ k}_{a}  c_{1133} = (224 ± 45) x 10^{22} m^{2}/V^{2} when the strainfree tensor components are used [5]. When elastic relaxation is allowed, the effective electrooptic and dielectric coefficients of Eqs. (5) and (6) must be used, resulting in c^{ k}_{a}  c_{1133} = (355 ± 70) x 10^{22} m^{2}/V^{2}. The errors are estimated on the basis of a ~10% uncertainty in the dielectric and electrooptic tensors [8]. The inclusion of piezoelectric effects does not lead to larger uncertainties because textemdash\ in this particular case in BaTiO_{3} textemdash\ only a shear deformation can be induced, and the effective tensors calculated from Eqs.(5)(6) correspond to the stressfree values, which can be taken directly from Ref. 9. The magnitude of the total piezoelectric contribution in this case is thus expected to be (130 ± 60) x 10^{22} m^{2}/V^{2} (the error takes into account that this value is calculated from a difference, but also that any errors in the strainfree and stressfree values of the material tensors are not necessarily uncorrelated). The contribution from the rectification grating with wavevector k_{b} is fixed at c^{ k}_{b} perp c_{1133} = (224 ± 45) x 10^{22} m^{2}/V^{2} because there any piezoelectric contribution vanishes by symmetry, as follows from Eq. (5). The experimental value of the effective DFWM thirdorder susceptibility should thus increase by approximately (25 ± 13)% over its (560 ± 80) x 10^{22} m^{2}/V^{2} strainfree value [5] when the wavevector of the rectified polarization increases to a magnitude that allows for elastic relaxation during the pulse duration. The beam intersection angle when this increase takes place is given by the modulation period L = 2 p/k_{a} of the rectified polarization grating for which the propagation time of a transverse acoustic wave is similar to the pulse duration of 7 ns. From the relevant speed of sound of a transverse wave, v=[C^{E}_{1313}/r]^{1/2}=3.18 mu m/ns (r is the mass density), one gets L ~ 20 mm. This corresponds to a beam intersection angle b between 50 and 60 mrad. An important fact that can be exploited for the measurements is that when the BaTiO_{3} crystal is oriented with its c axis perpendicular to the beam incidence plane, c^{ k}_{a} perp c_{3311} vanishes almost completely [5]. For this crystal orientation no dependence of the DFWM signal from the beam intersection angle b is expected. In addition, for both c  k_{a} and c perp k_{a}, c^{ k}_{b} perp c_{1133} is the same and independent of the angle between the interacting beams. It follows that the effects of possible variations of mirror reflectivity and/or beam overlap conditions when changing the beam intersection angle can be avoided by measuring the DFWM signal for both c  k_{a} and c perp k_{a} and taking the ratio between the two values. Because of the expected constancy of the signal for c perp k_{a}, the data measured in such a way should directly reflect the dependence of c^{ k}_{a}  c_{1133} on the magnitude of k_{a}, or the beam intersection angle b.
The results are shown in Fig. 3. The data points are averages over several measurements. The effect of acoustic phonons for longer wavevectors can clearly be seen as the increase of the thirdorder susceptibility when the angle between the beams increases. The transition from a lower to a higher susceptibility value does indeed take place for a period of the rectified polarization grating around 20 mum, as expected from the estimation above. The ratio between the two susceptibility values at short and long wavevectors is higher than the value calculated above, but is within the expected uncertainties. The value of c^{ k}_{a}  c_{1133}/\c^{ k}_{a} perp c_{3311} measured at short wavevectors is somewhat lower than what is expected from the calculated values of the second order contributions [5] but this small discrepancy is also within the experimental error. The piezoelectric effect whose signature is visible in Fig. 3 is peculiar to DFWM in noncentrosymmetric materials. In contrast to the example of thermal gratings, where textemdash\ given short enough pulses textemdash\ acoustic phonons can be induced and detected (as an oscillation of the signal intensity as a function of the delay time of the probe pulse) no matter what the magnitude of the wavevector k_{a} is, in the process observed here acoustic phonons only contribute to the signal if their characteristic frequency (given by the the speed of sound divided by the spatial period of the rectification grating) is higher than the frequency content of the laser pulses. Their effect can only be seen as a dependence of the DFWM signal from the pulse length, or from the angle between the beams. It should be noted that in this experiment the beams inducing the rectified polarization gratings are always polarized perpendicular to each other, so that there is no light intensity modulation with wavevector k_{a} or k_{b}. A generation of acoustic phonons by transient temperature gratings can therefore be excluded, together with any other effects relying on spatially inhomogeneous energy deposition. In conclusion, the measurements presented above give a clear demonstration of the existence of piezoelectric contributions to pulsed DFWM, which must be taken into account in order to determine reliable values of thirdorder susceptibilities in noncentrosymmetric materials, and to be able to compare DFWM experiments performed with different pulse lengths or experimental setups. References
