## Abstract

The effects of the axial field components of a focused beam under high NA on the second harmonic generation (SHG) in collagen was examined using a vectorial approach. We find that with high NA, the cross-component terms that are most likely to have an effect on SHG will be E_{x}E_{x}, E_{x}E_{y}, E_{x}E_{z} and E_{z}E_{z} as a result of tight focusing. By considering the tensor and the presence of the other electric field components the possibility of different polarization states of the generated second harmonic as a result of the nonlinear susceptibility tensor making it possible to generate radially polarized modes with linearly polarized beams.

©2006 Optical Society of America

## 1. Introduction

There has been an increasing interest in the use of nonlinear optical processes for microscopy among which second harmonic generation (SHG) is one of them. Some advantages that SHG microscopy enjoys over other conventional microscopy include the use of intrinsic optical properties for imaging, the absence of photodamage, a highly localized interaction region and structural identification through the resolution of the tensor components of SHG emitting structures [1,2]. Currently, SHG microscopy has been applied to studying surfaces [3,4], microtubules [5], collagen [1,6], tumours [7], muscle [8] and in-vivo developmental biology [9].

As a nonlinear optical process, SHG is dependent on the optical properties of the material. This property is normally described by a nonlinear susceptibility tensor. For SHG, the tensor that describes its nonlinear optical properties is a second-order tensor with 27 separate elements, which results in an electric polarization that is dependent on the various components of the electric field [10]. Currently, most papers dealing with SHG consider the effects of the focused beam have been approximated using an effectively scalar theory [11,12,13] or one that only focuses on the dominant component of the illuminating field [14,15]. In reality all three field components are present since one does not, in practice, encounter an infinite plane wave and also because microscopy utilizes high numerical aperture lenses to focus the beam. It is under such circumstances where a vectorial theory is appropriate [16]. In this paper, we consider the effects of the various electric field components on the SGH from a hypothetical collagen fibril using the Richards and Wolf approach. We also find that the axial components of the electric field play a role in SHG in the case of collagen through the susceptibility tensor.

## 2. Theory

SHG is dependent on the nonlinear susceptibility tensor of the sample. This tensor is a tensor with 27 elements similar to the piezoelectric tensor. Under certain conditions the number of non-zero elements in this tensor reduces according to Kleinmann’s symmetry [10]. In the case of collagen, a very good second harmonic generator, the tensor has been described as a C_{6} tensor [1,2,6]. When the axis of symmetry coincides with the laboratory *z* axis we can write:

Equation (1) represents the case when the axis of symmetry is parallel to the propagation of the excitation field. Alternatively, we may have the axis lying parallel to the polarization of the excitation field and assuming our field is polarized in the *x* direction, gives:

$${P}_{y}^{\mathrm{SHG}}={2d}_{\mathrm{yyx}}{E}_{y}{E}_{x},\phantom{\rule{8.2em}{0ex}}$$

$${P}_{z}^{\mathrm{SHG}}={2d}_{\mathrm{zzx}}{E}_{z}{E}_{x}.\phantom{\rule{8.5em}{0ex}}$$

Previous values of *d _{ijk}* been found experimentally [2] and take the values

*d*=

_{xzz}*d*=1,

_{xyy}*d*:=0.09 and

_{xxx}*d*=

_{yyx}*d*=1.15. The equations indicate that SHG polarization can be induced through ‘cross-components’ such as E

_{zzx}_{i}E

_{j}(

*i*≠

*j*). For low NA focusing, the component of the electric field in the direction of polarization is strongest and applying that to Eq. (2), induces only |P

_{x}

^{SHG}| or |P

_{y}

^{SHG}| or both depending on the orientation of the polarization with respect to the collagen axis. When the collagen axis is along

*z*the only SHG polarization induced will be |P

_{z}

^{SHG}| The resultant SHG signal from such a configuration is similar to a dipole radiating along the

*z*axis and will require a high NA collector for transmission type detection. This is perhaps why it has been reported that there is no apparent SHG signal when the orientation is such [17]. One interesting aspect of this tensor is the fact that it is possible to elicit various modes of polarization different from the excitation beam. For example, the P

_{z}

^{SHG}term in Eq. (2) implies a radially polarized mode obtainable with a mixing of the E

_{x}E

_{z}term. Likewise if we were to consider the case of

*P*=

_{z}^{SHG}*d*only we will have a radially polarized mode generated by a linearly polarized beam.

_{zxx}E_{x}E_{x}For high NA focusing, the electric field at the focus consists of transverse as well as axial components. Although weak, the transverse and axial components can have an effect when a tightly focused beam is used in nonlinear optical microscopy due to the tensorial nature of SHG. According to Richards and Wolf [18], the electric field components in the focal region are:

$${E}_{y}\left(u,v\right)=-i{I}_{2}\phantom{\rule{.2em}{0ex}}\mathrm{sin}2\varphi ,\phantom{\rule{2.8em}{0ex}}$$

$${E}_{z}\left(u,v\right)=-2{I}_{1}\phantom{\rule{.2em}{0ex}}\mathrm{cos}\varphi .\phantom{\rule{3em}{0ex}}$$

where

$${I}_{1}\left(u,v\right)={\int}_{0}^{\alpha}{\mathrm{cos}}^{1/2}\theta \phantom{\rule{.2em}{0ex}}{\mathrm{sin}}^{2}\theta {J}_{1}(\mathit{kr}\phantom{\rule{.2em}{0ex}}\mathrm{sin}\theta )\mathrm{exp}\left(\mathrm{ikz}\mathrm{cos}\theta \right)d\theta ,\phantom{\rule{4em}{0ex}}$$

$${I}_{2}\left(u,v\right)={\int}_{0}^{\alpha}{\mathrm{cos}}^{1/2}\theta \mathrm{sin}\theta \left(1-\mathrm{cos}\theta \right){J}_{2}(\mathit{kr}\phantom{\rule{.2em}{0ex}}\mathrm{sin}\theta )\mathrm{exp}\left(\mathrm{ikz}\phantom{\rule{.2em}{0ex}}\mathrm{cos}\theta \right)d\theta ,$$

where *r* = (*x*
^{2} + *y*
^{2})^{1/2}, α is the aperture half-angle and *J _{n}*( . ) is a Bessel function of the first type and order

*n*. In the high NA case (here we assume a NA 1.4 objective) the ratios of the electric field components |

*E*|:|

_{x}*E*|:|

_{y}*E*| are approximately 1:0.1:0.3 respectively. In Table 1 we find that the E

_{z}_{i}E

_{j}terms that are most likely to have an effect on SHG will be the cases where we have E

_{x}E

_{x}, E

_{x}E

_{y}, E

_{x}E

_{z}and E

_{z}E

_{z}. The other cross-component terms will have minimal effects since they are two orders of magnitude less than E

_{x}E

_{x}.

In order to solve for the far-field radiation of SHG, a Green’s function approach may be used [14] or a phased array approach [11,12]. We use the former approach and the interested reader is referred to the following references for further details [14,19]. We present only the end equation that is used for calculating the radiation pattern of SHG in the far-field (in spherical coordinates) [14]:

where **R̂** is the unit vector in the direction of the observer and *k* is now 2*k* due to frequency doubling, Θ is the polar angle of the observation point and Φ is the azimuthal angle of the observation point. For our calculations we took *λ*=1000nm, *n*=1.5, NA=1.4 and assumed that the refractive index of the lens and sample were the same.

In this work we assume, for simplicity, a hypothetical fibril of collagen that possesses the following properties. One is that the fibril is composed of many infinitesimally small subunits each possessing the C_{6} symmetry. These idealized subunits may then be arrayed in different orientations (with respect to their axis of symmetry) in ideal 1-dimensional lines or planes. Fig. 1 illustrates the geometry and orientation adopted in this theoretical study.

## 3. Results

We calculated the induced SHG polarization by assuming that the C_{6} structure was valid in the focal plane and applying Eqs. (1)–(4). Different orientations of the fibre axis result in different induced SHG polarization as can be seen from Fig. 2, which can be visualized as having the layout of Fig. 1(a) except that the linear array of subunits has been extended over the focal plane.

Although the field in the *x* direction is strongest for a focused beam linearly polarized along *x*, we find that it is not the |P_{x}
^{SHG}| that is dominant but rather the |P_{z}
^{SHG}|. The |P_{x}
^{SHG}| term is also significant as it is half as strong as the |P_{z}
^{SHG}| term when the axis is parallel to *z* axis. This |P_{x}
^{SHG}| component is generated by the product of the *x* and *z* components of the incident field. The |P_{x}
^{SHG}| component is nearly 5 times weaker than the |P_{z}
^{SHG}| term when the axis of the fibre is parallel to the *x* axis as seen in Fig. 3. This is explained by looking at Eqs. (1) and (2) and noting the values of the associated *d* coefficients. The dominant |P_{z}
^{SHG}| component is now generated by the product of the *x* and *z* components of the incident field.

Given that the induced SHG polarization is different from the applied field and exhibits different orientations, we calculated the SHG radiation patterns from a hypothetical collagen fibril of zero thickness (Fig. 1). We thus had four cases *viz*. orientation of the axis of symmetry in the *x* direction and “extending” in the *x* or *z* directions and orientation of the axis of symmetry in the *z* and ‘extending’ in the *x* or *z* directions. We used Eq. (5) to calculate the final radiation pattern in the Θ direction.

