## Abstract

We observe mixed two- and three-photon absorption in bulk rutile (TiO_{2}) around 800 nm using the open aperture *Z*-scan technique. We fit the data with an extended model that includes multiphoton absorption, beam quality, and ellipticity. The extracted two- and three-photon absorption coefficients are below 1 mm/GW and 2 mm^{3}/GW^{2}, respectively. We observe negligible two-photon absorption for 813-nm light polarized along the extraordinary axis. We measure the nonlinear index of refraction and obtain two-photon nonlinear figures of merit greater than 1.1 at 774 nm and greater than 12 at 813 nm. Similarly, we obtain three-photon figures of merit that allow operational intensities up to 0.57 GW/mm^{2}. We conclude that rutile is a promising material for all-optical switching applications around 800 nm.

© 2012 OSA

## 1. Introduction

Rutile (titanium dioxide) is a wide-bandgap semiconductor, advantageous for many optical applications due to its high refractive index, birefringence, and transparency for wavelengths ≥ 400 nm. In addition, it is naturally abundant, and highly stable. Rutile also has a large Kerr nonlinearity (*n*_{2} = 9×10^{−19} m^{2}/W) [1]. These properties make rutile a promising material for nonlinear all-optical applications such as ultrafast switching, logic [2] and wavelength conversion [3].

Even with a high Kerr nonlinearity, two- and three-photon absorption (2PA and 3PA, respectively) impose limitations for nonlinear optical devices that rely on high intensities [4]. For semiconductors, we expect 2PA (3PA) to occur for photon energies equal to or greater than the half- (third-) bandgap energy (*E _{g}*) [5, 6]. The half- (third-) bandgap energy corresponds to 800 nm (1200 nm) in rutile (

*E*= 3.101 eV) [7]. We expect device properties to be strongly wavelength dependent around 800 nm due to the transition from pure 3PA to mixed two- and three-photon absorption.

_{g}Studies have quantified pure 2PA in bulk rutile for visible wavelengths (532–750 nm) [8–10] and 3PA for near-infrared wavelengths (1054 μm) [11]. In bulk rutile, two-photon absorption was observed via weak luminescence for pump wavelengths near 800 nm, however two-photon absorption coefficients were not reported [12]. Meanwhile, the nonlinear properties of TiO_{2} films [13–20] and TiO_{2} composites [15, 16, 18, 19, 21, 22] have been extensively investigated. Of these thin film studies, those that address pure-phase TiO_{2} near 800 nm report effective two-photon absorption coefficients that vary in both magnitude and sign [14, 15, 17–20]. However, there are no reports on two-photon absorption in the bulk material using similar laser parameters for comparison. In addition, 3PA near 800 nm in bulk rutile has not been reported. Therefore, our understanding of multiphoton absorption in bulk rutile is incomplete and requires further study around 800 nm.

In this paper, we investigate multiphoton absorption near 800 nm in bulk rutile using the open aperture *Z*-scan technique [23]. We develop a theoretical model to fit the experimental data and extract two- and three-photon absorption coefficients for different wavelengths and crystal orientations. Using the extracted multiphoton absorption coefficients, we discuss nonlinear absorption processes near rutile’s half-bandgap and evaluate its applicability as a material for all-optical devices.

## 2. Experimental

We perform open-aperture *Z*-scan measurements using two commercially available 1-mm-thick single-crystal rutile samples grown using the floating zone method [24]. One sample is cleaved along the (001) plane and the second is cleaved along the (100) plane. The (001)-cut allows for measurements with the incident laser beam polarized along the ordinary axis of the crystal (*E$\perp $c*), while the (100)-cut allows for polarization along both the ordinary and extraordinary axes (*E$\left|\text{\hspace{0.17em}}\right|$c*).

Our experimental setup is similar to that described in [23]. We use a chirped-pulse amplified Ti:Sapphire system with pulse durations *τ* ≥ 50 fs, energies up to 3 μJ at a repetition rate of 250 kHz. We measure the pulse duration using frequency-resolved optical gating (FROG). The beam is spatially filtered to obtain a Gaussian shape then focused through a lens with a 0.15-m focal length. A scanning slit detector on a single axis stage measures the caustic (beam width versus *z*-position). The caustic determines the beam waist (*w*_{0}* _{x}*,

*w*

_{0}

*in the*

_{y}*x*and

*y*directions, respectively), quality factor (${M}^{2}$) and ellipticity (

*w*

_{0}

_{x}/w_{0}

*). A reference detector before the lens allows us account for small power fluctuations [25]. We adjust the central wavelengths,*

_{y}*λ*

_{0}, of the fs-pulses by spectrally filtering the seed using a physical aperture within the pulse-stretcher before the amplifier. Without filtering, the spectral width, Δ

*λ*(FWHM), is 30 nm at

*λ*

_{0}= 800 nm with a pulse duration of 50 fs. With filtering, we obtain narrow-band spectra centered at either 774 nm or 813 nm with spectral widths of 4 and 10 nm, respectively. Table 1 presents a summary of the measured laser parameters.

We carry out a set of open-aperture *Z*-scan measurements at the three different central wavelengths for each sample/orientation. Each measurement set consists of five different *Z*-scan curves for different incident powers. Measurements taken at low irradiances (< 0.06 GW/mm^{2}) characterize linear signal caused by sample imperfections. We use this signal to remove the linear response from the four higher-power measurements [23].

