Translator Disclaimer
1 July 2009 Theoretical simulation study of linearly polarized light on microscopic second-harmonic generation in collagen type I
Author Affiliations +
Abstract
A theoretical model is established for dealing with second-harmonic generation (SHG) in type I collagen excited by linearly polarized light focused by a microscope. With this model, the effects of the polarization angle α, numerical aperture (NA), as well as the ratio of hyperpolarizability ρ=β xxx/β xyy on SHG emission have been investigated. Simulation results reveal that SHG emission power changes periodically as α. The use of lower NA leads to weaker SHG emission but is more concentrated in two closer lobes, whereas more distributed emission in two detached lobes appear at higher NA. As the introduction of polarization direction, which is not along with the fiber axis (alpha ≠ 0 deg), one more element β xyy is valid in our case than β xxx alone, while their ratio ρ plays a very important role for collagen features characterization. SHG emission with ρ shows complicated modality that SHG emission is different at different α and not symmetric at ± ρ except at α=0 deg, suggesting the important impact of polarization working on ρ for SHG emission. Our theoretical simulation results provide useful clues for experimental study of microscopic SHG emission in collagen excited by linearly polarized beam.

1.

Introduction

Second-harmonic generation (SHG) microscopy, based on radiative SHG under a focused beam has attracted more and more interest because it has been proved to be a powerful tool for visualization of cell and tissue structures and functions with microresolution.1, 2

Collagen has the particular relevance of microscopic SHG applications onto biological imaging.3, 4 Collagen is the most abundant intrinsic structural protein in tissues. Collagen has a highly crystalline triple-helix structure, which is not centrosymmetric,5, 6 and the molecules are organized on the scale of the wavelength of light, which suggests that collagen is a promising biomedical material capable of SHG. Meanwhile, collagen possesses first hyperpolarizabilities β large enough for practical SHG.7

SHG in collagen fiber has been widely investigated, experimentally and theoretically. However, although many experiments on SHG in collagen have been accomplished under the microscope, theoretically, SHG in collagen under a focused beam has far from thorough investigation, such as the polarization effect. Thus far, polarization effect on SHG under focused beam has been investigated only on uniformly distributed dye dipoles on a planar cell membrane.8 In collagen, such studies have been concentrated on the excitation by collimated light rather than focused light.9, 10 Under a focused beam, only in a very special condition, such as the excitation beam, which polarizes parallel to the fibril axis, has been dealt with.11 SHG is a nonlinear optical phenomenon involving coherent photon scattering. Because of the coherent nature of SHG, the radiation of SHG is usually highly directional and depends critically on the emission source, such as its spatial extent, polarization of the excitation electric field.12, 13 Also, we realize that collagen fiber has its own physical properties compared to the uniformly distributed dipole on surface, and in practice, the polarization of the excitation light relative to the collagen fiber axis is indeterminable because of the unknown sample-laid position. It is therefore desirable for the theoretical exploration of SHG emission under the polarized focused beam on collagen fibrils.

Moreover, it is well known that the structure of collagen plays an important role in various pathological processes and diseases, such as wound healing, malignancy, or development,14 and SHG microscopy can be accomplished in vivo. Thus, this theoretical analysis exhibits promise for biomedical assessment of tissue structure, especially in pathological processes where the collagen is being actively remodeled.

Therefore, in this paper, we aim to explore the effects of polarized focused excitation beam on SHG in collagen. A single fibril is applied as the collagen model in this paper.

2.

Theory of SHG Field Propagation under a Linearly Polarized Focused Beam

We consider here an illumination beam focusing on to the collagen fiber through a microscopic objective with the numerical aperture of NA=nωsinϕ . The light is linearly polarized in the x-y plane (Eω) with the angle of α from the x direction and propagates in the z direction with an angular frequency of ω . Collagen fiber is assumed along the x direction. The fibrils within the fiber are generally aligned to the fiber axis, and the nonlinear scatterers (dipoles) are distributed uniformly within the collagen fibrils. The diagram of coordinate system for calculating the SHG field at an observation point (r,θ,φ) is shown in Fig. 1 .

Fig. 1

Schematic diagram defining SHG emission from collagen under linearly polarized focused beam. Collagen fiber is assumed along the x direction. The fibrils within the fiber are generally aligned to the fiber axis. The light is linearly polarized in x-y plane (Eω) with the angle of α from x direction and propagates in the z direction. SHG propagates in the emission plane, defined by (θ,φ) .

044016_1_011904jbo1.jpg

2.1.

Electric-field Eω Profile of the Linearly Polarized Focused Beam

The local field distribution of excitation beam within the focus is very important for the determination of the propagation direction of the SHG field. For nonlinear applications, a linearly polarized, z -directed beam at the focus can be well approximated by a three-dimensional Gaussian profile as follows:15

1.

Eq. 1a

Eω(x,y,z)=iEω(0)exp(x2+y2wxy2z2wz2+iξkωz)ε̂,

Eq. 1b

Eω(0)=Eω(0,0,0),
where Eω(x,y,z) is the complex amplitude of the driving field at the point (x,y,z) , ε̂ represents its polarization direction ( εz=0 in our case). kω=2πnωλω is the wave vector in a specimen with a refractive index of nω . ξ is the wave-vector reduction factor due to focusing that accounts for a reduction in axial momentum by conversion to lateral momentum components within the focus. wxy and wz are the 1e radii of the focal ellipse in the lateral and axial directions, respectively.1

2.

Eq. 2a

wxy={0.320λωNANA0.70.325λωNA0.91NA> 0.7},

Eq. 2b

wz=0.532λω[1nωnω2NA2].

2.2.

Induced Dipole Moment of SHG

The total dipole moment μ induced by the field Eω is often written as a Taylor series16

Eq. 3

