Hartmann-Shack wavefront sensing for nonlinear materials characterization

Hartmann-Shack wavefront sensing for nonlinear materials characterization D. Rativa, 1, R. E. de Araujo, 2 A. S. L. Gomes, 3 and B. Vohnsen 1 1 School of Physics, University College Dublin, Belfield, Dublin 4, Ireland 2 Department of Electronic and Systems, Federal University of Pernambuco, Recife, Brazil 3 Department of Physics, Federal University of Pernambuco, Recife, Brazil diego-jose.rativa-millan@ucd.ie Abstract: We present two new techniques exploiting a Hartmann-Shack wavefront sensor to characterize the optical self-focusing effect of nonlinear materials. We demonstrate that the defocus Zernike coefficient (C 5 ) can be used to quantify nonlinear optical properties of materials. In the first technique proposed, the wavefront of a collimated laser beam transmitted through a nonlinear sample is analyzed with different irradiance values. In the second technique,instead of conventional detectors, a Hartmann- Shack sensor is used in a Z-scan setup. The methods are demonstrated by measuring the nonlinear refractive indices of CS 2 and Quartz, using femtosecond Ti:sapphire lasers at 76MHz and 1KHz repetition rate. 2009 Optical Society of America OCIS codes: (190.4400) Nonlinear optics, materials; (280.4788) Optical sensing and sensors; (010.7350) Wave-front sensing. References and links 1. R. E. de Araujo and A. S. L. Gomes, Nonlinear optical Kerr coefficients of disordered media, Phys. Rev. A 57, 2037-2040 (1998). 2. H. Ma, A. S. L. Gomes, and C. B. de Araújo, Measurements of nondegenerate optical nonlinearity using a two-color single beam method, Appl. Phys. Lett. 59, 2666-2668 (1991). 3. M. Sheik-Bahae, A. A. Said, T. H. Wei, D. J. Hagan, and E. W. Van Stryland, Sensitive measurements of optical nonlinearities using a single beam, Quantum Electron. 26, 760-769 (1990). 4. T. Xia, D. J. Hagan, M. Sheik-Bahae, and E. W. Van Stryland, Eclipsing Z-scan measurement of λ/10 4 wavefront distortion, Opt. Lett. 19, 317-319 (1994). 5. J. Wang, M. Sheik-Bahae, A. A. Said, D. J. Hagan, and E. W. Van Stryland, Time resolved Z-scan measurements of optical nonlinearities, J. Opt. Soc. Am. B. 11, 1009-1017 (1994). 6. M. Falconieri and G. Salvetti, Simultaneous measurement of pure-optical and thermo-optical nonlinearities induced by high-repetition-rate, femtosecond laser pulses: application to CS2, Appl. Phys. B 69, 133-136 (1999). 7. A. S. L. Gomes, E. L. Filho, C. B. de Araújo, D. Rativa, and R. E. de Araujo, Thermally managed eclipse Z-scan, Opt. Express 15, 1712-1717 (2007). 8. R. E. de Araujo, A. S. L. Gomes, and C. B. de Araújo, Measurements of pk a of organic molecules using third-order nonlinear optics, Chem. Phys. Lett. 330, 347-353 (2000). 9. D. Rativa, R. E. de Araujo, C. B. de Araújo, A. S. L. Gomes, and L. Kassab, Femtosecond nonlinear optical properties of lead-germanium oxide amorphous films, Appl. Phys. Lett. 90, 231906-231909 (2007). 10. D. Rativa, R. E. de Araujo, and A. S. L. Gomes, One photon nonresonant high-order nonlinear optical properties of silver nanoparticles in aqueous solution, Opt. Express 16, 19244-19252 (2008). 11. B. Vohnsen and S. I. Bozhevolnyi, Near- and far-field second-harmonic imaging of quasi-phase-matching crystals, J. Microsc. 202, 244-249 (2001). 12. D. Rativa, A. S. L. Gomes, S. Wachsmann-Hogiu, D. L. Farkas, and R. E. de Araujo, Nonlinear Excitation of Tryptophan Emission Enhanced by Silver Nanoparticles, J. Fluoresc. 18, 1151-1155 (2008). 13. Q. Yang, J. Seo, S. Creekmore, D. Temple, A. Mott, N. Min, K. Yoo, S.Kim, and S.Jung, Distortions in Z-scan spectroscopy, Appl. Phys. Lett. 82, 19-21 (2003). 14. D. Malacara, Optical Shop testing 2nd Ed. (John Wiley & Sons, Inc., New York, 1992). (C) 2009 OSA 23 November 2009 / Vol. 17, No. 24 / OPTICS EXPRESS 22047

