Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2964859

Formats

Article sections

- Abstract
- I. Introduction
- II. Materials and Methods
- III. Results
- IV. Discussion and Conclusions
- References

Authors

Related links

IEEE J Sel Top Quantum Electron. Author manuscript; available in PMC 2011 January 1.

Published in final edited form as:

IEEE J Sel Top Quantum Electron. 2010; 16(3): 619–626.

doi: 10.1109/JSTQE.2009.2032666PMCID: PMC2964859

NIHMSID: NIHMS164622

Vladimir Turzhitsky, Department of Biomedical Engineering, Northwestern University, Evanston, IL 60208 USA (Phone: 847-491-7167; fax: 847-491-4928; Email: ude.nretsewhtron.u@tv).

See other articles in PMC that cite the published article.

Low-coherence enhanced backscattering (LEBS) is a technique that has recently shown promise for tissue characterization and the detection of early pre-cancer. Although several Monte Carlo models of LEBS have been described, these models have not been accurate enough to predict all of the experimentally observed LEBS features. We present an appropriate Monte Carlo model to simulate LEBS peak properties from polystyrene microsphere suspensions in water. Results show that the choice of the phase function greatly impacts the accuracy of the simulation when the transport mean free path (ls*) is much greater than the spatial coherence length (L_{SC}). When ls* < L_{SC}, a diffusion approximation based model of LEBS is sufficiently accurate. We also use the Monte Carlo model to validate that LEBS can be used to measure the radial scattering probability distribution (radial point spread function), *p(r)*, at small length scales and demonstrate LEBS measurements of *p(r)* from biological tissue. In particular, we show that pre-cancerous and benign mucosal tissues have different small length scale light transport properties.

Enhanced backscattering (EBS) is one of the two self-interference phenomena observed when a sample is illuminated with a coherent light source. The enhanced backscattering effect, also known as coherent backscattering, manifests as an enhancement of intensity in the backward direction of the incident light. This effect was initially tied with weak localization of electrons, and later extended to optical waves[1]. Further theoretical and experimental work uncovered that EBS can be used to measure the optical properties of diffusely scattering media[2, 3]. The width of the coherent backscattering peak was found to be proportional to λ/ls*, where λ is the wavelength and ls* is the transport mean free path. However, due to the inherent challenges in this method, the utilization of this phenomenon for characterizing biological tissue has been hampered until recent years.

The second self-interference phenomenon, known as speckle, is one of the primary impediments to measuring the enhanced backscattering in tissue. The large amplitude and ubiquity of the speckle pattern overwhelms the coherent backscattering signal in solid samples. This is why early studies primarily utilized liquid suspensions of particles such that the Brownian motion inside the sample provides a large number of ensembles, all of which are averaged within a single acquisition of the detector. The lack of Brownian motion in biological tissue requires the averaging of many ensembles by moving the entire sample. Furthermore, the narrow line shape of the coherent backscattering peak from tissue requires high angular resolution detection capable of resolving signals on the order of 0.001 degrees.

The extension of enhanced backscattering into the partial spatial coherence regime was first explored by Kim *et al* [4] and provides several advantages which make the technique amenable for tissue characterization. The technique overcomes the difficulties of measuring enhanced backscattering from tissue by broadening the coherent backscattering peak by several orders of magnitude and virtually eliminating the speckle noise.

The new method also offers further advantages that are essential for the biomedical application of enhanced backscattering. The reduced spatial coherence results in a depth selective signal by removing the contribution of light at longer scattered radii[5, 6]. The addition of spectrally-resolved detection allows for obtaining the interference phenomenon from a range of wavelengths, adding an extra dimension for more complete characterization of tissue properties[7]. Thus far, Low-coherence Enhanced Backscattering (LEBS) has shown promise for the early detection of colorectal cancer[8] and pancreatic cancer[9, 10] by the quantification of superficial tissue layer changes that are unapparent from conventional histology.

Earlier work in characterizing LEBS has revealed that the dependencies of the peak width and enhancement factor on sample properties are more complicated than the those of EBS[4]. For example, the LEBS peak width was found to have three ls* dependent regimes with distinctly different behaviors for varying values of the anisotropy factor (g) [7], this being very different from the λ/ls* relationship seen in EBS. The enhancement factor, defined as the relative height of the peak, was also found to be dependent on the mean free path of the medium. On the other hand, the enhancement factor of EBS has no dependence on the mean free path of the sample.

These experimental findings combined with the demonstrated diagnostic potential of the LEBS measurement have motivated several groups to develop Monte Carlo[11, 12] and theoretical[6, 12] models for the simulation and understanding of LEBS. Although these Monte Carlo models have given insight into important properties such as the penetration depth of the LEBS signal [6], a good correspondence between the experimental peak properties, including peak width and enhancement factor, with a simulated result has not been demonstrated until now.

In the following sections, this paper will discuss the experimental and modeling methodologies required to obtain agreement between simulated and experimental peak properties. We will present experimental and Monte Carlo results that illustrate the dependence of the LEBS peak width and enhancement factor on optical properties and the spatial coherence length. We will then present a diffusion based model for LEBS, which can be accurate in certain conditions. Finally, we will apply LEBS as a technique for measuring the scattering probability distribution, *p(r)*, at small length scales.

A random medium illuminated by a uniform plane wave will backscatter a fraction of the reflected power in the reverse direction of the incident beam. The principle of the reversibility of light paths states that there must also be an equivalent fraction of the incident beam that traverses an identical path in the opposite direction. The light from these time-reversed paths will have an identical phase and therefore constructively interfere in at the exact backward direction when the angular distribution of backscattered light is observed. For angles away from the backward direction, the interference depends on the phase difference which intern depends on the backscattering angle and location at which the partial waves exit the medium.

The resulting peak can be described as a summation of interference patterns observed from a series of Young's double pinhole experiments, where the distances between the pinholes are the exit locations of all of the possible time-reversed path pairs. The shape of the enhanced backscattering peak in a nonabsorbent medium therefore depends on the point spread function of light at the tissue surface in the 0° backward direction (i.e. 180° in respect to the direction of the incident light), termed *p(r)*, with a two dimensional Fourier transform relationship:

$$I\left(\theta \right)=\int \int p\left(\overrightarrow{r}\right)\phantom{\rule{thinmathspace}{0ex}}\mathrm{exp}\left(\mathit{ik}\overrightarrow{r}\phantom{\rule{thinmathspace}{0ex}}\mathrm{sin}\left(\theta \right)\right){d}^{2}r,$$

(1)

where *θ* is the angle between the incident wavevector and the backscattered wavevector and k is the wavenumber. The sinusoidal oscillations obtained in Youngs double pinhole experiment can be obtained by substituting *p(r)* with a pair of delta functions into (1). Thus, the width of the EBS cone is inversely related to the average backscattering exit radius of the beam. This average distance is proportional to ls*. When the diffusion approximation is used to model the shape of p(r), the EBS peak shape is approximated by [2]

$${I}_{\mathit{EBS}}\left(\theta \right)\approx \frac{3}{8\pi}\left[\begin{array}{c}\hfill 1+\frac{2{z}_{0}}{{\mathit{ls}}^{\ast}}+\frac{1}{{(1+k\theta {\mathit{ls}}^{\ast})}^{2}}\times \phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{thinmathspace}{0ex}}\hfill \\ \hfill \phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{thinmathspace}{0ex}}\left(1+\frac{1-\phantom{\rule{thinmathspace}{0ex}}\mathrm{exp}(-2k\theta {z}_{0})}{k\theta {\mathit{ls}}^{\ast}}\right)\hfill \end{array}\right],$$

(2)

where z_{0} is the location of the trapping plane.

The requirement in equations (1) and (2) is that all path pairs incident on the scattering medium are fully coherent. Under partial spatial coherence illumination, the interference depends on the separation distance of the path pairs and the degree of spatial coherence, *c(r)*. In this case:

$${I}_{\mathit{LEBS}}\left(\theta \right)=\int \int p\left(\overrightarrow{r}\right)c\left(\overrightarrow{r}\right)\phantom{\rule{thinmathspace}{0ex}}\mathrm{exp}\left(\mathit{ik}\overrightarrow{r}\phantom{\rule{thinmathspace}{0ex}}\mathrm{sin}\left(\theta \right)\right){d}^{2}r.$$

(3)

When *p(r)* is radially symmetric, the two dimensional Fourier transform can be written as

$${I}_{\mathit{LEBS}}\left(\theta \right)={\int}_{0}^{\infty}2\pi \mathit{rp}\left(r\right)c\left(r\right){J}_{0}\left(\mathit{kr}\phantom{\rule{thinmathspace}{0ex}}\mathrm{sin}\left(\theta \right)\right)\mathit{dr},$$

(4)

where *J _{0}* is the zero order Bessel function of the first kind. A symmetric p(r) will result when the illumination has normal incidence on the sample and either unpolarized or circularly polarized illumination is used.

When the degree of spatial coherence extends much farther than ls*, the diffusion approximation that conventionally used to model EBS can be extended to model the LEBS signal. This can be done by either substituting the diffusion approximation p(r) into equation (2) or by applying the convolution theorem:

$${I}_{\mathit{LEBS}}\left(\theta \right)=\mathit{FT}[p\left(r\right)\cdot c\left(r\right)]={I}_{\mathit{EBS}}\left(\theta \right)\ast S\left(\theta \right),$$

(5)

where S(θ) is the Fourier transform of c(r) and the angular profile of the source from the sample surface. By substituting equation (2) into equation (5), a diffusion approximation for the LEBS signal is obtained.

The properties of *p(r)* uniquely characterize the scattering media and determine the shape of the peak. The shape of *p(r)* at r<ls* has not been investigated in tissue due to the lack of an effective technique to measure this signal, despite the fact that scattering properties in this region can be sensitive to characteristics of the phase function. From equations (3) and (5), *p(r)* can be found from the inverse Fourier transform of the angular profile of an LEBS peak:

$$p\left(r\right)={\mathit{FT}}^{-1}\left[{I}_{\mathit{LEBS}}\left(\theta \right)\right]\u2215c\left(r\right).$$

(6)

As *c(r)* approaches 0, the measurement becomes dominated by noise and is no longer accurate. The use of LEBS allows the weight of the collected signal to be concentrated at small values of r, yielding a more accurate measurement of *p(r)* for small r (i.e. r<ls*) than what is possible with any other available method.

A schematic of the experimental setup used in earlier LEBS studies has been described in previous literature [4, 5, 7]. A modified experimental setup that simplifies alignment and improves sampling was implemented in the experimental studies discussed in this paper. A schematic of the setup is shown in Figure 1. A broadband 450 W Xe Oriel light source (Newport Corp., Stratford CT) relayed with a series of lenses onto a high power tungsten aperture that is positioned within an aperture wheel. The aperture acts as a secondary source of variable size that is then collimated by a lens (f=200mm). The degree of spatial coherence and the spatial coherence length (L_{SC}) is determined according to the Van Cittert-Zernike theorem as [13]:

$$\begin{array}{cc}\hfill c\left(r\right)& =\frac{2{J}_{1}(r\u2215{L}_{\mathit{SC}})}{r\u2215{L}_{\mathit{SC}}}\hfill \\ \hfill & \phantom{\rule{thickmathspace}{0ex}}{L}_{\mathit{SC}}=\lambda \u2215\pi \alpha \hfill \end{array}$$

(7)

Where *J _{1}* is the first order Bessel function of the first kind, λ is the wavelength, and α is the angular extent of the source. In the described setup, α is determined from the ratio of the diameter of the secondary source to the focal length of the collimating lens L

Schematic of LEBS experimental instrument. A broadbaned 450 W Xenon source (S) is imaged onto an aperture of variable size which serves as a secondary source (SS). The size of the aperture is selected by appropriately positioning the aperture wheel. The **...**

The setup in Figure 1 has two advantages over the earlier LEBS systems: i) the alignment is not disrupted when the spatial coherence is changed by rotating the aperture wheel to a different position. 2) The LEBS peak is not sampled with a spectrometer slit, allowing for higher resolution imaging of the entire backscattered cone. As an example of one of the advantages of enhanced backscattering measurement with partially coherent light, Figure 2 illustrates the speckle reduction capability of LEBS with backscattering measurement obtained from dried paint. In Figure 2a a coherent backscattering measurement from the solid sample is shown to be completely submerged in speckle. An EBS peak begins to be visible in Figure 2b, where 30 measurements at various positions in the sample are averaged. For comparison, only a single measurement is required to see the LEBS peak under partial coherence illumination, as shown in Figure 2c and Figure 2d for L_{SC} = 160μm and L_{SC} = 20μm respectively.

Demonstration of speckle reduction in LEBS from a dried paint sample (ls* ~ 4μm). Panel a: EBS measurement obtained with Helium Neon (HeNe) laser. Panel b: average of 30 ensembles of EBS measurements collected with HeNe laser. The EBS signal begins **...**

The LEBS peak was characterized for a variety of illumination conditions from measurements on polystyrene microsphere suspensions in water. Microspheres of sizes ranging from 0.36μm – 0.65μm (Duke scientific) were used to obtain a range of g values between 0.71 – 0.86. The samples were made to approximate a semi-infinite media by maintaining the thickness and the diameter of the sample more than 50·ls*. These suspensions were diluted to varying concentrations in order to achieve values of ls* in the range of 20μm – 1100μm. Measurements were obtained at λ = 680nm and all optical properties were calculated with Mie theory. The image obtained with the CCD camera was background subtracted and normalized by the an inensity measurement collected from a spectralon white reflectance standard (Ocean Optics, Dunedin FL). The incoherent baseline was subtracted by fitting a plane to the signal at the periphery of the peak. The enhancement factor of the LEBS peak was therefore bounded between 0 and 1. The LEBS peak width was measured as the full width at half maximum of the signal after averaging over the azimuth angle.

In order to illustrate the potential of *p(r)* measurements from tissue with LEBS, an animal study utilizing 8 rats was performed. Four of these rats were treated with Azoxymethane (AOM), a carcinogen that induces cancerous changes typically around 20-40 weeks after treatment. The remaining 4 rats were injected with saline and served as controls. The rats were sacrificed 10 weeks after the injection and LEBS measurements were obtained for L_{SC} = 87μm from 4-5 sites per rat. The LEBS peaks in each group were averaged and the inverse Fourier transform relationship (Eq. (6) was used to calculate *p(r)* for values of r up to the point were the degree of spatial coherence dropped below the half maximum (~ 2.2·Lsc).

The low-coherence enhanced backscattering phenomenon was modeled with polarized light Monte Carlo simulations. A publicly available Monte Carlo code [14, 15] was used to simulate the scattering properties of polystyrene microsphere suspensions in water. The code was modified to store small angle backscattered rays from 0-2 degrees following previously reported methodology[11]. A finite range of angles is necessary to accumulate backscattered rays simulated with the Monte Carlo method, while it is also important to collect as close to 0° in the backward direction as possible to obtain an accurate result[11]. The 2° range was chosen because, as reported in [11], p(r) changes negligibly for angles of backscattering from 0°-2°, while it differs for larger angles. A refractive index of 1.59 was used for the spheres and 1.33 for the medium. Single scattering events were excluded from the stored *p(r)* distribution. The polarization-dependent phase function was calculated with Mie theory at λ=680 and sampled with an angular resolution of π/500 radians. The Mie phase function is more accurate for modeling scattering from spherical particles than the often used Henyey Greenstein phase function and allows for the simulation of polarization effects. The probability distribution for backscattered light was collected and equation (3) was applied to calculate the two dimensional LEBS signal. The simulation modeled a thick slab geometry where less than 3% of the intensity was transmitted to the bottom of the slab. An infinitely narrow vertically polarized incident beam was propagated through the simulated random scattering medium with the Stokes vectors being tracked and stored once the beam re-emerged to the top surface. For each simulation run, between 10^{8} and 10^{9} rays were stored to obtain the p(r) for all of the four Stokes vectors. Mueller matrix multiplication was then used to obtain the co-polarized distribution of p(r). The grid resolution was maintained < L_{SC}/20 and the maximum collected r was maintained > 10L_{SC}. This ensured that the error due to under sampling was small.

A comparison of the experimentally observed enhancement factor with the enhancement factor calculated from Monte Carlo simulation for a variety of spatial coherence lengths is shown in Figure 3. A constant multiplicative coefficient of 0.5 was used to adjust the Monte Carlo enhancement factor for effects such as the contribution of the single scattering fraction and a finite spot size. Note that the selected polarization contribution to p(r) is already taken into account in the Monte Carlo model by tracking the Stokes vectors and Mueller matrix multiplication of the collected signal. As seen in Figure 3a, the enhancement factor from experiment and simulation is higher for increasing L_{SC}. When the enhancement factor is plotted with L_{SC}μ_{s}*, as shown in Figure 3b, and the graphs collapse to a single monotonic relationship. This relationship can be approximated with a straight line for values of L_{SC}μ_{s}* < 1, which agrees with previous observed trends from LEBS Monte Carlo simulations[12]. In contrast, the enhancement factor measured with EBS has no relationship with the scattering coefficient. This result can be explained from considerations of the Fourier transform relationship of LEBS with p(r)·c(r). The LEBS enhancement factor is the integral of p(r)·c(r),*I _{LEBS}* (0) = ∫∫

Enhancement factor as a function of μ_{s}* for varying L_{SC} (g=0.86). In panel a, the enhancement factor is shown to have three distinct dependencies on μ_{s}* for three values of L_{SC}. In panel b, the three dependencies coincide when enhancement **...**

The LEBS peak width, defined as the full width at half maximum, is another important measurable parameter that has previously been shown to have fascinating characteristics for varying optical properties. Earlier experimental results suggest that the LEBS peak width may be used to measure sub-diffusion optical properties such as the anisotropy factor g, the average cosine of the phase function. Indeed, by using a more accurate phase function to model the particle scattering, agreement between experiment and simulation is attained for the first time (Figure 4). Figure 4a shows the inverse LEBS peak width as function of ls* for Monte Carlo simulation and experiment (g=0.86) for three values of L_{SC}. The relationship between the peak width and L_{SC} is demonstrated in Figure 4b by appropriately scaling the axes as α/width vs. ls*/L_{SC}, where α is the angular extent of the source and can be calculated from equation (4). It is important to recognize that the shape of *p(r)* is uniquely characterized by the α/width vs. ls*/L_{SC} dependence. In Figure 5b, the appropriately scaled width is shown for a variety of phase functions (g=0.71 – g=0.86). Although, the LEBS peak width is related to λ/L_{SC} (i.e. α) at large ls*/L_{SC} as predicted by [12], the LEBS width is also highly dependent on sub-diffusion characteristics of the sample when ls*/L_{SC} >> 1. These characteristics are determined by the shape of *p(r)* at small values of r (r/ls* < 1).

LEBS width relationship with ls* and L_{SC}. Panel a shows the inverse width dependence on spatial coherence length (g=0.86) for Monte Carlo (lines) and experiment (symbols). Panel b shows the data from panel a when the width is scaled by the angular extent **...**

Figure 6a shows Monte Carlo simulation results compared with the diffusion approximation for LEBS, obtained from equations (2) and (5). The diffusion approximation for LEBS shows good agreement for ls* < 2L_{SC}, but becomes far less accurate for ls*/L_{SC} >> 1. This is because the diffusion approximation is a poor model of p(r) for r < ls*. Figure 6a also shows a comparison between the commonly used Henyey Greenstein phase function based Monte Carlo and the Mie phase function based Monte Carlo for the same value of g. This demonstrates that the shape of the phase function plays an important role in the accuracy of the simulation. Finally, figure 6b, shows the same Monte Carlo data observed in Figure 6a when the LEBS peak is sampled with a spectrometer based system that has been used to collect LEBS signals in prior work [4, 5, 7, 16]. This is simulated by averaging the peak in one angular dimension over a spectrometer slit width of 1.3α. Averaging over a slit width amplifies the dependence of the peak width on the shape of the phase function for ls* >> L_{SC}.

By applying equation (6), a measurement of *p(r)* is possible by obtaining the inverse Fourier transform of the LEBS intensity distribution and dividing by the degree of spatial coherence. This was accomplished with the discrete two-dimensional inverse Fast Fourier Transform in Matlab. The obtained p(r) is a sharply decaying function that is better visualized after integration over the azimuth angle, resulting in the more commonly shown dependence of 2*πrp(r)*. In principle, the calculation is valid from smallest value of r (determined by the collected angular range) up to the first zero of the spatial coherence function. Practically, the presence of noise limits the smallest and longest reliably measured r. The minimum reliable value of r is accurate from a resolution limit that is determined by the maximum angle of the I_{LEBS}(θ) that is above the noise level. The maximum reliable value of r is the point where the noise level overcomes the degree of spatial coherence. Figure 7 shows 2πrp(r) curves experimentally measured with LEBS compared with curves obtained from Monte Carlo simulations. The 2πrp(r) curves for two different particle sizes (D=0.36μm and D = 0.65μm) have different behaviors for r<ls* (Figure 7a). For larger values of r, the governing parameter is ls* and the curves are indistinguishable. Figure 7b shows the changes in p(r) when ls* is varied. As anticipated, an increase in ls* broadens p(r) and reduces the probability at any sampled radius.

Experimental measurements of p(r) compared with Monte Carlo Simulation. In panel a, a comparison for two different values of g and similar values of ls* is made between Monte Carlo and Experimentally collected data. In panel b, three values of ls* are **...**

Figure 8 shows measurements of 2πrp(r) obtained from pre-cancerous colonic rat tissue compared with healthy colonic rat tissue. The measurements were taken 10 weeks after either AOM or saline injection. In this animal model, colonic tumors normally develop 40 weeks after AOM injection. The colonic carcinogenesis process begins in the mucosa, a tissue layer that is <1mm thick. Alterations in p(r) at small radii suggest changes in the structure of the mucosal colonic tissue. Larger values of r correspond to more deeply penetrating light. In this case, the values of r are limited to small length-scales such that the majority of the signal is confined to the epithelium. This isolation of epithelial signal is of great value for the detection of early pre-cancer because most cancers originate due to the accumulation of mutations in epithelial tissue.

We have presented experimental data that validates the Monte Carlo model of LEBS for a broad range of optical properties. Polystyrene microspheres in water were used as samples, which where then simulated with Monte Carlo code that implements a phase function calculated with Mie theory and tracks the polarization state of each scattering event. The results obtained from the LEBS peak width, enhancement factor, and p(r) show unprecedented agreement between experiment and simulation. These results suggest that LEBS measurements are highly sensitive to the shape of the phase function, contrasting the case of EBS, which is well modeled by the diffusion approximation. In particular, Figure 5b illustrates the sensitivity of the LEBS width measurement to the anisotropy coefficient in the LEBS regime when ls*>>L_{SC}. As shown in Figure 7, when the LEBS signal is converted to p(r) for two different bead sizes (different phase functions), the difference in scattering at small values of r becomes apparent, while p(r) at r>ls* is insensitive to the shape of the phase function and simply determined by ls*. The agreement between the experimentally measured p(r) and that from Monte Carlo indicates that the measurement is accurate. In Figure 7, the p(r) measured using LEBS is shown from biological tissue measurements of control and carcinogen-treated rats. This tissue measurement of p(r) is the first validated measurement of p(r) using the enhanced backscattering technique.

A diffusion approximation for LEBS was presented and illustrated in Figure 7. The LEBS diffusion approximation was reliable at ls*/L_{SC} < 2, but quickly becomes inaccurate at longer values of ls*. Figure 7 also demonstrates the importance of the shape of the phase function. An identical value of g simulated with the Henyey Greenstein phase function results in different LEBS peak widths than those obtained from simulations with the Mie phase function. This underscores the importance of choosing an accurate phase function for LEBS simulation with Monte Carlo.

It is possible that the simulation results may be further improved by incorporating other considerations from the experimental setup. The Monte Carlo simulation models the propagation of light at a single wavelength, while in the experiment, the LCTF has a 7nm bandwidth. The interface between the air and water was not considered, and this may also improve the results [17]. Finally, it is possible that experimental errors such as the variation in shape or quality of the polystyrene microsphere suspensions limit the accuracy of the Mie phase function. In this case, the exact phase function must be measured experimentally in order to validate the Monte Carlo model of LEBS to a higher level of accuracy. One additional challenge in obtaining experimental LEBS measurements is an accurate fitting of the baseline. If angles that are too close to the peak are used to fit the baseline, over-subtraction may reduce the enhancement factor and narrow the peak width. On the other hand, noise considerations make baseline fitting at very large angles impractical. Finally, the assumption that the intensity throughout the secondary source is completely incoherent and perfectly uniform may be challenged. Any non-uniformity in the intensity distribution in the pinhole would result in an increased spatial coherence length. This effect would make the experimentally obtained peak widths have slightly narrower widths and higher enhancement factors and would be more problematic with larger apertures (shorter L_{SC}).

In summary, we have shown that experimental values of the LEBS enhancement factor and width are well modeled by the Monte Carlo method for a large range of optical properties. We have shown that the accuracy in the choice of the phase function has a large impact on the LEBS peak width when ls* >> L_{SC}. On the other hand, when ls* < L_{SC}, the choice of the phase function is not important and the use of the diffusion approximation to model LEBS behavior is sufficient. Because LEBS peak is directly related to p(r), this is consistent with a notion, which has been long held in the field of tissue optics: for r>ls*, light transport can be accurately modeled by the diffusion approximation (which depends on ls* only and is independent of the phase function). On the other hand, the shape of the phase function has a profound effect for short, sub-diffusion length scales of light propagation, i.e. r<<ls*. These results also suggest that the LEBS peak width may be used to measure parameters of the phase function such as g. Moreover, we have demonstrated that the calculation of *p(r)* at small length-scales is feasible by performing an inverse Fourier transform of the LEBS signal. The *p(r)* measurement from polystyrene microspheres showed good agreement between experiment and simulation. We then applied the p(r) calculation to an animal model demonstrating that there may be potential differences in the p(r) at small length-scales between pre-cancerous and healthy colonic rat tissue.

This work was supported in part by R01 EB003682, R01 CA112315, R01 CA128641.

**Vladimir Turzhitsky** received the B.S. in biomedical engineering from Rensselaer Polytechnic Institute, Troy, NY in 2004, and his M.S. and Ph.D. in biomedical engineering from Northwestern University, Evanston IL in 2006 and 2009, respectively. He is currently a postdoctoral fellow in the biomedical engineering department at Northwestern University. He has been involved in work leading to the development and characterization of Low-coherence Enhanced Backscattering Spectroscopy (LEBS), polarization gating spectroscopy for the early detection of colon and pancreatic cancer. His research interests include the understanding of light scattering phenomena under partial coherence illumination conditions and the application of these phenomena as cancer detection and screening techniques.

**Jeremy D. Rogers** received the B.S. degree in physics from Michigan Technological University in 1999, and the M.S. and Ph.D. degrees in Optical Sciences from the College of Optical Sciences, University of Arizona, Tucson, in 2003 and 2006, respectively.

He has previously worked at NanoSciences Corporation, Oxford, CT, and interned at Sandia National Laboratories, Albuquerque, NM. He is currently a post-doctoral Fellow in Biomedical Engineering at Northwestern University, Evanston. His research interests include theoretical and numerical modeling of light scattering as well as lens design and instrumentation for basic investigation and application of optical detection techniques.

**Nikhil N. Mutyal** received the B.Sc. and M.Sc. Degree in Chemistry from the Indian Institute of Technology, Kanpur. His master's thesis was in the use of endogenous NADH florescence markers for cancer diagnosis. He is currently pursuing his Ph.D. in biomedical engineering at Northwestern University. His research interests include finding new markers and therapeutics for the diagnosis and treatment of cancer. His current research focuses on finding the biological origins of changes in optical properties in the progression of cancer and developing instrumentation (fiber optic probes) to detect these changes using Low-coherence Enhanced Backscattering Spectroscopy.

**Hemant K. Roy** received the B.S. degree (summa cum laude) in molecular biology from Vanderbilt University, Nashville, TN, in 1985 and the M.D. degree with distinction from Northwestern University, Evanston, IL, in 1989. He completed his internal medicine training at Beth Israel Hospital/Harvard Medical School, Cambridge, MA, in 1992 and gastroenterology training at the University of Chicago, Chicago, IL, from 1992 to 1995. Dr. Roy's research interest focuses on colon cancer prevention. Specifically, his laboratory investigates non-cyclooxygenase mechanisms for apoptosis from nonsteroidal anti-inflammatory agents and has been one of the pioneers the use of polyethylene glycol in chemoprevention. His translational group works on development of novel colorectal cancer screening strategies. He currently is a clinical associate professor of medicine at University of Chicago Pritzker School of Medicine. He is at NorthShore University HealthSystems in Evanston Illinois where he serves Director of Research and Vice Chair of the section of Gastroenterology and Duckworth Chair of Cancer Research.

**Vadim Backman**, Ph.D. is a Professor of Biomedical Engineering at Northwestern University, Program Leader in Cancer Bioengineering, Nanotechnology and Chemistry at the Robert H. Lurie Comprehensive Cancer Institute and a Member of the Professional Staff in the Division of Gastroenterology at NorthShore University Health Systems. He received a Ph.D. in Medical Engineering from Harvard University and Massachusetts Institute of Technology in 2001. Dr. Backman received numerous awards including being selected as one of the top 100 young innovators in the world by the Technology Review Magazine and the National Science Foundation CAREER award. Dr. Backman's research is at the interface of physics, biology and medicine. His research interests include biomedical optics, spectroscopy, microscopy, development of analytical approaches to describe light transport in biological media, and optical microscopy for nanoscale cell analysis. Dr. Backman is also active in translational research focused on bridging these technological and biological innovations into clinical practice.

Vladimir Turzhitsky, Department of Biomedical Engineering, Northwestern University, Evanston, IL 60208 USA (Phone: 847-491-7167; fax: 847-491-4928; Email: ude.nretsewhtron.u@tv).

Jeremy D. Rogers, Department of Biomedical Engineering, Northwestern University, Evanston, IL 60208 (Email: ude.nretsewhtron@sregordj).

Nikhil N. Mutyal, Department of Biomedical Engineering, Northwestern University, Evanston, IL 60208 (Email: ude.nretsewhtron.u@2102laytumlihkin).

Hemant K. Roy, Department of Internal Medicine, Northshore University HealthSystems, Evanston, IL 60201 (Email: ude.nretsewhtron@ yor-h)..

Vadim Backman, Department of Biomedical Engineering, Northwestern University, Evanston, IL 60208 (Email: ude.nretsewhtron@namkcab-v).

1. Vanalbada MP, Lagendijk A. Observation of Weak Localization of Light in a Random Medium. Physical Review Letters. 1985;55:2692–2695. [PubMed]

2. Akkermans E, Wolf PE, Maynard R. Coherent Backscattering of Light by Disordered Media - Analysis of the Peak Line-Shape. Physical Review Letters. 1986 Apr 7;56:1471–1474. [PubMed]

3. Vandermark MB, Vanalbada MP, Lagendijk A. Light-Scattering in Strongly Scattering Media - Multiple-Scattering and Weak Localization. Physical Review B. 1988 Mar 1;37:3575–3592. [PubMed]

4. Kim YL, Liu Y, Turzhitsky VM, Roy HK, Wali RK, Backman V. Coherent backscattering spectroscopy. Optics Letters. 2004 Aug 15;29:1906–1908. [PubMed]

5. Kim YL, Liu Y, Turzhitsky VM, Wali RK, Roy HK, Backman V. Depth-resolved low-coherence enhanced backscattering. Optics Letters. 2005 Apr 1;30:741–743. [PubMed]

6. Subramanian H, Pradhan P, Kim YL, Backman V. Penetration depth of low-coherence enhanced backscattered light in subdiffusion regime. Physical Review E. 2007 Apr;75 [PubMed]

7. Kim YL, Liu Y, Wali RK, Roy HK, Backman V. Low-coherent backscattering spectroscopy for tissue characterization. Applied Optics. 2005 Jan 20;44:366–377. [PubMed]

8. Roy HK, Turzhitsky V, Kim Y, Goldberg MJ, Watson P, Rogers JD, Gomes AJ, Kromine A, Brand RE, Jameel M, Bogovejic A, Pradhan P, Backman V. Association between Rectal Optical Signatures and Colonic Neoplasia: Potential Applications for Screening. Cancer Research. 2009 May 15;69:4476–4483. [PMC free article] [PubMed]

9. Liu Y, Brand RE, Turzhitsky V, Kim YL, Roy HK, Hasabou N, Sturgis C, Shah D, Hall C, Backman V. Optical markers in duodenal mucosa predict the presence of pancreatic cancer. Clinical Cancer Research. 2007 Aug 1;13:4392–4399. [PubMed]

10. Turzhitsky V, Liu Y, Hasabou N, Goldberg M, Roy HK, Backman V, Brand R. Investigating population risk factors of pancreatic cancer by evaluation of optical markers in the duodenal mucosa. Disease Markers. 2008;25:313–321. [PMC free article] [PubMed]

11. Subramanian H, Pradhan P, Kim YL, Liu Y, Li X, Backman V. Modeling low-coherence enhanced backscattering using Monte Carlo simulation. Applied Optics. 2006 Aug 20;45:6292–6300. [PubMed]

12. Xu M. Low-coherence enhanced backscattering beyond diffusion. Optics Letters. 2008 Jun 1;33:1246–1248. [PubMed]

13. Born M, Wolf E, Bhatia AB. Principles of optics : electromagnetic theory of propagation, interference and diffraction of light. 7th (expanded) ed. Cambridge University Press; Cambridge [England] ; New York: 1999.

14. Ramella-Roman JC, Prahl SA, Jacques SL. Three Monte Carlo programs of polarized light transport into scattering media: part I. Optics Express. 2005 Jun 13;13:4420–4438. [PubMed]

15. Ramella-Roman JC, Prahl SA, Jacques SL. Three Monte Carlo programs of polarized light transport into scattering media: part II. Optics Express. 2005 Dec 12;13:10392–10405. [PubMed]

16. Kim YL, Turzhitsky VM, Liu Y, Roy HK, Wali RK, Subramanian H, Pradhan P, Backman V. Low-coherence enhanced backscattering: review of principles and applications for colon cancer screening. Journal of Biomedical Optics. 2006 Jul-Aug;11 [PubMed]

17. Ospeck M, Fraden S. Influence of Reflecting Boundaries and Finite Interfacial Thickness on the Coherent Backscattering Cone. Physical Review E. 1994 May;49:4578–4589. [PubMed]