μ=μ0+αEω+12!βEω2+13!γEω3+.
Here, μ0 is the permanent dipole moment, the frequency dependent α , β , and γ are the linear polarizability, and the first and second hyperpolarizability, respectively. SHG is related to the third term, and the total local dipole moment per unit scatterer concentration can be expressed in more detail as8

Eq. 4

μ2ω,i(x,y,z)=12Eω2(x,y,z)j,kβijkε̂jε̂k,
where β symbolizes a local ensemble average of the fibril hyperpolarizabilities. Under the assumption of Kleinman and cylindrical symmetry (please see the Appendix { label needed for app[@id='xA'] }), dipole moment introduced by electric field Eω along the x and y direction is respectively to be

5.

Eq. 5a

μ2ω,x(x,y,z)=12Eω2(x,y,z)(βxxxεxεx+βxyyεyεy),

Eq. 5b

μ2ω,y(x,y,z)=12Eω2(x,y,z)(βxyyεxεy+βxyyεxεy).
In particular, for a linearly polarized beam at an angle α with respect to the collagen fiber produces the SHG dipole moment as follows:

Eq. 6

μ2ω(x,y,z)=(μ2ω,x(x,y,z)μ2ω,y(x,y,z)μ2ω,z(x,y,z))=12Eω2(x,y,z)(βxxxεxεx+βxyyεyεyβxyyεxεy+βxyyεyεx0)=12Eω2(x,y,z)(βxxxcos2α+βxyysin2αβxyysin2α0)=12Eω2(x,y,z)β.

2.3.

Induced Electric Field E2ω of SHG

Following the achievement of the induced SHG dipole moment, we proceed to calculate the second-harmonic far field radiated in the direction (θ,φ) based on the coordinate system illustrated in Fig. 1. Far from the dipole, where it is called the wave or radiation zone, the configuration of E2ω is as follows:17

Eq. 7

E2ω(ψ)=μ2ωω2πε0c2rsin(ψ)exp(ik2ωr)ψ̂,
where ψ represents the angle between the x -axis and the emission direction r of SHG. sinψ=(sin2θsin2φ+cos2θ)12 illustrates the projection relationship between the direction of emission SHG dipole moment μ2ω and the excitation electric field E2ω of SHG. We define η=ω2πε0c2 , where ε0 is the spree-space permittivity, c is the speed of light, and ω is the frequency of fundamental beam.

The total radiated second-harmonic signal in the far field is the integrated one from all scatterers that have the dipole volume density defined by NV . Presented in polar coordinates, E2ω is described as follows:

Eq. 8

E2ω(θ,φ)=ηNVrsinψμ2ω(x,y,z)exp[ik2ω(zcosθ+ysinθsinφ+xsinθcosφ)]dxdydz=ηNVr(sin2θsin2φ+cos2θ)1212βEω(0)2exp(2x2+y2wxy22z2wz+2iξkωz)×exp[ik2ω(zcosθ+ysinθsinφ+xsinθcosφ)]dxdydz.
(Derivation from Eq. 8, 9, see Appendix { label needed for app[@id='xB'] }).

Eq. 9

E2ω(θ,φ)=(π2)3NVwxy2wzηr(sin2θsin2φ+cos2θ)12μ2ω(0)×exp{k2ω28[wxy2(sinθcosφ)2+wxy2(sinθsinφ)2+wz2(cosθξ)2]}=NE2ω(0)A(θ,φ),
where the introduced parameters are
μ2ω(0)=(μ2ω,x(0)μ2ω,y(0)μ2ω,z(0))=12Eω2(0,0,0)β=12Eω(0)2β,
which represents the contribution from the local induced polarization per unit volume density to the radiated electric field at the flocal center only [indicated by the superscript (0)].

Eq. 10

N=(π2)3wxy2wzNV

Eq. 11

E2ω(0)=ηr(sin2θsin2φ+cos2θ)12μ2ω(0)

Eq. 12

A=exp{k2ω28[wxy2(sinθcosφ)2+wxy2(sinθsinφ)2+wz2(cosθξ)2]}=exp{k2ω28[wxy2sin2θ+wz2(cosθξ)2]}.
A is the angular modulation term.

N=VNV , and V=(π2)3wxy2wz is the active SHG volume.18 N defines as the effective total number of dipoles that contribute to the generation of second-harmonic light.

2.4.

Emission Direction of Electric Field of SHG

Emission direction is confined by the emission plane (θ,φ) . Electric field of SHG is normal to the observed direction r . The radiated SHG with polarization components parallel (p) to and perpendicular (s) to this plane can be derived as follows:

Eq. 13

E2ω(0)(θ,φ)=(E2ω(0)p(θ,φ)E2ω(0)s(θ,φ))=ηrMμ2ω(0).
M is the projection matrix, which permutes coordinate (x,y,z) to (θ,φ) . It is defined by

Eq. 14

M(θφ)=(cosθcosφcosθsinφsinθsinφcosφ0).
On the basis of Eq. 9, E2ωp(θ,φ) and E2ωs(θ,φ) is then derived from E2ω(0)(θ,φ) . E2ωp is the amplitude of E2ω in the θ direction and E2ωs is the amplitude of E2ω in the φ direction.8

However, in practice, the measurements of components of p and s directions of SHG are difficult because they are not fixed as the emission plane changes. On the contrary, SHG are more easily to be detected along the parallel (∥) and perpendicular (⊥) direction to the excitation polarization (α) . Therefore, we further transform the (p,s) components to be (∥,⊥) components through the following relationship:

Eq. 15

(E2ωE2ω)=[cos(φα)sin(φα)sin(φα)cos(φα)](E2ωpE2ωs).

2.5.

Emission Power of SHG

On the basis of the above derivation, the relationship of total power distribution of SHG with the electric-field strength presented in Eq. 9 can be further described to be (see Appendix { label needed for app[@id='xC'] })