15. A. Roorda, F. Romero-Borja, W. Donnelly, H. Queener, T. Hebert, and M. Campbell, Adaptive optics scanning laser ophthalmoscopy, Opt. Express 10, 405-412 (2002). 16. B. Vohnsen, I. Iglesias, and P. Artal, Confocal scanning laser ophthalmoscope with adaptive optical wavefront correction, Proc. SPIE. 4964, 24-32 (2003). 17. J. Bueno, B. Vohnsen, L. Roso, and P. Artal, Temporal wavefront stability of an ultrafast high-power laser beam, Appl. Opt. 48, 770 777, (2009). 18. V. Mahajan, Optical Imaging and Aberrations: Part I. Ray Geometrical Optics (SPIE Press, Washington, 2001). 19. R. Boyd, Nonlinear Optics Third Edition. (Elsevier Science & Technology Books, New York 2008). 20. P. Chen, D. A. Oulianov, I. V. Tomov, and P. M. Rentzepis, Two-dimensional Z scan for arbitrary beam shape and sample thickness, J. Appl. Phys. 85, 7043 7050 (1999). 1. Introduction The field of Nonlinear Optics has provided many techniques to characterize photonic materials. The nonlinear refractive index can be determined by spectroscopic methods exploiting nonlinear interferometry [1], wave mixing [2] and beam distortion measurements [3]. One well established method, known as Z-scan (hereafter called conventional Z-scan), was introduced in 1989 by Sheik-Bahae and coworkers [3]. Since then, several variations of the method have been developed, including improvement of its sensitivity [4], allowing time-resolved measurements [5] and thermal and non-thermal analysis [6, 7]. The Z-scan technique [3] is the most widely used method to measure the nonlinear index of refraction and the nonlinear absorption of organic materials [8], optical materials [9], nanostructured materials [10], etc. The nonlinear characterization is important to describe effects such as second harmonic generation [11] and multiphoton fluorescence [12]. The Z-scan method exploits the wavefront distortion of a focused beam, analyzing the transmittance value of a sample measured through a finite aperture in the far field as the sample is moved along the propagation path (Z)[3]. Theoretically, it has a wavefront sensitivity of better than λ/300 to determine the nonlinear properties of optical materials [3]. However, for laser systems having pulsewidth instabilities, fluctuations of the irradiance or sample imperfections it may lead to a loss in resolution. As reported in Ref. [13] a 5% fluctuation of laser intensity is enough to decrease the wavefront resolution to approximately λ/10. Here we propose two different techniques to characterize the nonlinear refractive index of materials. The methods use a Hartmann-Shack wavefront sensor (HS) with near-ir light (800nm) to determine the nonlinear optical properties of materials in the femtosecond regime. The first technique uses a collimated setup, where the beam transmitted through a sample is analyzed with a HS sensor. The second technique is a new variation of the conventional Z- scan method, where the wavefront sensor measures simultaneously, the wavefront distortion and the transmittance changes along the focalized laser beam (hereafter called Hartmann Shack Z-scan or HS Z-scan). As the beam radius of curvature is not sensitive to intensity fluctuations, measuring the wavefront instead of the transmittance makes the HS Z-scan technique highly sensitive. The two techniques are demonstrated measuring the n 2 value of CS 2 and Quartz, transparent materials that are commonly used as reference materials in nonlinear optics [3]. The usefulness of the techniques is under evaluation with samples having strong (nonlinear) absorption. A HS wavefront sensor is commonly used for adaptive optics systems in astronomy [14], visual optics and retinal imaging [15, 16], optical components testing and laser beam characterization [17]. It consists of an array of microlenses of equal focal length, focusing an incident wavefront onto a CCD array. The local tilt of the wavefront across each lenslet can be calculated from the displacement of each focal spot on the sensor. By sampling an array of lenslets all of these local tilts can be measured and the whole wavefront determined. The wavefront can be expressed in terms of a set of Zernike polynomials as W(x,y) = 0 c nz n (x,y), where the (C) 2009 OSA 23 November 2009 / Vol. 17, No. 24 / OPTICS EXPRESS 22048

