OSA's Digital Library

Virtual Journal for Biomedical Optics

Virtual Journal for Biomedical Optics

| EXPLORING THE INTERFACE OF LIGHT AND BIOMEDICINE

  • Editor: Gregory W. Faris
  • Vol. 2, Iss. 9 — Sep. 26, 2007
« Show journal navigation

Focusing coherent light through opaque strongly scattering media

I. M. Vellekoop and A. P. Mosk  »View Author Affiliations


Optics Letters, Vol. 32, Issue 16, pp. 2309-2311 (2007)
http://dx.doi.org/10.1364/OL.32.002309


View Full Text Article

Acrobat PDF (391 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We report focusing of coherent light through opaque scattering materials by control of the incident wavefront. The multiply scattered light forms a focus with a brightness that is up to a factor of 1000 higher than the brightness of the normal diffuse transmission.

© 2007 Optical Society of America

Random scattering of light is what makes materials such as white paint, milk, or human tissue opaque. In these materials, repeated scattering and interference distort the incident wavefront so strongly that all spatial coherence is lost [1

1. P. Sebbah, ed., Waves and Imaging through Complex Media (Kluwer, 2001).

]. Incident coherent light diffuses through the medium and forms a volume speckle field that has no correlations on a distance larger than the wavelength of light. The complete scrambling of the field makes it impossible to control light propagation using the well-established wavefront correction methods of adaptive optics (see e.g., [2

2. R. K. Tyson, Principles of Adaptive Optics (Academic, 1998).

]).

We demonstrate focusing of coherent light through disordered scattering media by the construction of wavefronts that invert diffusion of light. Our method relies on interference and is universally applicable to scattering objects regardless of their constitution and scattering strength. We envision that, with such active control, random scattering will become beneficial, rather than detrimental, to imaging [1

1. P. Sebbah, ed., Waves and Imaging through Complex Media (Kluwer, 2001).

] and communication [3

3. A. Derode, A. Tourin, J. de Rosny, M. Tanter, S. Yon, and M. Fink, Phys. Rev. Lett. 90, 014301 (2003). [CrossRef] [PubMed]

, 4

4. S. H. Simon, A. L. Moustakas, M. Stoytchev, and H. Safar, Phys. Today 54(9), 38 (2001). [CrossRef]

, 5

5. G. Lerosey, J. de Rosny, A. Tourin, and M. Fink, Science 315, 1120 (2007). [CrossRef] [PubMed]

].

Figure 1 shows the principle of the experiment. Normally, incident light is scattered by the sample and forms a random speckle pattern [Fig. 1a]. The goal is to match the incident wavefront to the sample so that the scattered light is focused in a specified target area [Fig. 1b]. The experimental setup for constructing such wavefronts is shown in Fig. 2 . Light from a 632.8nm HeNe laser is spatially modulated by a liquid-crystal phase modulator and focused on an opaque, strongly scattering sample. The number of degrees of freedom of the modulator is reduced by grouping pixels into a variable number (N) of square segments. A CCD camera monitors the intensity in the target focus and provides feedback for an algorithm that programs the phase modulator.

We performed first tests of inverse wave diffusion using rutile (TiO2) pigment, which is one of the most strongly scattering materials known. The sample consists of an opaque, 10.1-μm-thick layer of rutile [6

6. R. H. J. Kop, P. de Vries, R. Sprik, and A. Lagendijk, Phys. Rev. Lett. 79, 4369 (1997). [CrossRef]

] with a transport mean free path of 0.55±0.10μm measured at λ=632.8nm. Since in this sample the transmitted light is scattered hundreds of times, there is no direct relation between the incident wavefront and the transmitted image [7

7. R. Pappu, B. Recht, J. Taylor, and N. Gershenfeld, Science 297, 2026 (2002). [CrossRef] [PubMed]

, 8

8. J. W. Goodman, Statistical Optics (Wiley, 2000).

].

In Fig. 3 we show the intensity pattern of the transmitted light. In Fig. 3a we see the pattern that was transmitted when a plane wave was focused onto the sample. The light formed a typical random speckle pattern with a low intensity. We then optimized the wavefront so that the transmitted light focused to a target area with the size of a single speckle. The result for a wavefront composed of 3228 individually controlled segments is seen in Fig. 3b, where a single bright spot stands out clearly against the diffuse background. The focus was over a factor of 1000 more intense than the nonoptimized speckle pattern. By adjusting the target function used as feedback it is also possible to optimize multiple foci simultaneously, as is shown in Fig. 3c where a pattern of five spots was optimized. Each of the spots has an intensity of approximately 200 times the original diffuse intensity. In Fig. 3d we show the phase of the incident wavefront corresponding to Fig. 3c. Neighboring segments are uncorrelated, which indicates that the sample fully scrambles the incident wavefront.

The algorithm that constructs the inverse diffusion wavefront uses the linearity of the scattering process. The transmitted field in the target, Em, is a linear combination of the fields coming from the N different segments of the modulator:
Em=n=1NtmnAneiϕn,
(1)
where An and ϕn are, respectively, the amplitude and phase of the light reflected from segment n. Scattering in the sample and propagation through the optical system is described by the elements tmn of the unknown transmission matrix. Clearly, the magnitude of Em will be the highest when all terms in Eq. (1) are in phase. We determine the optimal phase for a single segment at a time by cycling its phase from 0 to 2π. For each segment we store the phase at which the target intensity is the highest. At that point the contribution of the segment is in phase with the already present diffuse background. After the measurements have been performed for all segments, the phase of the segments is set to their stored values. Now the contributions from all segments interfere constructively and the target intensity is at the global maximum. A preoptimization with a small number of segments significantly improves the signal-to-noise ratio. This method is generally applicable to linear systems and does not rely on time reversal symmetry or absence of absorption. Although mathematically this algorithm is the most efficient, in noisy experimental conditions adaptive learning algorithms [9

9. D. E. Goldberg, Genetic Algorithms in Search, Optimization and Machine Learning (Addison-Wesley, 1989).

] might be more effective, and an investigation of such algorithms is on its way.

We tested the universal scaling behavior implied by Eq. (2) by changing N. Using the same TiO2 sample as before, the algorithm was targeted to construct a collimated beam. In Fig. 4 the enhancement is plotted as a function of the number of segments for different focusing conditions. The linear relation between the enhancement and the number of segments is evident until the enhancement saturates at η=1000. All measured enhancements were slightly below the theoretical maximum. This is understandable since all perturbations move the system away from the global maximum. The main reason for deviations from the optimal wavefront is residual amplitude modulation in the phase modulator, which introduced an uncontrolled bias in the field amounting to 14% of the total intensity.

The saturation of the enhancement is the result of slow changes in the speckle pattern. This instability effectively limited the number of segments for which the optimal phase could be measured. We estimate that the effective enhancement decreases to ηeff=η(1+NTTp), where T=1.2s is the time needed for one measurement and the persistence time Tp=5400s is the time scale at which the speckle pattern of the TiO2 sample remains stable. Depending on the environmental conditions, Tp can be consider- ably higher and enhancements of over 2000 have been measured overnight.

To verify the universal applicability of inversion of wave diffusion, we used a variety of materials of natural origin. Table 1 lists the intensity enhancement for different materials we used. Although the samples vary in thickness, composition, and scattering strength, they were all able to focus a properly prepared wavefront to a sharp spot. The intensity enhancement varies between 60 and 1000. The main reason for this variation is that the persistence time is not the same for all materials.

In summary, our results show that precise control of diffuse light is possible using an optimal, noniterative algorithm; light can be directed through opaque objects to form one or multiple foci. The brightness of the focal spot is explained by a model based on statistical optics. We expect inverse wave diffusion to have applications in imaging and light delivery in scattering media, possibly including metal nanostructures [13

13. M. I. Stockman, S. V. Faleev, and D. J. Bergman, Phys. Rev. Lett. 88, 067402 (2002). [CrossRef] [PubMed]

]. Dynamic measurements in biological tissue are possible when the time required for achieving a focus can be reduced to below 1ms per segment [14

14. B. J. Vakoc, S. H. Yun, J. F. de Boer, G. J. Tearney, and B. E. Bouma, Opt. Express 13, 5483 (2005). [CrossRef] [PubMed]

, 15

15. J. Li, G. Dietsche, D. Iftime, S. E. Skipetrov, G. Maret, T. Elbert, B. Rockstroh, and T. Gisler, J. Biomed. Opt. 10, 044002 (2005). [CrossRef]

]; we estimate that this time scale is technologically possible with the use of fast phase modulators [16

16. M. Hacker, G. Stobrawa, R. Sauerbrey, T. Buckup, M. Motzkus, M. Wildenhain, and A. Gehner, Appl. Phys. B: Lasers Opt. 76, 711 (2003).

]. Furthermore, the high degree of control over the scattered light should permit experimental verification of random matrix theories for the transport of light [11

11. C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). [CrossRef]

, 12

12. J. B. Pendry, A. MacKinnon, and P. J. Roberts, Proc. R. Soc. London, Ser. A 437, 67 (1990).

].

