Third-harmonic generation microscopy with focus-engineered beams : a numerical study

We use a vector field model to analyze third-harmonic generation (THG) from model geometries (interfaces, slabs, periodic structures) illuminated by Hermite-Gaussian (HG) and Laguerre-Gaussian (LG) beams focused by a high NA lens. Calculations show that phase matching conditions are significantly affected by the tailoring of the field distribution near focus. In the case of an interface parallel to the optical axis illuminated by an odd HG mode, the emission patterns and signal level reflect the relative orientation of the interface and the focal field structure. In the case of slabs and periodic structures, the emission patterns reflect the interplay between focal field distribution (amplitude and phase) and sample structure. Forward-to-backward emission ratios using different beam shapes provide sub-wavelength information about sample spatial frequencies. © 2008 Optical Society of America OCIS codes: (180.6900) Three-dimensional microscopy; (190.4160) Multiharmonic generation; (170.3880) Medical and biomedical imaging. References and links 1. E. Yew and C. Sheppard, ”Second harmonic generation microscopy with tightly focused linearly and radially polarized beams,” Opt. Commun. 275, 453-457 (2007). 2. K. Yoshiki, R. Kanamaru, M. Hashimoto, N. Hashimoto, and T. Araki , ”Second-harmonic-generation microscope using eight-segment polarization-mode converter to observe three-dimensional molecular orientation,” Opt. Lett. 32, 1680-1682 (2007). 3. V. V. Krishnamachari and E. O. Potma, ”Focus-engineered coherent anti-Stokes Raman scattering microscopy: a numerical investigation,” J. Opt. Soc. Am. A 24, 1138-1147 (2007). 4. V. V. Krishnamachari and E. O. Potma, ”Imaging chemical interfaces perpendicular to the optical axis with focus-engineered coherent anti-Stokes Raman scattering microscopy,” Chem. Phys. 341, 81-88 (2007). 5. Y. Barad, H. Eisenberg, M. Horowitz, and Y. Silberberg, ”Nonlinear scanning laser microscopy by third harmonic generation,” Appl. Phys. Lett. 70, 922-924 (1997). 6. M. Müller, J. Squier, K. R. Wilson, and G. J. Brakenhoff, ”3D-microscopy of transparent objects using thirdharmonic generation,” J. Microsc. 191, 266-274 (1998). 7. D. Débarre and E. Beaurepaire, ”Quantitative characterization of biological liquids for third-harmonic generation microscopy,” Biophys. J. 92, 603-612 (2007). 8. D. Yelin, and Y. Silberberg, ”Laser scanning third-harmonic generation microscopy in biology,” Opt. Express 5 (1999). 9. D. Oron, D. Yelin, E. Tal, S. Raz, R. Fachima, and Y. Silberberg, ”Depth-resolved structural imaging by thirdharmonic generation microscopy,” J. Struct. Biol. 147, 3-11 (2004). 10. D. Débarre, W. Supatto, E. Farge, B. Moulia, M.-C. Schanne-Klein, and E. Beaurepaire, ”Velocimetric thirdharmonic generation microscopy: micrometer-scale quantification of morphogenetic movements in unstained embryos,” Opt. Lett. 29, 2881-2883 (2004). 11. C.-K. Sun, S.-W. Chu, S.-Y. Chen, T.-H. Tsai, T.-M. Liu, C.-Y. Lin, and H.-J. Tsai, ”Higher harmonic generation microscopy for developmental biology,” J. Struct. Biol. 147, 19-30 (2004). #98503 $15.00 USD Received 8 Jul 2008; revised 21 Aug 2008; accepted 21 Aug 2008; published 3 Sep 2008 (C) 2008 OSA 15 September 2008 / Vol. 16, No. 19 / OPTICS EXPRESS 14703 12. W. Supatto, D. Débarre, B. Moulia, E. Brouzés, J.-L. Martin, E. Farge, and E. Beaurepaire, ”In vivo modulation of morphogenetic movements in Drosophila embryos with femtosecond laser pulses,” Proc. Nat. Acad. Sci. USA 102, 1047-1052 (2005). 13. D. Débarre, W. Supatto, A.-M. Pena, A. Fabre, T. Tordjmann, L. Combettes, M.-C. Schanne-Klein, and E. Beaurepaire, ”Imaging lipid bodies in cells and tissues using third-harmonic generation microscopy,” Nat. Methods 3, 47-53 (2006). 14. D. Débarre, W. Supatto, and E. Beaurepaire, ”Structure sensitivity in third-harmonic generation microscopy,” Opt. Lett. 30, 2134-2136 (2005). 15. J.-X. Cheng and X. S. Xie, ”Green’s function formulation for third harmonic generation microscopy,” J. Opt. Soc. Am. B 19, 1604-1610 (2002). 16. B. Richards and E. Wolf, ”Electromagnetic diffraction in optical systems II. Structure of the image field in an aplanetic system.,” Proc. Royal Soc. A 253, 358-379 (1959). 17. L. Novotny and B. Hecht, Principles of nano-optics (Cambridge Univ Press, 2006). 18. Boyd, R. W. Nonlinear optics, 2nd edition, (Academic Press 2003). 19. H. Kogelnik, and T. Li, ”Laser beams and resonators,” Appl. Opt. 5, 1550 (1966). 20. K. Youngworth and T. Brown, ”Focusing of high numerical aperture cylindrical-vector beams,” Opt. Express 7, 77-87 (2000). 21. D. Débarre, N. Olivier, and E. Beaurepaire, ”Signal epidetection in third-harmonic generation microscopy of turbid media,” Opt. Express 15, 8913-8924 (2007). 22. E. Y. S. Yew, and C. J. R. Sheppard, ”Fractional Gouy phase,” Opt. Lett. 33, 1363-1365 (2008). 23. J. Mertz, and L. Moreaux, ”Second-harmonic generation by focused excitation of inhomogeneously distributed scatterers,” Opt. Commun. 196, 325-330 (2001). 24. C. J. R. Sheppard, ”High-aperture beams,” J. Opt. Soc. Am. A 18, 1579-1587 (2001). 25. S. Quabis, R. Dorn, M. Eberler, O. Glöckl, and G. Leuchs, ”Focusing light to a tighter spot,” Opt. Commun. 179, 1-7 (2000). 26. S. Carrasco, B. E. A. Saleh, M. C. Teich, and J. T. Fourkas, ”Secondand third-harmonic generation with vector Gaussian beams,” J. Opt. Soc. Am. B 23, 2134-2141 (2006). 27. S. S. Sherif, M. R. Foreman, and P. Török, ”Eigenfunction expansion of the electric fields in the focal region of a high numerical aperture focusing system,” Opt. Express 16, 3397-3407 (2008). 28. M. R. Foreman, S. S. Sherif, P. R. T. Munro, and P. Török, ”Inversion of the Debye-Wolf diffraction integral using an eigenfunction representation of the electric fields in the focal region,” Opt. Express 16, 4901-4917 (2008).


