# Deterministic control of broadband light through a multiply scattering medium via the multispectral transmission matrix

Daria Andreoli, Giorgio Volpe, Sébastien Popoff, Ori Katz, Samuel Grésillon, Sylvain Gigan

► To cite this version:

Daria Andreoli, Giorgio Volpe, Sébastien Popoff, Ori Katz, Samuel Grésillon, et al., Deterministic control of broadband light through a multiply scattering medium via the multispectral transmission matrix. Scientific Reports, 2015, 5, pp.10347. 10.1038/srep10347. hal-01223433

HAL Id: hal-01223433

https://hal.sorbonne-universite.fr/hal-01223433

Submitted on 2 Nov 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

OPEN

# Deterministic control of broadband light through a multiply scattering medium via the multispectral transmission matrix

Received: 17 November 2014

Accepted: 08 April 2015

Published: 12 May 2015

Daria Andreoli<fnref n="1" />,², Giorgio Volpe<fnref n="1" />,³, Sébastien Popoff², Ori Katz<fnref n="1" />,², Samuel Grésillon² & Sylvain Gigan<fnref n="1" />

We present a method to measure the spectrally-resolved transmission matrix of a multiply scattering medium, thus allowing for the deterministic spatospectral control of a broadband light source by means of wavefront shaping. As a demonstration, we show how the medium can be used to selectively focus one or many spectral components of a femtosecond pulse, and how it can be turned into a controllable dispersive optical element to spatially separate different spectral components to arbitrary positions.

Wave propagation through a multiply scattering medium is an important problem for different disciplines, from optics to acoustics and microwaves<fnref n="1" />⁻³, since it finds important applications from fundamental physics to biomedical imaging⁴, ⁵.

For a monochromatic coherent input beam, the medium typically generates a speckle pattern⁶ - the result of complex interference between the scattered waves - that is usually considered a major drawback for imaging and focusing. In optics, wavefront shaping by spatial light modulators (SLMs) has emerged as a powerful tool to control the speckle and thus light propagation through such complex media⁷. Different approaches have been proposed, in particular via optimization methods⁸, to focus light through⁹ or into¹⁰ a complex medium. However, a more general way to understand light propagation is to describe it via its transmission matrix (TM)¹¹, which in essence describes the fact that propagation through a disordered medium, however complex, remains a linear process, and therefore can be described by a linear operator H linking input to output modes, i.e. it can be described as $E^{out} = H \cdot E^{in}$. Recently, we introduced a method to measure the monochromatic transmission matrix of a complex medium and to subsequently demonstrate focusing and imaging through it¹², ¹³. Matricial methods are now widely used to characterize disordered media¹⁴, ¹⁵ or multimode fibers¹⁶⁻¹⁹, and even focusing inside complex media²⁰, ²¹.

Generally speaking, the medium has a spectrally dependent response, with a spectral correlation bandwidth, $\Delta\lambda_c$, that is inversely related to the confinement time of light in the medium $\Delta t$ ($\Delta\lambda_c = (\lambda_c^2)/(c\Delta t)^{22-25}$). The monochromatic treatment is, therefore, only valid as long as the bandwidth of the source is contained within this spectral correlation bandwidth. Consequently, when light from a broadband source passes through a multiply-scattering medium, its different spectral components generates different speckle patterns. This is the case where a short pulse - shorter than the average confinement time through the medium - is scattered in space and spread in time. Thanks to the inherent coupling between spatial and spectral degrees of freedom in a scattering medium²⁶, it is still possible to achieve spatiotemporal focusing through it by optimising the spatial input wavefront²⁷⁻²⁹.