We thank Ad Lagendijk for valuable discussions, Willem Vos and Vinod Subramaniam for a critical reading of the manuscript, and the Photon Scattering group of the Institute for Atomic and Molecular Physics (AMOLF) for providing samples. This work is part of the research program of the Stichting voor Fundamenteel Onderzoek der Materie (FOM), which is financially supported by the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO).

Table 1. Measured Intensity Enhancement for Different Materialsa

table-icon
View This Table
Fig. 1 Design of the experiment. (a) A plane wave is focused on a disordered medium, and a speckle pattern is transmitted. (b) The wavefront of the incident light is shaped so that scattering makes the light focus at a predefined target.
Fig. 2 Schematic of the apparatus. A 632.8nm HeNe laser beam is expanded and reflected off a Holoeye LR-2500 liquid crystal spatial light modulator (SLM). Polarization optics select a phase mostly modulation mode. The SLM is imaged onto the entrance pupil of the objective with a 1:3 demagnifying lens system (not shown). The objective is overfilled; we use only segments that fall inside the pupil. The shaped wavefront is focused on the strongly scattering sample (S), and a CCD camera images the transmitted intensity pattern. λ4, quarter-wave plate; λ2, half-wave plate; M, mirror; BS, 50% nonpolarizing beam splitter; P, polarizer.
Fig. 3 Transmission through a strongly scattering sample consisting of TiO2 pigment. (a) Transmission micrograph with an unshaped incident beam. (b) Transmission after optimization for focusing at a single target. The scattered light is focused to a spot that is 1000 times brighter than the original speckle pattern. (c) Multibeam optimization. The disordered medium generates five sharp foci at the defined positions. (a)–(c) are presented on the same logarithmic color scale that is normalized to the average transmission before optimization. (d) Phase of the incident wavefront used to form (c).
Fig. 4 Measured intensity enhancement as a function of the number of segments. Squares, sample in focus; triangles, sample 100μm behind focus; solid curve, ideal enhancement [Eq. (2)], dotted curve, corrected for residual amplitude modulation and finite persistence time of Tp=5400s. The experimental uncertainty is of the order of the symbol size.
1.