The intensity dependent attenuation coefficient is approximately given by *α*(*I*, *λ*) = *α*_{0} *+ α*_{2} *I + α*_{3} *I ^{2}*, where

*I*is the intensity,

*α*

_{0}is the linear attenuation coefficient (

*α*

_{0}≈0 for

*λ*= 800 nm in rutile),

*α*

_{2}is the 2PA coefficient, and

*α*

_{3}is the 3PA coefficient. These coefficients are wavelength dependent. We extract multiphoton absorption coefficients by numerically fitting the measured data to a theoretical model that takes into account the experimental beam and sample parameters (see Appendix A). Our mixed multiphoton absorption model has four degrees of freedom (

*α*

_{2},

*α*

_{3}, waist position, and astigmatism). Using this model, we fit the multiple

*Z*-scan curves taken at different peak irradiances simultaneously to extract a single value of

*α*

_{2}and

*α*

_{3}for a given wavelength and polarization.

## 3. Results and discussion

#### 3.1 Open aperture Z-scan measurements

We show open aperture *Z*-scan measurements of bulk rutile for different incident wavelengths and pulse durations in Figs. 1
through 4
. Figure 1 (left) shows a typical experimental open aperture *Z*-scan curve in rutile for 800-nm, 50-fs pulses with *E$\left|\text{\hspace{0.17em}}\right|$c*. Figure 1 (right) shows the transmittance change (Δ*T* = 1–*T*, at the focus) as a function of the peak on-axis irradiance for the same sample and pulse duration (referred to here as a Δ*T*-plot). Figure 2
shows data collected using *λ*_{0} = 774 nm and an input beam polarized with *E$\perp $c* and *E$\left|\text{\hspace{0.17em}}\right|$c* (Figs. 2, left and right, respectively). Similarly, Figs. 3
and 4 show measurements and fits for *λ*_{0} = 813 nm and *λ*_{0} = 800, respectively.

The primary sources of uncertainty in our *Z*-scan measurements are the noise of the laser and sample imperfections. We measure laser noise to be ±1–2%. Using low intensity and by adjusting the sample position, we reduce the linear background variations caused by sample imperfections to less than 1%. To reduce the overall noise, we use a reference photodiode, drift detection, temporal averaging, and background subtraction. These methods suppress the normalized transmittance variation to ±0.0015. Consecutive scans resulted in measurements that were within this transmittance variation. To characterize the repeatability of our technique, we analyzed four measurement sets using *E$\perp $*c, taken over non-consecutive days for both the (001)- and (100)-cut. We could not compare these measurements on a point-for-point basis, consequently, we analyzed these measurement sets using our mixed two- and three-photon model (see Appendix A) and determined the standard deviation of *α*_{2} and *α*_{3}. We obtain a standard deviation of 5% for *α*_{2} and 18% for *α*_{3}. In addition, we have tested our system using CS_{2} with a 50-fs pulse at 800 nm and we obtain a 2PA coefficient of 0.06 mm/GW using standard analysis. This value falls within the range of values obtained for CS_{2} in previous studies using 110-fs pulses [26, 27] and thermal methods [28, 29].

Sample damage [30] and thermal effects [28, 29] can mask electronic nonlinearities. We verified that no damage or alterations occurred to the sample by comparing low irradiance scans taken directly before and after each high intensity measurement set. We performed thermally managed *Z*-scan [28, 29] using our bulk rutile samples in both open- and closed-aperture configurations. By utilizing an optical chopper and fast photodiodes to resolve pulse-to-pulse transmittance changes, we did not observe cumulative (i.e. thermal) changes in the signal. Similarly, measurements taken at 10 kHz produced equivalent measurements to those taken at 250 kHz, confirming that no cumulative effects were present.

#### 3.2 Justification and fitting of the mixed-multiphoton absorption model

To establish the optimum fitting procedure, we first fit our data with standard open-aperture *Z*-scan models and compare to our mixed-multiphoton absorption model. Figure 1 (left) presents a typical experimental *Z*-scan curve that is fit using a pure 2PA-model [23], a pure 3PA-model [31] and our mixed-multiphoton model (dashed, dotted, and solid lines, respectively). In Fig. 1 (right), we show an experimental Δ*T*-plot which is fit using a pure 2PA-, 3PA- and our mixed-multiphoton model (dashed, dotted, and solid lines, respectively).

Comparing a pure 2PA to a pure 3PA model, we find that fitting individual *Z*-scan curves appears to favor the 2PA model. However, a 2PA model predicts a very different Δ*T* behavior with increasing irradiance. Three-photon absorption better fits the high intensity Δ*T* data, yet is insufficient for low intensity data and individual *Z*-scan curves. We see that our mixed-multiphoton model provides an excellent fit to both *Z*-scan curves and Δ*T*-plots. These observations highlight two important points. First, neither a pure 2PA- nor a pure 3PA-model can describe the data sufficiently, necessitating a mixed-multiphoton absorption model. Second, both curve fitting and the Δ*T*-plot method provide complimentary information. Curve fitting is advantageous because signal away from the focus always contains the lowest order absorption, such as 2PA. However, curve fitting can mask higher order absorption and is susceptible to beam distortions. In contrast, the Δ*T*-plot provides a fast measure on higher order absorption processes. Yet, this method may overlook lower order absorption if sufficient low irradiance data is not included.

