|Home | About | Journals | Submit | Contact Us | Français|
Quantitative determination of fluorophore content from fluorescence measurements in turbid media, such as tissue, is complicated by the influence of scattering properties on the collected signal. This study utilizes a Monte Carlo model to characterize the relationship between the fluorescence intensity collected by a single fiber optic probe (FSF) and the scattering properties. Simulations investigate a wide range of biologically relevant scattering properties specified independently at excitation (λx) and emission (λm) wavelengths, including reduced scattering coefficients in the range μ′s(λx) [0.1 – 8]mm−1 and μ′s(λm) [0.25 – 1] × μ′s(λx). Investigated scattering phase functions (P(θ)) include both Henyey-Greenstein and Modified Henyey-Greenstein forms, and a wide range of fiber diameters (df [0.2 – 1.0] mm) was simulated. A semi-empirical model is developed to estimate the collected FSF as the product of an effective sampling volume, and the effective excitation fluence and the effective escape probability within the effective sampling volume. The model accurately estimates FSF intensities (r=0.999) over the investigated range of μ′s(λx) and μ′s(λm), is insensitive to the form of the P(θ), and provides novel insight into a dimensionless relationship linking FSF measured by different df.
Detection and quantitation of fluorescence is important for many biomedical and clinical applications. The optical detection of fluorescent endogenous compounds  such as collagen and NADH, or exogenous compounds that include labelled markers, can be used for diagnostic purposes [2, 3]. The measurement of therapeutic compounds, such as photosensitizers used in photodynamic therapy [4, 5], may provide insight into the pharmacokinetic distribution and pharmacodynamic activity in tissues of interest and may play a role in monitoring administered therapies . However, quantitation of fluorescence in tissue in vivo is complicated by the influence of the tissue optical properties on the collected fluorescence signal . Absorption by chromophores within the tissue causes attenuation that is (non-linearly) proportional to the absorption coefficient at the excitation and emission wavelengths. Scattering within tissue is known to have a complicated effect on fluorescence measurements: the properties at the excitation wavelength (λx) affect the delivered excitation light profile and the properties at the emission wavelength (λm) determine the likelihood that fluorescent emission photons propagate to the detector used in the measurement. In order to quantitatively analyze fluorescence in tissue, it is important to obtain an intrinsic fluorescence signal that is independent of the optical property effects [8–10]. This approach would yield a quantity that is proportional to the product of the concentration and quantum yield of the fluorophore within the optically sampled volume, and would be comparable between measurements of samples with different background tissue optical properties.
Previously developed methods to extract intrinsic fluorescence spectra involve the acquisition of a paired measurement of fluorescence and white-light reflectance, where the latter is used to inform a correction of the influence of optical properties on fluorescence. This general approach has been extensively investigated for multi-fiber fluorescence probes, with separate source(s) and detectors [8–16]. These probes collect multiply scattered, or diffuse, light and sample volumes of tissue on the orders of several mm3. An alternative approach for fluorescence measurements is to use small fiber optic probes that utilize a single optical fiber to both deliver excitation light and collect emitted fluorescence [5, 17–21]; such a measurement results in a localized sampling volume, with the majority of the collected signal originating very close to the probe face . Single fiber fluorescence (SFFL) measurements collect photons that have undergone few scattering events, and in turn, have a very small light propagation path, making the collected intensity less sensitive to tissue chromophores and scattering properties than diffuse measurements. The influence of scattering on collected SFFL intensity has been previously investigated and was observed to be nonlinear and fiber-diameter specific . Furthermore, the SFFL intensity was observed to be insensitive to variations in the scattering phase function (PF) . The underlying mechanism of these factors was not fully elucidated. These and other previous studies accounted for the influence of scattering on SFFL by characterizing ranges of fiber diameters and optical property combinations where the SFFL signal was insensitive to optical properties [17, 19, 20]. While this approach may be useful for specific applications with well-known ranges of optical properties, it does not return a quantitative description of tissue fluorescence that is independent of optical properties, and therefore, does not provide a reliable comparison of measurements performed on different tissue locations or with different fiber diameters.
To the best of the authors’ knowledge, there is currently no analytical or empirical description of the influence of scattering properties on the fluorescence intensity sampled by a single fiber. The present study investigates the detailed mechanisms associated with the influence of scattering properties on the SFFL intensity measured in a turbid medium, and develops a mathematical model to correct for these influences. This represents a first step towards a full correction of collected SFFL intensities for the influence of optical properties (i.e. both scattering and absorption). Monte Carlo (MC) simulations are used to investigate SFFL measurement of a wide range of scattering properties that are independently specified at excitation and emission wavelengths; simulations also included a wide range of fiber diameters. Simulated data are used to identify and characterize a semi-empirical model that expresses SFFL intensity as a function of a dimensionless scattering property (given as the product of scattering coefficient and fiber diameter). The resulting model is applicable to all investigated fiber diameters and provides insight into the physics underlying the SFFL measurement.
The Monte Carlo (MC) code utilized in this study is a customized version of the MCML program  that is modified to emulate single fiber fluorescence measurements of a homogeneous turbid medium. The code allows independent specification of both the scattering coefficient (μs) and scattering phase function (P(θ)) at excitation (λx) and emission (λm) wavelengths. Excitation photons were initialized by selecting a location on the fiber face, which is modeled in contact with the turbid medium at the air/medium interface z = 0, and were launched into a direction within the fiber cone of acceptance, where the acceptance angle was given as ; both the location and the direction were sampled from uniform distributions. The index of refraction (n) of the medium and fiber were specified at 1.37 and 1.45, respectively, and were held constant between λx and λm. The numerical aperture (NA) of the fiber was set as 0.22. Reflection and refraction due to the index of refraction mismatch at the medium/fiber and the surrounding medium/air interface were calculated using the Fresnel equations and Snell’s law. This code simulated propagation of excitations photons by stochastically selecting step sizes (sn) from an exponential distribution weighted by μs(λx), and each scattering angle was selected from P(θ)(λx). At discrete points along each individual step, excitation photons were stochastically checked for a fluorescence event, with the probability given by , where is the specific absorption coefficient of the fluorophore. Stochastic absorption by the fluorophore resulted in an isotropic scattering event, and propagation of the emission photon was continued at the scattering properties at λm. Emission photons propagating within the turbid medium that cross the medium interfacial boundary at z = 0, were checked for contact with the fiber face; those in contact and traveling at an angle within the fiber cone of acceptance were collected, the rest were terminated. Excitation photons contacting the fiber face at any angle were terminated and did not contribute to the collected fluorescence intensity. This calculation returned the fraction of the number of collected fluorescence photons and the number of excitation photons for each simulation, calculated as:
where TXPL is the total number of excitation photons launched and TMPC the total number of emission photons collected. Excitation and emission photons propagating within the medium far from the fiber face do not contribute to the collected fluorescence intensity and were terminated at a hemispherical limit from the fiber face of ; a limit that was confirmed to not influence model outputs for the range of optical properties investigated in this study. Model outputs of were validated by comparison with previously published fluorescence intensities over a range of background optical properties .
During photon propagation, the photon positions were tracked in a discrete voxel grid to yield individual 2D(r,z) probability density profiles for all incident excitation photons, for all fluorescence emission photons, and a separate profile for all collected fluorescence photons. Specifically, the code generated 2D maps of the relative excitation light fluence (Φx(r,z) [m−2]), which is calculated as previously described , and of the photon probability density of fluorescence collected by the fiber (Fcol(r,z) [m−3]), which represents the spatial location of origin for all collected fluorescence photons . Note that these quantities involve ratio calculations and do not depend on the number of launched excitation photons. From these maps, the dimensionless escape probability density profile of emission photons , which is defined as the probability of emission photon collection per fluorescence photon generated, was calculated. Note that the fluorescence generated at a location (r,z) is proportional to the product of and Φx(r,z). These 2-D spatial profiles were used to calculate effective values for the volume sampled and the excitation fluence and escape probability within the sampled volume, by properly weighting each respective quantity by the collected fluorescence that originated at the corresponding location. A scalar effective optical sampling depth (ZMC [m]) is calculated as the weighted average depth of the collected emission photons, given as
where Δaj is the area of a voxel at position rj and Δz is the z-dimensional length of each voxel . A scalar effective excitation fluence within the optically sampled volume was calculated from the weighted average of Φx(r,z), with the collected fluorescence photon probability density Fcol(r,z) as weight factors, as
Similarly, a scalar for the effective escape probability within the optically sampled volume was calculated from the weighted average of the escape probability density distribution Hm(r,z), with the collected fluorescence photon probability density Fcol(r,z) as weight factors, as
MC simulations were performed over a broad range of biologically relevant  reduced scattering coefficient (μ′s) values that were individually specified at λx and λm, with: μ′s(λx) = [0.1,0.25,0.5,1,2,4,8] mm−1 and μ′s(λm) = [0.25,0.5,0.75,1.0] × μ′s(λx). This series of simulations was performed at all specified μ′s combinations using the Modified Henyey-Greenstein (MHG) PF  with the anisotropy specified as g1 = 0.9 and γ, which characterizes the first two moments of the phase function and is given as , was set as γ = 1.4.
A subset of simulations further investigated the influence of PF over a selected range of reduced scattering values, μ′s(λx) = [0.5,1,2] mm−1 and μ′s(λm) = [0.5,1.0] × μ′s(λx), using the Henyey-Greenstein (HG) PF with g1 = [0.5,0.9] and γ = [1.5,1.9] and the MHG PF with combinations of g1 = [0.8,0.9,0.95] and γ = [1.4,1.5,1.6,1.7,1.8,1.9].
Additionally, simulations investigated variations in NA from the baseline value of 0.22 over the range [0.1 – 0.4]. This subset of simulations was performed using the same scattering properties as the subset of simulations used to investigate the influence of PF.
Simulations of each possible combination of scattering properties were performed for a range of fiber diameters, with df = [0.2,0.4,0.6,1.0] mm. The absorption of the fluorophore was given as in all simulations; this study did not consider absorption due to background chromophores. In total, the data presented in this study include 616 MC simulations, each launching at least 20 million photons.
The fluorescence signal F (in units of Joules [J]) collected by a fiber optic probe is given by the integral 
where Φx(r) [J m−2] is the excitation fluence, Hm(r) [−] is the escape probability of emission photons, is the absorption coefficient of the fluorophore at the excitation wavelength, and Qf [−] is the fluorescence quantum yield. The ratio (λx/λm) accounts for the difference in photon energy between the emission and excitation wavelengths, and will be omitted in the remainder of the equations since in the Monte Carlo simulations this ratio is set to unity. The volume integral of ΦxHm is dependent on optical properties; however the intrinsic fluorescence, given by the product of is independent of optical properties and has dimensions [m−1].
This study develops an approximate solution to Eq. (5) for a SFFL measurement by representing the volume integral of ΦxHm as the product of an effective optically sampled volume and the effective Φx and Hm values within that volume, thus redefining Eq. (5) as,
where V is the effective sampling volume, and and are the effective excitation fluence and effective escape probability within the sampled volume, respectively. These quantities can be related to the scalar outputs from the MC simulations defined in Section 2.1 by approximating the effective sampling volume as
with ZMC the effective sampling depth and A1 a proportionality factor that depends on the geometrical shape of the effective volume. Furthermore, the effective excitation fluence and escape probability within the sampled volume can be written as
where Px is the total power output from the fiber, which in these simulations is proportional to the number of launched photons TXPL. Eq. (8) properly accounts for differences in the incident excitation intensity emitted from the fiber face for different numbers of launched photons TXPL, but does not correct for differences in incident excitation light intensity due to differences in fiber diameter. Since the incident excitation intensity is inversely proportional to the fiber area, it is expected that scales with . Substituting Eqs. (7–9) in Eq. (6), dividing by Px, which is then replaced by TXPL yields SFFL expressed as a percentage of incident excitation photons, as
As described in Section 2.1, the MC simulations used in this study were used to return information about how SFFL intensity and the effective terms presented in Eq. (10) are influenced by scattering properties at the excitation and emission wavelengths. Inspection of the simulated data led to the identification of candidate empirical expressions to describe each quantity; from these a set of equations was selected on the basis of fit quality and model simplicity, and is given as
where [A1,2,3, B1,2, C1,2,3] in Eqs. (10–13) are fitted parameters. The effective sampling depth ZMC was observed to follow an exponential decay with respect to the product of μ′s,avg(λx,λm)df, where μ′s,avg(λx,λm) is calculated as the average of μ′s(λx) and μ′s(λm). The effective relative excitation fluence was observed to scale with as expected, and to follow an exponential expression that depended on the product μ′s(λx)df. The effective escape probability was observed to follow an exponential function that dependent on the product μ′s(λm)df. Substituting Eqs. (11–13) into Eq. (10) results in
where [ζ1,ζ2,ζ3] are fitted parameters. This represents a reduction from the parameter set specified in Eqs. (10–13). Here, ζ1 represents the product of A1, A2, B1 and C1. Fitted parameters were estimated using a Levenberg-Marquardt algorithm coded into a Matlab script (version 2009a, MathWorks). Confidence intervals of the estimated parameters were calculated from the square root of the diagonal of the covariance matrix . During the model fit analysis, the estimated values for A3, B2 and C2 were observed to have overlapping 95% confidence intervals , which led to the reduction of these terms to a single fitted parameter, ζ2. This substitution did not result in a significant increase in model residual error.
Continuing the description of the terms in Eq. (14), νn represents the influence of the index of refraction mismatch at z = 0 (between fiber/medium and the annular air/medium interfaces). This parameter was found to be dependent on df, and to follow the form: , with = 0.17mm−1. This form was identified from comparing simulations of the fiber surrounded by air with the fiber surrounded by a refractive index matching the fiber; this factor is analogous to offset factors described previously .
Equation (14) represents a fiber diameter dependent expression that relates fluorescence collected by a single fiber with diameter df that has been distorted by scattering at excitation and emission wavelengths, to the intrinsic fluorescence within the sampled turbid medium. For brevity, the quantity will be used throughout this manuscript to refer to the expression
MC simulations investigated the relationship between single fiber fluorescence and variations in μ′s, initially specified as equivalent at λx and λm, and varied over the range [0.1 – 8.0] mm−1. Figure 1 A and B shows collected by single fiber probes with df [0.2 – 1.0] mm. These data show a fiber-diameter specific nonlinear relationship between and μ′s. Inspection of data sampled by the df = 0.2 mm fiber shows a 60% decrease in intensity as μ′s increases across the investigated range. However, the df = 1.0 mm fiber shows an initial decrease in of 25% as μ′s increases from 0.1 to 0.5 mm−1, and then doubles in intensity as μ′s increases from 0.5 to 8 mm−1.
Inspection of the fiber-diameter specific vs. μ′s profiles led to the identification of two dimensionless transformations that are important for interpretation of the data. First, transformation of the abscissa to dimensionless reduced scattering, given as the product μ′sdf, shifted the data on the x-axis such that the minimum values for each fiber specific profile aligned at the μ′sdf value of 0.5; the effect of this transformation is clearly shown in Figures 1 C and D. Second, expression of the ordinate as the dimensionless ratio of brought measurements from different fiber diameters onto an overlapping profile; this observed proportionality between fiber diameter and collected fluorescence is consistent with previous analysis of SFFL . Figures 1 E and F show the resulting dimensionless relationship between and μ′sdf that is observed for measurements from all investigated fiber diameters; there exists more than a factor of 2 variation in the observed magnitude of across the investigated μ′sdf range. These data exhibit a distinct U-shaped profile characterized by two phases: (1) for small μ′sdf values (μ′sdf < 0.5), decreases in response to increases in μ′sdf, and (2) for larger μ′sdf values (μ′sdf > 0.5), increases in response to increases in μ′sdf. This bi-phasic behavior is consistent with previous observations of the influence of scattering on fluorescence collected at or near the source [20, 21]; the underlying mechanisms of these phases and their respective dependence on scattering parameters is described in detail in Section 4.
The data investigated in Figure 1 are for the case μ′s(λx) = μ′s(λm); however, in tissue, μ′s(λ) is understood to follow a wavelength-dependent expression (e.g. Mie or Rayleigh approximations) such that μ′s(λx) > μ′s(λm). MC simulations were used to investigate for the case of independent variation of μ′s(λx) (range: [0.1 – 8.0] mm) and μ′s(λm) (specified as μ′s(λm) = [0.25,0.5,0.75,1.0]×μ′s(λx)). Figures 2 A and B show linear and log representations of the full data set plotted vs. μ′s(λx)df. Here, stratification of measurements at μ′s(λx)df values are attributable to the influence of μ′s(λm) on the collected intensity. These data show clear deviation of from the smooth curve displayed in Figures 1 E and F due to the independent influence of both μ′s(λx) and μ′s(λm) on SFFL intensity.
In tissue the exact form and wavelength-dependence of the PF is not well characterized. This study utilized a subset of MC simulations to investigate in detail the influence of PF on , as described in Section 2.2. The showed minimal influence from variation among different phase functions, with < 3% variation between values returned from the 19 simulated PFs at each of the dimensionless reduced scattering values (data not shown). For simulations specifying different PFs at λx and λm, the simulated values showed no observable difference if the PF were interchanged between the wavelengths. These results demonstrate that SFFL is insensitive to the form of the PF for all investigated scattering properties and fiber diameters.
This study utilized a subset of MC simulations to investigate in the influence of fiber NA on , as described in Section 2.2. Simulated data showed that the effect of fiber NA on is well approximated by an NA2 proportionality, with < 5% mean residual error between estimates of FSF measured by fibers of NA= [0.22] and NA= [0.1,0.4] in the investigated scattering range (data not shown), with increasing deviations associated with decreasing dimensionless reduced scattering values.
MC simulations were used to investigate the dependence of optical sampling depth, excitation fluence, and emission escape probability within the sampled volume on μ′s(λx) and μ′s(λm) within the sampled medium; these quantities were calculated as described in Section 2.1. Figure 3A shows a dimensionless description of effective optical sampling depth, given here as ZMC/df, plotted vs. μ′s,avgdf, with μ′s,avg calculated as the average of μ′s(λx) and μ′s(λm) for each measurement. These ZMC/df data exhibit a power law that shows a decreasing relationship with increasing μ′s,avgdf, resulting in a 10-fold decrease over the investigated μ′s,avgdf range. This relationship is well-characterized by Eq. (11); fitting this equation to these data yielded estimated values for A2 = 0.71 ± 0.01 and A3 = 0.36 ± 0.01,and resulted in accurate estimates of ZMC/df over the full range of investigated μ′s,avgdf values (r = 0.996); model predictions are visualized by the solid black line on the plot.
MC simulations also returned scalar metrics representative of effective excitation fluence and effective emission escape probability within the optically sampled volume. Figure 3B displays vs. μ′s(λx); these data show that nonlinearly depends on μ′s(λx)df, with an observed 2.5-fold increase across the investigated range. The form of the observed relationship is empirically described by Eq. (12); fitting the data to this model yielded estimated parameters of B1 = 0.88 ± 0.01 and B2 = 0.27 ± 0.03 and resulted in an accurate description of the simulated data (r = 0.977); fit quality is visualized by the model estimated black line. Figure 3C displays vs. μ′s(λm)df ; these data show a nonlinear dependence on μ′s(λm)df, with a 2.9-fold increase in the likelihood of collection associated with increasing μ′s(λm)df. This relationship is described by Eq. (13); the model fit to these data provides estimated parameter values of C1 = 0.12 ± 0.07, C2 = 0.10 ± 0.04, and C3 = 2.22 ± 0.61 and returns accurate estimates of (r = 0.967) over the investigated range of μ′s(λm). These results indicate that the Eqs. (12) and (13) describe both the magnitude and dynamic trends of the respective dependencies of vs. μ′s(λx)df and vs. μ′s(λm)df.
Figure 4 shows simulated by the MC model vs. estimated by the fit of Eq. (14). Here the estimated parameter values of ζ1 = 0.0935 ± 0.003, ζ2 = 0.31 ± 0.01, and ζ3 = 1.61 ± 0.05 resulted in the minimum weighted residual error between simulated and model-estimated values. The model estimates were strongly correlated with simulated outputs, with the quality of the fit given by the Pearson correlation coefficient of r = 0.991 and displayed by the proximity of the data points to the plotted line of unity. The mean absolute residual between simulated and model estimated values is < 3% and all data points have a mean residual error that is < 10% of the simulated value. Figures 5A and B show simulated and model estimated vs. μ′s(λx)df ; this plot visualizes the capability of the model to describe the influence of both μ′s(λx) and μ′s(λm) on the collected fluorescence intensity. These results indicate that Eq. (14) provides an accurate description of the SFFL intensity over a wide range of μ′s(λx), μ′s(λm), and df, and is valid for all investigated forms of the PF.
This study utilizes a Monte Carlo model to characterize the relationship between the fluorescence intensity collected by a single fiber (FSF) and the scattering properties within an optically sampled turbid medium. Simulated data were used to identify a relationship between dimensionless fluorescence intensity, , and dimensionless reduced scattering. We found that the collected fluorescence does not scale exclusively with dimensionless reduced scattering at the excitation wavelength, nor with dimensionless reduced scattering at the emission wavelength; rather it shows a more-complicated dependence on the reduced scattering coefficients at both wavelengths. These data were used to develop a semi-empirical model that expresses as the product of an effective sampling volume, and the effective excitation fluence and the effective escape probability within the effective sampling volume. The influence of scattering properties on each of these components was identified and mathematically described using simulation outputs. The semi-empirical model of accurately describes simulated fluorescence intensities over a wide range of biologically relevant scattering properties.
The fluorescence model, given in Eq. (14), utilizes empirical functions to represent the individual components of the SFFL measurement, including ZMC/df, , and . This approach provides insight into the mechanisms underlying the bi-phasic relationship observed between and the dimensionless reduced scattering coefficient, as visualized in Figure 1F. For μ′sdf < 0.5, denoted as phase (1), increases in μ′s at either λx and λm result in a decrease in . In this scattering region, the average depth of origin for collected fluorescence (and in turn the sampling depth) follows a similar trend, while the effective excitation fluence and effective emission probability are relatively insensitive to changes in this dimensionless scattering region; these trends are visualized in Figure 3. Figure 6A shows data in phase (1) following a smooth and continuous dependence on μ′s,avgdf ; this scattering dependence is shared by ZMC/df. These observations suggest that the left hand side of the profile is dominated by volume effects. Here, collected fluorescent photons originate from relatively deep locations in the medium, and an increase in μ′s,avg represents an impediment to light transport (either for excitation or emission photons), resulting in a decreased tissue volume optically sampled, and a reduced collected intensity. Conversely, for μ′sdf > 0.5, denoted as phase (2), increases in μ′s at excitation or emission wavelengths result in an increase in . In this scattering region, increases as μ′s(λx) increases, as shown in Figure 3B. This observation is attributable to the fluence ’build-up’ within the turbid medium near the fiber-tip for increasing scattering . Also in this scattering region, increases in response to an increase in μ′s(λm), as shown in Figure 3C. This phenomenon can be understood as follows. Scattering at the emission wavelength has two counteracting effects on fluorescence collection. First, fluorescent photons traveling towards the detecting fiber may be scattered away from the fibertip, decreasing the SFFL signal. This attenuating effect will be more pronounced for fluorescent photons that are emitted from relatively deep locations within the sample, while for photons originating close to the fibertip the attenuation due to scattering is expected to be small due to the small path traveled to the fibertip. Second, fluorescent photons traveling away from the detecting fiber may be backscattered towards the fibertip, increasing the SFFL signal. The balance of these counteracting effects will be depth dependent; photons originating from large depths inside the medium are expected to suffer more from attenuation due to scattering than benefit from fluorescence backscattering, while the opposite is true for photons originating close to the fibertip. Since for high scattering coefficients the effective sampling depth is relatively small, the net effect is that the benefit from fluorescence backscattering outweighs the attenuation of fluorescence due to scattering, resulting in an increase in effective escape probability with increasing scattering coefficient. data in phase (2) were observed to smoothly follow a dimensionless reduced scattering parameter dependent on the harmonic average of μ′s(λx) and μ′s(λm); this relationship was gained from inspection of the dependence of and on the reduced scattering coefficient. Specifically, the harmonic average reduced scattering coefficient is given as
These observations suggest that the right hand side of the profile (phase (2)) is dominated by a combination of excitation fluence build-up and increased fluorescence escape probability close to the fibertip for increasing reduced scattering coefficients at excitation and emission wavelengths, respectively. These observations are consistent with mechanisms that were previously proposed, but not explicitly investigated, in studies of localized  or single fiber measurements  of fluorescence.
The semi-empirical model developed in this study provides a method to return scattering-independent FSF quantities provided that μ′s(λx) and μ′s(λm) are determined, e.g. from a white-light reflectance measurement. This approach is in contrast to other techniques that utilize raw reflectance to correct raw fluorescence for the influence of scattering properties. Such an approach is not appropriate for single fiber measurements, because reflectance intensities collected by single fibers (RSF) are not only sensitive to μ′s, but (in contrast to SFFL) are also heavily influenced by the PF [24, 31, 32]. Due to this difference in PF dependence of FSF and RSF, the ratio of these two quantities will also be PF dependent. The magnitude of this dependence can best be appreciated by considering RSF measurements of two (hypothetical) turbid media with μ′s values of 0.5 and 2.0 mm−1, both with the same intrinsic fluorescence, and measured by a fiber with df = 1.0 mm. If the PF within the two media were varied from γ = 1.9 to γ = 1.4 (a change that would increase the likelihood of large-angle scattering events), the resulting RSF would increase by a factor of 2.3 for μ′s = 0.5 mm−1 and a factor of 1.4 for μ′s = 2.0 mm−1 [24, 25]. For a smaller fiber of df = 0.2 mm, the effects are amplified to factors of 3.1 and 1.5 for each respective case. Importantly, the variation in PF would have a negligible effect on the raw FSF ; such a difference in sensitivity to PF is attributable to the isotropic release of emission photons during propagation of fluorescent light. In contrast to RSF, which relies on the likelihood of forward directed incident light to undergo a large-angle scattering event (defined by the PF), the isotropic release of a fluorescent photon greatly reduces the sensitivity of FSF to PF. Therefore, for single fiber measurements (and likely other geometries which collect light close to the source fiber), a fluorescence correction algorithm that utilized a ratio of FSF and RSF could result in inaccurate estimation of intrinsic fluorescence by up to a factor of > 3 for small dimensionless scattering values.
The PF-specific analysis presented in this study indicates that quantitative analysis of SFFL requires determination of μ′s(λx) and μ′s(λm) independent of PF. This could be achieved using a multi-diameter SFR measurement, as described recently by our group [25,26]. The MDSFR approach utilizes the γ-specific RSF vs. μ′sdf relationship for measurements using multiple fibers at each investigated wavelength. By specification of a background scattering model within the sampled tissue (e.g. Mie and or Rayleigh scattering) it is possible to determine μ′s and γ across the a range of wavelengths. Moreover, this calculation can be made in the presence of absorption from tissue chromophores, requiring only specification of the basis set of absorbing constituents and their respective specific absorption coefficients. This multi-fiber approach can be executed using as few as two optical fibers with different diameters ; moreover, such a device can easily be developed to sample both RSF and FSF. The combined multi-diameter SFR and SFFL would return paired local measurements of fluorescence and tissue optical properties within the same (shallow) sampling volume. Such a technique has the potential to provide clinically useful information for tissue diagnostics and monitoring of administered therapies. The localized measurement volume would allow quantitative characterization of heterogeneities in the spatial distribution of an administered fluorescent compound; this may be advantageous compared with a volume-averaged metric gained from diffuse optical measurements. Moreover, the measurement volume can be selected at a specific area of interest (e.g. in the center of an identified malignant area, or on the border between suspicious and normal tissue). This multi-fiber approach faces challenges that must be properly assessed, including proper identification of background scattering models for determination of μ′s(λ), the influence of μa on FSF, and the influence of heterogeneities on both RSF and FSF ; ongoing studies are investigating these issues.
In order to appropriately utilize the semi-empirical model of SFFL presented in this study, it is important to consider the assumptions and approximations utilized in its development. The mathematical modeling approach utilized in this study represents the collected fluorescence intensity in terms of the product of three factors contributing to fluorescence that were extracted from Monte Carlo models outputs; these relationships are presented in the transition from Equation 6 to 10. A critical assumption of this modeling approach is that the effective scalar values for these components are representative of the more complicated 2-D maps of these properties. The empirical models of each of the components expresses a high quality of fit, providing evidence that this assumption is reasonable. Another important point of this study is the specific investigation of a single optical fiber in contact with a turbid medium; the exact form of the expressions governing light transport have been defined for this geometry. While the approach to modeling SFFL utilized here is extensible to modifications in measurement geometry, it is important to note that changes to the geometry will result in changes to the excitation and emission light distributions, and will require assessment of the appropriateness and accuracy of the specified model structures. Such modifications include interstitial placement of the fiber optic in the sampled medium, or placement of the fiber optic into a probe face surrounded by epoxy, metal, or other optical fibers; ongoing work is investigating these influences. Another important consideration is that this study characterized the scattering dependence of FSF, and did not consider background absorption effects. Absorption within the sampled medium, at both excitation and emission wavelengths, is expected to have a substantial influence on the raw fluorescence intensity collected and the volume probed during measurement. Further complicating matters, the magnitude of the absorption attenuation is expected to be heavily influenced by the paired scattering properties at excitation and emission wavelengths. An ongoing study will characterize the influence of absorption on the individual components of the SFFL model. Additionally, the MC model utilized in this study was validated by comparison with model returned outputs reported in the literature; future work will conduct experimental validation in optical phantoms.
In summary, the current study utilized MC simulations to investigate the influence of scattering properties on fluorescence intensity collected by a single fiber probe. Simulated data were used to identify an underlying dimensionless relationship between fluorescence intensity and dimensionless reduced scattering. Results indicate that the mathematical model of FSF is valid over a wide range of reduced scattering coefficients, in the range μ′s(λx) [0.1 – 8] mm−1 and μ′s(λm) [0.25 – 1] × μ′s(λx), and scattering phase functions (P(θ)), with both Henyey-Greenstein and Modified Henyey-Greenstein forms with anisotropy in the range 0.5 – 0.95 and γ [1.4 – 1.9], and a wide range of fiber diameters (df [0.2 – 1.0] mm). The model accurately estimates FSF given μ′s(λx) and μ′s(λm), and is insensitive to the anisotropy and higher order moments of the PF. Results indicate that correction for the influence of scattering on FSF requires estimation of scattering optical properties from a paired measurement of white-light reflectance.