wavefront W(x,y) is derived from a summation of Zernike polynomials Z n (x,y) weighted by the Zernike coefficients c n [18]. Different wavefront aberrations such as Tilt, Astigmatism and Defocus can be described using the Zernike functions. Common low-order Zernike polynomials in the Malacara index scheme are shown in Table 1 as a function of polar coordinates (r,θ) [14]. Table 1. Malacara Zernike index scheme Zernike Polynomial Name Zernike Polynomial Name 1 1 Piston 4 r 2 sin(2θ) Astigmatism ±45 2 rsin(θ) Tilt y 5 2r 2 1 Defocus 3 rcos(θ) Tilt x 6 r 2 cos(2θ) Astigmatism 0 /90 If a sufficiently energetic laser pulse propagates through a nonlinear medium, it will induce an intensity-dependent change in the index of refraction given by n(r,t)=n 0 + n 2 I(r,t), where I is the irradiance of the laser beam, and n 0 and n 2 are the linear and nonlinear indices of refraction, respectively [19]. When n 2 > 0orn 2 < 0 the induced wavefront distortion results in a self-focusing or self-divergence effect of the beam. Assuming a Gaussian beam travelling in the z direction we can write the complex electric field exiting the sample as E out (z,r,t) =E in (z,r,t)e αl/2 e iδφ L(z,r,t) where E in is the incident electric field, L is the sample length and α is the linear absorption. The phase change inside of the sample is exclusively due to linear refraction given by Δφ L (z,r,t)=(2π/λ)n i L where n i is the linear refractive index of the sample. For a high intensity beam a new phase change Δφ NL (t) is introduced due to self-phase modulation effects, being directly proportional to n 2 and the (ideally) Gaussian beam intensity I(r,t)=I 0 (t)e [ 2r2 /w 2]. Therefore, the nonlinear phase change is given by [3] : Δφ NL (t)=(2π/λ)n 2 I(r,t)L ef f (1) where L ef f is given by L ef f = ( 1 e αl) /α. Therefore, the total wavefront at the end of the sample is given by W total = W (z=0) +W L(z=L) +W NL(z=L), where W (z=0) contains the wavefront profile including the aberrations of the experimental setup, W L(z=L) the aberration induced by the linear refraction and sample imperfections, and W NL(z=L) the nonlinear phase change accumulated in the sample. The nonlinear phase change Δφ NL (t) can be expressed in terms of Zernike polynomials. 2. System configurations The two methods were performed using an infrared (800nm) femtosecond laser, Coherent Libra (mj, 80fs) delivering optical pulses at 1kHz repetition rate. In the first technique (see Fig. 1), we also used a Coherent Mira (nj, 150 fs) at 76MHz, ideal to create thermo-optical effects inducing nonlinear negative refraction effects in the sample [7]. The wavefront was acquired using a H-S wavefront analyzer (Thorlabs WFS150C) with a 15Hz frame rate (i.e. 66ms). 2.1. Beam collimated configuration As shown schematically in Fig. 1, the laser intensity was controlled with a polarizer and a half-wave plate (λ/2) and the beam diameter was collimated and reduced to 1.1mm using a keplerian telescope (T1) such that the beam Rayleigh length (L R 5m) remained much larger than the sample thickness (L 1cm). At the exit of the sample the laser beam was magnified (C) 2009 OSA 23 November 2009 / Vol. 17, No. 24 / OPTICS EXPRESS 22049