Quantifying mixed two- and three-photon absorption requires a balance between the curve fitting and the Δ*T*-plot methods. It also requires further constraint to account for the additional fitting coefficients. Consequently, we use a more robust method whereby we fit an entire measurement set, consisting of several full *Z*-scan curves taken at different peak irradiances, using a single set of parameters. This technique produces equivalent results to the Δ*T*-plot method within the uncertainty.

We tested the optimization algorithm by using various starting points while keeping all other parameters at experimentally determined values. These variations included starting both *α*_{2} and *α*_{3} at zero as well as other magnitudes and signs of *α*_{2}, *α*_{3} and the astigmatism. Other combinations, such as two- and four-photon absorption (with astigmatism) could produce reasonable fits in some, but not all cases (notably the 813-nm data set). We found that including two- and three-photon absorption in our model resulted in good fits for all experimental data while minimizing the number of fitted parameters.

We quantitatively compare the fit quality using standard 2PA [23] and 3PA [31] models to our mixed-multiphoton model. Fitting a full data set (for example, a 50 fs pulse, *E$\left|\text{\hspace{0.17em}}\right|$c*) using a 2PA or 3PA model produces *R*^{2} values of 0.974 and 0.925 for ± 4 Rayleigh distance (respectively). We fit the same data set with our mixed-multiphoton model and observe an *R*^{2} of 0.991 (solid lines, Fig. 4 right). We performed a similar analysis with the 813-nm data set (*E$\left|\text{\hspace{0.17em}}\right|$c*) and obtain *R*^{2} values of 0.841, 0.944 and 0.950 (2PA, 3PA and mixed-multiphoton models, respectively). From this analysis, the *R*^{2} correlation shows that the mixed-multiphoton model fits the data better than either a pure two- or three-photon absorption model in all cases. Therefore, we conclude that our mixed-multiphoton model is adequate for fitting our experimental data and use this model to fit the *Z*-scan data in Figs. 2–4 (shown in solid lines).

Two additional fitting parameters were included in the two- and three-photon absorption fitting method. We fit the beam waist location, which aligns the model to the experimental waist position. We found that fitting the astigmatism resulted in a significant improvement in fit quality that we could not achieve by other means (i.e. using additional multiphoton absorption terms or fitting the waist). However, the fit astigmatism is relatively large, ranging from 25% to 100% of the measured Rayleigh distance. We find the largest fit astigmatism using 50-fs pulses. Meanwhile, the measured astigmatism is only up to 12% of the Rayleigh distance. This discrepancy represents the largest source of uncertainty within our analysis. By including the effects of fitting with and without the astigmatism, we estimate the uncertainty to be ± 15% in *α*_{2} and a factor of 2 for *α*_{3}.

#### 3.3 Two- and three-photon absorption coefficients

Table 2
summarizes bulk rutile two- and three-photon absorption coefficients. Overall, *α*_{2} ranges from less than 10^{−7} mm/GW to 0.9 mm/GW ( ± 15%). Comparing 774 nm to 813 nm measurements, we find that the 2PA coefficients at these two wavelengths differ considerably. At 774 nm, *α*_{2} is 0.54 mm/GW (*E$\perp $c*) and 0.89 mm/GW (*E$\left|\text{\hspace{0.17em}}\right|$c*). At 813 nm and *E$\left|\text{\hspace{0.17em}}\right|$c*, *α*_{2} is below the measureable limit. For *E$\perp $c*, we extract a value of 0.08 mm/GW. Lastly, two-photon absorption coefficients for 800 nm are between the 774 nm and 813 nm coefficients. Rutile has typical values of *α*_{3} on the order of 1 mm^{3}/GW^{2} as shown in Table 2. The values of *α*_{3} at 800 nm are 0.2 mm^{3}/GW^{2}. We observe 4 to 9 times higher three-photon absorption coefficients for 774 and 813 nm wavelengths. The difference in *α*_{3} between 774-nm and 813-nm measurements is less than a factor of two and is therefore within the experimental uncertainty (considering both orientations).

These values of *α*_{2} and *α*_{3} in bulk rutile are consistent with values obtained under similar conditions [11]. Our degenerate 2PA coefficients around 774 nm are in good agreement with other nondegenerate two-photon absorption measurements (using a 1054-nm pump and a 612-nm probe producing an effective *α*_{2} coefficient of 1–5 mm/GW) [11]. Similarly, our values of *α*_{3} (0.2–2 mm^{3}/GW^{2}) are within the range of measurements taken at 1054 nm (0.5 mm^{3}/GW^{2}) in the same study.

Recently, both negative and positive 2PA coefficients of significantly higher magnitude have been measured in TiO_{2} thin films in the 750–830 nm range [14, 15, 17–19]. However, the thin films investigated vary widely in growth methods, composition, and crystalline phase, thus making it difficult to compare directly to the bulk measurements reported here. Furthermore, while our values are consistent with other bulk measurements, there are inconsistencies with reported thin film values in both sign and magnitude for polycrystalline rutile [14] and anatase TiO_{2} [14, 15, 18, 19]. Thin films are more challenging to measure than bulk samples, and often require intensities approaching the damage threshold to achieve sufficient signal-to-noise, which may lead to linear effects that appear as nonlinear signal [30]. Therefore, further investigation of TiO_{2} thin film versus bulk nonlinearities is required.

