Abstract
The recent interpretation of experiments on the nonlinear non-resonant birefringence induced in a weak probe beam by a high intensity pump beam in air and its constituents has stimulated interest in the non-resonant birefringence due to higher-order Kerr nonlinearities. Here a simple formalism is invoked to determine the non-resonant birefringence for higher-order Kerr coefficients. Some general relations between nonlinear coefficients with arbitrary frequency inputs are also derived for isotropic media. It is shown that the previous linear extrapolations for higher-order birefringence (based on literature values of n 2 and n 4) are not strictly valid, although the errors introduced in the values of the reported higher- order Kerr coefficients are a few percent.
©2011 Optical Society of America
1. Introduction
There has been a growing interest in higher-order nonlinear coefficients χ (2 m + 1 ) of odd order where m = 2, 3, 4 in Kerr media [1–5]. Earlier work on higher-order nonlinearities in semiconductors were related to charge carrier excitation either due to cascaded effects, saturation effects such as band filling etc. [6,7]. To the best of our knowledge, the first direct measurement of the fifth order Kerr nonlinearity is the work of Arabat and Etchepare who measured the non-resonant χ (5) for a WG630 Schott glass at intensities of 100’s GW/cm2 [3]. More recently, fifth order nonlinearities have been measured in a number of glasses and organic materials [1,4,5]. Chen and associates at Cornell also verified that there is resonant enhancement of the fifth order nonlinearity for wavelengths approaching the absorption edge of a glass and were even able to use this to estimate the seventh order susceptibility in a chalcogenide glass [4].
The fact that an intense beam induces a nonlinear birefringence, (I – local intensity), in any medium is well-known since the early days of nonlinear optics [8]. Such a birefringence is usually formulated in terms of the third order susceptibility χ (3) for Kerr nonlinearities involving electronic states in a medium. For isotropic Kerr media, and the proportionality constant depends on how many unique eigenmodes are present. For example, the numerical factor is 1/3 for a single intense beam where-as it is 2/3 for a strong pump, weak probe geometry.
The general formulation of the nonlinear birefringence problem requires calculating the nonlinear index changes produced by a strong pump beam either for the pump itself, or for a second beam, usually a weak probe beam, with different frequency, propagation direction and/or polarization properties from those of the pump. In isotropic media this normally requires knowledge of the ratio of at least two nonlinear susceptibilities and their dispersion with frequency [8]. This can be a daunting problem since the number of different susceptibility terms increases rapidly with the order of the nonlinearity, i.e. with “m” in The situation simplifies considerably for isotropic media in the non-resonant regime for the susceptibilities since there is only one independent nonlinear susceptibility for each value of m [3,9,10]. The formulation of the nonlinear birefringence problem described here relies strongly on this fact and a formula is derived for arbitrary order nonlinearities.
The most recent interest in nonlinear birefringence due to higher-order Kerr coefficients was stimulated by experiments at ~800nm on filaments which form in air at high (>10 TW/cm2) laser intensities [1,11]. In order to explain their birefringence measurements, Loriot et al. assumed non-resonant nonlinear index coefficients up to n 10 (involving χ (11)). There is some controversy in the filamentation community concerning the interpretation of the measured birefringence but in this paper we simply focus on their analytical expression for the birefringence due to higher-order Kerr effects [12]. They obtained their contributions to the birefringence from the well-known relations between the tensor coefficients for χ (3) and that obtained for χ (5) by Arabat and Etchepare based on an anharmonic oscillator model, and then linear extrapolation to higher-orders [3]. Although this nonlinear oscillator model fails to reproduce accurately the frequency dispersion of the third (and presumably higher-order) nonlinearities obtained from quantum mechanics, it does give non-resonant (ω→0) results for χ (3), χ (5) and χ (7), albeit not in terms of physically measurable parameters [2,9,10]. We are not aware of any extension to yet higher-order nonlinearities.
The most frequently cited nonlinear index coefficient n 2 in isotropic media is defined for a single intense beam (x-polarized, for example) as [8]:
Note that we have introduced a notation for n 2 in which the beam which induces the nonlinear index change is the second ( + )ω argument and the first argument (-ω) identifies the eigenmode in which the index change occurs. If another eigenmode is present such as a weak “probe” beam of the same or different frequency ωp, also x-polarized but travelling at a small angle to the “pump” beam, the appropriate nonlinearity in this case is defined asIn this paper we derive from first principles the nonlinear birefringence introduced by higher-order Kerr coefficients in the non-resonant limit in an isotropic medium. We start by reformulating the well-known relations for χ (3) in terms of combinatorial expressions which then provide a simple procedure for dealing with the higher-order Kerr nonlinearities. We find that the extrapolation used by Loriot et al. is not correct for the pump-probe geometry they considered [1].The procedure followed here is a three-step process:
- 1. The nonlinear polarizations and are calculated in terms of the nonlinear susceptibilities by permuting the input eigenmodes via their frequencies for an isotropic medium.
- 2. The relation between thesusceptibilities is found for an isotropic medium by permuting the polarizations. This utilizes the concept that the nonlinear polarization in an isotropic medium must be independent of the choice of axes.
- 3. The square of the refractive indices, i.e. are calculated from the respective polarizations and the square root of each is taken to give the nonlinear birefringence.
2. Pump-probe geometry
Here we consider the specific geometry of the Loriot et al. experiment shown in Fig. 1 [1]. An intense plane wave of the form
is assumed to propagate along the z-axis in an isotropic material, i.e. the x-axis is chosen parallel to the polarization of the intense beam. A second probe beam (subscript “p”) of frequency ωp = ω is also present but propagating at a small angle from the z-axis in the y-z plane (making it a different eigenmode from the pump beam). Its polarization has equal x and y-components written asand the nonlinear polarization induced in the probe beam is written as3. Nonlinear polarizabilities
The third order nonlinear polarization induced by the pump beam in the molecules of the air, as experienced by the probe beam, is
The susceptibilitiesare values of the coefficient averaged over the constituent air molecules, i.e.here wq is the fraction of the number density corresponding to species q, i.e. nitrogen, oxygen etc. In the non-resonant limit (identified by the superscript ~), the imaginary part of the susceptibility is negligibly small, zero for ω = 0, and [8]are all real. Therefore,A different way to arrive at this result is to note that there are three separate input positions for frequency in the expression for χ (3) giving 3!( = 3x2x1) different possibilities when they are permuted over the three input fields. (In nonlinear optics can be treated as separate eigenmodes because they have different frequencies in a mixing process, i.e. + ω and –ω.) Thus there are three separate NLO (nonlinear optics) eigenmodes, each of which appears just once, so that the total number of unique terms is given by 3!/1!1!1!, i.e.Similarly, for the polarization nonlinearly induced along the y-axis by the strong x-polarized field,In the non-resonant limit all sixare equal so that with 3 independent eigenmodes each of which appears only once,The nonlinear susceptibilities are now abbreviated so that and are written as and respectively. (This will also subsequently be extended to higher-order susceptibilities.) Furthermore, since the labeling of the axes in isotropic media is arbitrary, . Applying the same arguments as for the probe case but with two equal co-polarized fields at + ω for but not for, the nonlinear polarizations experienced by the pump beam are Generalizing these results to the 2m + 1 case [4], Here the co-ordinate subscripts (2m + 2)x mean that there is a total of 2m + 2 “x co-ordinates,” one referring to the output polarization of the probe, and one of the remaining 2m + 1 refers to the input probe polarization, interspersed amongst the 2m others associated with ± ω of the pump beam. The (2)y,(2m)x means that there are 2 “y” co-ordinates, one always being the first co-ordinate which refers to the output probe polarization, and the second to the input probe beam polarization interspersed amongst the 2m ± ω x-polarizations associated with the pump beam.Therefore the total polarization for the probe beam is given by
For the pump beam, the ωp on the input side is replaced by another ω and hence there are m + 1 + ω ’s but still m -ω’s so thatFor, the first 2 y’s belong to the probe beam and there are still m ω’s and m -ω’s, just like in so that4. Total nonlinear birefringence
It is clear from Eqs. (15)–(18), that in order to find the birefringence, the relationship between the nonlinear susceptibilitiesandmust be found. This depends on the symmetry properties of the medium. Even for isotropic media these are relatively complicated calculations and hence they are summarized in the Appendix along with some general results valid for all frequencies. Making the results specific to the non-resonant, isotropic medium case, Eq. (A17) is
For the pump-probe geometry in the non-resonant limit, Eq. (19) is inserted into Eqs. (16) and (18) to give so that both the x- and y-components of the nonlinear polarization are given in terms of the same susceptibilities. Noting that from Eq. (13) and combining Eqs. (15), (17), (20), and (21) now leads directly toin which the coefficient is given byThis form was chosen so that for the individual nonlinearities m In order to make contact with the experimental data in reference 1 we focus on the nonlinear refractive indices for the pump-probe case so that the nonlinear birefringence is given byThe expansion offor small b is well known from textbooks [13], to be:ThereforeThe leading term (s = 1), expanded up to n 10 (largest term reported in reference 1), isTerms with s≥2 contain products of the nonlinear coefficients. Including all of the terms up to I 5,Note that all the numerical pre-factors in this case are all less than 2.5. The products of different nonlinear coefficients are limited to 2 here. However, products of more than two nonlinear coefficients occur for higher-orders in intensity, the first one being From Eq. (29) it is evident that in a strict mathematical sense the nonlinear birefringence cannot be used as a means to measure the nonlinear coefficients higher than.There is no direct correlation between the coefficientand the corresponding power of the intensity Im for m>1 due to the existence of the product terms. However, it makes sense to use the simplified notation of Eq. (29) if the relation holds.5. Comparison with experiments on air
Reference 1 contains data measured in air and its constituents for n 2(-ωp;ω) → n 8(-ωp;ω) and also n 10(-ωp;ω) for argon. Based on their values, with m = rq + vu and m≤5 is always satisfied in air. Assuming that the only nonlinear mechanism present is the Kerr effect, the nonlinear birefringence is given by the leading term, Eq. (28), which can be expressed as the series
This result should be compared with the expansion used by Loriot et al. [1]. Based on a linear extrapolation from the first two terms which Loriot et al. obtained from the literature [3,8] they assumed the seriesin their analysis of their data. Note that in both series the numerical pre-factors 2m/(2m + 1) and 2m/(2m + 1) respectively converge to unity for large m. A graphical comparison of the two expansions is given in Fig. 2 . In Fig. 2(a) we compare the expansion terms as deduced from Eq. (30), , to the ones derived by Loriot et al. . As m is increased their difference is maximized for m = 11. The relative deviation of Loriot et al.’s expansion terms as compared to the analytically derived factors is depicted in Fig. 2(b). For m = 11 the relative error peaks at 6.25%. Furthermore, the Loriot et al. formulation systematically underestimates the expansion term coefficients and thus leads to an overestimation of the correspondingcoefficient for m>2.
Fig. 2 Comparison between the expansion coefficients estimated by the two models. (a) Coefficients corresponding to χ(m) terms. (□) analytical model, (•) Loriot et al. estimation, dotted/dashed lines are a guide to the eye. (b) Relative error for the various coefficients of the χ(m) terms. (Dotted lines are guides to the eye).
6. Conclusions
Expressions for the non-resonant, nonlinear birefringence induced in a probe beam (frequencyωp) by a strong pump beam of the same frequency in an isotropic medium have been derived for nonlinear Kerr indices n 2 m(-ωp;ω) for arbitrary m. This was made possible by using combinatorial approaches and by assuming that in isotropic media there is only one unique value forfor each value of m which was verified previously in the literature for m = 1, 2. Some general relations for arbitrary frequency inputs were also derived.
Because the polarization, linear and nonlinear, induced in a material depends on the square of the refractive index, the nonlinear birefringence was found to depend not only on the intensity-dependent refractive index coefficients n 2 m(-ωp;ω) but also on the products of the various nonlinear index coefficients. Comparison with existing experiments in air and its constituents showed that the product terms were negligible in that case.
An analytical series was found to describe the nonlinear birefringence. This series was different from that assumed by Loriot et. al based on a linear extrapolation of two points. Since in both cases the individual numerical factors for n 2 m(-ωp;ω) converged to unity for increasing m, the errors introduced into the analysis of the data were relatively small.
Appendix A. Relationships between the nonlinear susceptibilities
In this Appendix the relations between the and are derived, some for arbitrary frequency inputs. Isotropy requires that each coordinate (x and y) comes in pairs. It also requires that the nonlinear polarization should be independent of the orientation of any axis system used. Consider first the general case (unrelated to the previous discussion) of three, parallel, co-polarized (along the x-axis) input fields E1, E2 and E3 with different frequencies ω 1, ω 2 and ω 3 producing the field ω 4 via. The third order nonlinear polarization (along the x-axis) is
Now consider the axis system (x', y') rotated 45° from the original x-axis [8]. The three fields have the following components along the x'-axis and y'-axis
For isotropic media, etc., and hence the nonlinear polarization induced along the x'-axis is given by
The nonlinear polarization in Eq. (A3) can also be obtained by projecting the nonlinear polarization given by Eq. (A1) onto the x'-axis to give
Since Eqs. (A3) and (A4) must give the same result which is valid for any frequencies,
Note that any isotropic material, for example a mature electron plasma, which exhibits third order effects such as third harmonic generation [14,15] must have all of these coefficients non-zero and related as given by Eq. (A5). In the non-resonant limit it can easily be shown that
The same result holds for pump beam, i.e. Although this result is valid for a single medium, extension to multi-component air is trivial giving
An alternate and more compact approach for arriving at the same result is to again resort to combinatorial mathematics. Since there are three input polarization components, two y'-polarized and one x'-polarized, which can be permuted among the three input eigenmodes (frequencies), there are 3! possibilities for permuting the corresponding polarization components in. Because there must be two identical polarization components (y') and only one x′, there are 3!/2!1!unique possibilities and Eq. (A5) can be re-written in the non-resonant limit as
The evaluation of the relation between (and subsequently the yet higher-order susceptibilities) has additional aspects (relative to the χ (3) case) associated with the etc. terms. Again assuming the general case of five, parallel, co-polarized (along the x-axis) input fields namely E1, E2, E3, E4, and E5 with different frequencies ω 1, ω 2, ω 3, ω 4 and ω 5 producing the field ω 6 via . This produces the nonlinear polarization (along the x-axis)
Now consider again the axis system (x', y') rotated 45° from the original x-axis. The five input x-polarized fields again have components along the x'-axis and y'-axis. Note that both mixed polarization terms likeas well as contribute to the nonlinear polarization induced along the x'-axis, For the first one, there are 5! input slots for the polarization of which 3 are identical (x') and the two others are also identical (y') and, for the second one, there are 4 (y') identical slots and only the x' is a single slot. Hence the number of unique combinations are 5!/3!2! and 5!/4!1! respectively for and for . There are further simplifications in the non-resonant limit etc. so that
The nonlinear polarization in Eq. A10 can also be obtained by projecting the nonlinear polarization given by Eq. (A9) onto the x'-axis to give
Again Eqs. (A10) and (A11) must yield identical results and noting again from references 2 and 8 that =etc. yields for the cases of interest here in the non-resonant limit
Consider briefly the 7’th and 9’th order susceptibilities. The same procedures as for the 3′rd and 5′th order cases are used. In order to derive the relationship between the different etc. seven co-polarized input fields are considered, first in the x, y, z coordinate system and then in an axis system rotated 45° in the x-y plane. In this case, the mixed polarization terms
all contribute to the nonlinear polarizations induced along the x'-axis, . The number of unique combinations are 7!/5!2!, 7!/4!3!and 7!/6!1!respectively for the three cases. Thus again in the non-resonant limit
Based on the preceding results, only one, unique, nonlinear susceptibility is expected for an isotropic material in the non-resonant limit for each order “2m+1” of. Therefore all the mixed polarization susceptibilities are equal which gives
Again using the same approach, for the 9’th order susceptibility,
In the non-resonant limit
These results suggest simple relations governing the relationship between the susceptibilities, namely
For frequency inputs ω 1, ω 2, ω 3,.. ω 2 m +1 giving an output frequency ω 2 m +2 for isotropic media, the above formulas suggest the following general result:
which gives
Acknowledgments
This work was partially supported by the European Union Marie Curie Excellence Grant “MULTIRAD” MEXT-CT-2006–042683.
References and links
1. V. Loriot, E. Hertz, O. Faucher, and B. Lavorel, “Measurement of high order Kerr refractive index of major air components,” Opt. Express 17(16), 13429–13434 (2009). [CrossRef] [PubMed]
2. V. Loriot, E. Hertz, O. Faucher, and B. Lavorel, “Measurement of high order Kerr refractive index of major air components: erratum,” Opt. Express 18(3), 3011–3012 (2010). [CrossRef]
3. J. Arabat and J. Etchepare, “Nonresonant fifth-order nonlinearities induced by ultrashort intense pulses,” J. Opt. Soc. Am. B 10(12), 2377–2382 (1993). [CrossRef]
4. Y.-F. Chen, K. Beckwitt, F. W. Wise, B. G. Aitken, J. S. Sanghera, and I. D. Aggarwal, “Measurement of fifth- and seventh-order nonlinearities in glasses,” J. Opt. Soc. Am. 23(2), 347–352 (2006). [CrossRef]
5. For example,R. A. Ganeev, M. Baba, M. Morita, A. I. Ryasnyansky, M. Suzuki, M. Turu, and H. Kuroda, “Fifth-order optical nonlinearity of pseudoisocyanine solution at 529 nm,” J. Opt. A, Pure Appl. Opt. 6(2), 282–287 (2004). [CrossRef]
6. For example,S. Wu, X.-C. Zhang, and R. L. Fork, “Direct experimental observation of interactive third and fifth order nonlinearities in a time- and space-resolved four-wave mixing experiment,” Appl. Phys. Lett. 61, 1919–1921 (1992). [CrossRef]
7. E. J. Canto-Said, D. J. Hagan, J. Young, and E. W. Van Stryland, “Degenerate four-wave mixing measurements of high order nonlinearities in semiconductors,” IEEE J. Quantum Electron. 27(10), 2274–2280 (1991). [CrossRef]
8. R. W. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, Amsterdam, 2008).
9. V. I. Zavelishko, V. A. Martynov, S. M. Saltiel, and V. G. Tunkin, “Optical nonlinear fourth- and fifth-order susceptibilities,” Sov. J. Quantum Electron. 5(11), 1392–1393 (1975). [CrossRef]
10. S. V. Popov, Y. P. Svirko, and N. I. Zheludev, Susceptibility Tensors for Nonlinear Optics (Taylor and Francis, 1995)
11. A. Couairon and A. Mysyrowicz, “Femtosecond filamentation in transparent media,” Phys. Rep. 441(2–4), 47–189 (2007). [CrossRef]
12. M. Kolesik, D. Mirell, J.-C. Diels, and J. V. Moloney, “On the higher-order Kerr effect in femtosecond filaments,” Opt. Lett. 35(21), 3685–3687 (2010). [CrossRef] [PubMed]
13. M. Abramowitz, and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th ed. (Dover, 1972).
14. S. Suntsov, D. Abdollahpour, D. G. Papazoglou, and S. Tzortzakis, “Efficient third-harmonic generation through tailored IR femtosecond laser pulse filamentation in air,” Opt. Express 17(5), 3190–3195 (2009). [CrossRef] [PubMed]
15. S. Suntsov, D. Abdollahpour, D. G. Papazoglou, and S. Tzortzakis, “Filamentation-induced third-harmonic generation in air via plasma-enhanced third-order susceptibility,” Phys. Rev. A 81(3), 033817 (2010). [CrossRef]