with a second telescope (T2) to obtain a beam waist about 3mm and imaged on a HS wavefront analyzer. The wavefront data was acquired using the Thorlabs software and the National Instruments DataSocket server integrated with the LabView platform. Fig. 1. Experimental set-up for HS sensing. A half-wave plate (λ/2) and a polarizer are used to adjust the laser intensity. T1 and T2 are Keplerian telescopes. 2.2. HS Z-scan configuration The experimental setup for HS Z-scan is shown schematically in Fig. 2. The optical beam was focused using a 15cm focal length lens (L1), and after exiting the sample the beam were captured by the HS sensor. The HS sensor is located at a distance such that the laser beam overfills its sensitive area (diameter=4mm). A computer-controlled motor translation stage is used to scan the sample along the beam while the Zernike coefficients and the CCD array intensities are recorded. The beam waist measured at the sample position was 20μm. When a convergent Gaussian beam passes through a self-focusing absorbing medium, the radius of curvature of the wavefront and the beam intensity changes. Due to the HS sensor characteristics, it is possible to measure simultaneously the sample transmittance on the CCD array and the wavefront distortion. Therefore, just one branch is necessary to measure the wavefront distortions due to the nonlinear self-focusing effect as well as transmittance changes due to the nonlinear absorption during the axial scan [20]. Fig. 2. HS Z-scan experimental set-up. 3. Experimental results and discussion 3.1. Beam collimated configuration A wavefront reference was taken using lower intensities with the sample placed in the collimated beam, W re f = W (z=0) +W L(z=L). In that case the measured wavefront changes are based on this reference to allow accurate measurements of further wavefront distortions W NL(z=L) induced by higher intensities. We have used the initial 6 Zernike modes indicated in Table 1 to fit the measured wavefront within the pupil area (the piston term, not measurable with the HS sensor has been excluded). Figure 3 shows the Zernike coefficients used to describe the wavefront aberrations induced by the CS 2 sample to different intensities (I a = 3.0GW /cm 2, I b = 4.5GW /cm 2, I c = 5.5GW /cm 2 and I d = 6.0GW /cm 2 ). While Zernike coefficients corresponding to tilt and astigmatism (C 2, C 3, C 4, C 6 ) remain approximately constant there is a dependence of the defocus coefficient C 5 (C) 2009 OSA 23 November 2009 / Vol. 17, No. 24 / OPTICS EXPRESS 22050

with the laser intensity. In accordance with the Zernike orthogonality properties [14], the defocus Zernike is thus enough to describe the wavefront distortions caused by nonlinear effects. Fig. 3. Zernike coefficients used to describe the wavefront aberrations induced for different intensities: (a) I a = 3.0GW /cm 2, (b) I b = 4.5GW /cm 2, (c) I c = 5.5GW /cm 2 and (d) I d = 6.0GW /cm 2. Due to the self-focusing effect the initial wavefront reference becomes defocused with a negative or positive radius of curvature according to the sign of n 2. Neglecting the evolution of the beam diameter within the sample, the nonlinear index of refraction can be estimated as n 2 (r,t)= λ 2πI(r,t)L ef f C 5 (2) where C 5 can be described in waves (wavelength) or radians. Figure 4 shows the defocus coefficient C 5 dependence on laser intensity for CS 2 (L=2mm) and Quartz (L=10mm) samples. In Fig. 4(a), the evolution of C 5 for different irradiances is plotted with a positive linear behavior for a low rate repetition laser excitation (1KHz). Figure. 4(b) shows the positive and negative dependence for 1kHz and 76MHz rate repetition lasers in a CS 2 sample. From the measured C 5 coefficient and the laser beam power at the sample, the n 2 value of the given sample could be calculated from Eq. 2. The measured values for CS 2 and Quartz are: n CS 2 2 (2.7 ± 0.2) 10 15 cm 2 /W and n Quartz 2 (2.4 ± 0.2) 10 16 cm 2 /W respectively, in agreement with previously reported values [3, 7]. To demonstrate that the method measures not only the magnitude but also the sign of n 2, a high rate repetition laser was used to induce negative nonlinearities in the CS 2 sample due to thermal effects [7]. The thermal nonlinearity value measured for CS 2 was n CS 2 2(thermal) ( 4.4 ± 0.2) 10 15 cm 2 /W, in accordance with the values obtained under similar conditions using the TMEZ-scan technique [7]. In conclusion, the collimated beam technique is suitable to measure the n 2 value and sign of a transparent material. Though the wavefront sensitivity is an order of magnitude less than reported values using the Z-scan technique [3] larger areas can be evaluated (focalization is not necessary), and sample thickness limitations are negligible. These properties make the method suitable for many applications. (C) 2009 OSA 23 November 2009 / Vol. 17, No. 24 / OPTICS EXPRESS 22051