Introduction
Coherent nonlinear microscopies based on parametric processes such as coherent anti-Stokes Raman scattering (CARS), second-harmonic generation (SHG) or third-harmonic generation (THG) are receiving considerable attention.All these imaging modalities are compatible with two-photon excited fluorescence microscopy and provide different information on biological and non biological media with micrometer 3D resolution.A remarkable property of coherent nonlinear imaging techniques is that they are very sensitive to both the sub-micrometer sample structure and to the focal field structure.Indeed, the far field signal results from the coherent superposition in the detection plane of waves emitted at different locations near focus, and interference phenomena define the visibility of a particular distribution of scatterers.
Engineering the focal field structure is therefore a logical step in coherent nonlinear imaging.Intensity, phase and polarization may be modified by controlling the wavefront at the pupil of the objective, resulting in a modulation of phase-matching conditions and far-field emission patterns.This concept has been explored recently for SHG microscopy, where focused beams with strong axial components were used to enhance signal from fibers parallel to the optical axis [1,2], and in CARS microscopy where focus engineering was used to highlight interfaces [3,4].
In this article, we study the use of engineered beams in THG microscopy.THG microscopy relies on the third-order nonlinear susceptibility χ (3) of the sample to provide contrast [5,6,7], and has proved effective for imaging biological samples [8,9,10,11,12,13].The imaging properties of THG microscopy strongly depend on the field distribution near focus.The most salient characteristic of THG microscopy with Gaussian beams is that no signal is obtained from a homogeneous normally dispersive sample [5].Signal is obtained around χ (3) inhomo-geneities, with an efficiency depending on the relative sizes of the inhomogeneity and of the focal volume [14].Since most materials have a non-negligible χ (3) , THG signal creation for a particular sample geometry is essentially determined by interference effects.Focus engineered THG microscopy is therefore expected to give access to sub-wavelength structural information about the sample.We here present a numerical study of vectorial and phase-matching aspects of THG by tightly focused Gaussian, Hermite-Gaussian (HG), and Laguerre-Gaussian (LG) beams incident on slabs, interfaces, and axially periodic samples.These calculations provide insight on the interplay between field and sample structure in THG microscopy with focused complex beams, and should more generally prove useful for designing coherent nonlinear microscopy (SHG, THG, CARS) experiments with engineered beams.Our strategy for simulations follows the framework described in [15], but integrates a complete vector field model to account for vectorial effects with arbitrary excitation beam profiles.