We find that the 2PA coefficients decrease significantly between 774 and 813 nm. These wavelengths correspond to two-photon energies of 2*ћω* = 3.20 eV and 3.05 eV, respectively. These two-photon energies lie on opposite sides of rutile’s indirect bandgap (3.101 eV) [7]. Although we observe no effective 2PA for *E$\left|\text{\hspace{0.17em}}\right|$c* at 813 nm, we still measure a small 2PA-signal for *E$\perp $c*. Other studies attribute weak absorption at 3.062 eV, which is only present for *E$\perp $c*, to a direct forbidden bandgap [7, 12, 32]. This energy is within the two-photon bandwidth of the 813-nm pulse and is likely the cause of the small 2PA-signal observed. Between these two effects, the resonance at rutile’s indirect bandgap is the largest source of 2PA. This analysis implies that using wavelengths further red-shifted from the two-photon resonance (longer than 800 nm) should result in reduced 2PA.

The measurements at 800 nm using a 50-fs pulse require additional interpretation. We find the 2PA coefficients using 800-nm light lie between the values obtained at 774 nm and 813 nm. As shown in Table 2, as the wavelength increases, *α*_{2} decreases monotonically as it passes through the two-photon resonance at the band-edge. This behavior at 800 nm may be partially due to the pulse’s extended spectrum (785–815 nm, FWHM). This bandwidth includes spectral components experiencing both higher (774 nm, Δ*λ* = 4 nm) and lower (813 nm, Δ*λ* = 10 nm) 2PA. This implies that the extracted values of *α*_{2} are a convolution between the pulse spectrum and the two-photon absorption spectrum.

We also observe that the three-photon absorption coefficients measured using a 50-fs pulse are 4–10 times smaller than for either 774 or 813 nm. This reduction is unexpected given that the 3PA coefficients do not change significantly between these wavelengths. This reduction is likely due to the model not taking dispersive and nonlinear pulse broadening into account. In the normal group velocity dispersion-regime, with rutile’s positive *n*_{2}, a pulse broadens temporally during propagation. This broadening decreases the peak intensity and leads to reduced multiphoton absorption. Dispersive broadening is significant if the dispersion length is shorter than or comparable to the sample length. Using the group velocity dispersion for rutile’s extraordinary axis (*β*_{2} = 1250 fs^{2}/mm) and a pulse duration of *τ* = 50 fs, we calculate a dispersion length of *τ*^{2}/*β*_{2} = 2 mm, reducing the peak intensity by 40% over the sample length. In contrast, pulse broadening is negligible for a 174-fs (290-fs) pulse where the dispersion length is 2.4 cm (6.7 cm), resulting in a peak reduction of less than 1%. This intensity reduction has a more pronounced effect on 3PA than on 2PA owing to the *I*^{3}- versus *I*^{2}-dependence. Taking dispersion into account, the uncertainty in the 800-nm measurements increases to ± 80% in *α*_{2} and a factor of 6 for *α*_{3}. Because of the large uncertainty introduced by both dispersive and spectral convolution effects at shorter pulse durations, it is difficult to draw strong conclusions on the time-dependence of nonlinear processes in rutile using the present analysis.

#### 3.4 Nonlinear figures of merit

Multiphoton absorption has important implications for nonlinear optical devices. Material performance is quantified by the nonlinear figures of merit [4]. The 2PA figure of merit (FOM) is given by *n*_{2}/(*α*_{2}*λ*) and the 3PA FOM is given by *n*_{2}/(*α*_{3}*λI*_{max}), where *I*_{max} is the peak operational intensity. Both FOMs should be greater than unity for all-optical applications. By setting the 3PA FOM equal to one, we calculate the intensity at which 3PA becomes a limit, given by ${I}_{\mathrm{max}}^{3PA}$ = *n*_{2}/(*α*_{3}*λ*). We measured the nonlinear index of refraction around 800-nm (50-fs pulse) using closed-aperture *Z*-scan and obtain a value of 7.9 × 10^{−19} m^{2}/W for both polarizations using standard analysis methods [33]. Although theory predicts a resonant enhancement of the nonlinear index near the half-bandgap [5], we did not measure a significant deviation from values reported at 1064 nm [1], which may be due to dispersive broadening as explained in the previous section. Therefore, the calculated nonlinear figures of merit should be considered a lower limit. Using the nonlinear index and multiphoton absorption coefficients, we evaluate rutile as a candidate material for all-optical applications around 800 nm by calculating its nonlinear figures of merit.

Table 2 shows the calculated 2PA FOMs and ${I}_{\mathrm{max}}^{3PA}$ for each measurement. The 2PA FOMs at 774 nm are the lowest observed (1.9 and 1.1 for *E$\perp $c* and *E$\left|\text{\hspace{0.17em}}\right|$c*, respectively). At 813 nm, we obtain a FOM of 12.1 for *E$\perp $c* and greater than 10^{6} for *E$\left|\text{\hspace{0.17em}}\right|$c*. Calculated ${I}_{\mathrm{max}}^{3PA}$ is lowest for *E$\left|\text{\hspace{0.17em}}\right|$c* using 774 nm, with a value of 0.57 GW/mm^{2}, and greater than 0.93 GW/mm^{2} for all other measurements.

All 2PA figures of merit calculated are greater than 1.1, which implies that rutile is compatible with all-optical applications. Measurements at 813 nm show extremely large FOMs, which is favorable for devices. We observe that 2PA can be neglected at 813 nm (*E$\left|\text{\hspace{0.17em}}\right|$c*) and expect similar results for *E$\perp $c* at longer wavelengths [34]. These FOMs imply compatibility with popular femtosecond pulsed sources such as Ti:Sapphire and fiber lasers.

Although 2PA may be negligible for certain wavelengths, we must still avoid significant 3PA by operating below${I}_{\mathrm{max}}^{3PA}$. For the largest 3PA measured, ${I}_{\mathrm{max}}^{3PA}$ is 0.57 GW/mm^{2} (*E$\left|\text{\hspace{0.17em}}\right|$c*, 774 nm). This limit is far below thin-film damage thresholds for TiO_{2} [35, 36] and provides a reasonable upper limit for integrated photonic devices [37–40].

## 4. Conclusion

We observe mixed two- and three-photon absorption in bulk rutile around 800 nm. By fitting open aperture *Z*-scan data to a numeric model, we separate the effects of mixed multiphoton absorption. Two- and three-photon absorption coefficients are below 1 mm/GW and 2 mm^{3}/GW^{2}, respectively. Two-photon absorption falls below our measurement capabilities for 174-fs pulse with a central wavelength of 813 nm (*E$\left|\text{\hspace{0.17em}}\right|$c*). We determine 2PA FOMs greater than 1 above and below rutile’s half-bandgap. We expect superior performance for *λ*_{0} ≥ 813 nm, where the 2PA FOMs are greater than 12.1. Additionally, 3PA is sufficiently small to enable all-optical applications for intensities below 0.57 GW/mm^{2}. From the low two- and three-photon absorption, we conclude that rutile TiO_{2} is a promising material for all-optical applications such as switching and logic near 800 nm.

## Appendix A: theoretical model

Several reports of modeling multiphoton absorption using the *Z*-scan method exist in the literature [23, 31, 41–44]. Models including single nonlinear absorption mechanisms such as pure 2PA [23], higher-order nonlinear absorption [31], and mixed multiphoton absorption [41–43] have been presented. In contrast to single absorption processes, mixed multiphoton absorption is only occasionally observed [45] and measured [46–48]. Non-idealities, such as ellipticity [44] and laser beam quality [49, 50], can obscure interpretation of mixed multiphoton data by changing the shape of the *Z*-scan curve. For example, beam ellipticity reduces the total signal and can cause asymmetry [44]. Meanwhile, the beam quality factor, ${M}^{2}$, reduces the apparent Rayleigh distance by $1/{M}^{2}$ while maintaining the waist [49–51]. Mixed multiphoton absorption is particularly sensitive to such shape changes [41–43]. These effects are not included in the previous models. Thus, we develop a mixed multiphoton absorption model that includes these non-idealities here.

We apply an open aperture *Z*-scan derivation using an elliptical beam, similar to that presented in Ref [44]. and highlight the differences. We start with Equation 10 for an elliptical Gaussian beam from reference [44]:

Here, *z* is the sample position, *x* and *y* are the transverse spatial coordinates, *P*(*t*) is the instantaneous power as a function of time, and *w _{x}*(

*z*) and

*w*(

_{y}*z*) are the beam widths as a function of the sample position,

*z*. We include the beam quality factor through the width as follows:

In this equation, *w _{x}*

_{0}is the beam waist,

*M*

_{x}_{2}is the beam quality factor [51],

*z*

_{x}_{0}is the location of the waist and

*λ*

_{0}is the center wavelength. We use a similar expression in the

*y*-direction. We define the ellipticity by

*w*

_{x}_{0}/

*w*

_{y}_{0}and the astigmatism by

*z*

_{x}_{0}–

*z*

_{y}_{0}. We include 3PA in the attenuation coefficient, defined for a single wavelength by

*α*(I) =

*α*

_{0}+

*α*

_{2}

*I*+

*α*

_{3}

*I*

_{2}. Here

*I*is the intensity, and

*α*

_{0},

*α*

_{2}, and

*α*

_{3}are the one-, two-, and three-photon absorption coefficients, respectively. Rutile is highly transparent for 800 nm, so we set

*α*

_{0}to zero. We solve the following ordinary different equation using standard numerical methods along the sample length,

*L*:

Here, *z´* is the position inside the sample. We solve for each point in *x* and *y* at a fixed sample position (*z*). Then, we integrate over *x* and *y* at the entrance and exit of the sample to calculate the transmittance (power out normalized to power in). It should be noted that by following conventional theory [23], we do not include the effects of self-phase modulation, self-focusing (e.g. filamentation), or dispersion. We model pulsed transmittance using a Gaussian shape, given by *P*(*t*), integrated over time using:

In this equation, *T*(*z*) is the observed energy transmittance, *T* [*z*, *P*(*t*)] is the instantaneous power transmittance, and *E*_{in} and *E*_{out} are the input and output energies, respectively. We take into account the reflection from the front surface. Numerical methods require that we discretize our parameters including the *x*- and *y*-position and time. Consequently, we test for the convergence of all discretization parameters. Additionally, we find excellent agreement between this model and analytical solutions for two- [23] and three-photon absorption [31].

Unless otherwise stated, we fit each data set (multiple *Z*-scan traces at different powers) using this model to a single set of fit parameters, keeping all others fixed to experimental values. We fit our data using the Levenberg-Marquardt method. Our fit parameters are *α*_{2}, *α*_{3}, the average *z*-position of the waist (*z*_{0}) and the astigmatism. From this fit, we determine a single value for *α*_{2} and *α*_{3} that describes the *Z*-scan traces for all irradiances measured.

## Acknowledgments

Several people contributed to the work described in this paper. C.E. and E.M-P. conceived of the basic idea for this work. C.E., E.M-P., and J.B. planned the experiments. C.E. designed and carried out the experiments, implemented the model/fitting, and analyzed the results. E.M., E.M-P. and J.B. supervised the research and the development of the manuscript. C.E. wrote the first draft of the manuscript; all authors subsequently took part in the revision process and approved the final copy of the manuscript. Benjamin Franta, Phillip Munoz and Orad Reshef provided feedback on the manuscript throughout its development. The research described in this paper was supported by the NSF under contract ECCS-0901469. One of the authors, EM-P, thanks CONACYT Mexico for the support through the project 120346, as well as the David Rockefeller Center for Latin American Studies and the Mazur Group for additional financial support.

## References and links

**1. **R. Adair, L. L. Chase, and S. A. Payne, “Nonlinear refractive index of optical crystals,” Phys. Rev. B Condens. Matter **39**(5), 3337–3350 (1989). [CrossRef]

**2. **M. Jinno and T. Matsumoto, “Nonlinear Sagnac interferometer switch and its applications,” IEEE J. Quantum Electron. **28**(4), 875–882 (1992). [CrossRef]

**3. **M. A. Foster, A. C. Turner, R. Salem, M. Lipson, and A. L. Gaeta, “Broad-band continuous-wave parametric wavelength conversion in silicon nanowaveguides,” Opt. Express **15**(20), 12949–12958 (2007). [CrossRef]

**4. **G. I. Stegeman and W. E. Torruellas, “Nonlinear materials for information processing and communications,” Philos. Trans. R. Soc. Lond. A **354**(1708), 745–756 (1996). [CrossRef]

**5. **M. Sheik-Bahae, D. J. Hagan, and E. W. Van Stryland, “Dispersion and band-gap scaling of the electronic Kerr effect in solids associated with two-photon absorption,” Phys. Rev. Lett. **65**(1), 96–99 (1990). [CrossRef]

**6. **E. W. Van Stryland, M. A. Woodall, H. Vanherzeele, and M. J. Soileau, “Energy band-gap dependence of two-photon absorption,” Opt. Lett. **10**(10), 490–492 (1985). [CrossRef]

**7. **J. Pascual, J. Camassel, and H. Mathieu, “Fine structure in the intrinsic absorption edge of TiO_{2},” Phys. Rev. B **18**(10), 5606–5614 (1978). [CrossRef]

**8. **Y. Watanabe, M. Ohnishi, and T. Tsuchiya, “Measurement of nonlinear absorption and refraction in titanium dioxide single crystal by using a phase distortion method,” Appl. Phys. Lett. **66**(25), 3431–3432 (1995). [CrossRef]

**9. **A. Penzkofer and W. Falkenstein, “Direct determination of the intensity of picosecond light pulses by two-photon absorption,” Opt. Commun. **17**(1), 1–5 (1976). [CrossRef]

**10. **H. S. Waff and K. Park, “Structure in the two-photon absorption spectrum of TiO_{2} (Rutile),” Phys. Lett. A **32**(2), 109–110 (1970). [CrossRef]

**11. **P. Sathy and A. Penzkofer, “Three-photon absorption and its limitation of third-order nonlinear optical effects in rutile,” Appl. Phys. B **61**(2), 127–134 (1995). [CrossRef]

**12. **K. Watanabe, K. Inoue, and F. Minami, “Resonant phenomena of hyper-Raman-scattering of optic phonons in a TiO_{2} crystal,” Phys. Rev. B Condens. Matter **46**(4), 2024–2033 (1992). [CrossRef]

**13. **T. Hashimoto, T. Yoko, and S. Sakka, “Sol-gel preparation and third-order nonlinear optical properties of TiO_{2} thin films,” Bull. Chem. Soc. Jpn. **67**(3), 653–660 (1994). [CrossRef]

**14. **H. Long, A. Chen, G. Yang, Y. Li, and P. Lu, “Third-order optical nonlinearities in anatase and rutile TiO_{2} thin films,” Thin Solid Films **517**(19), 5601–5604 (2009). [CrossRef]

**15. **D. Torres-Torres, M. Trejo-Valdez, L. Castañeda, C. Torres-Torres, L. Tamayo-Rivera, R. C. Fernández-Hernández, J. A. Reyes-Esqueda, J. Muñoz-Saldaña, R. Rangel-Rojo, and A. Oliver, “Inhibition of the two-photon absorption response exhibited by a bilayer TiO_{2} film with embedded Au nanoparticles,” Opt. Express **18**(16), 16406–16417 (2010). [CrossRef]

**16. **M. Trejo-Valdez, R. Torres-Martínez, N. Peréa-López, P. Santiago-Jacinto, and C. Torres-Torres, “Contribution of the two-photon absorption to the third order nonlinearity of au nanoparticles embedded in TiO_{2} films and in ethanol suspension,” J. Phys. Chem. C **114**(22), 10108–10113 (2010). [CrossRef]

**17. **H. Long, G. Yang, A. Chen, Y. Li, and P. Lu, “Femtosecond Z-scan measurement of third-order optical nonlinearities in anatase TiO_{2} thin films,” Opt. Commun. **282**(9), 1815–1818 (2009). [CrossRef]

**18. **M. Kyoung and M. Lee, “Z-scan studies on the third-order optical nonlinearity of Au nanoparticles embedded in TiO_{2},” Bull. Korean Chem. Soc. **21**, 26–28 (2000).

**19. **H. I. Elim, W. Ji, A. H. Yuwono, J. M. Xue, and J. Wang, “Ultrafast optical nonlinearity in poly(methylmethacrylate)-TiO_{2} nanocomposites,” Appl. Phys. Lett. **82**(16), 2691–2693 (2003). [CrossRef]

**20. **L. Irimpan, B. Krishnan, V. P. N. Nampoori, and P. Radhakrishnan, “Luminescence tuning and enhanced nonlinear optical properties of nanocomposites of ZnO-TiO_{2.},” J. Colloid Interface Sci. **324**(1-2), 99–104 (2008). [CrossRef]

**21. **A. H. Yuwono, B. Liu, J. Xue, J. Wang, H. I. Elim, W. Ji, Y. Li, and T. J. White, “Controlling the crystallinity and nonlinear optical properties of transparent TiO_{2}-PMMA nanohybrids,” J. Mater. Chem. **14**(20), 2978–2987 (2004). [CrossRef]

**22. **Q. F. Zhou, Q. Q. Zhang, J. X. Zhang, L. Y. Zhang, and X. Yao, “Preparation and optical properties of TiO_{2} nanocrystalline particles dispersed in SiO_{2} nano-composites,” Mater. Lett. **31**(1-2), 39–42 (1997). [CrossRef]

**23. **M. Sheik-Bahae, A. A. Said, T. H. Wei, D. J. Hagan, and E. W. Van Stryland, “Sensitive measurement of optical nonlinearities using a single beam,” IEEE J. Quantum Electron. **26**(4), 760–769 (1990). [CrossRef]

**24. **M. Higuchi, T. Hosokawa, and S. Kimura, “Growth of rutile single crystals by floating zone method,” J. Cryst. Growth **112**(2-3), 354–358 (1991). [CrossRef]

**25. **E. W. V. Stryland and M. Sheik-Bahae, “Z-Scan measurements of optical nonlinearities,” in *Characterization techniques and tabulations for organic nonlinear optical materials*, M. G. Kuzyk and C. W. Dirk, eds. (Marcel Dekker, 1998).

**26. **S. Couris, M. Renard, O. Faucher, B. Lavorel, R. Chaux, E. Koudoumas, and X. Michaut, “An experimental investigation of the nonlinear refractive index (n_{2}) of carbon disulfide and toluene by spectral shearing interferometry and z-scan techniques,” Chem. Phys. Lett. **369**(3-4), 318–324 (2003). [CrossRef]

**27. **R. A. Ganeev, A. I. Ryasnyansky, M. Baba, M. Suzuki, N. Ishizawa, M. Turu, S. Sakakibara, and H. Kuroda, “Nonlinear refraction in CS_{2},” Appl. Phys. B **78**(3-4), 433–438 (2004). [CrossRef]

**28. **M. Falconieri and G. Salvetti, “Simultaneous measurement of pure-optical and thermo-optical nonlinearities induced by high-repetition-rate, femtosecond laser pulses: application to CS_{2},” Appl. Phys. B **69**(2), 133–136 (1999). [CrossRef]

**29. **A. Gnoli, L. Razzari, and M. Righini, “Z-scan measurements using high repetition rate lasers: how to manage thermal effects,” Opt. Express **13**(20), 7976–7981 (2005). [CrossRef]

**30. **B. M. Patterson, W. R. White, T. A. Robbins, and R. J. Knize, “Linear optical effects in z-scan measurements of thin films,” Appl. Opt. **37**(10), 1854–1857 (1998). [CrossRef]

**31. **D. S. Corrêa, L. De Boni, L. Misoguti, I. Cohanoschi, F. E. Hernandez, and C. R. Mendonça, “Z-scan theoretical analysis for three-, four- and five-photon absorption,” Opt. Commun. **277**(2), 440–445 (2007). [CrossRef]

**32. **K. Watanabe and K. Inoue, “Two-photon resonant effect of hyper-Raman scattering in the vicinity of the direct forbidden gap in a rutile crystal,” Phys. Rev. B Condens. Matter **41**(11), 7957–7960 (1990). [CrossRef]

**33. **C. C. Evans, J. D. B. Bradley, F. Parsy, K. C. Phillips, R. Senaratne, E. A. Martí-Panameño, and E. Mazur, “Thermally managed Z-scan measurements of titanium dioxide thin films,” presented at Photonics West, San Francisco, CA, USA, 27 Jan. 2011.

**34. **B. Imangholi, M. P. Hasselbeck, and M. Sheik-Bahae, “Absorption spectra of wide-gap semiconductors in their transparency region,” Opt. Commun. **227**(4-6), 337–341 (2003). [CrossRef]

**35. **J. Yao, Z. Fan, Y. Jin, Y. Zhao, H. He, and J. Shao, “Investigation of damage threshold to TiO_{2} coatings at different laser wavelength and pulse duration,” Thin Solid Films **516**(6), 1237–1241 (2008). [CrossRef]

**36. **M. Mero, J. Liu, W. Rudolph, D. Ristau, and K. Starke, “Scaling laws of femtosecond laser pulse induced breakdown in oxide films,” Phys. Rev. B **71**(11), 115109 (2005). [CrossRef]

**37. **V. G. Ta’eed, N. J. Baker, L. B. Fu, K. Finsterbusch, M. R. E. Lamont, D. J. Moss, H. C. Nguyen, B. J. Eggleton, D. Y. Choi, S. Madden, and B. Luther-Davies, “Ultrafast all-optical chalcogenide glass photonic circuits,” Opt. Express **15**(15), 9205–9221 (2007). [CrossRef]

**38. **S. J. Madden, D. Y. Choi, D. A. Bulla, A. V. Rode, B. Luther-Davies, V. G. Ta’eed, M. D. Pelusi, and B. J. Eggleton, “Long, low loss etched As(_{2)}S(_{3}) chalcogenide waveguides for all-optical signal regeneration,” Opt. Express **15**(22), 14414–14421 (2007). [CrossRef]

**39. **I. D. Rukhlenko, M. Premaratne, and G. P. Agrawal, “Analytical study of optical bistability in silicon ring resonators,” Opt. Lett. **35**(1), 55–57 (2010). [CrossRef]

**40. **P. Koonath, D. R. Solli, and B. Jalali, “Limiting nature of continuum generation in silicon,” Appl. Phys. Lett. **93**(9), 091114 (2008). [CrossRef]

**41. **B. Gu, J. Wang, J. Chen, Y.-X. Fan, J. Ding, and H.-T. Wang, “Z-scan theory for material with two- and three-photon absorption,” Opt. Express **13**(23), 9230–9234 (2005). [CrossRef]

**42. **B. Gu, X. Q. Huang, S. Q. Tan, M. Wang, and W. Ji, “Z-scan analytical theories for characterizing multiphoton absorbers,” Appl. Phys. B **95**(2), 375–381 (2009). [CrossRef]

**43. **J. Wang, B. Gu, X.-W. Ni, and H.-T. Wang, “Z-scan theory with simultaneous two- and three-photon absorption saturation,” Opt. Laser Technol. **44**(2), 390–393 (2012). [CrossRef]

**44. **S. M. Mian, B. Taheri, and J. P. Wicksted, “Effects of beam ellipticity on Z-scan measurements,” J. Opt. Soc. Am. B **13**(5), 856–863 (1996). [CrossRef]

**45. **M. G. Vivas, T. Shih, T. Voss, E. Mazur, and C. R. Mendonca, “Nonlinear spectra of ZnO: reverse saturable, two- and three-photon absorption,” Opt. Express **18**(9), 9628–9633 (2010). [CrossRef]

**46. **G. Boudebs, S. Cherukulappurath, M. Guignard, J. Troles, F. Smektala, and F. Sanchez, “Experimental observation of higher order nonlinear absorption in tellurium based chalcogenide glasses,” Opt. Commun. **232**(1-6), 417–423 (2004). [CrossRef]

**47. **M. Chattopadhyay, P. Kumbhakar, C. S. Tiwary, A. K. Mitra, U. Chatterjee, and T. Kobayashi, “Three-photon-induced four-photon absorption and nonlinear refraction in ZnO quantum dots,” Opt. Lett. **34**(23), 3644–3646 (2009). [CrossRef]

**48. **T. Sarma, P. K. Panda, P. T. Anusha, and S. V. Rao, “Dinaphthoporphycenes: synthesis and nonlinear optical studies,” Org. Lett. **13**(2), 188–191 (2011). [CrossRef]

**49. **P. B. Chapple, J. Staromlynska, J. A. Hermann, T. J. McKay, and R. G. McDuff, “Single-beam Z-scan: measurement techniques and analysis,” J. Nonlinear Opt. Phys. Mater. **6**(3), 251–293 (1997). [CrossRef]

**50. **A. Eriksson, M. Lindgren, S. Svensson, and P.-O. Arntzen, “Numerical analysis of Z-scan experiments by use of a mode expansion,” J. Opt. Soc. Am. B **15**(2), 810–816 (1998). [CrossRef]

**51. **I. O. f. Standardization, ISO 11146–1:2005 Lasers and laser-related equipment - Test methods for laser beam widths, divergence angles and beam propagation ratios (Geneva, Switzerland, 2005).