Fig. 4. Evolution of the defocus Zernike coefficient: (a) induced in a 2mm-thick cuvette CS 2 sample (squares) and a 10mm-thick quartz sample (circles) at 1KHz, (b) in a 2mm CS 2 sample using 1KHz (Non-thermal) and 76MHz (Thermo-optical) pulse repetition. 3.2. HS Z-scan. In the HS Z-scan sensor configuration the wavefront of the transmitted beam with the sample far away from the focal point was used as reference wavefront with low beam intensity and negligible nonlinear refraction. The measurements were based on this reference and allow wavefront distortions at the pre- and post-focal positions due to self-focalization to be determined. Figure 5 shows the results of the HS Z-scan compared to Z-scan performed at 800nm in CS 2 (2mm thick cuvette). Figure 5(a) shows the HS Z-scan curves, obtained for different intensities (I b = 6.8GW /cm 2, I c = 0.5GW /cm 2 and I d = 0.2GW /cm 2 )at a repetition rate of 1KHz. The changes in radius of curvature can be described using the defocus Zernike coefficient C 5. To demonstrate the method, we measured under the same conditions the Z-scan curves shown in Fig. 5(b). The n CS 2 2 is positive, therefore the self-lensing effect prior to focus will tend to increase the beam divergence, causing an increase in the curvature radius (Zernike C 5 ) and a decrease of the transmittance detected by the conventional Z-scan setup. As the scan continues the same self-lensing increase the beam divergence, collimating the beam and decreasing the radius of curvature. As a consequence, the HS Z-scan curve looks inverted compared with the Z-scan curve. The experimental relation between the normalized peak and valley transmittance ΔT p v and the phase distortion Δφ (using Z-scan and HS Z-scan, respectively) is ΔT p v 0.4Δφ, in agreement with the relation obtained using a Gaussian decomposition [3]. The HS Z-scan method has a wavefront distortion resolution of approximately λ/50. Although the HS Z-scan sensitivity is lower than reported by Bahae [3], the fact that it measures the wavefront slope of the entire laser beam instead of the transmittance laser intensity makes it potentially more sensitive. This could be explored with a higher precision wavefront sensor than used in the present work. In turn, a higher resolution with the conventional Z-scan technique depends strongly on the aperture size and would eventually become signal limited. As shown in the inset figure in Fig. 5(a), using a lower irradiance I d = 0.2GW /cm 2, we measured the wavefront change using the HS Z-scan for CS 2. Nevertheless, for intensities lowers than I c = 0.5GW /cm 2, it was not possible to distinguish the transmittance variance from laser fluctuations with the Z-scan technique as shown in Fig. 5(b). Naturally, there are other experimental ways to overcome power fluctuations which can be implemented in a laboratory set-up. (C) 2009 OSA 23 November 2009 / Vol. 17, No. 24 / OPTICS EXPRESS 22052

Fig. 5. (a) HS Z-scan and (b) Z-scan signature for a CS 2 sample (2mm thick cuvette) using different intensities: I a = 120GW /cm 2 ( ), I b = 6.1GW /cm 2 ( ), I c = 0.5GW /cm 2 (+) and I d = 0.2GW /cm 2 ( ). The case of I a has been excluded from plot (a) since the beam diverges beyond the entrance pupil of the HS. 4. Conclusions We have demonstrated that the HS wavefront sensor can be used to characterize nonlinear properties of materials using two techniques. In the first one, no scanning is required, the method is insensitive to misalignments, linear scattering and sample imperfections. A second technique, referred to as HS Z-scan is insensitive to fluctuations during the measurement. The HS based methods proposed here are useful low-cost techniques to characterize nonlinear optical properties of biological samples, films and others materials. Acknowledgments Financial support from Science Foundation Ireland (07/SK/B1239a) and Enterprise Ireland (PC/2008/125) is gratefully acknowledged. The authors would like to thank the Brazilian Agencies CNPq and CAPES for financial support, which was also carried out under auspices of the Photonics National Institute of Science and Technology. (C) 2009 OSA 23 November 2009 / Vol. 17, No. 24 / OPTICS EXPRESS 22053