#### 3.1 Extension in the x direction

When the axis of symmetry is parallel to *z* and the extension is in *x* we would expect short lengths of the fibril to radiate in a manner similar to that of a dipole along the *z* axis. This is because for ranges very close to zero the only induced SHG polarization is in the axial direction (Fig. 2(a) and Fig. 2(c)). A very short fibril length will thus approximate a dipole along the *z* axis and will radiate as shown in Fig. 4(a). Similarly for the case when the axis of symmetry is parallel to *x*, Fig. 3 indicates that a very short length of such a layout will result in a radiation pattern similar to a dipole along the *x* axis as seen in Fig. 5(a). Fig. 5(a) is, however, not rotationally symmetric. This is due to the fact that in Eq. (5), the P_{z}
^{SHG} term is dependent only on Θ and not on Φ.

We expect the radiation pattern to change as the length varies simply by looking at the relative strengths of the induced SHG polarization in Fig. 2(a) and Fig. 2(c). Figures 4 and Fig. 5 give an idea of how the SHG radiation pattern ‘evolves’ as the fibril gets longer in the *x* direction. The need for phase matching also creates a radiation pattern that has lobes radiating away from the *z* axis [12,14]. Since the fibril is only along the *x* axis, it is possible for a b-SHG (backward propagating SHG) as seen in Fig. 4 and Fig. 5 (where dotted lines denote the *z* = 0 level) illustrate the relative strengths of the b-SHG to the f-SHG.

#### 3.2 Extension in the z direction

For the case when the extension of the collagen is in the *z* direction the radiation pattern for a short length is determined by the induced SHG polarization that is non-zero near the origin as previously. More interestingly is the lack of b-SHG as seen in Fig. 6 and Fig. 7. This is primarily due to coherence effects. In the axial direction, the phase changes and a Gouy phase shift occurs. Each SHG-producing element along the axial direction thus radiates equally backwards and forwards. The relative phase shift between each element implies that the b-SHG will interfere destructively and thus bulk objects typically have a strong f-SHG signal only. Figure 7 has also been calculated by calculated previously based on inhomgenous scatterers in membranes [11] and hemispherical interfaces [14].

Thus there exists a certain amount of observable SHG in reflection for a hypothetical line of SHG emitting elements along the transverse axes. The presence of a ‘thickness’ or extension in the axial direction means that the SHG radiates in a forward direction due to coherence effects. Phase matching considerations also indicate that SHG will rarely propagate along the *z* axis but will, instead, propagate at an angle away from the *z* axis. In Fig. 4, the angle of propagation away from the *z* axis is approximately 20°, 29° for figure 5, 47° for Fig. 6 and 40° for Fig. 7. A high NA condenser is thus necessary to collect transmitted SHG light more efficiently.

## 4. Discussion

We have examined the far-field SHG radiation patterns for a variety of cases using a full vectorial approach to describe a linearly polarized excitation field. We find that for high NA focusing, the axial component of the electric field plays a role in the generation of SHG through the susceptibility tensor in the case of collagen. We find that the E_{i}E_{j} terms that are most likely to have an effect on SHG will be the cases where we have E_{x}E_{x}, E_{x}E_{y}, E_{x}E_{z} and E_{z}E_{z}. A further point of interest is the possibility of different polarization states of the generated second harmonic as a result of the nonlinear susceptibility tensor. This possibility of generating radially polarized modes with a linearly polarized beams is of interest. There have also been several reports on aperturing on harmonic generation [20,21]. Effects of aperturing could be made independent of the order of the harmonic as well as the gas species and that the harmonic efficiencies in terms of aperture size was peaked indicating a certain aperture size which is optimum [20]. It was also demonstrated that higher NA objectives improved the resolution and varied the intensity of the SHG [21]. The reason assigned to this was that in certain cases, the orientation of the crystal was such that it favoured SHG with the excitation beam parallel to the axis of propagation. With an increasing NA, the spectrum of waves was spread out over a wider angle and hence less energy was available for SHG in the direction optimal for phase matching. Similarly, cases where the intensity was observed to increase with NA could be explained in a similar fashion. It is clear that as the NA influences the strengths of the various components of the electric field, further study of the effects of NA on harmonic generation incorporating a vectorial theory of the electric field at the focus is advantageous for nonlinear microscopy techniques. Applications for such techniques will include identifying the three-dimensional orientation of molecules and tracking changes in the susceptibility tensor which will have potential interest in cancer and various diseases caused by modifications to the tissue/cellular structure.

## References and Links

**1. **P. Stoller, B. M. Kim, A. M. Rubenchik, K. M. Reiser, and L. B. Da Silva, “Polarization-dependent optical second-harmonic imaging of rat-tail tendon,” J. Biomed. Opt. **7**, 205–214 (2002) [CrossRef] [PubMed]