Theory and numerical implementation
A general method for analyzing signal generation in nonlinear microscopy can be described as follows (see Figure 1 for notations).First, the focal field distribution is calculated using a Debye-Wolf diffraction integral [16].Then the induced nonlinear polarization in the focal volume is calculated for a given sample geometry.Finally, the resulting nonlinear field is propagated using Green's functions into the far field [15], where signal level and radiation patterns are analyzed.This theoretical description is summarized below.

Excitation field near focus
The field distribution near the focus of a high numerical aperture (NA) and anti-reflexion-coated objective lens given an arbitrary field at the back pupil can be calculated using the angular spectrum representation method [16,17]: with: where E 0 (θ , Φ) describes the field phase and intensity distribution at the back aperture of the objective, k = k ω = 2πω/n ω is the wavenumber, f is the focal length of the objective, n ω is the refractive index at frequency ω, (ρ, θ , z) are cylindrical coordinates near focus, and θ max = sin −1 (NA/n) is the maximum focusing angle of the objective.
If we expand E 0 (θ , Φ) as a polynomial expansion of cos(Φ) and sin(Φ) functions, this 2D integral can be reduced to a 1D integral involving Bessel functions J n .We can then use the following abbreviations to express the focal fields of the various beam modes considered in this study: where f w (θ ) = exp(−(sin(θ )/( f 0 sinθ max )) 2 ) is a filling factor that takes into account the ratio ( f 0 ) between the beam size (related to w 0 ) and the back aperture of the objective ( f sinθ max ).
We introduce the following shorthand notations: Expressions for focused Hermite-Gaussian and Laguerre-Gaussian modes can then be derived as [17] : Focused x− polarized HG 00 (Gaussian) mode: Focused x− polarized HG 10 mode: Focused x− polarized HG 01 mode: Focused x− polarized HG 20 mode: Focused linearly polarized LG lin 01 ('donut') mode: LG lin 01 = HG 10 + iHG 01 (9) Focused azimuthally polarized LG az 01 mode: Focused radially polarized LG rad 01 mode: For future reference in this article, we summarize the calculated field distributions near focus for these various cases in Figure 2. Phase distributions are presented without the propagation term exp(−ik ω z) in order to highlight the differences between the modes.

Calculation of the induced third-order non-linear polarization
In a medium characterized by its third-order nonlinear tensor χ (3) i jkl (r), the excitation field induces a polarization density described by: The χ (3) tensor of a homogeneous isotropic medium can be expressed as [18]: We can then express the nonlinear polarization induced by the exciting field E in cartesian coordinates as:

Propagation of the harmonic field
Finally, the harmonic field originating from all positions r in the focal region and propagated to a position R in the collection optics aperture can be expressed as [17,15]: where V spans the excitation volume and G FF is the far field Green's function: where R is the coordinate of a point in the far field (see Fig. 1) and I is the third-order identity tensor.
Emission diagrams can be analyzed from these equations by calculating the squared harmonic field |E FF (R)| 2 at different positions R. Alternatively, total THG power emitted in the forward (F-THG) or backward (B-THG) directions can be estimated by integrating |E FF (R)| 2 over the front aperture of an epicollecting or trans-collecting objective.
To simplify the analysis of the results, we assume no linear index mismatch and we neglect temporal aspects such as group velocity mismatch for ultrashort pulses.However we assume that the samples consist of normally dispersive media (which is usually the case in biological THG imaging), since dispersion plays a significant role in THG contrast formation: for example, negative dispersion in homogeneous gas samples can result in bulk emission [18].

