|
1.IntroductionNonlinear optics has been proved to be a powerful tool for biomedical studies.1 Since its first demonstration by Denk 2 in 1990, two-photon fluorescence microscopic imaging has been widely applied to a variety of biological imaging tasks in aspects of physiology, morphology, and cell-cell interactions, etc.1, 2, 3 More recently, nonlinear microscopies based on radiative second-harmonic generation (SHG) have dramatically emerged as a viable microscope imaging contrast mechanism for visualization of cell and tissue structures and functions.4, 5, 6, 7, 8, 9 The nonlinear optical effect known as SHG was recognized in the earliest days of laser physics. Two decades ago, it was successfully demonstrated through a microscope by Gannaway and Sheppard10 and was first applied11 in biological tissue in 1986. SHG is commonly called frequency doubling and the second-harmonic light emerging from the material is at precisely half the wavelength of the light entering the material. SHG has been widely investigated both experimentally and theoretically. Second-harmonic signals from collagen fibrils have been quantitatively explored based on theoretical physical models12, 13, 14 and the biophysical features characteristic of collagen structures can be resolved from the experimental second-harmonic microscopic images on the basis of the SHG phase-matching fundamental principles.15 SHG from both single and arrays of spherical particles have been theoretically explored or experimentally observed.16, 17, 18, 19, 20 Meanwhile, the second-harmonic (SH) cross section, as well as the average cosine of the SH scattering angle is theoretically calculated in a dilute suspension of spherical particles confined within a very thin slab21 (the thickness of the slab scattering length of the turbid medium). The angle-resolved pattern of the SH light scattering from suspensions of different sizes of centrosymmetric microsize polystyrene spheres with surface-absorbed dye (malachite green) has been experimentally measured.22 However, as we have noticed, all the researchers have focused on the generation mechanism study of SH signals in either an isolate or a suspension of particles. The scattering effect of the turbid medium itself in which the SH signals generate and go through has not been paid much more attention. As our intension is directed to the SH microscopic imaging studies in the biological tissue, a kind of highly scattered media, the studies of the SH signals through a turbid medium becomes absolutely necessary. Because biological tissue is usually composed of scatters with a size variation from nanometers to a few micrometers,23 the dominant effect caused by these scatterers is Mie scattering. To understand the effect of multiple scattering in multiphoton fluorescence microscopy, one usually uses the Monte Carlo simulation method based on Mie scattering because the conventional image theory based on Fourier optics is not applicable.24, 25, 26 So far the established Monte Carlo method in multiphoton microscopy is in fluorescence applications.27, 28, 29 This paper establishes a Monte Carlo simulation model for SH microscopic imaging through turbid media. The paper is organized as follows. Section 2 describes the simulation model of the SH microscopic imaging through a turbid medium. In Sec. 3, we demonstrate the simulation results of the angle-resolved SH signal distribution and the signal level under different constraints of the scattered numbers . The effect of the parameters such as the thickness of the turbid media, the numerical aperture (NA) of the objective, and the particles size that form the turbid medium are investigated. Discussion and conclusions are presented in Sec. 4. 2.Simulation Model for SH Microscopic Imaging through a Turbid MediumFigure 1 shows the schematic diagram for SH microscopic imaging through turbid media in a transmission mode. It is supposed that the excitation light is focused through an imaging objective with on to a focal plane. An SH signal of a certain strength and direction within the focal spot are generated and then propagate through a thickness in the turbid medium. The focal depth is the distance between the turbid medium’s surface and the focal plane. In this paper, we choose (imaging near the surface) to avoid the influence coming from the polarization change of the excitation light due to scattering. The turbid medium is supposed to be a slab formed by microsize scattering particles . Fig. 1Schematic diagram of imaging through a turbid medium in a transmission-mode SH microscope, where is the focal depth, and is the thickness of the turbid medium. The coordinate system defining the SHG emission direction is also shown. The shaded part demonstrates the membrane surface, which can be approximated to be planar at the length scales considered. The focused excitation beam propagates in the direction. The resultant SHG signal is mostly confined to an interaction area schematically depicted as a dotted ellipsoid and radiates in the directions defined by and . ![]() The Monte Carlo simulation method for multiphoton fluorescence microscopic imaging through the turbid medium is adopted and the basic principle is the same as that described elsewhere.27, 28, 29 However, for SH microscopic imaging, the implementation is much more complicated than that in fluorescence imaging. Unlike the process of two-photon-excited fluorescence (TPF), an SHG process is a coherent scattering process, in which the scattered photons must satisfy a phase-matching constraint, and thus produces highly directed radiation rather than isotropic emission. According to a simplified biophysical model, the total power distribution of SHG emitted from the cell membrane labeled by Di-6ASPBS styryl dye is (the details can be found in the Appendix30): Therefore, three steps are involved to implement it. First, the average photon intensity distribution at focal depth is calculated using the Monte Carlo simulation. Since in this paper, it is supposed that a pulse laser with pulse duration and a repetition frequency of is used as the excitation light to approach reality, with the pulse laser, the two-photon excitation depends on the average squared intensity rather than the squared average intensity , which yields1 where the dimensionless factor is used for a Gaussian pulse shape.Second, the simulation starts from the beginning. This time, for each excitation photon, a direction is produced according to , (where is a uniformly distributed random number) and SH signals will generate with the weighting factor and the efficiency factor at that direction . Once it is generated, the propagation of the SH light follows the same way as that detailed in Ref. 29 for TPF. The excitation wavelength is assumed to be and the SH emitted is therefore assumed to be . The excitation power after the objective is . The collective NA of the objective in this transmission model is supposed to be always the same as that of the excitation objective. According to the Mie scattering theory,31 the corresponding parameters of anisotropy value , the scattering cross section and the scattering mean-free-path length of the scatterers (at a concentration of ) in the simulated turbid medium are shown in Table 1 . Table 1Scattering parameters of the turbid medium under 2p excitation.
The concentration of the molecular surface density is assumed to be and the first hyperpolarizability of the dye molecule is . Two indicators, the angle-resolved SH signal distribution on the emitted plane along the direction and the signal level collected immediately after the emitted plane, are adopted in this paper. The signal level is defined as the number of SH photons collected by the detector under different maximum scattering number limits and normalized by the number when no scattering exists (or the thickness ). Here, we should clarify that we use the limit of to implicate that the SH signals are dominated by ballistic photons, which means the basic features of the photon, such as the polarization, the coherence, etc., are kept at a maximum degree; the limit of includes the effect from both ballistic and snakes SH photons; while means all the effects including the multiple scattered photons are considered. 3.Angle-Resolved SH Signals and Signal Level through Turbid Media3.1.Effect of the Thickness of the Turbid MediaThe thickness of the turbid medium is important when we decide to explore the signal changes through a turbid medium as it directly influences the number of the scattering events, thus inducing the direction change of the photons and the collection efficiency. Therefore, in this paper, the effect of the thickness is first investigated. Figures 2a and 2b show the angle-resolved distributions of the total SH signals under in a turbid medium formed by the scattering particle size of under and , respectively. From the comparison of the Figs. 2a and 2b, we can see that as the increases, the angle distribution of the SH signals becomes disperse along both the and directions. However, the lobe effect induced by the Gouy shift31 becomes unobvious because of the multiple scattering. Fig. 2Angle-resolved distributions of the SH signals under in the turbid medium formed by scattering particles of a size of at (a) and (b) . ![]() Figure 3 shows the change of the signal level as the increase of the thickness under different maximum scattering number limits , respectively. It is shown that as the increases, the collected SH signals decrease under all three conditions. In practice, we are much more interested in the SH signals that are composed of ballistic photons (the solid line in Fig. 3 presents for the ballistic SH signal level versus ), since the fewer scattering events they experience, the less carried information they lose. It is shown that as the thickness reaches , the SH signals made up of ballistic photons almost disappear. In another words, the penetration limit of ballistic SH signals is under the condition described. After an approximate penetration depth of , the signal level drops to half the height of the total signals. 3.2.Effect of the Objective NAThe size of the objective NA is proved to play a role in multiphoton fluorescence microscopic imaging.32 In this paper, its effect on the SH microscopic imaging is thus investigated. Figures 4a and 4b demonstrate the angle distribution of the total collected SH signal under under thicknesses and , respectively. It is clear from a comparison of Fig. 4 with Fig. 2 that NA has a significant effect on both the angle distribution and the signal strength of SH signals. As NA becomes smaller, the angle distribution becomes narrower, while many more SH signals are collected. According to a relationship31 of to NA, where , in our situation, when and are used, the peak angle of the lobes of SH emission will locate at and and and , respectively. This provides the explanation for why at a low NA, the distribution of the ballistic SH signals are narrow (Fig. 4). Meanwhile, because and are1 inversely proportional to NA, as the NA becomes low, the effective total number of molecules that contribute to the generation of SH light becomes large. Therefore, the power of SH light becomes strong (Fig. 4). Figure 5 shows the corresponding signal level versus under three scattering number constraints. We notice that NA has no obvious effect on the penetration limit regarding the ballistic SH signals if the same turbid medium is provided (see solid lines in Figs. 3 and 5). However, it does influence the collected snake and multiple scattered SH signals, which can be realized from a comparison of the dotted and the dashed lines in Figs. 3 and 5. For a small NA (Fig. 5), the collected snake and multiple scattered SH signals have a greater rate of decrease as increases than that for a large NA, which means a small NA has a suppression function to snake and multiple scattered SH signals. 3.3.Size Effect of the Spherical Particles in the Turbid MediumThe size effect of the spherical particles that form the turbid medium is also part our investigation scope in this paper, since different scatterer sizes have different scattering cross sections and anisotropy values according to Mie scattering, which may directly influence the behavior of SH light during the turbid medium. In this paper, two sizes of particles are chosen as our objects. One is the larger one with a diameter of and the other is the smaller one with a diameter of . Figures 6a and 6b demonstrate the total angle-distribution of SH signals under thicknesses of and , respectively. It is shown that the angle-distribution has no obvious difference under thicknesses of compared with that for . Figure 7 is the signal level versus under three scattering number constraints. We noticed that the signal level under three scattering number constraints also has a slight difference along the whole thickness and the drop of the signal level is very low. This is because the small scattering particles have a much smaller cross section than that of large particles (Table 1). Under the same concentration, the mean-free path length is much longer than that of large particles, which leads to fewer scattering events at a given thickness . Thus, the penetration limit is much deeper than that in turbid medium with large scattering particles. 4.Discussion and ConclusionThe fluorescence signals to the excitation intensity have the relationship of ( corresponding to , , and excitation). In our previous investigation, the model of multiphoton fluorescence microscopic imaging through the turbid medium based on the Monte Carlo simulation method was successfully established,27, 28, 29 where fluorescence photons are generated by a uniform random generator with the weighting factor of and the fluorescence efficiency factor . In the case of an SH microscopic imaging condition, the phase matching between the SHG and the excitation fields largely prevents SHG propagation from the forward direction, and forces it to propagate off-axis to form two well-defined lobes. In this case, the SH photons generated are no longer as uniformly randomized as these of fluorescence photons, instead, the probability of the generation of the SH photons in each direction are determined based31 on the angular structure parameter . Therefore, at direction as the efficiency factor is implemented by Monte Carlo simulation to get the SH photons. To efficiently demonstrate our model, in this paper, we simplified our condition of . As imaging at a nonzero focal depth in the turbid medium, the effect of the polarization change of the excitation beam due to the scattering must be considered. This may greatly increase the complication of the problem, and it becomes another interesting topic calling for further investigation. To summarize, we demonstrated the simulation model based on the Monte Carlo method for dealing with SH signals through a turbid medium for the first time. The angle-resolved distribution and signal level under various conditions were explored. Simulation results reveal that the use of the small objective NA results in a narrow angle distribution and strong SH signals. The turbid medium composed of large scattering particles has a strong influence on the angle distribution and the signal level of SH light, which results in a low penetration limit in the turbid medium for ballistic SH signals, which is approximate . 5.AppendixIt is assumed that SHG is produced from Di-6ASPBS styryl dye molecules labeled on the cell’s membrane surface. A tightly focused driving field under a microscope is polarized in the direction and the molecules in the membrane are perfectly oriented in the direction and strictly uniaxial. Under this condition, we can assume that the amplitude of the excitation light field with frequency (wavelength ) has a Gaussian profile in both the axial and the lateral directions about the focal center and that its phase near the focal center progresses linearly. That is,30 where is the local wave vector ( are the indices of refraction at ); is a parameter that characterizes the phase shift, commonly referred to as a Gouy shift or a phase anomaly, experienced by a Gaussian beam in the vicinity of a focal center; and and are the axial and the radial field waists at , respectively, given by1According to Ref. 30, the total power distribution of SHG is where is the second-harmonic wave vector, is the excitation intensity in units of photons per second per area, is the molecular surface density, are the indices of refraction at , is Planck’s constant, and denotes the first hyperpolarizability of the dye molecular.AcknowledgmentThe authors gratefully thank the National natural Science Foundation of China (Grant No. 30470495) for its support and the authors also thank Prof. Min Gu of the Centre for Microphotonics, Swinburne University of Technology, for his kind help in modification. ReferencesW. R. Zipfel,
R. M. Williams, and
W. W. Webb,
“Nonlinear magic: multiphoton microscopy in the biosciences,”
Nat. Biotechnol., 21 1369
–1377
(2003). https://doi.org/10.1038/nbt899 1087-0156 Google Scholar
W. Denk,
J. H. Stricker, and
W. W. Webb,
“Two-photon laser scanning fluorescence microscopy,”
Science, 248 73
–76
(1990). 0036-8075 Google Scholar
K. König,
“Multiphoton microscopy in life sciences,”
J. Microsc., 200 83
–104
(2000). https://doi.org/10.1046/j.1365-2818.2000.00738.x 0022-2720 Google Scholar
X. Deng,
E. D. Williams,
E. W. Thompson,
X. Gan, and
M. Gu,
“Second-harmonic generation from biological tissues: effect of excitation wavelength,”
Scanning, 24 175
–178
(2002). 0161-0457 Google Scholar
P. J. Campagnola and
L. M. Loew,
“Second-harmonic imaging microscopy for visualizing biomolecular arrays in cells, tissues and organisms,”
Nat. Biotechnol., 21 1356
–1360
(2003). https://doi.org/10.1038/nbt894 1087-0156 Google Scholar
A. Zoumi,
A. Yeh, and
B. J. Tromberg,
“Imaging cells and extracellular matrix in vivo by using second-harmonic generation and two-photon excited fluorescence,”
Proc. Natl. Acad. Sci. U.S.A., 99 11014
–11019
(2002). https://doi.org/10.1073/pnas.172368799 0027-8424 Google Scholar
E. Brown,
T. Mckee,
“Dynamic imaging of collagen and its modulation in tumors in vivo using second-harmonic generation,”
Nature Med., 9 796
–800
(2003). 1078-0432 Google Scholar
D. A. Dombeck,
W. W. Webb,
“Uniform polarity microtubule assemblies imaged in native brain tissue by second-harmonic generation microscopy,”
Proc. Natl. Acad. Sci. U.S.A., 100 7081
–7086
(2003). https://doi.org/10.1073/pnas.0731953100 0027-8424 Google Scholar
P. J. Capagnola,
W. A. Mohler,
“Three-dimensional high-resolution second-harmonic generation imaging of endogenous structural proteins in biological tissues,”
Biophys. J., 81 493
–508
(2002). 0006-3495 Google Scholar
J. N. Gannaway and
C. J. R. Sheppard,
“Second-harmonic imaging in the optical scanning microscope,”
Opt. Quantum Electron., 10 435
–439
(1978). https://doi.org/10.1007/BF00620308 0306-8919 Google Scholar
I. Freund and
M. Deutsch,
“2nd-harmonic microscopy of biological tissue,”
Opt. Lett., 11 94
–96
(1986). 0146-9592 Google Scholar
P. Stoller,
A. M. Rubenchik,
“Quantitative second-harmonic generation microscopy in collagen,”
Appl. Opt., 42 5209
–5219
(2003). 0003-6935 Google Scholar
P. Stoller,
A. M. Rubenchik,
“Polarization-modulated second-harmonic generation in collagen,”
Biophys. J., 82 3330
–3342
(2002). 0006-3495 Google Scholar
P. Stoller and
L. B. Da Silva,
“Polarization-dependent optical second-harmonic imaging of a rat-tail tendon,”
J. Biomed. Opt., 7 205
–214
(2002). https://doi.org/10.1117/1.1431967 1083-3668 Google Scholar
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). 0006-3495 Google Scholar
V. L. Brudny,
B. S. Mendoza, and
W. L. Mochán,
“Second-harmonic generation from spherical particles,”
Phys. Rev. B, 62 11152
–11162
(2000). https://doi.org/10.1103/PhysRevB.62.11152 0163-1829 Google Scholar
J. I. Dadap,
J. Shan,
K. B. Eisenthal, and
T. F. Heinz,
“Second-harmonic Rayleigh scattering from a sphere of centrosymmetric material,”
83 4045
–4048
(1999). Google Scholar
N. Yang,
W. E. Angerer, and
A. G. Yodh,
“Second-harmonic microscopy of single micrometer-size particles on a substrate,”
Phys. Rev. A, 64 045801
(2001). https://doi.org/10.1103/PhysRevA.64.045801 1050-2947 Google Scholar
J. I. Dadap,
J. Shan, and
T. F. Heinz,
“Theory of optical second-harmonic generation from a sphere of centrosymmetric material: small-particle limit,”
J. Opt. Soc. Am. B, 21 1328
–1347
(2004). https://doi.org/10.1364/JOSAB.21.001328 0740-3224 Google Scholar
W. Luis Mochán,
J. A. Maytorena,
B. S. Mendoza, and
V. L. Brudny,
“Second harmonic generation in arrays of spherical particles,”
Reys. Rev. B., 68 085318
(2003). Google Scholar
E. V. Makeev and
S. E. Skipetrov,
“Second-harmonic generation in suspensions of spherical particles,”
Opt. Commun., 224 139
–147
(2003). https://doi.org/10.1016/S0030-4018(03)01756-5 0030-4018 Google Scholar
N. Yang,
W. E. Angerer, and
A. G. Yodh,
“Angle-resolved second-harmonic light scattering from colloidal particles,”
Phys. Rev. Lett., 87 103902
(2001). https://doi.org/10.1103/PhysRevLett.87.103902 0031-9007 Google Scholar
A. G. Loewy and
P. Siekevitz, Cell Structure and Function, 2nd ed.Holt, New York
(1971). Google Scholar
L. H. Wang and
G. Liang,
“Absorption distribution of an optical beam focused into a turbid medium,”
Appl. Opt., 38 4951
–4958
(1999). 0003-6935 Google Scholar
C. M. Blanca and
C. Saloma,
“Monte Carlo analysis of two-photon fluorescence imaging through a scattering medium,”
Appl. Opt., 37 8092
–8102
(1998). 0003-6935 Google Scholar
A. K. Dunn,
V. P. Wallace,
M. Coleno,
M. W. Berns, and
B. Tromberg,
“Influence of optical properties on two-photon fluorescence imaging in turbid samples,”
Appl. Opt., 39 1194
–1201
(2000). 0003-6935 Google Scholar
X. Deng,
X. Gan, and
M. Gu,
“Multiphoton fluorescence microscopic imaging through double-layer turbid media,”
J. Appl. Phys., 91 4659
–4665
(2002). https://doi.org/10.1063/1.1459107 0021-8979 Google Scholar
X. Deng,
X. Gan, and
M. Gu,
“Monte Carlo simulation of multiphoton fluorescence microscopic imaging through inhomogeneous tissue like turbid media,”
J. Biomed. Opt., 8 440
–449
(2003). https://doi.org/10.1117/1.1577116 1083-3668 Google Scholar
X. Gan and
M. Gu,
“Fluorescence microscopic imaging through tissue-like turbid media,”
J. Appl. Phys., 87 3214
–3221
(2000). https://doi.org/10.1063/1.372326 0021-8979 Google Scholar
L. Moreaux,
O. Sandre, and
J. Mertz,
“Membrane imaging by second-harmonic generation microscopy,”
J. Opt. Soc. Am. B, 17 1685
–1694
(2000). 0740-3224 Google Scholar
C. F. Bohren and
D. R. Huffman, Absorption and Scattering of Light by Small Particles, Wiley, New York
(1983). Google Scholar
X. Deng and
M. Gu,
“Penetration depth of single-, two- and three-photon fluorescence microscopic imaging through human cortex structures: Monte Carlo simulation,”
Appl. Opt., 42 3321
–3329
(2003). 0003-6935 Google Scholar
|