Characterization of the angular memory effect of scattered light in biological tissues
High resolution optical microscopy is essential in neuroscience but suffers from scattering in biological tissues. It therefore grants access to superficial layers only. Recently developed techniques use scattered photons for imaging by exploiting angular correlations in transmitted light and could potentially increase imaging depths. But those correlations (‘angular memory effect’) are of very short range and, in theory, only present behind and not inside scattering media. From measurements on neural tissues and complementary simulations, we find that strong forward scattering in biological tissues can enhance the memory effect range (and thus the possible field-of-view) by more than an order of magnitude compared to isotropic scattering for 1 mm thick tissue layers.
Institut Langevin, ESPCI ParisTech, CNRS UMR 7587, ESPCI, 1 rue Jussieu, 75005 Paris, France.
Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK.
Physics and Astronomy Department, University of Exeter, Stocker Road, Exeter EX4 4QL, UK.
Ecole Normale Supérieure, Institut de Biologie de l’ENS, IBENS, 46 rue d’Ulm, Paris, F-75005 France.
Inserm, U1024, Paris, F-75005 France.
CNRS, UMR 8197, Paris, F-75005 France.
Laboratoire Kastler Brossel, Université Pierre et Marie Curie, Ecole Normale Supérieure, CNRS, Collège de France, 24 rue Lhomond, 75005 Paris, France.
Advances in deep tissue optical microscopy have been a driving force in biology and neuroscience for the last twenty years. Two-photon excited fluorescence microscopy has emerged as a powerful tool to study synaptic, cellular and integrative processes. However, scattering in neural tissues has limited optical imaging to superficial cortex layers only – even after surgical removal of the skull . Deeper cortical layers (layers 4 to 6) are of substantial biological interest due to their interactions with other brain regions such as the thalamus . Together with subcortical structures, e.g, the hippocampus, they still remain inaccessible without removing the overlying tissue. Thus, even small improvements in imaging depth and resolution are bound to reveal novel insights on the functioning of mammalian brains.
Efforts to increase the penetration depth of optical techniques have aimed to improve the quantity and focusing of ballistic photons by using regenerative amplifiers , longer wavelength lasers  and adaptive optics  as well as enhancing the collection of fluorescence photons . Those efforts commonly sacrifice intensity in favor of resolution as the number of ballistic photons decreases exponentially with imaging depth. In an alternative approach, multiply scattered light itself can be used for imaging [7, 8]. Recently, the ‘optical memory effect’  has emerged as a powerful tool to either enable the reconstruction of images from scattered light or to accelerate the imaging process[10, 11, 12, 13, 14]. In essence, the memory effect (ME) describes an intrinsic isoplanatism of the scattering process for small angles even in strongly scattering media. By tilting the wavefront incident on a scattering medium, the emerging speckle patterns shifts accordingly and does not immediately decorrelate.
This isoplanatism can be exploited for adaptive optics imaging techniques. Spatial modulation of the incident wavefront permits focusing behind or inside turbid media if the transmission matrix is accessible (e.g., by placing a point source or detector behind the scattering layer) [12, 15]. Since this comes at the cost of determining the transmission matrix for each image point, the ME plays a key role: if a focus is achieved, it can then be scanned over the ME range by simply introducing a linear phase tilt in the modulated wavefront .
In a more recent approach, a light-emitting object hidden behind an opaque layer is numerically reconstructed from its transmitted speckle pattern . Due to the ME, the autocorrelation of the speckle pattern is essentially equal to the object’s autocorrelation and numerical inversion yields the object itself. This method is both non-invasive and does not require sophisticated imaging equipment but its field of view is intrinsically limited by the ME range.
Albeit important for imaging, the ME is in theory very limited and scales inversely with the thickness of the medium. Furthermore, in theory, it can not be observed inside a scattering medium where the concept of isoplanatism breaks down, but only in transmission at a distance from the output plane. This can be understood when considering a tilt of the incident wavefront, i.e., the introduction of a linear phase gradient. The ME effect states that such a gradient, if sufficiently small, is preserved during the scattering process and results in the same linear phase gradient imposed on the emerging distorted wavefront. This gradient at the output plane then becomes a shift of the speckle pattern only after propagation through space. We thus expect to see the same decorrelation with the tilt of incident angle of the intensity pattern at the output plane itself but without the shift that is useful for imaging.
Recent work in biomedical imaging with adaptive optics [18, 19] suggests otherwise: phase corrections compensating for random scattering remain valid for a field-of-view of several microns inside brain tissues. This paper investigates how strong forward scattering inside biological media provides a much larger ME range that predicted by multiply scattering theory, in particular at intermediary depths of 1 mm. Complementing angular correlation measurements in neural tissues, we simulate the impact of anisotropic scattering on wave propagation through multiple forward scattering layers and qualitatively reproduce the experimental results.
2 Memory effect in anisotropic media
In the limit of weak disorder where is the wavelength and the mean free path (MFP), the ME is well described by a first order approximation to ensemble averaged
when following a tilt in both the incident wavefront and the emerging speckle pattern ( denotes the wave number). As a figure of merit, we chose the angular difference at which the correlation drops to 1/5, the value at which visual resemblance to the reference image vanishes. This ME range scales linearly with the wavelength and is inversely proportional to the sample thickness:
Note that Eq. (1) does not depend on the specific realization of the material such as the location of scatters or the transport mean free path so long as photons are scattered multiple times before reaching the output plane . For a sample thickness of mm and a wavelength of 532 nm, the speckle correlation already drops to 1/5 at an angular shift of 11.5 mdeg.
|Slice, old rat||800||47||(470)|||
The situation is expected to improve somewhat for biological tissues where the above conditions are no longer valid since waves are scattered preferentially in the forward direction in biological samples. This anisotropy is characterized by the average scattering angle
Here, denotes the probability of a photon being scattered into the angle relative to its incident direction with and the solid angle . It is easy to see that isotropic scattering yields while complete forward scattering, i.e., no scattering at all, gives . In the presence of such anisotropy, the relevant transport parameter is no longer but the transport mean free path (TMFP) with . While the is the average distance between scattering events, can be interpreted as the distance at which the direction of wave propagation has become independent from the initial direction. The condition for Eq. (1) now reads and is easily broken for large -factors and intermediate sample thicknesses. This applies in particular to biological tissues that typically show large g-factors between 0.8 and 0.98 . Scattering parameters of the rodent cortex from literature are given in Table 1 for in-vivo and ex-vivo measurements. In that regime, the ME range is presumed to be larger than predicted by Eq. (1) and possibly present inside scattering samples as well where some directionality is preserved for distances inferior to . However, the ME in this forward scattering regime is not well understood and we are not aware of any experimental studies to date linking speckle correlations with scattering anisotropy.
Here, we present a systematic study of the ME in chicken breast and rat cortex samples. We identify the anisotropy of scattering as a key factor for the extended correlation range by matching experimental results with a simulation of wave propagation through multiple forward scattering layers. In order to get a qualitative understanding of this effect, we introduce an effective thickness that corresponds to the equivalent thickness of an diffusive medium producing the same memory effect. In essence, a larger memory effect corresponds to a thinner effective sample thickness.
A sketch of the experimental setup for ME measurements is shown in Fig. 1. To measure the decorrelation of speckle patterns with changing incident angle, we focus a continuous wave laser beam (wavelength, 532 nm) on the rotation axis of a mirror mounted on a motorized rotation stage and image that spot onto the scattering sample. The lensesâ focal lengths were chosen to provide 100/9 fold magnification and a corresponding increase in angular resolution. A CCD camera placed in the far field behind the sample then captures the transmitted light. By taking a series of images while rotating the mirror at constant speed, we capture speckle patterns for a series of incident angles with a fixed point of illumination. Such speckle patterns are shown for three incident angles in the case of brain slices in Fig. 2a-c. The correlation between two images is then calculated as a function of possible offsets between patterns, giving the actual correlation to its reference from the maximum value (Fig. 2d).
Measurements on a ground-glass diffuser and zinc oxide (ZnO) scattering sample were performed to confirm alignment and accuracy of the experimental setup. Albeit being completely opaque, ground glass diffusers effectively consist of a single scattering layer and thus have a very large ME range. In contrast, ZnO is an ideal multiply scattering sample: it does not absorb visible light but scatters it quite effectively with a refractive index contrast, . Correlation curves of a ZnO diffuser are well reproduced by multiple scattering theory (Fig. 3a) and the fit parameter is close to the actual layer thickness of . This gives us the expected small ME range of mdeg.
Additionally, we confirmed that the angular correlation curves are independent of illumination spot size and thus speckle size for both isotropic and forward scattering (Fig. S1).
As model biological systems, we chose chicken breast muscle tissue and rat cortex slices, the first because it is readily available and the second since it is an actively researched system in neuroscience. To obtain an imaging geometry similar to in vivo experiments, three month old Wistar were transcardially perfused with 4% paraformaldeide (PFA) in phosphate buffer. The cortices of extracted brains were then flattened and stored overnight in a solution of 4% PFA between two glass slides separated by 1.5 mm thick spacers. Tangential slices of different thicknesses were cut from the flattened cortices. We mounted rat cortex samples in water and performed all measurements within 1-2 days after preparation while chicken breast samples were mounted as is. Respective correlation curves are shown in Fig. 3 and both effective thickness and ME range values are given in Table 2.
The most striking observation is that the fit parameter is around one order of magnitude smaller than the actual thickness , the ME range is thus much larger than for a multiple scattering sample of the same thickness.
For chicken muscle tissue, a TMFP of (MFP, ) and an anisotropy factor of have been measured . Our sample thicknesses were chosen to range from less than one to above . But due to the inherent inhomogeneity of biological tissues, those values are likely to fluctuate both within and between samples. This becomes evident when comparing the and slices: although the thickness increases, we observe that the ME ranges slightly decreases. Keeping this variability in mind, there is still a clear trend emerging from the chicken breast measurements. The ME range is up to 35 times larger than expected for isotropic scattering but this difference decreases to a factor of 5 for the thickest sample. Eventually, it is expected to vanish entirely for samples much thicker than , i.e., in the cm range. In addition to the extended ME range, we also observe a deviation from the theoretical bell shape of angular correlation curves towards a more exponential shape. For the thickest sample however, the ideal theory shape is recovered.
A similar picture emerges from the rat cortex samples. Two series of cortex slices were prepared to account for brain layers with different properties. The first series starts at the cortex surface (layer 1) and is labeled ‘top’ in Table 2 while the second starts at the cortex middle (label ‘mid’) and the thickest sample consists of an entire flattened cortex. With a smaller TMFP of [29, 21] in rat cortices, all of our samples are thicker than . Nevertheless, we observe the same qualitative differences in ME shape and range as for chicken tissue. We recover the ideal theory shape for angular correlation curves for the two thicker samples and at a thickness of , the experimental ME range is only twice as big as the theory value.
The change in ME curve shape indicates that the additional correlations are of a different nature than in the conventional ME. In multiply scattering media, the ME emerges from correlations between the input and output positions for light: a pencil beam incident on a thin slab spreads through the medium and emerges as a diffuse spot at the output plane, or a size roughly given by the thickness, but with a total loss of the incident direction. It corresponds to a macroscopic structure in the transmission matrix, with larger transmission amplitudes close to the diagonal in position space. But in case of strong forward scattering and a sample thickness close to the TMFP, the direction of wave propagation is not completely randomized when passing through the sample. Instead, as Judkewitz et al. argue , one input mode in k-space is transmitted only to a narrow cone of output directions around the incident k, thereby giving the transmission matrix a similar structure with large near-diagonal amplitudes in k-space. They derive translational (instead of angular) correlations arising from that k-space structure, even for a negligible traditional ME, i.e., for very thick slabs. In practice for our sample, we expect those correlations to contribute to the angular ME in forward scattering tissues as well.
In the above discussion, we have neglected the effects of absorption on the ME. In fact, absorption does increase angular correlations by narrowing the diffusive cone that would emerge from a pencil beam (larger near-diagonal amplitudes in the real space transmission matrix). Van Rossum and Nieuwenhuizen  derived a modified version of Eq. (1) depending on the absorption coefficient but resulting increases in correlation are very small even for strongly absorbing materials with and retain the shape of Eq. (1). Therefore, we can safely exclude absorption effects as the origin of larger correlations.
Another question is whether correlations observed in the far field are retained close to the output plane where only a fraction of the scattered amplitudes contribute to the interference pattern. This is of special interest for imaging inside biological tissues where we ideally would like to scan a focus or exploit the memory effect inside a medium. We mimicked a ânear fieldâ situation by placing a diaphragm immediately behind the sample output plane. Resulting correlations in the far field speckles thus correspond to angular correlations at the aperture. The obtained correlation curves were indeed the same as without diaphragm (Fig. S2). However, in a strongly multiple scattering sample, we do not expect an actual shift of the output plane intensity pattern since the ME relies on well-defined output modes in k-space. Instead, we would expect a mere decorrelation that only becomes a shift in the far field. While this control experiment is not a definite confirmation that ME is present inside tissues, we believe the combination of lower effective thicknesses together with the presence of translational invariance as measured in  is in agreement with the experimental observation of the ME inside tissues exploited in [18, 19] for imaging.
4 Phase mask simulation
To confirm the scattering anisotropy as cause for a longer ME range, we aim to reproduce the impact of larger g-values on correlation curves in a multiple scattering simulation. Again, we are especially interested in optically thin samples with predominant forward scattering and therefore model the scattering process with a number of consecutive phase masks (Fig. 4a). Each of those phase masks is created from a random matrix by multiplying its Fourier spectrum with a 2D Gaussian function and thus determining the frequencies of its spatial fluctuations. By changing the Gaussian’s width, we change the scattering angles and thus the -factor via (Fig. 4b)
At a phase mask , the scalar electric field accumulates a phase of and we use the Rayleigh-Sommerfeld equation , an extract solution to the scalar Helmholtz equation, for light propagation. There are two obvious shortcomings of this simplified model: First, phase masks allow for forward scattering only and all back scattering is thus assumed to be negligible and second, scattering events take place only at certain planes unlike the reality of randomly distributed scatters with anisotropies on a large range of length scales . We trade accuracy for shorter computation times but nevertheless are able to control MFP (distance between phase masks), -factor and TMFP of the scattering process.
Keeping the above in mind, we still see good qualitative agreement in both ME range and curve shape. We simulated correlation curves for anisotropy factors from 0.70-0.98 while keeping all remaining parameters fixed. As expected, we see the strongest impact of g on correlation curves in case of only a few scattering events and strong forward scattering (Fig. 4 c,d). We qualitatively reproduce the curve shape for high -values and retrieve the bell-like shape from multiple theory at lower -values. In a second series of simulations, we kept the MFP (distance between phase masks) fixed at while varying the total sample thickness from to . The resulting ME ranges are shown in Fig. 3d together with experimental values. As expected, the ratio between simulated forward scattering and theoretical multiple scattering ME range again increases with thinner samples. However, the simulation somewhat underestimates the impact of anisotropy and is uniformly bigger than the simulated values. This comes as a surprise since we expect the exclusion of backward scattering to produce stronger correlations. The discrepancy might originate from the discrete scattering planes in our model. Usually, distances between scattering events are distributed around the MFP and some photons undergo less scattering events than others (âsnake photonsâ). The phase mask simulation assumes the same number of scattering events for all photons and might thus reduce speckle correlations.
Even though angular correlations are underestimated, we can still qualitatively reproduce the cusped curve shape for weak forward scattering in thin slices together with the transition to the theory shape of Eq. (1) for stronger scattering in thicker samples.
Several emerging imaging techniques for scattering media rely on angular correlations of the transmitted light [10, 11, 17]. We have shown here that anisotropic scattering can extend the range of this optical memory effect by more than an order of magnitude when propagation through tissues.
Both the experiments and our phase mask simulation show that the preservation of directionality during scattering results in correlations that become visible in the distinct exponential like shape of angular correlation curves. From our measurements on chicken and rat tissues, we expect that diffraction-limited imaging with a field-of-view of 50-100 mdeg through a tissue layer of 1 mm should be possible at wavelength of 533 nm. The ME range scales linearly with and the TMFP is known to increase by factor of 5-10 in the near infrared and for in-vivo compared to in-vitro slices (see Table 1). At common wavelengths used in two-photon fluorescence microscopy, the ME range should increase accordingly.
Further improvements might be achieved through several methods. First, translational correlations that arise from anisotropic scattering  could well be exploited together with temporal and spectral correlation alongside the traditional angular ME to further increase the field-of-view. Second, we have seen that the ME is strongest for photons that have undergone only a limited number of scattering events. Combining traditional gating techniques that only retain snake photons  with imaging processes that use scattered light should yield further improvements for both field-of-view and imaging depth.
With an emerging better understanding of speckle correlations in anisotropic media, we expect imaging techniques that utilize scattered photons to become far more powerful for applications in biological imaging or photostimulation.
This work was funded by European Research Council Grant 278025 and the Agence Nationale de la Recherche (Investissements d’Avenir ANR-10-LABX-54 MEMO LIFE, ANR-11-IDEX-0001-02 PSL* Research University). We thank Prof. Georg Maret for enabling Sam Schott’s stay at institut Langevin and his support of the project and David Martina for technical help in the development of the experimental setup.
- Ensemble average denotes the average over possible realizations of disorder (here: locations of scatterers).
- V. Ntziachristos, “Going deeper than microscopy: the optical imaging frontier in biology.” Nat. Methods 7, 603–14 (2010).
- D. a. Dombeck, C. D. Harvey, L. Tian, L. L. Looger, and D. W. Tank, “Functional imaging of hippocampal place cells at cellular resolution during virtual navigation.” Nat. Neurosci. 13, 1433–1440 (2010).
- E. Beaurepaire, M. Oheim, and J. Mertz, “Ultra-deep two-photon fluorescence excitation in turbid media,” Opt. Commun. 188, 25–29 (2001).
- N. G. Horton, K. Wang, D. Kobat, C. G. Clark, F. W. Wise, C. B. Schaffer, and C. Xu, “In vivo three-photon microscopy of subcortical structures within an intact mouse brain.” Nat. Photonics 7, 205–209 (2013).
- N. Ji, T. R. Sato, and E. Betzig, “Characterization and adaptive optical correction of aberrations during in vivo imaging in the mouse cortex,” Proc. Natl. Acad. Sci. 109, 22–27 (2012).
- M. Oheim, E. Beaurepaire, E. Chaigneau, J. Mertz, and S. Charpak, “Two-photon microscopy in brain tissue: parameters influencing the imaging depth.” J. Neurosci. Methods 111, 29–37 (2001).
- I. Vellekoop, I. Vellekoop, A. Mosk, and A. Mosk, “Universal optimal transmission of light through disordered materials,” Phys. Rev. Lett. 101, 120601 (2008).
- S. M. Popoff, G. Lerosey, R. Carminati, M. Fink, a. C. Boccara, and S. Gigan, “Measuring the transmission matrix in optics: An approach to the study and control of light propagation in disordered media,” Phys. Rev. Lett. 104, 1–5 (2010).
- S. Feng, C. Kane, P. Lee, and A. Stone, “Correlations and Fluctuations of Coherent Wave Transmission through Disordered Media,” Phys. Rev. Lett. 61, 834–837 (1988).
- J. Bertolotti, E. G. van Putten, C. Blum, A. Lagendijk, W. L. Vos, and A. P. Mosk, “Non-invasive imaging through opaque scattering layers.” Nature 491, 232–4 (2012).
- O. Katz, E. Small, and Y. Silberberg, “Looking around corners and through thin turbid layers in real time with scattered incoherent light,” Nat. Photonics 6, 549–553 (2012).
- M. Nixon, O. Katz, E. Small, Y. Bromberg, A. A. Friesem, Y. Silberberg, and N. Davidson, “Real-time wavefront shaping through scattering media by all-optical feedback,” Nat. Photonics 7, 919–924 (2013).
- X. Yang, Y. Pu, and D. Psaltis, “Imaging blood cells through scattering biological tissue using speckle scanning microscopy.” Opt. Express 22, 3405–13 (2014).
- K. T. Takasaki and J. W. Fleischer, “Phase-space measurement for depth-resolved memory-effect imaging,” Opt. Express 22, 784–790 (2014).
- T. Chaigne, O. Katz, A. C. Boccara, M. Fink, E. Bossy, and S. Gigan, “Controlling light in scattering media non-invasively using the photoacoustic transmission matrix,” Nat. Photonics 8, 58–64 (2013).
- I. M. Vellekoop and C. M. Aegerter, “Scattered light fluorescence microscopy: imaging through turbid layers.” Opt. Lett. 35, 1245–1247 (2010).
- O. Katz, P. Heidmann, M. Fink, and S. Gigan, “Non-invasive real-time imaging through scattering layers and around corners via speckle correlations,” Nat. Photonics 8, 784–790 (2014).
- J. Tang, R. N. Germain, and M. Cui, “Superpenetration optical microscopy by iterative multiphoton adaptive compensation technique,” Proc. Natl. Acad. Sci. 109, 8434–8439 (2012).
- L. Kong and M. Cui, “In vivo fluorescence microscopy via iterative multi-photon adaptive compensation technique,” Opt. Express 22, 2649–2655 (2014).
- R. Berkovits, M. Kaveh, and S. Feng, “Memory effect of waves in disordered systems: A real-space approach,” Phys. Rev. B 40, 737–740 (1989).
- M. Mesradi, A. Genoux, V. Cuplov, D. Abi Haidar, S. Jan, I. Buvat, and F. Pain, “Experimental and analytical comparative study of optical coefficient of fresh and frozen rat tissues.” J. Biomed. Opt. 18, 117010 (2013).
- D. Kobat, M. E. Durst, N. Nishimura, A. W. Wong, C. B. Schaffer, and C. Xu, “Deep tissue multiphoton microscopy using longer wavelength excitation.” Opt. Express 17, 13354–13364 (2009).
- D. Kleinfeld, P. P. Mitra, F. Helmchen, and W. Denk, “Fluctuations and stimulus-induced changes in blood flow observed in individual capillaries in layers 2 through 4 of rat neocortex.” Proc. Natl. Acad. Sci. U. S. A. 95, 15741–15746 (1998).
- D. Kobat, N. G. Horton, and C. Xu, ‘‘In vivo two-photon microscopy to 1.6-mm depth in mouse cortex,” J. Biomed. Opt. 16, 106014 (2011).
- a. N. Yaroslavsky, P. C. Schulze, I. V. Yaroslavsky, R. Schober, F. Ulrich, and H. J. Schwarzmaier, “Optical properties of selected native and coagulated human brain tissues in vitro in the visible and near infrared spectral range.” Phys. Med. Biol. 47, 2059–2073 (2002).
- E. Akkermans and G. Montambaux, Mesoscopic Physics of Electrons and Photons (Cambridge University Press, 2007).
- W. F. W. F. Cheong, S. A. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron. 26, 2166–2185 (1990).
- T. J. Coutts, X. Li, T. M. Barnes, B. M. Keyes, C. L. Perkins, S. E. Asher, S. Zhang, S.-H. Wei, and S. Limpijumnong, Zinc Oxide Bulk, Thin Films and Nanostructures (Elsevier Science Ltd, 2006).
- I. Nishidate, S. Kawauchi, S. Sato, M. Ishihara, M. Kikuchi, and M. Sato, ‘‘In vivo determination of absorption and scattering properties in rat cerebral cortex using single reflectance fiber probe with two source- collector geometries Izumi,” Proc. SPIE 8087, 80872P (2011).
- B. Judkewitz, R. Horstmeyer, I. M. Vellekoop, and C. Yang, “Translation correlations in anisotropically scattering media,” arXiv:1411.7157 (2014).
- M. van Rossum and T. Nieuwenhuizen, “Influence of skin layers on speckle correlations of light transmitted through disordered media,” Phys. Lett. A 177, 452–458 (1993).
- A. S. Marathay and J. E. McCalmont, “On the usual approximation used in the Rayleigh-Sommerfeld diffraction theory.” J. Opt. Soc. Am. A. Opt. Image Sci. Vis. 21, 510–516 (2004).
- J. R. Mourant, M. Canpolat, C. Brocker, O. Esponda-Ramos, T. M. Johnson, A. Matanock, K. Stetter, and J. P. Freyer, “Light scattering from cells: the contribution of the nucleus and the effects of proliferative status.” J. Biomed. Opt. 5, 131–137 (2000).
- C. Dunsby and P. M. W. French, “Techniques for depth-resolved imaging through turbid media including coherence-gated imaging,” J. Phys. D. Appl. Phys. 36, R207–R227 (2003).