<footnote n="1">¹Laboratoire Kastler Brossel, ENS-PSL Research University, CNRS, UPMC-Sorbonne universites, College de France; 24 rue Lhomond, F-75005 Paris, France. ²Institut Langevin, UPMC, ESPCI ParisTech, CNRS UMR 7587, 1 rue Jussieu, 75005 Paris. ³Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK. Correspondence and requests for materials should be addressed to S.G. (email: sylvain.gigan@lkb.ens.fr)</footnote>
In this work, we show that we are able to measure the multispectral transmission matrix (MSTM) of a complex medium - i.e. the tridimensional matrix that is essentially a set of monochromatic TMs - for several chosen wavelengths. First, we show that we are able to retrieve the spectral correlation bandwidth of the medium from the spectral correlations in the MSTM. Then, exploiting the MSTM, we can achieve efficient and deterministic spatiospectral control of a broadband and spatially coherent light source, here a femtosecond pulse, by spatial-only wavefront shaping. The MSTM gives more flexibility than both the optimization<sup>27, 28</sup> and the pulse shaping approach<sup>30</sup> to control separately and individually spatial and spectral components of the output fields. We demonstrate selective spatial focusing of a given wavelength as well as multispectral focusing. Finally, we show that the medium can be turned into a controllable dispersive optical element to focus different spectral components of the pulse at different spatial positions, similarly to a grating<sup>31, 32</sup>.

## Methods

In essence, in the monochromatic TM measurement approach, a set of input wavefronts, that describes a complete basis of the input modes, is sent through the medium. Part of the beam remains unmodulated and generates a stationary speckle that serves as a reference<sup>12</sup>. Using phase-shifting interferometry, the amplitude and phase of the output signal can then be retrieved from a set of speckle images. In order to measure the full MSTM, the laser is operated in CW mode at a given wavelength followed by the procedure for measuring the monochromatic TM of the medium<sup>12</sup>; this step is then reiterated at every desired wavelength by changing the laser wavelength. In this 3D matrix, the output field at the i-th pixel on the CCD camera is connected to the input field by

$$E_i^{out}(\omega_k) = \sum_{j=1}^{N} h_{i,j,k} E_j^{in}(\omega_k), \quad (1)$$

where $h_{i,j,k}$ is the element of the MSTM linking the i-th output CCD pixel to the j-th input SLM pixel at frequency $\omega_k$. Just as in<sup>12</sup>, the reference is an unknown speckle whose amplitude and phase information is encoded in the $h_{i,j,k}$ coefficients. Since this reference is frequency dependent, the information on the relative local phase between monochromatic TMs cannot be accessed in this type of measurement. This measurement, therefore, precludes exact temporal control but still allows spatiospectral control of a broadband light source.

Fig. 1(a) shows a simplified scheme of the experimental setup used to measure the MSTM. A titanium-sapphire laser (Ti:Sa, Spectra-Physics Mai-Tai) illuminates a two-dimensional phase-only spatial light modulator (SLM, Hamamatsu X10468-04). The laser is meant to work in mode-locked pulsed regime (central wavelength tunable between 710 nm and 920 nm, pulse duration ~100 fs and corresponding bandwidth $\Delta\lambda_L \sim 10$ nm) but it can also be unlocked and operated as a continuous wave (CW), tunable laser between 710 nm and 920 nm with a coherence length of >10 cm. With such a coherence length, the laser in CW can be considered monochromatic for the scattering media studied in this work. A spectrometer (HR400 Ocean Optics) is used to monitor the spectrum and the mode of operation of the laser, with a spectral resolution of 0.2 nm. The illuminated SLM pixels are grouped into macro-pixels, each of which induces a controllable phase shift $\phi_j$ between 0 and $2\pi$. The light is reflected by the SLM and passes through a strongly scattering layer of ZnO or TiO<sub>2</sub> powder of different thicknesses<sup>33</sup>. Figs. 1(b-d) show examples of measured speckle patterns transmitted through the same zone of one sample at different wavelengths but for a same input wavefront: these speckles are strongly wavelength-dependent. By correlating speckle images taken while tuning the wavelength of the laser in CW mode by steps of 0.2 nm, the spectral correlation bandwidth of the medium can be quantified<sup>23</sup>. Fig. 1(e) shows these correlation curves for samples of different thicknesses and materials. Their full-width-half-maximum defines the spectral correlation bandwidth $\Delta\lambda_c$ of the corresponding medium.

Fig. 2(a) shows that the correlation between different spectral components of the MSTM is consistent with the spectral correlation bandwidth of the medium as measured in Fig. 1(e). In order to maximise the information and minimise the amount of measurements required to control a polychromatic pulse of bandwidth $\Delta\lambda_L$, the spectral sampling of the MSTM can be chosen to match the spectral correlation bandwidth $\Delta\lambda_c$ of the medium. The elements of the MSTM are sufficiently uncorrelated and yet contain enough information to describe accurately the scattering process. Here, therefore, the medium response is sampled by a MSTM composed of $N_\lambda$ monochromatic TMs with $N_\lambda = (\Delta\lambda_L)/(\Delta\lambda_c)$ (Fig. 2(b)).

In the following, all measurements are performed on a sample with spectral correlation bandwidth of 1.6 nm, corresponding to an average confinement time of 625 fs. The input light is a broadband pulse with a bandwidth of 10 nm, corresponding to a Fourier-limited duration of 100 fs. Therefore, an experimental $N_\lambda = 6$ is chosen. The number of input macro-pixels is set to 256, corresponding to 16 × 16 macro-pixels, and, on the CCD camera, an observation region, approximately corresponding to 10 × 10 speckle grains, is selected.

## Results and discussion

**Spatial focusing.** The first most straightforward use of the MSTM is as a set of monochromatic TMs to focus light at a specific wavelength through a strongly scattering medium by phase conjugation, as
<figure><img src="image_2.png" /><figcaption>**Figure 1.** (a) Experimental setup. A tunable femtosecond Ti:Sapphire laser (Ti:Sa), working either in CW or pulsed mode, is expanded by a telescope (T) to illuminate a two-dimensional phase-only spatial light modulator (SLM) operated in reflection. A 4f system ($f_1$ = 75 mm) conjugates the SLM with the back focal plane of a low numerical aperture objective (O₁, NA = 0.3) that focuses the laser beam on a strongly scattering medium (M). The transmitted light is then imaged on a CCD camera after polarization selection by a linear polariser (P). A spectrometer (S) is used to monitor the output wavelength of the laser. (b-d) Examples of the speckle patterns at the output of a ZnO sample with spectral correlation bandwidth $\Delta\lambda_c$ = 1.6 nm for three totally uncorrelated wavelengths: (b) $\lambda$ = 797 nm, (c) $\lambda$ = 800 nm, and (d) $\lambda$ = 803 nm. The white circles highlight the decorrelation of one bright speckle grain at the central wavelength $\lambda$ = 800 nm. (e) Normalised spectral correlation functions $C(\Delta\lambda)$ of different strongly scattering media obtained by correlating different speckle images acquired while scanning the laser wavelength from $\lambda$ = 794 nm to $\lambda$ = 806 nm with a step of $\Delta\lambda$ = 0.2 nm. The full width half maximum (FWHM) of these curves gives the spectral correlation bandwidth $\Delta\lambda_c$ of the different samples: TiO₂ with $\Delta\lambda_c$ = 7.2 nm (circles), ZnO with $\Delta\lambda_c$ = 2.6 nm (squares), ZnO with $\Delta\lambda_c$ = 1.4 nm (diamonds), and ZnO with $\Delta\lambda_c$ = 0.4 nm (triangles). For the ZnO samples the correlation bandwidth was lowered by increasing the thickness of the scattering layer. The shaded areas represent one standard deviation around the average values.</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Measurement of the multispectral transmission matrix (MSTM). (a) The cyan dots show the spectral correlation C(Δλ) among different monochromatic transmission matrices (TMs) for a ZnO medium (Δλc = 1.6 nm) measured in the wavelength range between λ = 794 nm and λ = 806 nm with a step of Δλ = 0.2 nm. This correlation curve is consistent with the spectral correlation curve C(Δλ) calculated by directly correlating speckle intensity images (Fig.1(e)) and here represented by the red dashed line. The shaded areas represent one standard deviation around the average values. (b) In order to maximise the information about the scattering process and minimise the measurement time, the MSTM is a 3D matrix composed of N monochromatic matrices recorded at regular intervals within the bandwidth ΔλL of the polychromatic pulse propagating through the scattering medium.</figcaption></figure>