16.

P2ω(θ,φ)=18n2ωε0cη2A2N2Eω(0)4,
(1sin2θcos2φ)(βxxxcos2α+βxyysin2αβxyysin2α0)2.
By integrating Eq. 16 in solid angle, total power of SHG is derived, which only depends on the polarized direction α

Eq. 17

P2ω=12n2ωε0cη2N2Θ(μ2ω,x(0)2+μ2ω,y(0)2+μ2ω,z(0)2)=18n2ωε0cη2N2ΘEω(0)4[(βxxxcos2α+βxyysin2α)2+(βxyysin2α)2],
where
Θ=A2(θ,φ)(1sin2θcos2φ)sinθdθdφ,

Accordingly, the power distribution in parallel (∥) and perpendicular (⊥) components are as follows (see Appendix { label needed for app[@id='xD'] }):

18.

Eq. 18a

P2ω(θ,φ)=12n2ωε0cr2E2ω(θ,φ)2=18n2ωε0cη2N2A2(θ,φ)Eω(0)4{cos(φα)[cosθcosφ(βxxxcos2α+βxyysin2α)+βxyycosθsinφsin2α]sin(φα)[sinφ(βxxxcos2α+βxyysin2α)+βxyycosφsin2α]}2,

Eq. 18b

P2ω(θ,φ)=12n2ωε0cr2E2ω(θ,φ)2=18n2ωε0cη2N2A2(θ,φ)Eω(0)4{sin(φα)[cosθcosφ(βxxxcos2α+βxyysin2α)+βxyycosθsinφsin2α]+cos(φα)[sinφ(βxxxcos2α+βxyysin2α)+βxyycosφsin2α]}2.

SHG power in these two directions (∥,⊥), respectively, can be calculated by integrating in solid angle as follows:

19.

Eq. 19a

P2ω=P2ω(θ,φ)sinθdθdφ,

Eq. 19b

P2ω=P2ω(θ,φ)sinθdθdφ.

3.

Simulation Study of Parameters of Linearly Polarized Focused Beam on SHG Emission

3.1.

Effects of Illumination Polarization on SHG Emission Power under Focused Beam

In this section, we aim to explore the effects of illumination polarization on SHG emission power under focused beam. From the theoretical expression of total, parallel, and perpendicular SHG power in Eqs. 17, 19a, 19b, we note that there remain two valid hyperpolarizability elements βxxx and βxyy . βxyy is added in this situation due to the introduce of the illumination polarization α , which is not consistent with the fiber axis (α=0deg) . ρ=βxxxβxyy has been evaluated to vary between 3 to 3,19 in this section, we take ρ to be 2.6 (type I collagen in rat liver) and further assume βxyy=1 to simplify the model. The effect of ρ=βxxxβxyy is to be investigated in more detail in Sec. 3.3.

Fig. 2 shows the total P2ω , parallel P2ω , and perpendicular P2ω SHG emission power from collagen fiber as the change of polarization angle α with respect to fiber axis. All data are normalized by (18)n2ωε0cη2N2Eω(0)4=1 . We can see from Fig. 2 that P2ω , P2ω , and P2ω SHG emission power change as the polarization angle α periodically. When polarization direction locates at α=0 , 180, 360deg , indicating that the polarization of the excited beam has the same direction as the fibrils axis, the total power P2ω , and parallel P2ω power component achieve their maximum value, while the perpendicular P2ω power component has its lowest value. On the contrary, we find that when the polarization of excited beam is perpendicular on the fibril axis ( α=90 , 270deg ) the total power P2ω and parallel P2ω power component have their lowest value and perpendicular P2ω component is the reverse. P2ω and P2ω power component keep the similar configuration in amplitude and P2ω stands predominance in contribution to the SHG total power.

Fig. 2

Effects of polarization angle α on total P2ω , parallel P2ω and perpendicular P2ω components of SHG emission power.

044016_1_011904jbo2.jpg

Parallel P2ω and perpendicular P2ω power distribution along (θ,φ) , using Eqs. 18a, 18b, at polarization angles α=0 , α=30 , α=60 , and α=90deg , respectively have been demonstrated in three-dimension in Fig. 3 . Basically, the parallel and perpendicular distribution of SHG emission possess two types of patterns. Under α=0 , α=30 , and α=60deg conditions, these three parallel distribution has similar patterns, and the perpendicular one has another similar patterns; however, under the α=90deg situation, the distribution shows reverse patterns in parallel and perpendicular components. As the polarization angle increases from α=0 to α=90deg , SHG parallel emission power in the z -axis (θ=0deg) decreases, while the perpendicular one increases. Although the emission distribution pattern of the parallel component at α=90deg and that of perpendicular components at α=0deg are very similar, their amplitude is still slight different, which indicates that although the relative direction of measurement direction and fibril axis is the same as the above two conditions (which is perpendicular to the fiber axis), the polarization direction relative to the fiber axis has an important influence on the SHG emission power.

Fig. 3

Effects of polarization angle α on SHG emission power distribution (θ,φ) of parallel P2ω and perpendicular P2ω components.

044016_1_011904jbo3.jpg

3.2.

Effects of the Numerical Aperture (NA) on SHG Emission

The most important parameter to characterize the focused beam is the numerical aperture, NA=nωsinϕ , where ϕ represents the incident angle of the illumination light (please refer to Fig. 1).

The focus of the light causes the phase retardation near the focus center along the axial direction. This phase retardation is called Guoy phase shift and can be represented by a wave-vector reduction factor ξ . The change of NA will induce the change of ξ as approximately by ξcos(ϕ2) (when NA<1.2 at most8).