**2. **S. W. Chu, S. Y. Chen, G. W. Chern, T. H. Tsai, Y. C. Chen, B. L. Lin, and C. K. Sun, “Studies of χ^{(2)}/χ^{(3)} tensors in submicron-scaled bio-tissues by polarization harmonics optical microscopy,” Biophys. J. **86**, 3914–3922 (2004) [CrossRef] [PubMed]

**3. **R. Gauderon, P. B. Lukins, and C. J. R. Sheppard, “Three-dimensional second-harmonic generation imaging with femtosecond laser pulses,” Opt. Lett. **23**, 1209–1211 (1998) [CrossRef]

**4. **C.K. Sun, S.W. Chu, S.P. Tai, S. Keller, A. Abare, U. K. Mishra, and S. P. DenBaars, “Mapping piezoelectric-field distribution in Gallium Nitride with scanning second-harmonic generation microscopy,” Scanning **23**, 182–192 (2001) [CrossRef] [PubMed]

**5. **P. J. Campagnola, A. C. Millard, M. Terasaki, P.E. Hoppe, C. J. Malone, and W. A. Mohler, “Three-dimensional high-resolution second-harmonic generation imaging of endogenous structural proteins in biological tissues,” Biophys. J. **81**, 493–508 (2002) [CrossRef]

**6. **S. Roth and I. Freund, “Second harmonic generation in collagen,” J. Chem. Phy. **70**, 1637–1643 (1979) [CrossRef]

**7. **E. Brown, T. McKee, E. diTomaso, A. Pluen, B. Seed, Y. Boucher, and R. K. Jain, “Dynamic imaging of collagen and its modulation in tumours in vivo using second-harmonic generation,” Nature Med. **9**, 796–800 (2003) [CrossRef] [PubMed]

**8. **M. Both, M. Vogel, O. Friedrich, F. von Wegner, T. Kunsting, R. H. A. Fink, and D. Uttenweiler, “Second harmonic imaging of intrinsic signals in muscle fibres in situ,” J. Biomed. Opt. **9**, 882–892 (2004) [CrossRef] [PubMed]

**9. **S. W. Chu, S. Y. Chen, T. H. Tsai, T. M. Liu, C. Y. Lin, H. J. Tsai, and C. K. Sun, “In vivo developmental biology study using non-invasive multi-harmonic generation microscopy,” Opt. Express **11**, 3093–3099 (2003). [CrossRef] [PubMed]

**10. **D. A. Kleinmann, “Nonlinear dielectric polarization in optical media,” Phys. Rev. **126**, 1977–1979 (1962) [CrossRef]

**11. **L. Moreaux, O. Sandre, and J. Mertz, “Membrane imaging by second-harmonic generation microscopy,” J. Opt. Soc. Am. B **17**, 1685–1694 (2000) [CrossRef]

**12. **J. Mertz and L. Moreaux, “Second-harmonic generation by focused excitation of inhomogeneously distributed scatterers,” Opt. Commun. **196**, 325–330 (2001) [CrossRef]

**13. **R. W. Boyd, “*Nonlinear Optics*,” 2nd Ed(Academic Press, Amsterdam, 2003)

**14. **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) [CrossRef]

**15. **R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophys. J. **88**, 1377–1386 (2005) [CrossRef]

**16. **A.A. Asatryan, C. J. R. Sheppard, and C. M. de Sterke, “Vector treatment of second harmonic generation produced by tightly focused vignetted Gaussian beams,” J. Opt. Soc. Am. B **21**, 2206–2212 (2004) [CrossRef]

**17. **T. Yasui, Y. Tohno, and T. Araki, “Characterization of collagen orientation in human dermis by two-dimensional second-harmonic-generation polarimetry,” J. Biomed. Opt. **9**, 259–264 (2004) [CrossRef] [PubMed]

**18. **B. Richards and E. Wolf, “Electromagnetic diffraction in optical systems. II. Structure of the image field in an aplanatic system,” Proc. Roy. Soc. Lond. Ser. A , **253**, 358–379 (1959) [CrossRef]

**19. **C. T. Tai, *Dyadic Green’s Functions in Electromagnetic Theory*, (Intext Educational Publishers, 1971)

**20. **S. Kazamias, F. Weihe, D. Douillet, C. Valentin, T. Planchin, S. Sebban, G. Grillon, F. Auge, D. Huilin, and Ph. Baclou, “High order harmonic generation optimization with an apertured laser beam,” Eur. Phys. J. D **21**, 353–359 (2002) [CrossRef]

**21. **R. Gauderon, P. B. Lukins, and C. J. R. Sheppard, “Optimization of second-harmonic generation,” Micron **32**, 691–700 (2001) [CrossRef] [PubMed]