shown in Fig. 3. This corresponds to selecting a vertical section of the MSTM (in the CCD-SLM pixels plane) in Fig. 2(b). In this case, only the fraction of the input beam at the chosen wavelength (within the spectral correlation bandwidth) will be effectively controlled<sup>23</sup>. Figs. 3(a-i) show focusing images of monochromatic waves at various wavelengths (columns) using the phase mask given by a vertical section of the MSTM recorded at a different wavelength (rows). As expected, a focus appears only when the incident wavelength matches the wavelength at which the monochromatic TM was recorded.

**Multispectral focusing.** Just like it is possible to focus in the monochromatic case in several points simultaneously by adding coherently the corresponding phase masks on the SLM<sup>12</sup>, it is also possible to focus several wavelengths at the same position by displaying a coherent sum of different phase-conjugate wavefronts for different wavelengths extracted from the MSTM. When one considers a single output position, an elegant way to formalize the above solution is to take a horizontal section of the MSTM as shown in Fig. 4(a). This section is a two-dimensional matrix $\tilde{H}$ that links the spectral components of the input SLM pixels to the spectral components at a given output spatial position $i_t$. In a quasi-monochromatic case, when the bandwidth of the polychromatic source is much smaller than the central wavelength, the dependency in $\omega_k$ at the input can be neglected. In this configuration, Eq. (1) reduces to

$$E_{i_t}^{\text{out}}(\omega_k) = \sum_{j=1}^{N} \tilde{h}_{i_t, j, k} E_j^{\text{in}} \qquad (2)$$

This problem is equivalent to the one treated in<sup>12</sup> but now the output dimension is spectral rather than spatial. For a given desired output spectrum vector $E^{\text{target}}$, the input vector can be obtained by phase conjugation as $E^{\text{in}} = \tilde{\tilde{H}}^{\dagger} \cdot E^{\text{target}}$.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Selective monochromatic focusing through a strongly scattering medium with the MSTM. (a-i) Focusing images of a monochromatic beam at (a,d,g) $\lambda = 797$ nm, (b,e,h) $\lambda = 800$ nm, and (c,f,i) $\lambda = 803$ nm, using the phase mask given by a vertical section of the MSTM, in the CCD-SLM pixels plane, (Fig.2(b)) recorded at (a-c) $\lambda = 797$ nm, (d-f) $\lambda = 800$ nm, and (g-i) $\lambda = 803$ nm. A focus appears only when the TM recorded at the corresponding wavelength of the beam is selected in (a) for $\lambda = 797$ nm, in (e) for $\lambda = 800$ nm, and in (i) for $\lambda = 803$ nm, respectively.</figcaption></figure>

In this case, with homogeneous illumination and phase-only control, the phase of the calculated $E^{in}$ needs to be displayed on the SLM. In order to focus spatially one spectral component $\omega_k$, $E^{target}(\omega_k)$ is set to 1 while all the other components are set to zero: the resulting wavefront is identical to what is obtained by phase conjugation of the monochromatic TM at $\omega_k$. In order to focus all spectral components with the same weight, the imposed condition is $E^{target} = 1$. Using a phase masks that focuses several wavelengths simultaneously, and scanning the wavelength in CW mode, Figs. 4(b-d) show that all wavelengths are indeed focused at the same spatial position.

This approach becomes even more relevant for broadband illumination. Fig. 5 shows, for an ultrashort pulse, how the intensity at the focal point increases over the background when the number $N_\lambda$ of different spectral components set to one in $E^{target}$, all other components being set to zero. In Fig. 5(c), all 6 independent spectral components of the broadband pulse measured with the MSTM are controlled, showing clear spatial focusing. The signal to backroung ratio (SBR), defined as the ratio between the average value of all the pixels inside the FWHM and the average value of the speckle image before focusing, is 30.7. Since these are linear images, they are not sensitive to possible fluctuations of the temporal duration of the pulse at this position. When we use the speckle as reference, the relative phase between the different spectral components is unknown, and no temporal compression of the pulse occurs.

**Turning the medium into a generalized grating.** As a last example, Fig. 6 demonstrates another possible use of the MSTM by focusing different wavelengths to different positions. In this way the scattering medium is used as a deterministic dispersive optical element ³¹. In Fig. 6(a), three different wavelengths of a broadband pulse are focused at three arbitrary positions along a line. Scanning the wavelength of a monochromatic laser reveals that each focus corresponds indeed to a different wavelength within the
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Simultaneous focusing of different wavelengths at the same position (a) A transversal cut of the MSTM (grey area) allows one to spatially focus on a given CCD pixel all the independent wavelengths within the bandwidth Δλ<sub>L</sub>. (b-d) Focusing images for a monochromatic beam at various frequency for this optimal phase mask: (b) λ = 797 nm, (c) for λ = 800 nm, and (d) λ = 803 nm. The focusing efficiency is reduced by a factor roughly N<sub>λ</sub> compared to Fig 3(a)(e) and (i) since the phase mask tries to focus many wavelengths simultaneously.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Polychromatic focusing through a strongly scattering medium with the MSTM. (a-c) The contrast of the spatial focus increases with N<sub>λ</sub>, number of controlled wavelength: (a) SBR = 7.2 for N<sub>λ</sub> = 1, (b) SBR = 26.1 for N<sub>λ</sub> = 3, and (c) SBR = 30.7 for N<sub>λ</sub> = 6.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Implementation of a generalized grating through a strongly scattering medium with the MSTM. (a) Simultaneous spatial focusing in different output points of three independent wavelengths (λ = 797 nm, λ = 800 nm and λ = 803 nm) within the bandwidth Δλ~L~ of a polychromatic pulse using the MSTM. (b-d) Focusing images of a monochromatic beam at (b) λ = 797 nm, (c) λ = 800 nm, and (d) λ = 803 nm, using the same phase mask of (a). The three beams are focused in independent foci corresponding to the three output points in (a), as highlighted by the dashed lines.</figcaption></figure>

pulse (797 nm in Fig. 6(b), 800 nm in Fig. 6(c) and 803 nm in Fig. 6(d)). These results can be interpreted as a generalized grating. In terms of resolution and efficiency, the main difference between our approach and a conventional dispersive element is that, in our case, scattering by the medium itself determines both spectral resolution and spatial localization. Spectral resolution is given by the spectral correlation bandwidth of the medium, linked to the confinement time of the light in the medium, whereas the latter is diffraction-limited, i.e. limited by the effective numerical aperture of the medium. On the other hand, the efficiency of the spatio-spectral focusing, given by the SBR, is ultimately limited by how many degrees of freedom are controlled, i.e. how many SLM pixels and how many speckle grains are controlled simultaneously. In practice, the limitation in number of degrees of control is due to the measurement time, which is limited by the stability of the experimental setup and the medium. Here the measurement time is of the order of 15 minutes for 256 input modes and 6 wavelengths. This can easily be scaled up by increasing the stability and the modulation speed. The spectral control is, in this case, deterministic³. Interestingly, unlike a grating, the spatial dispersion of the different spectral foci can be arbitrary as demonstrated in Fig. 6, since each independent wavelength can be focused at any spatial position.

**Conclusion**

In conclusion, we introduced a method to measure the spatially and spectrally resolved transmission matrix, or MSTM, of a multiply scattering medium. Compared to other approaches based on optimization, the main advantage of our technique is that it is deterministic. Exploiting the information contained in the MSTM, we are able to selectively focus one or many spectral components of a broadband light source to one or many arbitrary spatial positions. As a proof of principle, we demonstrate multispectral focusing of a short pulse, and show the medium can be turned into a deterministic dispersive optical element, thus generalizing the concept of grating. Here we demonstrate spatiospectral control, and our results are valid not only for short laser pulse, but also for temporally incoherent sources such as SLEDs or supercontinuum sources. Further control of the spectral phase with a plane wave as reference should enable deterministic spatiotemporal control of an ultrashort pulse.

**References**

1. Sebbah, P. *Waves and imaging through complex media*. (Springer, 2001).

2. Ishimaru, A. *Wave propagation and scattering in random media*. (Wiley-IEEE Press, 1999).

3. Stoytchev, M. & Genack, A. Z. Microwave transmission through a periodic three-dimensional metalwire network containing random scatterers. *Phys. Rev. B*. **55**, 14 (1997).

4. Sheng, P. *Introduction to Wave Scattering, Localization, and Mesoscopic Phenomena*. (Academic Press, New York, 1995).

5. Wang, L. V. & Wu, H. *Biomedical Optics: Principles and Imaging*. (Wiley, 2007).

6. Goodman, J. *Speckle Phenomena in Optics*. (Roberts and Company Publishers, 2010).

7. Mosk, A. P., Lagendijk, A., Lerosey, G. & Fink, M. Controlling waves in space and time for imaging and focusing in complex media. *Nat. Phot.* **6**, 283292 (2012).

8. Vellekoop, I. & Mosk, A. Focusing coherent light through opaque strongly scattering media. *Opt. Lett.* **32**, 12 (2007).

9. Vellekoop, I., Lagendijk, A. & Mosk, A. Exploiting disorder for perfect focusing. *Nat. Phot.* **4**, 320–322 (2010).
10. Vellekoop, I., Van Putten, E., Lagendijk, A. & Mosk, A. Demixing light paths inside disordered metamaterials. *Opt. Exp.* **16**, 67–80 (2008).

11. Beenakker, C. W. J. Random-matrix theory of quantum transport. *Rev. Mod. Phys.* **69**, 731 (1997).

12. Popoff, S. M. *et al.* Measuring the transmission matrix in optics: an approach to the study and control of light propagation in disordered media. *Phys. Rev. Lett.* **104**, 100601 (2010).

13. Popoff, S. M., Lerosey, G., Fink, M., Boccara, A. C. & Gigan, S. Image transmission through an opaque material. *Nat. Commun.* **1**, 81 (2010).

14. Choi, Y. *et al.* Overcoming the diffraction limit using multiple light scattering in a highly disordered medium. *Phys. Rev. Lett.* **107**, 2 023902 (2011).

15. Kim, M. *et al.* Maximal energy transport through disordered media with the implementation of transmission eigenchannels. *Nat. Phot.* **6**, 581–585 (2012).

16. Čžimár, T. & Dholakia, K. Shaping the light transmission through a multimode optical fibre: complex transformation analysis and applications in biophotonics. *Opt. Exp.* **19**, 18871–18884 (2011).

17. Yoon, C. *et al.* Experimental measurement of the number of modes for a multimode optical fiber. *Opt. Lett.* **37**, 21 (2012).

18. Papadopoulos, I. N., Farahi, S., Moser, C. & Psaltis, D. Dynamic bending compensation while focusing through a multimode fiber. *Biomed. Opt. Exp.* **4**, 260270 (2013).