In our single-fibril model, the phase match term, including the Guoy phase shift and the dispersion due to the refraction index change between the incident (nω) and second harmonic frequency (n2ω) , can be expressed as

Eq. 20

Δk=k2ω2ξkω(kω=2πnωλωandk2ω=2πn2ωλ2ω).

ξ appears in the axial direction of the focal ellipse (wz) exclusively. Momentum conservation along the axial direction then forces the SHG (k2ω) to propagate at an off-axis angle θ . This angle θ is θ=cos1[ξ(nωn2ω)] to satisfy the phase-match condition, k2ωcosθ2ξkω=0 , thus achieves the max SHG emission intensity. Figure 4a illustrates the relationship between the phase match terms. The solid arrow stands for the condition of phase match, and the dash arrow represents for condition of the phase mismatch. ξ=ξnωn2ω has been introduced in Eq. 33c in Appendix { label needed for app[@id='xB'] } thus θ=cos1(ξ) . If nω=n2ω , where dispersion effect has been ignored, then θ=cos1(ξ) .

Fig. 4

Phase-match condition makes SHG propagate in the direction of k2ω , which deviates an angle of θ with respect to incident beam kω due to the Guoy phase shift (ξ) and dispersion. (b) Effects of NA on SHG electrical radiation pattern (on the right is the enlargement part from the square center of the image on the left).

044016_1_011904jbo4.jpg

In this paper, the refractive index of SHG and fundamental excitation light in collagen are correspondingly taken to be n2ω=1.503 , nω=1.453 based on the research results from Bashkatov.20 To demonstrate the influence of NA under the consideration of dispersion, SHG emissions under NA=0.8 and NA=0.2 has been compared as those shown in Fig. 4b. We note that due to the influence of dispersion, no matter at high NA or low NA, SHG emission presents two off-axis lobes. However, at higher NA, the SHG emission lobes are much more distributed and detached, the max value of SHG emission appears at the angle of θ=27deg . On the other hand, at lower NA, the SHG emission concentrates along small angle θ=16deg and has much weaker intensity.

To investigate how NA under polarized beam influences the total power, parallel, and perpendicular components, P2ω , P2ω , and P2ω [all data are normalized by (18)n2ωε0cη2N2Eω(0)4=1 ] along NA, which vary from 0–1.0 under polarization angle α=0,30,60,90deg is demonstrated in Fig. 5 . In general, we note that the total emission power of SHG increases as the increase of NA, but the change does not show linearity. As the increase of NA, SHG parallel emission at α=0,30,60deg has the similar trend of change but different degree and at α=90deg , there is no change. Whereas SHG perpendicular emission also increases at different degree when α=30,60,90deg and the slight increase at α=0deg can be ignored. As the polarization angle α changes from 0to90deg , the total and parallel components of SHG emission power decrease, while the perpendicular one increases.

Fig. 5

Effects of NA on total P2ω , parallel P2ω , and perpendicular P2ω components of SHG emission power under different polarization angle α .

044016_1_011904jbo5.jpg

3.3.

Effects of the Ratio ρ=βxxxβxyy under Different Polarization Direction

β is the hyperpolarizability of the collagen fibrils which is related to the electronic transition in the material, thus reflects its biological structure and chemical features.21, 22 In our specific collagen model with the polarized illumination beam, two elements βxxx , βxyy remain and have the contribution to SHG. The ratio ρ=βxxxβxyy plays very important role in characterizing the collagen features. It has been proved that ρ increases with the increasing of age.9, 11 As mentioned above, ρ=βxxxβxyy has been predicted to vary between 3 to 3.19 The minus value of ρ suggests the polarity of the fibrils reverse from the excited electrical field. The study of the effect of ρ on SHG emission power may provide us significant information for tissue characterization.

How the value of ρ=βxxxβxyy varies from 3 to 3 influence the total power, parallel and perpendicular components at different polarization angle α has been three-dimensionally presented in Fig. 6 . The overall configuration of the SHG emission pattern indicates that the influence of minus value and positive value of ρ is not symmetric when α0deg . For more details, two-dimensional figures under four particular polarization angles α=0,30,60,90deg are selectively demonstrated in Fig. 7 . We note that no matter the total, parallel or the perpendicular components, at different polarization conditions, the change patterns as ρ are quite different. And from minus value to positive value of ρ , the change patterns are not symmetric except under the excitation condition of α=0deg . At polarization angle α=90deg , different ρ makes P2ω , P2ω , and P2ω have no change. At polarization α=0deg , as the absolute value of ρ increases, all SHG emission power (total, parallel, and perpendicular components) increases; however, the increase; of ρ has minor contribution to the perpendicular one. At polarization α=30deg , as ρ changes from 3 to 3, the total, parallel, and perpendicular components of the SHG emission power change from high to low value and then increase again, but the lowest turning point does not locate at ρ=0 . At polarization α=60deg , the curves increase slowing and have no obvious turning point along ρ except in the perpendicular one, which slightly shows a turning.

Fig. 6

Three-dimensional demonstration of effects of ratio ρ of hyperpolarizability elements of βxxx to βxyy on total, parallel, and perpendicular SHG power emission at different polarization angle α .

044016_1_011904jbo6.jpg

Fig. 7

Two-dimensional demonstration of effects of ratio ρ of hyperpolarizability elements of βxxx to βxyy on total, parallel, and perpendicular SHG power emission at different polarization angle α .

044016_1_011904jbo7.jpg

4.

Discussion and Conclusions