Numerical implementation
Calculations are performed using Matlab.We typically discretize the focal volume over a 200 × 140 × 140 × (λ /40) grid, and evaluate the excitation field using quadrature algorithms.Unless otherwise stated, we use the following parameters: λ = 1.2μm,NA = 1.4 or 1.2, f 0 = 2, n ω = 1.5, n 3ω = 1.52.We note that incorporating positive dispersion in the model is numerically advantageous because smaller focal volumes can be considered, and calculations are generally less noise-sensitive than in the limit case of zero-dispersion.For a given sample/focal field combination, we calculate the projection of the forward-and backward-emission patterns on planes perpendicular to the optical axis located at Z = ±10cm.We choose to present projected far-field patterns rather than angular emission diagrams because they appeared to be more readable in the case of complex emission profiles.For the interface and slab sample geometries, we assume that the focal volume encompasses two homogeneous isotropic media with thirdorder nonlinear susceptibilities χ (3) 1 = 1 and χ (3) 2 = 0.This choice is motivated by the fact that, for excitation geometries where bulk THG emission is canceled by destructive interference, THG from an interface scales as |χ For periodic samples, we assume a sine-like variation χ (3) = 1 + sin(2πz/δ e)/2 along the optical axis.We then iterate for each beam shape and for various sample positions the calculation of emission patterns, F-THG and B-THG powers.Normalization is done by considering the same total intensity in the focal volume for every mode.

Vectorial aspect of THG microscopy with tightly focused beams
Although studies of THG by focused Gaussian beams generally neglect vectorial aspects, a general analysis of THG microscopy requires a priori a vector field model because high NA focusing does not preserve linear polarization.Furthermore it is seen from Eq.14 that the induced .Conversely if E z does not overlap with E x , only the E 3 x term will contribute to the TH signal.In particular, in the case of a tightly focused Gaussian beam with initial linear polarization the axial component near focus is important (see Fig. 3): Max(E z ) ≈ Max(E x )/3 for NA = 1.4.However in this case there is little overlap between E x and E z , so that E z contributes little to THG.Thus, a scalar approximation will usually work well for THG from simple interfaces excited by a focused linearly polarized Gaussian beam.However it will typically not be accurate for higher-order beam shapes or other input polarization patterns.Recalling that the phase distribution (including the Gouy shift) is generally different for the various field components [17,19], cross-terms may define different coherence lengths within the focal volume and affect the imaging properties.Furthermore, focused radially polarized beams typically exhibit strong axial components (see e.g.LG rad 01 mode in Fig. 3) [20, 1, 2] which give them original imaging properties.
We also point out that a well-described vectorial aspect of THG microscopy is the absence of THG signal when a HG 00 beam with circular polarization is focused on an interface between isotropic media [9].We begin our study of THG microscopy with non-standard beams by considering the simplest sample geometry, namely an interface between two media with different nonlinear susceptibil- (3)

THG imaging of XY interfaces with HG and LG beams
1 and χ (3) 2 .A well-known fact in the case of a XY-interface (i.e.perpendicular to the optical axis) excited by a tightly focused HG 00 beam is that a z-scan produces a symmetric Gaussian (or Lorentzian if the objective aperture is overfilled) curve peaking when the interface is at the focus.Our simulations predict similar behaviors for THG emission with all the higher-order HG and LG modes considered in this study (see Fig. 4 for characteristic examples), whether or not they possess cylindrical symmetry.The main difference lies in the widths of the z-scan responses, which are related to the axial extents of the corresponding excitation field distributions.The single-peaked nature of the z-scans reflects the fact that the modes considered here exhibit a single axial maximum.We note that emission is essentially forward-directed, as in the case of Gaussian excitation [21].However different modes produce different emission patterns.For example, on-axis harmonic emission is prevented with focused HG 01 modes due to the laterally antisymmetric nature of the focal phase distribution, as predicted for CARS emission from bulk media [3].
We also point out that the HG 00 case is qualitatively well-described by the paraxial approximation (not shown) even at high NA (in contrast with [15] where an inappropriate value of the confocal parameter was used).Paraxial approximation provides meaningful results for nonlinear processes involving Gaussian beams and simple geometries because over the region of highest intensity where most signal creation occurs, both the intensity and the Gouy phase shift (which varies almost linearly with z) are accurately approximated.Of course, the situation can be quite different with complex field distributions.More interesting is the case of an interface parallel to the optical axis excited by an asymmetric field distribution such as a focused HG 01 or HG 10 beam.When a Y Z interface is x−scanned across a focused HG 01 beam (Fig. 5(a)), the F-THG response exhibits a double peak reflecting the field distribution in the focal plane (see Fig. 5(c)), contrasting with the case of a focused HG 00 or HG 10 .Even more striking is the case of a XZ interface being y−scanned across a focused HG 10 beam (Fig. 5(b)).In this case the THG response exhibits a triple peak.The central peak results from the presence of a significant axially polarized component in the strongly focused x-polarized HG 10 field (I 030 term in Eq. 6 which is not present in the HG 01 case, see also Fig. 5(c)).This vectorial interpretation is corroborated by the double-peaked shape of the THG y−scan obtained when the axial component is omitted in the simulation (Fig. 5(b)).

THG imaging of XZ/YZ interfaces with focused HG beams
Since this situation is equivalent to that of a single half-space with susceptibility |χ (3) 2 |, maximum emission is obtained when one of the two main excitation peaks is incident on the interface.We now seek to take advantage of the laterally asymmetric nature of odd HG beams to probe sample orientation.We analyze the THG response obtained from HG 01 and HG 10 beams focused on an interface parallel to the optical axis as a function of the angle φ that it makes relative to the X axis (see Fig. 6(a)).As anticipated from the shape of the focal fields, THG emission is strongly modulated (>50%) as a function of the interface angle.This is because the geometry is roughly equivalent to that of a single-peaked excitation where the distance between the focal spot and the interface is proportional to the sine (or cosine) of the interface angle.Therefore, anti-correlated behaviors are predicted with HG 01 and HG 10 beams (Fig. 6(b)), and a combination of measurements with two such beam shapes provides sub-μm information about sample orientation.The corresponding far-field emission patterns exhibit even more subtle variations (Fig. 6(c)).When the two main peaks of the focal field lie on both sides of the interface (i.e.φ = 0 with HG 10 excitation or φ = 90 with HG 01 excitation), F-THG emission occurs along two off-axis lobes.These emission patterns are reminiscent of the case of Gaussian excitation of a XY interface [15], where the Gouy shift prevents on-axis phase matching and deflects the emission.However when both excitation peaks are incident on the interface (i.e.φ ≈ 90 with HG 01 or φ ≈ 0 with HG 10 ), the emission originates from two sources and interferences structure it into a non-trivial 8-lobes pattern.The relative intensities of these lobes are tightly related to the sample orientation (Fig. 6).See (Media1).HG 01 excitation (red disc) produce a double-peaked response as a function of thickness, corresponding to distinct emission patterns.The double peak behavior is blurred for HG 10 excitation (see text).The HG 20 case without dispersion is also presented for comparison (empty triangles).For all the modes considered here, the peak TH signal intensity is between 1.5 and 2 times higher than that obtained from a semi-infinite slab.focused higher-order beams are an active area of research [24,25,22], and so are the vectorial aspects of nonlinear microscopy [2,26].Studying focus-engineered THG is informative because third-harmonic generation can be obtained from simple (isotropic) sample geometries and is highly sensitive to the focal phase distribution.Therefore it is a convenient means to study the impact of focus engineering on phase-matching conditions.Moreover, THG microscopy is usually a non-spectroscopic, structure-sensitive [14] imaging technique that provides morphological information about unstained samples.The results presented here show the potential of focus-engineered THG microscopy to provide sub-resolution information about complex samples: angles and characteristic lengths in the 130 − 250nm range are reflected in the emission patterns and can be probed using simple ratiometric measurements.A perspective is to design pupil functions producing a targeted field distribution [27,28], axial field engineering (bottle beams [3], etc) being of particular relevance.Finally we note that the ideas explored in this article are generally transposable to other imaging modalities such as CARS and SHG.For example, CARS emission from an axial interface excited with a laterally asymmetric field distributions should exhibit a behavior similar that shown in Fig. 6(b).Also, focus engineering can be used to alter the spatial resonances resulting in efficient SHG and CARS emission in a particular direction, as in Fig. 8.It is anticipated that additional studies will explore the potential of focal field engineering for coherent nonlinear microscopy.

Fig. 2 .
Fig. 2. Distributions in the xy and xz planes of the focal field intensity and of the phase of the x-polarized component when relevant, for the modes described by Eqs.(5-11).Arrows indicate the direction of polarization in the xy plane for focused LG 01 beams.Intensity plots are normalized to their maximum values.Phase color table ranges from white (−π rad) to black (π rad).NA = 1.4,x, y ∈ [−1 1]μm, z ∈ [−2 2]μm.

Fig. 3 .
Fig. 3. Distribution of the different field polarization components and total intensity in the transverse focal plane for focused HG 00 and LG rad 01 beams.x, y ∈ [−1 1]μm.

Fig. 4 .
Fig. 4. F-THG during an axial scan through a xy interface with HG and LG beams.Curves are normalized by the factors indicated in the inset.

Fig. 5 .
Fig. 5. F-THG during lateral scans through interfaces parallel to the optical axis.(a) x−scan through a YZ interface.(b) y−scan through a XZ interface.The HG 10x curve (empty green triangles) is the behavior predicted when the z component of HG 10 is omitted.(c) Excitation field and intensity distribution in the focal plane for focused HG 01 and HG 10 beams.

Fig. 6 .
Fig. 6.Sensitivity to the orientation of an interface parallel to the optical axis using asymmetric excitation (HG 01 and HG 10 ).(a) Geometry of the sample and distribution of the excitation intensity in the focal plane.(b) Normalized F-THG signal as a function of interface angle φ for HG 00 (black squares), HG 01 (red discs) and HG 10 (green triangles).Normalization factors are indicated in the inset.The HG 10 /HG 01 signal ratio (purple stars) probes the interface orientation within the focal volume with good contrast.(c) Projected far-field emission patterns as a function of interface angle for HG 01 and HG 10 excitation.(Media1): TH emission patterns for HG 00 , HG 01 , and HG 10 excitation, as a function of interface orientation.Patterns are evaluated at z=+10cm over a 15×15cm area transverse to the optical axis, which corresponds to a detection NA of approximately 0.5.

Fig. 7 .
Fig.7.F-THG and B-THG from slabs of varying thicknesses using different beam shapes.(a) F-THG as a function of slab thickness indicates the forward coherence length associated with a particular field profile.The inset depicts the corresponding geometry.HG 20 excitation (blue triangles) results in larger forward coherence length than HG 00 (black squares).HG 01 excitation (red disc) produce a double-peaked response as a function of thickness, corresponding to distinct emission patterns.The double peak behavior is blurred for HG 10 excitation (see text).The HG 20 case without dispersion is also presented for comparison (empty triangles).For all the modes considered here, the peak TH signal intensity is between 1.5 and 2 times higher than that obtained from a semi-infinite slab.(b) Far-field emission patterns using HG 01 and HG 10 excitation, for different slab thicknesses.(c) B-THG as a function of slab thickness, according to the geometry depicted in the inset.Oscillation period indicates the backward coherence length.(d) On-axis phase distribution (without propagation term) for HG 00 and HG 20 modes with different NAs.
Fig.7.F-THG and B-THG from slabs of varying thicknesses using different beam shapes.(a) F-THG as a function of slab thickness indicates the forward coherence length associated with a particular field profile.The inset depicts the corresponding geometry.HG 20 excitation (blue triangles) results in larger forward coherence length than HG 00 (black squares).HG 01 excitation (red disc) produce a double-peaked response as a function of thickness, corresponding to distinct emission patterns.The double peak behavior is blurred for HG 10 excitation (see text).The HG 20 case without dispersion is also presented for comparison (empty triangles).For all the modes considered here, the peak TH signal intensity is between 1.5 and 2 times higher than that obtained from a semi-infinite slab.(b) Far-field emission patterns using HG 01 and HG 10 excitation, for different slab thicknesses.(c) B-THG as a function of slab thickness, according to the geometry depicted in the inset.Oscillation period indicates the backward coherence length.(d) On-axis phase distribution (without propagation term) for HG 00 and HG 20 modes with different NAs.
(3ω)can linearly depend on a particular field component.For example if E z is strong at a particular location near focus and spatially overlaps with E x , a cross-term proportional to E z E 2x will significantly contribute to P #98503 -$15.00USD Received 8 Jul 2008; revised 21 Aug 2008; accepted 21 Aug 2008; published 3 Sep 2008 (C) 2008 OSA 15 September 2008 / Vol.16, No. 19 / OPTICS EXPRESS 14708 nonlinear polarization P