P. Sebbah, ed., Waves and Imaging through Complex Media (Kluwer, 2001).

2.

R. K. Tyson, Principles of Adaptive Optics (Academic, 1998).

3.

A. Derode, A. Tourin, J. de Rosny, M. Tanter, S. Yon, and M. Fink, Phys. Rev. Lett. 90, 014301 (2003). [CrossRef] [PubMed]

4.

S. H. Simon, A. L. Moustakas, M. Stoytchev, and H. Safar, Phys. Today 54(9), 38 (2001). [CrossRef]

5.

G. Lerosey, J. de Rosny, A. Tourin, and M. Fink, Science 315, 1120 (2007). [CrossRef] [PubMed]

6.

R. H. J. Kop, P. de Vries, R. Sprik, and A. Lagendijk, Phys. Rev. Lett. 79, 4369 (1997). [CrossRef]

7.

R. Pappu, B. Recht, J. Taylor, and N. Gershenfeld, Science 297, 2026 (2002). [CrossRef] [PubMed]

8.

J. W. Goodman, Statistical Optics (Wiley, 2000).

9.

D. E. Goldberg, Genetic Algorithms in Search, Optimization and Machine Learning (Addison-Wesley, 1989).

10.

N. Garcia and A. Z. Genack, Phys. Rev. Lett. 63, 1678 (1989). [CrossRef] [PubMed]