In this paper, a theoretical model of SHG on type I collagen under polarized focused beam has been established. On the basis of this model, we have investigated the total, parallel, perpendicular SHG emission power and their distribution along (θ,φ) , induced by different polarization angles α of the illumination light. The total, parallel, and perpendicular SHG emission power periodically change as the polarization angle α , similar to the pattern produced from the uniformly distributed dye dipoles on the cell membrane. The distribution of parallel, perpendicular components along (θ,φ) is complicated and possesses multiple peaks, but they symmetrically distribute around the plane, where fibrils locate ( θ=0deg , φ=180deg ). NA is another important factor that influences the illumination source. Therefore, the effects of NA on SHG emission are also investigated. At lower NA (for example, NA=0.2 ), the SHG emission concentrates in two closer lobes (like two lines), while at relative higher NA, such as NA=0.8 , SHG emission pattern has changed to be two obvious distributed lobes. At lower NA, where the light is loosely focused, it acts like a collimated light on the samples; thus, the influence from the Guoy phase shift can be ignored and dispersion alone plays role on SHG emission. According to θ=cos1[ξ(nωn2ω)] , when ξ=1 (unfocused light by very low NA), the emission angle θ is smaller than that at ξ<1 (focused light by high NA). However, at relative higher NA, where the light is tightly focused, because of the Guoy shift (phase shift) induced by focusing as well as the dispersion effect, to keep the momentum equilibrium, SHG emission strays far away from the direction of illumination light. As NA increases, we realizes that the waists of wxy and wz decreases. On the basis of expression A in Eq. 12, we understand that NA not only affects the emission angle θ , but also the distribution of the emission. At higher NA, SHG emission is much more distributed. The combination effect of NA under different polarization angles α is also studied. We find the influence of NA on SHG emission under different excitation conditions of polarization angle α is different, and the contribution of different NA to SHG emission is not linearity. However, our theoretical analysis results of NA are not applicable for linearly polarized light being focused with an objective with very high NA, because under such conditions, the polarization is no longer constant within the focus volume. At last, the ratio ρ of two existing components of β , which are βxxx and βxyy in our case is explored. Twenty-seven items of β are reduced to seven nonzero elements (two independent) through cylindrical symmetry and rule of Kleinman symmetry. For the light linearly polarized in x-y plane (εz=0) , those leave only two elements βxxx and βxyy . Here, the inconsistency of the polarization angle α with respect to the fiber axis keeps one more valid hyperpolarizability element βxyy of collagen in the model rather than βxxx alone at α=0deg . As we have mentioned in Sec. 3.3, the significance of ρ to collagen is obvious; therefore, SHG emission study of the changes of ρ may offers us important clues of the relationship between SHG emission with collagen features, such as its maturity, polarity, as well as the pathologically induced characterization changes.

In conclusion, our theoretical simulation results provide useful clues for experimental study of microscopic SHG in collagen type I under linearly polarized beam.

Appendices

{ label needed for app[@id='xA'] }

Appendix A: Simplification of Nonlinear Optical Hyperpolarizability

μ2ω,i(x,y,z) has three components, which are in the x , y , and z directions respectively,

Eq. 21

μ2ω,i(x,y,z)=12Eω2(x,y,z)j,kβijkε̂jε̂k.
Its subscript i presents these three directions. On the basis of the definition of the second-order tensor, hyperpolarizabilities β owns 27 items

Eq. 22

(βxxxβxyyβxzzβxyzβxzyβxzxβxxzβxxyβxyxβyxxβyyyβyzzβyyzβyzyβyzxβyxzβyxyβyyxβzxxβzyyβzzzβzyzβzzyβzzxβzxzβzxyβzyx).

Kleinman Symmetry

The assumption of Kleinman symmetry condition is applied to lossless media.23 It is valid if the frequency ω of excitation light is far from any natural frequencies of material. Because near-infrared light is usually used for SHG in biological tissue, in our case, if light with 800-nm wavelength is assumed to be used, then it is far from the wavelength of the first electronic transition in collagen (310nm) ,24 and therefore, Kleinman symmetry applies. Within in the validity of Kleinman symmetry,25, 26, 27, 28

Eq. 23

βijk=βjki=βkij=βikj=βjik=βkji.
As a result, only ten independent elements of β remain, which are as follows:

24.

xxx,yyy,zzz,xyy=yxy=yyx,xzz=zxz=zzx,xxy=xyx=yxx,
xxz=xzx=zxx,yzz=zyz=zzy,zyy=yzy=yyz,xyz=xzy=yxz=yzx=zxy=zyx.

Cylindrical Symmetry

Collagen fiber possesses cylindrically symmetrical feature in structure29, 30, 31 (as shown in Fig. 8 ). Under such condition, the hyperpolarizability tensor can be further simplified. The cylinder symmetry indicates that by taking the changes of xx , yy , zz on the subscripts of matrix βijk in Eq. 22, the derived new matrix of βijk in Eq. 5 should be the same as those in Eq. 22

Eq. 25

(xxxxyyxzzxyzxzyxzxxxzxxyxyxyxxyyyyzzyyzyzyyzxyxzyxyyyxzxxzyyzzzzyzzzyzzxzxzzxyzyx).
Therefore, the elements in Eq. 25 with minus sign should be zero and the elements of β now become

Eq. 26

(xxxxyyxzzxyzxzy000000000yzxyxzyxyyyx00000zzxzxzzxyzyx).

Fig. 8

Schematic diagram of cylindrical symmetry of collagen fiber along x -axis.

044016_1_011904jbo8.jpg

Also, from Fig. 8, we get the relation between y and z because the cylinder symmetry by applying xx , yy , zz

Eq. 27

x=x,y=ycosϕzsinϕ,z=ysinϕ+zcosϕ.
A project matrix is used to present the relationship.

Eq. 28

(xyz)=(1000cosϕsinϕ0sinϕcosϕ)(xyz)=(lxxlxylxzlyxlyylyzlzxlzylzz)(xyz).
Or, further to be

Eq. 29

βijk=liiljjlkkβijk.
According to the above expression, we have the following relationship:
βxxx=lxxlxxlxxβxxx,βxyy=lxxlyzlyzβxzz=sin2ϕβxzz,
βxyz=lxxlyzlzyβxzy=sin2ϕβxzy,βyyx=lyzlyzlxxχzzx=sin2ϕβzzx
βyxy=lyzlxxlyzβzxz=sin2ϕβzxz,βyzx=lyzlzylxxχzyx=sin2ϕβzyx
βyxz=lyzlxxlzyχzxy=sin2ϕβzxy
The rotation between y and z , ϕ=90 or 270deg , makes sin2ϕ=1 . Accoding to Eq. 4, we get

30.

Eq. 30a

βxxx,

Eq. 30b

βxyy=βxzz=βyyx=βzzx=βyxy=βzxz,

Eq. 30c

βxyz=βxzy=βyxz=βyzx=βzxy=βzyx=0.
As a result, there are only seven nonzero elements, in particular, two independent elements.

In our case, for the light linearly polarized in xy plane (Eω) , εz=0 . It remains only two (βxxx,βxyy) of the above seven nonzero elements has the contribution to the dipole moment μ , thus to the SHG emission.

{ label needed for app[@id='xB'] }

Appendix B: Integration of SHG emission electric field

In this part, we demonstrate the integration of the SHG emission electric field in x , y , z direction from negative infinite to positive infinite.

31.

Eq. 31a

E2ω(θ,φ)=ηNVrsinψμ2ω(x,y,z)exp[ik2ω(zcosθ+ysinθsinφ+xsinθcosφ)]dxdydz=ηNVr(sin2θsin2φ+cos2θ)1212βEω(0)2exp(2x2+y2wxy22z2wz+2iξkωz)×exp[ik2ω(zcosθ+ysinθsinφ+xsinθcosφ)]dxdydz=NVηr(sin2θsin2φ+cos2θ)12μ2ω(0)

Eq. 31b

×exp(2x2wxy2ik2ωxsinθcosφ)dx

Eq. 31c

×exp(2y2wxy2ik2ωysinθsinφ)dy

Eq. 31d

×exp(2z2wz2+2iξkωzik2ωzcosθ)dz.

According to the following integral equation,

Eq. 32

+exp(p2x2±qx)dx=exp(q24p2)πp(Rep2> 0)
Equations 31b, 31c, 31d have the corresponding parameters of p and q , such as

33a.

p2=2wxy2,p=2wxy,q=ik2ωsinθcosφ,q2=k2ω2(sinθcosφ)2
exp(2x2wxy2ik2ωxsinθcosφ)dx=exp[k2ω28wxy2(sinθcosφ)2]π2wxy

33b.

p2=2wxy2,p=2wxy,q=ik2ωsinθsinφ,q2=k2ω2(sinθsinφ)2
exp(2y2wxy2ik2ωysinθsinφ)dy=exp[k2ω28wxy2(sinθsinφ)2]π2wxy

33c.

p2=2wz2,p=2wz,q=ik2ωcosθ2iξkω
kω=2πnωλω,k2ω=2πn2ωλ2ω,ξ=ξnωn2ω,
q2=k2ω2(cosθξ)2
exp(2z2wz2+2iξkωzik2ωzcosθ)dz=exp[k2ω28wz2(cosθξ)2]π2wz.

{ label needed for app[@id='xC'] }

Appendix C: Total Power of SHG

The relation between total power distribution of SHG and electric field is18

Eq. 34

P2ω(θ,φ)=12n2ωε0cr2E2ω2(θ,φ)=12n2ωε0cr2N2A2(θ,φ)E2ω(0)2(θ,φ)=12n2ωε0cr2N2A2(θ,φ)η2r2(sin2θsin2φ+cos2θ)μ2ω(0)2=18n2ωε0cη2N2A2(θ,φ)Eω(0)4(sin2θsin2φ+cos2θ)β2=18n2ωε0cη2N2Eω(0)4A2(θ,φ)(1sin2θcos2φ)(βxxxcos2α+βxyysin2αβxyysin2α0)2.
The total power of SHG then is derived from double integral on θ and φ in a solid angle,

35.

P2ω=P2ω(θ,φ)sinθdθdφ=12n2ωε0cr2N2A2(θ,φ)η2r2(sin2θsin2φ+cos2θ)μ2ω(0)2sinθdθdφ=12n2ωε0cη2N2A2(θ,φ)(sin2θsin2φ+cos2θ)μ2ω(0)2sinθdθdφ=12n2ωε0cη2N2Θ(μ2ω,x(0)2+μ2ω,y(0)2+μ2ω,z(0)2)=18n2ωε0cη2N2ΘEω(0)4[(βxxxcos2α+βxyysin2α)2+(βxyysin2α)2].
Θ=A2(θ,φ)(1sin2θcos2φ)sinθdθdφ
μ2ω(0)=(μ2ω,x(0)μ2ω,y(0)μ2ω,z(0))=12Eω(0)2(βxxxcos2α+βxyysin2αβxyysin2α0).

{ label needed for app[@id='xD'] }

Appendix D: Parallel and Perpendicular Components of SHG Emission Power

(E2ωE2ω)=[cos(φα)sin(φα)sin(φα)cos(φα)](E2ωpE2ωs),
E2ω=E2ωpcos(φα)E2ωssin(φα),
E2ω=E2ωpsin(φα)+E2ωscos(φα).

The parallel and perpendicular SHG emission power distributions are

36.

Eq. 36a

P2ω(θ,φ)=12n2ωε0cr2E2ω(θ,φ)2=12n2ωε0cr2[cos(φα)E2ωpsin(φα)E2ωs]2,

Eq. 36b

P2ω(θ,φ)=12n2ωε0cr2E2ω(θ,φ)2=12n2ωε0cr2[sin(φα)E2ωp+cos(φα)E2ωs]2.
While
E2ω(0)(θ,φ)=(E2ω(0),p(θ,φ)E2ω(0),s(θ,φ))=ηrMμ2ω(0)=ηrM12Eω(0)2β=η2rEω(0)2(cosθcosφcosθsinφsinθsinφcosφ0)(βxxxcos2α+βxyysin2αβxyysin2α0)=η2rEω(0)2(cosθcosφ(βxxxcos2α+βxyysin2α)+βxyycosθsinφsin2αsinφ(βxxxcos2α+βxyysin2α)+βxyycosφsin2α).
Therefore,
E2ω(0),p(θ,φ)=η2rEω(0)2[cosθcosφ(βxxxcos2α+βxyysin2α)+βxyycosθsinφsin2α],
E2ω(0),s(θ,φ)=η2rEω(0)2[sinφ(βxxxcos2α+βxyysin2α)+βxyycosφsin2α],
E2ω(θ,φ)=NE2ω(0)A(θ,φ)=ηNA(θ,φ)2rEω(0)2(cosθcosφ(βxxxcos2α+βxyysin2α)+βxyycosθsinφsin2αsinφ(βxxxcos2α+βxyysin2α)+βxyycosφsin2α),
E2ωp(θ,φ)=ηNA(θ,φ)2rEω(0)2[cosθcosφ(βxxxcos2α+βxyysin2α)+βxyycosθsinφsin2α],
E2ωs(θ,φ)=ηNA(θ,φ)2rEω(0)2[sinφ(βxxxcos2α+βxyysin2α)+βxyycosφsin2α].
At last,

37.

Eq. 37a

P2ω(θ,φ)=18n2ωε0cη2N2A2(θ,φ)Eω(0)4{cos(φα)[cosθcosφ(βxxxcos2α+βxyysin2α)+βxyycosθsinφsin2α]sin(φα)[sinφ(βxxxcos2α+βxyysin2α)+βxyycosφsin2α]}2.

Eq. 37b

P2ω(θ,φ)=18n2ωε0cη2N2A2(θ,φ)Eω(0)4{sin(φα)[cosθcosφ(βxxxcos2α+βxyysin2α)+βxyycosθsinφsin2α]+cos(φα)[sinφ(βxxxcos2α+βxyysin2α)+βxyycosφsin2α]}2.

Acknowledgment

The authors gratefully thank the National Natural Science Foundation of China (Grant No: 30470495) for its support.

References

1. 

W. R. Zipfel, R. M. Williams, and W. W. Webb, “Nonlinear magic: multiphoton microscopy in the biosciences,” Nat. Biotechnol., 21 (11), 1368 –1376 (2003). https://doi.org/10.1038/nbt899 1087-0156 Google Scholar

2. 

F. Helmchen and W. Denk, “Deep tissue two-photon microscopy,” Nat. Methods, 2 (12), 932 –940 (2005). https://doi.org/10.1038/nmeth818 1548-7091 Google Scholar

3. 

E. Brown, T. McKee, E. diTomaso, A. Pluen, B. Seed, Y. Boucher, and R. K. Jain, “Dynamic imaging of collagen and its modulation in tumors in vivo using second-harmonic generation,” Nat. Methods, 9 (6), 796 –800 (2003). https://doi.org/10.1038/nm879 1548-7091 Google Scholar

4. 

G. Cox, E. Kable, A. Jones, I. Fraser, F. Manconi, and M. D. Gorrell, “3-Dimensional imaging of collagen using second harmonic generation,” J. Struct. Biol., 141 53 –62 (2003). https://doi.org/10.1016/S1047-8477(02)00576-2 1047-8477 Google Scholar

5. 

D. J. Prockop and A. Fertala, “The collagen fibril: the almost crystalline structure,” J. Struct. Biol., 122 111 –118 (1998). https://doi.org/10.1006/jsbi.1998.3976 1047-8477 Google Scholar

6. 

K. Beck and B. Brodsky, “Supercoiled protein motifs: the collagen triple-helix and the helical coiled coil,” J. Struct. Biol., 122 17 –29 (1998). https://doi.org/10.1006/jsbi.1998.3965 1047-8477 Google Scholar

7. 

I. Freund and M. Deutsch, “Second-harmonic microscopy of biological tissue,” Opt. Lett., 11 94 –96 (1986). https://doi.org/10.1364/OL.11.000094 0146-9592 Google Scholar

8. 

L. Moreaux, O. Sandre, and J. Mertz, “Membrane imaging by second-harmonic generation microscopy,” J. Opt. Soc. Am. B, 17 168 –1694 (2000). https://doi.org/10.1364/JOSAB.17.001685 0740-3224 Google Scholar

9. 

S. Roth and I. Freund, “Second harmonic generation in collagen,” J. Chem. Phys., 70 (4), 1637 –1643 (1979). https://doi.org/10.1063/1.437677 0021-9606 Google Scholar

10. 

I. Freund, M. Deutsch, and A. Sprecher, “Connective tissue polarity optical second-harmonic microscopy, crossed-beam summation, and small-angle scattering in rat-tail tendon,” Biophys. J., 50 693 –712 (1986). https://doi.org/10.1016/S0006-3495(86)83510-X 0006-3495 Google Scholar

11. 

R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second harmonic generation images of collagen I fibrils,” Biophys. J., 80 (2), 1377 –1386 (2005). https://doi.org/10.1529/biophysj.104.047308 0006-3495 Google Scholar

12. 

P. J. Campagnola, M.-d. Wei, A. Lewis, and L. M. Loew, “High-resolution nonlinear optical imaging of live cells by second harmonic generation,” Biophys. J., 77 3341 –3349 (1999). https://doi.org/10.1016/S0006-3495(99)77165-1 0006-3495 Google Scholar

13. 

