Two-step four-photon absorption

We report a theoretical investigation of effective four-ph oton absorption (4PA) process and propose a two-step 4PA model of threephoton-induced excited-state absorption (ESA). Based on t hree-level rate-equation theory, we find an analytical result for the ef fective 4PA coefficient that depends on the three-photon absorption (3P A) cross-section, excited-state photophysical properties, and laser pulse d uration. We present the analytical theory of the z-dependent nonlinear transmission for straightforwardly yet unambiguously evaluating the 3PA an d effective 4PA coefficients simultaneously. © 2008 Optical Society of America OCIS codes: (190.4180) Multiphoton processes; (320.7100) Ultrafast me asurements References and links 1. D. R. Larson, W. R. Zipfel, R. M. Williams, S. W. Clark, M. P. B ruchez, F. W. Wise, and W. W. Webb, “Watersoluble quantum dots for multiphoton fluorescence imaging in v ivo,” Science300, 1434 (2003). 2. G. S. He, P. P. Markowicz, T. C. Lin, and P. N. Prasad, “Obser vation of stimulated emission by direct three-photon excitation,” Nature (London) 415, 767 (2002). 3. S. Kawata, H. B. Sun, T. Tanaka, and K. Takada, “Finer featu r s for functional microdevices-Micromachines can be created with higher resolution using two-photon absorpt i n,” Nature (London)412, 697 (2001). 4. G. S. He, K. T. Yong, Q. Zheng, Y. Sahoo, A. Baev, A. I. Ryasny anskiy, and P. N. Prasad, “Multi-photon excitation properties of CdSe quantum dots solutions and optical limiti ng behavior in infrared range,” Opt. Express 15, 12818 (2007). http://www.opticsinfobase.org/abstract.cfm?URI=oe-15 -20-21818 5. D. A. Parthenopoulos and P. M. Rentzepis, “Three-dimensio nal ptical storage memory,” Science 245, 843 (1989). 6. J. D. Bhawalkar, G. S. He, and P. N. Prasad, “Nonlinear multi photon processes in organic and polymeric materials,” Rep. Prog. Phys. 59, 1041 (1996). 7. D. S. Corr̂ ea, L. De Boni, D. T. Balogh, and C. R. Mendonça, “Threeand f our-photon excitation of poly(2methoxy-5-(2’-ethylhexyloxy)-1,4-phenylenevinylene)( MEH-PPV),” Adv. Mater.19, 2653 (2007). 8. F. Yoshino, S. Polyakov, M. Liu, and G. Stegeman, “Observat ion of three-photon enhanced four-photon absorption,” Phys. Rev. Lett. 91, 063902 (2003). 9. H. Matsuda, Y. Fujimoto, S. Ito, Y. Nagasawa, H. Miyasaka, T . Asahi, and H. Masuhara, “Development of nearinfrared 35 fs laser microscope and its application to the det ction of threeand four-photon fluorescence of organic microcrystals,” J. Phys. Chem. B 110, 1091 (2006). 10. S. Delysse, P. Filloux, V. Dumarcher, C. Fiorini, and J. M. Nunzi, “Multiphoton absorption in organic dye solutions,” Opt. Mater.9, 347 (1998). 11. C. G. Lu, Y. P. Cui, W. Huang, B. F. Yun, Z. Y. Wang, G. H. Hu, J . Cui, Z. F. Lu, and Y. Qian, “Vibrational resonance enhanced broadband multiphoton absorption in a tr iphenylamine derivatives,” Appl. Phys. Lett. 91, 121111 (2007). 12. G. S. He, L. S. Tan, Q. D. Zheng, and P. N. Prasad, “Multipho on absorbing materials: molecular designs, characterizations, and applications,” Chem. Rev. 108, 1245 (2008). 13. A. Penzkofer and W. Falkenstein, “Three-photon absorpt ion and subsequent excited-state absorption in CdS,” Opt. Commun.16, 247 (1976). #96204 $15.00 USD Received 14 May 2008; revised 18 Jun 2008; accepted 19 Jun 2008; published 24 Jun 2008 (C) 2008 OSA 7 July 2008 / Vol. 16, No. 14 / OPTICS EXPRESS 10208 14. E. Parilov and M. J. Potasek, “Generalized theoretical t reatment and numerical method of time-resolved radially dependent laser pulses interacting with multiphoton absorbers,” J. Opt. Soc. Am B23, 1894 (2006). http://www.opticsinfobase.org/abstract.cfm?URI=josab -23-9-1894 15. R. L. Sutherland, M. C. Brant, J. Heinrichs, J. E. Slagle, D. G. McLean, and P. A. Fleitz, “Excitedstate characterization and effective three-photon absorp tion model of two-photon-induced excited-state absorption in organic push-pull charge-transfer chromopho res,” J. Opt. Soc. Am. B22, 1939 (2005). http://www.opticsinfobase.org/abstract.cfm?URI=josab -22-9-1939. 16. Y. Gao and M. J. Potasek, “Effects of excited-state absor ption on two-photon absorbing materials,” Appl. Opt. 45, 2521 (2006).http://www.opticsinfobase.org/abstract.cfm?URI=ao-45 -11-2521 17. M. Fakis, G. Tsigaridas, I. Polyzos, V. Giannetas, P. Per sphonis, I. Spiliopoulos, J. Mikroyannidis, “Intensity dependent nonlinear absorption of pyrylium chromophores,” Chem. Phys. Lett. 342, 155 (2001). 18. B. Gu, W. Ji, P. S. Patil, S. M. Dharmaprakash, and H. T. Wang , “Two-photon-induced excited-state absorption: Theory and experiment,” Appl. Phys. Lett. 92, 091118 (2008). 19. B. Gu, W. Ji, P. S. Patil, and S. M. Dharmaprakash, “Ultrafa st optical nonlinearities and figures of merit in acceptor-substituted 3,4,5-trimethoxy chalcone derivati ves: Structure-property relationships,” J. Appl. Phys. 103, 103511 (2008). 20. M. I. Demchuk, V. S. Konevskii, N. V. Kuleshov, V. P. Mikhai lov, P. V. Prokoshin, and K. V. Yumashev, “Nonlinear transmission, ultrafast relaxation and three-photo n absorption in reduced SrTiO 3,” Opt. Commun.82, 273 (1991). 21. D. S. Corr̂ ea, L. De Boni, L. Misoguti, I. Cohanoschi, F. E. Hernandez, a nd C. R. Mendonça, “Z-scan theoretical analysis for three-, fourand five-photon absorption,” Opt . Commun.277, 440 (2007).


Introduction
The proliferation of laser systems with ultrafast pulse durations gives the opportunity to investigate high-order nonlinear optical effects.Apart from fundamental academic interest, such investigations are motivated by potential applications for ultrafast all-optical devices and ultrafast information processing.One of the investigations is multiphoton absorption (MPA), which is of direct interest to fluorescence imaging [1], frequency-up-conversion lasing [2], optical microfabrication [3], optical power limiting [4], optical data storage [5], and photodynamic therapy [6].Recently, studies of higher-order MPA related processes, in particular, three-and four-photon absorption (3PA and 4PA), have received considerable attention [7][8][9][10][11].From the viewpoint of both fundamental research and technological applications, unambiguous characterization is crucial to full understanding of MPA processes.
As well known, MPA conventionally refers to the simultaneous absorption of n identical photons, promoting an electron from the ground state of a system, S 0 , to an excited state, S 1 , by virtual intermediate states [12].This process is called a one-step n-photon absorption.For instance, Fig. 1(a) shows the one-step 4PA process.However, many of the strongest MPA molecules are hybrid chromophores where the MPA is accompanied by excited-state absorption (ESA).Under the excitation of high-power ultrafast-pulsed laser, the system simultaneously absorbs n identical photons, promoting an electron from S 0 to S 1 by virtual states.Subsequently, the electron is excited to higher-lying state S h by absorbing another single photon, resulting in MPA-induced ESA.We name this process as a two-step (n + 1)-photon absorption.As an example, Fig. 1(b) illustrates an effective 4PA process of 3PA-induced ESA process.Many investigations have revealed that the MPA-induced ESA response of a given sample is different in the magnitude -depending on the laser characteristics, in particular, on the laser pulse duration [12][13][14][15][16]. Two-photon-induced ESA has been experimentally and theoretically investigated [14][15][16][17][18][19].For a 3PA-induced ESA process, observations have been experimentally reported for CdS [13] and SrTiO 3 [20].However, there is no report in the literature on the unambiguous characterization of 3PA-induced ESA.
In the work presented here, we theoretically investigate a three-level two-step 4PA process of 3PA-induced ESA.The effective 4PA coefficient, which depends on the 3PA cross-section, excited-state photophysical properties, and laser pulse duration, is analytically presented.The analytical expression allows ones to unambiguously yet conveniently separate the contribution of 3PA from the effective 4PA in the case that both processes occur simultaneously.

Theory
Under the excitation of ultrafast laser pulses, a simplified three-level model can adequately explain nonlinear absorption processes in polyatomic molecules because long-lived states such as triplet-triplet state can be justifiably neglected [15].Figure 1(b) illustrates the proposed mechanism for three-level two-step 4PA process.This model is similar to the model of two-photoninduced ESA, in which 2PA is replaced by 3PA from the ground state [18].The photodynamics of the system is described as follows.Initial absorption of the laser pulse promotes electrons from S 0 to S 1 , by simultaneous absorption of three photons.From S 1 , the electrons relax to S 0 or are promoted to S h , by absorbing single photons.The propagation equation of the pulse intensity is given by and the time derivatives of populations of states S 0 , S 1 , and S h are governed by where z ′ is the propagation length inside the sample; σ 3 and σ e are the molecular 3PA and ESA cross-sections, respectively; τ e and τ h are the lifetimes of S 1 and S h , respectively; and hω is the incident photon energy.Usually the relaxation time τ h from S h to S 1 is very short, so the electrons pumped onto S h jump back promptly, and the population on S h is close to zero.In this analysis, we assume that for the laser intensities used, the ground state S 0 is not sufficiently depleted, therefore N S 1 ≪ N and N S h ≪ N, we set N S 0 ≃ N.Under these conditions, Eq. ( 2) becomes Solving this equation, the population density of S 1 is given by The above expression can be written as where To characterize the MPA properties of a material, time averaging techniques have been extensively exploited in both Z-scan and nonlinear transmission experiments [15,21].It should be noted that integrating Eq. ( 1) over time yields an equation that describes the attenuation of the pulse fluence.The last term on the right-hand side will contain an integral, +∞ −∞ F(t)I 4 (t)dt.The equation is unchanged if we replace this integral with F +∞ −∞ I 4 (t)dt, where is the temporal average of F(t) weighted by intensity distribution I 4 (t).
In solving Eq. ( 1), the transmitted intensity is numerically calculated and then it is integrated over time and space to give the transmitted energy.We assume that a reasonable approximation to the transmitted energy can be yielded by replacing F(t) with its corresponding time average F in Eq. ( 7), solving Eq. ( 1) to determine the transmitted intensity, and then integrating temporally and spatially to obtain the transmitted energy.We also assume a temporal Gaussian profile for the incident pulses.With these approximations, Eq. ( 1) can be revised as dI(r, z,t) dz ′ = −[α 3 I 2 (r, z,t) + α 4 I 3 (r, z,t)]I(r, z,t), ( where Here α 3 = σ 3 N is the 3PA coefficient [3PA process is usually quantified by 3PA cross-section σ 3PA through the conversion formula σ 3PA = α 3 (hω) 2 /N = σ 3 (hω) 2 ], α 4 is the effective 4PA coefficient, and τ is the half-width at e −1 of the maximum for the pulse duration.Obviously, Eq. ( 9) predicts that the effective 4PA coefficient depends on 3PA coefficient, ESA crosssection, lifetime of the first excited state, and laser pulse duration.In the limit of τ e ≪ τ, Eq. ( 9) becomes α 4 = σ e α 3 τ e /3hω.Under the extreme condition of τ e ≫ τ, Eq. ( 9) gives α 4 = √ πσ e α 3 τ/(3 √ 3hω).In particular, we find an empirical analytic expression as follows: Figure 2 shows the effective 4PA coefficient α 4 as a function of the pulse-duration τ for τ e = 1 ps obtained by numerical simulation (circles) and our analytic solution (solid line), respectively.Clearly, our analytical solution agrees well with the numerical simulation.For the 3PA-induced ESA process, the ESA lifetimes and ESA cross-section are of a few picoseconds and usually close to 10 −17 cm 2 , respectively [10,13,20].It is interesting to note that the difference of effective 3PA coefficients measured in the nanosecond regime and in the femtosecond regime for the same sample could be too large in order of magnitude [12,13].This fact implies that, as illustrated in Fig. 2, in the nanosecond regime, the 3PA-induced ESA could play a more important role than the pure 3PA process.With the excitation of femtosecond laser pulses, however, 3PA is the dominant mechanism of whole nonlinear absorption process.At low intensity, 3PA process is predominant because the N S 1 population is not sufficient and no electron could be excited from S 1 to S h .At excitation intensity exceeding a critical value, the competition between 3PA and ESA occurs.That is, Eqs.(8) and Eq. ( 12) are valid only above a critical intensity.Using the steady-state condition in Eq. (3), we estimate the critical intensity where N c S 1 is the critical population of S 1 .Under the excitation of low intensities (I < I c ), the sample only exhibits 3PA effect.When the intensity is larger than the critical intensity, however, the simultaneous appearance of 3PA and 3PA-induced-ESA effects is observable.
The popular yet effective techniques to characterize MPA related studies are Z-scan and nonlinear transmittance measurements [12,21].In solving Eq. ( 8), we assume that the spatiotemporal profiles of incident pulses are Gaussian.To simulate the z-dependent nonlinear transmission, the transmitted intensity is numerically calculated and then it is integrated over time and space to give the transmitted energy.We determine the same Z-scan traces using Eqs.( 1) and (3), and ( 8) and ( 9), respectively.It indicates that the effective 4PA theory of 3PA-induced ESA we presented correctly describes the Z-scan measurements.
For a thin sample with simultaneous 3PA and effective 4PA, we obtain the normalized energy transmittance T (z, p, h), in polynomial series form, as follows Here p = p 0 /(1 + z 2 /z 2 0 ) and h = h 0 /(1 + z 2 /z 2 0 ), where p 0 = (2α 3 I 2 00 L) 1/2 and h 0 = (3α 4 I 3 00 L) 1/3 are the on-axis peak phase shifts caused by the 3PA and 4PA processes, respectively; I 00 is the on-axis peak intensity at the beam waist; z 0 is the Rayleigh length of the Gaussian beam; and L is the sample physical length.The coefficient a mn in Eq. ( 12) is listed in Table 1.The above empirical expression has an utmost error of 0.3 % with the exact numerical results when 0 ≤ p 0 ≤ π and/or 0 ≤ h 0 ≤ π.For the pure 3PA (α 4 = 0) and 4PA (α 3 = 0), Eq. ( 12) degenerates into the following forms, respectively: In Eq. ( 12), there are two free parameters, namely α 3 and α 4 .Nevertheless, the α 3 value is unambiguously obtained by Eq. ( 13) from the measurements at low intensities where ESA is negligible.With the measured α 3 , the value of α 4 can be extracted from the z-dependent nonlinear transmission measurements at high intensity by the use of Eq. (12).
It should be emphasized that Eq. ( 12) is applicable to the following three cases: (i) simultaneous occurrence of 3PA and 3PA-induced ESA processes in polyatomic molecules, (ii) appearance of both 3PA and 3PA-generated free-carrier absorption in semiconductors such as in CdS [13], and (iii) interplay between 3PA and intrinsic 4PA, a situation that exists, for example, in a conjugated polymer at 1700 nm [7].

Conclusion
In summary, we have proposed a three-level two-step 4PA model of 3PA-induced ESA.We have derived an analytical expression for the effective 4PA coefficient that depends on the 3PA cross-section, excited-state photophysical properties, and laser pulse duration.We have also presented the analytical theory of the z-dependent nonlinear transmission for straightforwardly yet unambiguously characterizing 3PA-induced ESA.

Fig. 2 .
Fig.2.Pulse-width dependence of effective 4PA coefficient for τ e = 1 ps.The circles are numerical simulations, while the solid line is the analytical result.The effective 4PA coefficient is normalized to σ e α 3 τ e /3hω.