19. Bianchi, S. *et al.* Focusing and imaging with increased numerical apertures through multimode fibers with micro-fabricated optics. *Opt. Lett.* **38**, 4935–4938 (2013).

20. T. Chaigne, *et al.* Controlling light in scattering media non-invasively using the photoacoustic transmission matrix. *Nat. Phot.* **8**, 5864 (2014).

21. Choi, Y. *et al.* Measurement of the time-resolved reflection matrix for enhancing light energy delivery into a scattering medium. *Phys. Rev. Lett.* **111**, 243901 (2013).

22. Small, E., Katz, O., Guan, Y. & Silberberg, Y. Spectral control of broadband light through random media by wavefront shaping. *Opt. Lett.* **37**, 16 (2012).

23. Van Beijnum, F., Van Putten, E. G., Lagendijk, A. & Mosk, A. P. Frequency bandwidth of light focused through turbid media. *Opt. Lett.* **36**, 3 (2011).

24. Faez, S., Johnson, P. M. & Lagendijk, A. Varying the effective refractive index to measure optical transport in random media. *Phys. Rev. Lett.* **103**, 053903 (2009).

25. Curry, N. *et al.* Direct determination of diffusion properties of random media from speckle contrast. *Opt. Lett.*, **36**, 17 (2011).

26. Lemoult, F., Lerosey, G., De Rosny, J. & Fink, M. Manipulating spatiotemporal degrees of freedom of waves in random media. *Phys. Rev. Lett.* **103**, 173902 (2009).

27. Aulbach, J., Gjonaj, B., Johnson, P., Mosk, A. & Lagendijk, A. Control of light transmission through opaque scattering media in space and time. *Phys. Rev. Lett.* **106**, 103901 (2011).

28. Katz, O., Small, E., Bromberg, Y. & Silberberg, Y. Focusing and compression of ultrashort pulses through scattering media. *Nat. Phot.* **5**, 372 (2011).

29. Paudel, H. P., Stockbridge, C., Mertz, J. & Bifano, T. Focusing polychromatic light through strongly scattering media. *Opt. Exp.* **21**, 14 (2013).

30. McCabe, D. J. *et al.* Spatio-temporal focusing of an ultrafast pulse through a multiply scattering medium. *Nat. Commun.* **2**, 447 (2011).

31. Freund, I. *Phys. A*, Looking through walls and around corners. **168**, 4965 (1990).

32. Park, J. H., Park, C., Yu, H., Cho, Y. H. & Park, Y. K. Active spectral filtering through turbid media. *Opt. Lett.* **37**, 15 (2012).

33. García, P., Sapienza, R., Blanco, Á. & López, C. Photonic glass: a novel random material for light. *Adv. Mater.* **19**, 2597 (2007).

34. Redding, B., Liew, S. F., Sarma, R. & Cao, H. Compact spectrometer based on a disordered photonic chip. *Nat. Phot.* **7**, 746751 (2013).

## Acknowledgements

This work was funded by the European Research Council (grant no. 278025), ANR ROCOCO and Programme Emergence(s) from the city of Paris. O.K. was supported by the Marie Curie Intra-European fellowship for career development (IEF). We thank David Martina and Mickael Mounaix for helpful discussions.

## Author Contributions

S.G. (Gigan) conceived the project. D.A., S.G. and S.G. planned the experiment. D.A. performed the experiment and analyzed the data. S.P. proposed the analytic formalism for spatiospectral control. D.A., G.V. and S.G. (Gigan) wrote the manuscript. S.G., D.A., G.V., O.K., S.P. and S.G. discussed the results and revised the manuscript.

## Additional Information

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Andreoli, D. *et al.* Deterministic control of broadband light through a multiply scattering medium via the multispectral transmission matrix. *Sci. Rep.* **5**, 10347; doi: 10.1038/srep10347 (2015).

<figure><img src="image_8.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/