R. Gauderon, P. B. Lukins, and C. J. R. Sheppard, “Simultaneous multichannel nonlinear imaging: combined two-photon excited fluorescence and second-harmonic generation microscopy,” Micron, 32 685 –689 (2001). https://doi.org/10.1016/S0968-4328(00)00067-6 0968-4328 Google Scholar

14. 

P. J. Campagnola and L. M. Loew, “Second-harmonic imaging microscopy for visualizing biomolecular arrays in cells, tissues and organisms,” Nature (London), 21 (11), 1356 –1360 (2003). https://doi.org/10.1038/nbt894 0028-0836 Google Scholar

15. 

J. E. Bjorkholm, “Optical second-harmonic generation using a focused Gaussian laser beam,” Phys. Rev., 142 (1), 126 –136 (1966). https://doi.org/10.1103/PhysRev.142.126 0031-899X Google Scholar

16. 

A. Willetts, J. E. Rice, and D. M. Burland, “Problems in the comparison of theoretical and experimental hyperpolarizabilities,” J. Chem. Phys., 97 (10), 7590 –7599 (1992). https://doi.org/10.1063/1.463479 0021-9606 Google Scholar

17. 

E. Hecht, Optics, Higher Education Press, Beijing (2005). Google Scholar

18. 

J. Mertz and L. Moreaux, “Second-harmonic generation by focused excitation of inhomogeneously distributed scatterers,” Opt. Commun., 196 325 –330 (2001). https://doi.org/10.1016/S0030-4018(01)01403-1 0030-4018 Google Scholar

19. 

C. Odin, Y. L. Grand, A. Renault, L. Gailhouste, and G. Baffet, “Orientation fields of nonlinear biological fibrils by second harmonic generation microscopy,” J. Microsc., 229 32 –38 (2008). https://doi.org/10.1111/j.1365-2818.2007.01868.x 0022-2720 Google Scholar

20. 

A. N. Bashkatov, E. A. Genina, V. I. Kochubey, and V. V. Tuchin, “Estimation of wavelength dependence of refractive index of collagen fibers of scleral tissue,” Proc. SPIE, 4162 265 –268 (2000). https://doi.org/10.1117/12.405952 0277-786X Google Scholar

21. 

S.-W. Chu, S.-Y. Chen, G.-W. Cher, T.-H. Tsai, Y.-C. Chen, B.-L. Lin, and C.-K. Sun, “Studies of x(2)x(3) tensors in submicron-scaled bio-tissues by polarization harmonics optical microscopy,” Biophys. J., 86 3914 –3922 (2004). https://doi.org/10.1529/biophysj.103.034595 0006-3495 Google Scholar

22. 

Y. Fu, H. Wang, R. Shi, and J. X. Cheng, “Second harmonic and sum frequency generation imaging of fibrous astroglial filaments in ex vivo spinal tissues,” Biophys. J., 92 3251 –3259 (2007). https://doi.org/10.1529/biophysj.106.089011 0006-3495 Google Scholar

23. 

R. W. Boyd, Nonlinear optics, (1992) Google Scholar

24. 

P. Stoller, K. M. Reiser, P. M. Celliers, and A. M. Rubenchik., “Polarization-modulated second harmonic generation in collagen. Biophys,” Biophys. J., 82 3330 –3342 (2002). https://doi.org/10.1016/S0006-3495(02)75673-7 0006-3495 Google Scholar

25. 

K. A. Kleinman, “Nonlinear dielectric polarization in optical media,” Phys. Rev., 126 (6), 1977 –1979 (1926). https://doi.org/10.1103/PhysRev.126.1977 0031-899X Google Scholar

26. 

T. L. Mazely and W. M. Hetherington, “Second-order susceptibility tensors of partially ordered molecules on surfaces,” J. Chem. Phys., 86 (6), 3640 –3647 (1987). https://doi.org/10.1063/1.451967 0021-9606 Google Scholar

27. 

C. A. Dailey, B. J. Burke, and G. J. Simpson, “The general failure of Kleinman symmetry in practical nonlinear optical applications,” Chem. Phys. Lett., 390 8 –13 (2004). https://doi.org/10.1016/j.cplett.2004.03.109 0009-2614 Google Scholar

28. 

J. F. Ward and P. A. Franken, “Structure of nonlinear optical phenomena in potassium dihydrogen phosphate,” Phys. Rev., 133 A183 –190 (1964). https://doi.org/10.1103/PhysRev.133.A183 0031-899X Google Scholar

29. 

P. Stoller, P. M. Celliers, K. M. Reiser, and A. M. Rubenchik, “Quantitative second-harmonic generation microscopy in collagen,” Appl. Opt., 42 (25), 5209 –5219 (2003). https://doi.org/10.1364/AO.42.005209 0003-6935 Google Scholar

30. 

S. V. Plotnikov, A. C. Millard, P. J. Campagnola, and W. A. Mohler, “Characterization of the Myosin-based source for second-harmonic generation from muscle sarcomeres,” Biophys. J., 90 693 –703 (2006). https://doi.org/10.1529/biophysj.105.071555 0006-3495 Google Scholar

31. 

M. Born and E. Wolf, Principles of Optics, 6th ed.Oxford, Pergamon Press, London (1993). Google Scholar
©(2009) Society of Photo-Optical Instrumentation Engineers (SPIE)
Ying Chang, Changshui Chen, Jianxin Chen, Ying Jin, and Xiao-Yuan Deng "Theoretical simulation study of linearly polarized light on microscopic second-harmonic generation in collagen type I," Journal of Biomedical Optics 14(4), 044016 (1 July 2009). https://doi.org/10.1117/1.3174427
Published: 1 July 2009
JOURNAL ARTICLE
12 PAGES


SHARE
Advertisement
Advertisement
Back to Top