11.

C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). [CrossRef]

12.

J. B. Pendry, A. MacKinnon, and P. J. Roberts, Proc. R. Soc. London, Ser. A 437, 67 (1990).

13.

M. I. Stockman, S. V. Faleev, and D. J. Bergman, Phys. Rev. Lett. 88, 067402 (2002). [CrossRef] [PubMed]

14.

B. J. Vakoc, S. H. Yun, J. F. de Boer, G. J. Tearney, and B. E. Bouma, Opt. Express 13, 5483 (2005). [CrossRef] [PubMed]

15.

J. Li, G. Dietsche, D. Iftime, S. E. Skipetrov, G. Maret, T. Elbert, B. Rockstroh, and T. Gisler, J. Biomed. Opt. 10, 044002 (2005). [CrossRef]

16.

M. Hacker, G. Stobrawa, R. Sauerbrey, T. Buckup, M. Motzkus, M. Wildenhain, and A. Gehner, Appl. Phys. B: Lasers Opt. 76, 711 (2003).

OCIS Codes
(030.6600) Coherence and statistical optics : Statistical optics
(290.1990) Scattering : Diffusion
(290.4210) Scattering : Multiple scattering

ToC Category:
Scattering

History
Original Manuscript: March 6, 2007
Revised Manuscript: June 14, 2007
Manuscript Accepted: June 21, 2007
Published: August 2, 2007

Virtual Issues
Vol. 2, Iss. 9 Virtual Journal for Biomedical Optics

Citation
I. M. Vellekoop and A. P. Mosk, "Focusing coherent light through opaque strongly scattering media," Opt. Lett. 32, 2309-2311 (2007)
http://www.opticsinfobase.org/vjbo/abstract.cfm?URI=ol-32-16-2309


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. P.Sebbah, ed., Waves and Imaging through Complex Media (Kluwer, 2001).
  2. R. K. Tyson, Principles of Adaptive Optics (Academic, 1998).
  3. A. Derode, A. Tourin, J. de Rosny, M. Tanter, S. Yon, and M. Fink, Phys. Rev. Lett. 90, 014301 (2003). [CrossRef] [PubMed]
  4. S. H. Simon, A. L. Moustakas, M. Stoytchev, and H. Safar, Phys. Today 54(9), 38 (2001). [CrossRef]
  5. G. Lerosey, J. de Rosny, A. Tourin, and M. Fink, Science 315, 1120 (2007). [CrossRef] [PubMed]
  6. R. H. J. Kop, P. de Vries, R. Sprik, and A. Lagendijk, Phys. Rev. Lett. 79, 4369 (1997). [CrossRef]
  7. R. Pappu, B. Recht, J. Taylor, and N. Gershenfeld, Science 297, 2026 (2002). [CrossRef] [PubMed]
  8. J. W. Goodman, Statistical Optics (Wiley, 2000).
  9. D. E. Goldberg, Genetic Algorithms in Search, Optimization and Machine Learning (Addison-Wesley, 1989).
  10. N. Garcia and A. Z. Genack, Phys. Rev. Lett. 63, 1678 (1989). [CrossRef] [PubMed]
  11. C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). [CrossRef]
  12. J. B. Pendry, A. MacKinnon, and P. J. Roberts, Proc. R. Soc. London, Ser. A 437, 67 (1990).
  13. M. I. Stockman, S. V. Faleev, and D. J. Bergman, Phys. Rev. Lett. 88, 067402 (2002). [CrossRef] [PubMed]
  14. B. J. Vakoc, S. H. Yun, J. F. de Boer, G. J. Tearney, and B. E. Bouma, Opt. Express 13, 5483 (2005). [CrossRef] [PubMed]
  15. J. Li, G. Dietsche, D. Iftime, S. E. Skipetrov, G. Maret, T. Elbert, B. Rockstroh, and T. Gisler, J. Biomed. Opt. 10, 044002 (2005). [CrossRef]
  16. M. Hacker, G. Stobrawa, R. Sauerbrey, T. Buckup, M. Motzkus, M. Wildenhain, and A. Gehner, Appl. Phys. B: Lasers Opt. 76, 711 (2003).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1 Fig. 2 Fig. 3
 
Fig. 4
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited