###### Abstract

We study theoretically how multiple scattering of light in a disordered medium can spontaneously generate quantum correlations. In particular we focus on the case where the input state is Gaussian and characterize the correlations between two arbitrary output modes. As there is not a single all-inclusive measure of correlation, we characterise the output correlations with three measures: intensity fluctuations, entanglement, and quantum discord. We found that, while a single mode coherent state input can not produce quantum correlations, any other Gaussian input will produce them in one form or another. This includes input states that are usually regarded as more classical than coherent ones, such as thermal states, which will produce a non zero quantum discord.

Quantum correlation of light scattered by disordered media

Ilya Starshynov, Jacopo Bertolotti, Janet Anders

University of Exeter, Stocker Road, Exeter EX4 4QL, United Kingdom

is283@exeter.ac.uk

OCIS codes: (290.4210) Multiple scattering; (290.7050) Turbid media; (270.0270) Quantum optics.

## References and links

- [1] P. Sheng, Introduction to Wave Scattering, Localization and Mesoscopic Phenomena (Springer, 2006).
- [2] J.W. Goodman, Speckle Phenomena in Optics (Roberts & Company, 2009).
- [3] E. Akkermans and G. Montambaux, Mesoscopic Physics of Electrons and Photons (Cambridge University Press, 2011).
- [4] S. Feng, C. Kane, P.A. Lee, and A.D. Stone, “Correlations and Fluctuations of Coherent Wave Transmission through Disordered Media,” Phys. Rev. Lett. 61, 834 (1988).
- [5] R. Berkovits, and S. Fen, “Correlations in Coherent Multiple Scattering,” Phys. Rep. 238, 135 (1994).
- [6] B.J. Berne, Dynamic Light Scattering: With Applications to Chemistry, Biology, and Physics (Dover, 2003).
- [7] D.A. Boas, and A.K. Dunn, “Laser speckle contrast imaging in biomedical optics,” J. Biomed. Opt. 15, 011109 (2010).
- [8] J.C. Dainty, ed. Laser Speckle and Related Phenomena (Springer, 1984).
- [9] I.M. Vellekoop, and A.P. Mosk, “Focusing coherent light through opaque strongly scattering media,” Opt. Lett. 32, 2309 (2007).
- [10] J.Bertolotti, E.G. van Putten, C. Blum, A. Lagendijk, W.L. Vos, A.P. Mosk, “Non-invasive imaging through opaque scattering layers,” Nature 491, 232 (2012).
- [11] O. Katz, P. Heidmann, M. Fink, and S. Gigan, “Non-invasive single-shot imaging through scattering layers and around corners via speckle correlations,” Nat. Phot. 8, 784 (2014).
- [12] P. Lodahl and A. Lagendijk, “Transport of Quantum Noise through Random Media,” Phys. Rev. Lett. 94, 153905 (2005).
- [13] J.R. Ott, N.A. Mortensen, and P. Lodahl, “Quantum Interference and Entanglement Induced by Multiple Scattering of Light,” Phys. Rev. Lett. 105, 090501 (2010).
- [14] S.R. Huisman, T.J. Huisman, T.A.W. Wolterink, A.P. Mosk, and P.W.H. Pinkse, “Programmable multiport optical circuits in opaque scattering materials,” Opt. Expr. 23, 229054 (2015).
- [15] H. Defienne, M. Barbieri, B. Chalopin, B. Chatel, I.A. Walmsley, B.J. Smith, S. Gigan, “Nonclassical light manipulation in a multiple-scattering medium,” Opt. Lett. 39, 6090 (2014).
- [16] R. Glauber, “Quantum theory of coherence,” in Quantum optics, S. Kay and A. Maitland, ed. (Academic Press, 1970).
- [17] P. Lodahl, A. P. Mosk, A. Lagendijk, “Spatial quantum correlations in multiple scattered light,” Phys. Rev. Lett. 95, 173901 (2005).
- [18] R. Loudon, The Quantum Theory of Light, (OUP, 2000).
- [19] J. Laurat, G. Keller, J. A. Oliveira-Huguenin, C. Fabre, T. Coudreau, A.Serafini, G. Adesso, F. Illuminati, “Entanglement of two-mode Gaussian states: characterization and experimental production and manipulation,” J. Opt. B 7, S577–S587 (2005).
- [20] C. W. Beenakker, “Random-matrix theory of quantum transport,” Rev. Mod. Phys. 69, 731–808 (1997).
- [21] G. Adesso, A. Datta, “Quantum versus Classical Correlations in Gaussian States,” Phys. Rev. Lett. 105, 030501 (2004).
- [22] A. Peres, “Separability Criterion for Density Matrices,” Phys. Rev. Lett. 77, 1413-1415 (1996).
- [23] M. Horodecki, P. Horodecki, R. Horodecki, “Separability of mixed states: necessary and sufficient conditions,” Phys. Rev. A 223, 1-8 (1996).
- [24] R. Simon, N. Mukunda, B. Dutta, “Quantum-noise matrix for multimode systems: U(n) invariance, squeezing, and normal forms,” Phys. Rev. A 49, 1567-1583 (1994).
- [25] R. Simon, “Peres-Horodecki separability criterion for continuous variable systems,” Phys. Rev. Lett. 84, 2726 (1999).
- [26] Lu-Ming Duan, G. Giedke, J. I. Cirac, P. Zoller, “Inseparability Criterion for Continuous Variable Systems,” Phys. Rev. Lett. 84, 2722-2752 (2000).
- [27] A. Datta, S. T. Flammia, C. M. Caves, “Entanglement and the power of one qubit,” Phys. Rev. A 72, 042316 (2005).
- [28] H. Ollivier, W. H. Zurek, “Quantum Discord: A Measure of the Quantumness of Correlations,” Phys. Rev. Lett. 88, 017901 (2001).
- [29] S. Rahimi-Keshari, C. M. Caves, and T. C. Ralph, “Measurement-based method for verifying quantum discord,” Phys. Rev. A 87, 012119 (2013).
- [30] S. Hosseini, S. Rahimi-Keshari, J. Y. Haw, S. M. Assad, H. M. Chrzanowski, J. Janousek, T. Symul, T. C. Ralph, and P. K. Lam, “Experimental verification of quantum discord in continuous-variable states,” J. Phys. B At. Mol. Opt. Phys. 47, 025503 (2014).
- [31] A. S. Holevo and R. F.Werner, “Evaluating capacities of bosonic Gaussian channels,” Phys. Rev. A 63, 032312 (2001).
- [32] P. Giorda, M. Paris, “Gaussian quantum discord,” Phys. Rev. Lett. 105, 020503 (2010).
- [33] L. Mandel, E. Wolf. Optical Coherence and Quantum optics (Cambridge University Press, 1995).
- [34] M. C. Teich, B. E. Saleh, “Squeezed states of light,” Quantum Opt. 1, 153â191 (1989).
- [35] M. Paris, “Entanglement and visibility at the output of a Mach-Zehnder interferometer,” Phys. Rev. Lett. 59, 1615-1621 (1998).
- [36] A. Gatti, E. Brambilla, M. Bache, L. A. Lugiato, “Ghost imaging with thermal light: Comparing entanglement and classical correlation,” Phys. Rev. Lett. 93, 093602 (2004).
- [37] Y. Lahini, Y. Bromberg, D. N. Christodoulides, and Y. Silberberg, “Quantum correlations in two-particle Anderson localization,” Phys. Rev. Lett. 105, 163905 (2010).
- [38] A. Ferraro and M. G. a Paris, “Nonclassicality criteria from phase-space representations and information-theoretical constraints are maximally inequivalent,” Phys. Rev. Lett. 108, 260403 (2012).
- [39] B. E. A. Saleh, A. F. Abouraddy, A. V. Sergienko, and M. C. Teich, “Duality between partial coherence and partial entanglement,” Phys. Rev. A 62, 043816 (2000).
- [40] T. Wellens and B. Gremaud, “Nonlinear Coherent Transport of Waves in Disordered Media,” Phys. Rev. Lett. 100, 033902 (2008).

## Appendix A Introduction

Light propagation in a scattering medium can often be described as a diffusive process [1], where any memory of the initial state is lost almost immediately, and transport is represented by the incoherent sum over many Brownian random walks. However it was early realized that temporal coherence survives multiple elastic scattering, and thus interference is still possible even after passing through a diffusive medium [2]. As a consequence of interference the diffusive picture has to be modified, and the light scattered into different channels develop correlations [3]. Classical correlations between the intensity scattered in one direction and the intensity scattered into another direction were studied since the ’60s [4, 5] and found many applications, especially in imaging, e.g. in dynamic light scattering [6], speckle contrast imaging [7], stellar speckle interferometry [8], wavefront shaping [9], and speckle scanning microscopy [10, 11]. The study of quantum correlations of multiply scattered light started much later, probably due to the widespread idea that quantum features are frail and thus unlikely to play a significant role in the presence of strong disorder. Nevertheless it was recently shown that not only certain quantum features can survive multiple scattering [12], but that quantum correlations can even be spontaneously created in the scattering process [13]. The quantum properties of scattered light proved to be so robust that control over the transport of single photons in a disordered medium via wavefront shaping was demonstrated [14, 15].

In this paper we consider the case of a generic Gaussian input state and study theoretically the necessary conditions for the output state to present quantum correlations. In particular we show that, while entanglement requires specific conditions to be produced, quantum discord is always present in the output state when a thermal input state is used but, surprisingly, not when the input mode is in a coherent state.

## Appendix B Quantum correlations between two modes

Classical correlations of scattered light are commonly described by a correlation function that measures the correlations between intensity fluctuations of two different modes and [3, 4, 5],

(1) |

where is the light intensity and represents either a time or an ensemble average (if the system is ergodic the two are equivalent). can be used beyond the classical case to study certain classes of quantum correlations by substituting with the mode’s number operator [16, 17], and since classical light can never lead to , a correlation value below 1 is considered a clear signature of quantumness [18].

In order to study the effect of multiple scattering on a Gaussian input state we employ the formalism of density operators and covariance matrices [19]. The covariance matrix of two output modes is defined as

(2) |

where is the vector of quadrature operators, that satisfies the commutation relation , which imposes constrains on :

(3) |

where is a direct sum. Importantly, for Gaussian states the covariance matrix captures all the correlation properties between the modes. As will be discussed below, different measures of quantum correlations, such as entanglement and discord, can be formulated as conditions on the elements of the covariance matrix.

When monochromatic light propagates through a disordered medium it is multiply scattered, and its wavefront becomes completely irregular. As a result the output modes take the form of a speckle pattern [2]. Here we consider a model of linear scattering, such that light propagation can be described by the scattering matrix [20], which couples the fields of the output modes, , with the input modes, [17]. The electromagnetic field of a mode with polarization , wave vector and frequency can be expressed as:

(4) |

where and are the creation and annihilation operators of the input mode, is the vacuum permittivity and is the mode volume [18]. The scattering matrix links the input fields to the output ones as . Likewise for the ladder operators we have:

(5) |

Here and are the annihilation operators related to, respectively, the input and the output modes as shown in Fig. 1. The elements of the matrix are the complex transmission coefficients from the -th input mode to the -th output mode.

In this paper we focus on the correlations between two of the scattered modes, whose covariance matrix has the structure:

(6) |

where , , and are matrices. The matrices and refer to individual properties of the two modes, and the matrix describes the correlations between them. The determinants of these matrices, together with the determinant of the whole covariance matrix, are invariant under local transformations (i.e. those acting only on one of the modes). Because of this local invariance the determinants characterize entanglement and other non-local correlation properties of the state [21, 22, 23, 24]. Such two-mode covariance matrix can be written in the form

(7) |

keeping the invariants of the state: and the same [25].

The intensity correlation contains fourth order moments of the field distribution. For Gaussian states these can always be written as a function of the second order moments, contained in the covariance matrix [19]. This leads to the expression for the correlation function:

(8) |

From Eq. (8) we can see that Gaussian states always have , as as a consequence of the commutation relations in Eq. (3).

While is useful to describe classical correlations and certain features of quantum light, like photon anti-bunching [16], it does not capture other quantum correlations, and thus other measures have been introduced, notably entanglement and discord. The characterization of entanglement for general mixed states of multiple modes is a challenging task. Necessary and sufficient conditions for entanglement have been identified for bipartite discrete quantum systems of and dimensions [22, 23] and for continuous two-mode Gaussian states [25, 26]. In the latter case, these mathematical conditions can be understood as inverting time in one of the modes of the system and checking if the resulting quantum state is still a valid physical state, i.e. hermitian and positive [22, 23, 25, 26]. These conditions on the quantum state can be recast as conditions on the matrix , obtained from the covariance matrix Eq. (7), by changing the sign to the momentum of one of the two modes, i.e. . If obeys the commutation relation stated in Eq. (3) then the two modes are separable, otherwise they are entangled. This separability condition for the two modes can be rewritten as [21],

(9) |

where refers to the sub-matrices of the covariance matrix and (see Eq. (6) and Eq. (7)). are the symplectic eigenvalues of the matrix [19].

When the parameters and (i.e. the number of photons in each of the modes) are fixed, Eq. (3) and Eq. (9) define closed regions in the and parameter space (see Fig. 2). The condition in Eq. (3) defines an area enclosed by the thick solid line in Fig. 2, which corresponds to all valid states. The entangled states lie within this areas but outside the region defined by the condition in Eq. (9) (gray areas in Fig. 2), which can be graphically obtained as the first region rotated by (dotted lines in Fig. 2).

For a long time entanglement was considered to be the quintessential non-classical ingredient for quantum computation and communication. However, it was shown that even using separable (i.e. non-entangled) states it is possible to perform certain computational tasks exponentially faster than any known classical algorithm [27]. This means that just the fact of the system being quantum may lead to applications impossible for classical systems, irrespectively of separability. This fact has been formalized by introducing quantum discord [28], a measure which broadens the concept of non-classical correlations to states without entanglement.

For a pure bipartite entangled state the measurement of one of the parties completely determines the result of the measurement on the second, while for classical states measurement on one party does not effect the other. On the other hand there are separable states where the measurement of one party influences the other in a probabilistic sense [29, 30], making them non-classical, but not entangled either. Quantum discord quantifies this influence and is defined as the discrepancy of two measures of mutual information for a joint quantum state of systems and :

(10) |

Here and are the reduced states of and is the von Neumann entropy: . The second equation depends on the measurement choice , and is the state of the system A conditioned on a certain outcome of this measurement performed on the subsystem B.

Quantum discord is quantified by the difference between the two expressions in Eq. (10) minimized over all possible sets of measurement operators,

(11) |

When the input state is Gaussian and its covariance matrix is , the expression for the entropy of this state reduces to: , where are the symplectic eigenvalues of and [31]. Here we limit ourselves to Gaussian measurements, i.e. those that preserve the Gaussian nature of the state they are applied to. Under this assumptions the minimization of Eq. (11), written as a function of the invariants of the covariance matrix from Eq. (7), gives:

(12) |

which is known as Gaussian discord [32]. The Gaussian discord vanishes if and only if [21, 29, 32].

## Appendix C Correlations between two scattered modes

To characterize the correlation of two modes of the scattered light we explicitly calculate the elements of the corresponding covariance matrix. We will consider the experimentally common situation where the input light is in a single mode , and all the other input modes are in a vacuum state, as depicted in Fig. 1. Replacing the quadrature operators in the vector with the ladder operators: and , and substituting Eq. (5) in, the general expression for the elements of the covariance matrix of two output modes is

(13) | ||||

where

This expression allows us to analyse the correlation properties between two output modes, e.g. entanglement and quantum discord, for different possible inputs.

### Coherent state:

If the input mode is in a coherent state Eq. (8) gives the value of as expected, which means that there is no correlation between the intensity fluctuations of the two output modes according to this measure [33]. Since all expectation values of the operators and are 0, the covariance matrix of the output modes will be . Substituting the elements of into Eq. (9), the allowed region shrinks to a point . This means that for any coherent state as an input, any two output modes will simply be a product of two coherent states and no quantum correlations are present.

### Squeezed state:

If the input mode is in a squeezed state with a squeezing parameter and phase , the expectation values for the operators in Eq. (13) are: , and [34]. If we set , for which we expect maximal entanglement [19], we can express the entanglement criterion as , which is always true if and both scattering matrix elements are non-zero. This means that we will get entanglement for any non-zero degree of squeezing [35]. It is remarkable that the degree of entanglement does not depend on the phases of transmission coefficients of the scattering matrix, but only on their moduli.

Although quantum features of intensity fluctuation correlations are often linked to , implying a reduction of coincidences in simultaneous detections of photons in the two modes, for squeezed states entanglement leads to positive correlation, which can lead to high . The maximal possible value of which could be achieved with squeezed entangled states (the top left or bottom right points of the region of allowed states in Fig. 2) is:

(14) |

where is the average number of photons in the input mode. When , approaches 2, which corresponds to the value expected for thermal states. The presence of entanglement allows to reach values of inaccessible for thermal states, and this is exploited in quantum imaging where it can allow faster recovery of information, especially in the low photon number regime [36]. However, this needs to be treated with some care. In fact, as it can be seen from Fig. 2b, with the increase of and the circle and the boundary of the gray area (corresponding to the entangled states) cross, and therefore it is possible to find non-entangled states with higher intensity correlations than some entangled states.

### Thermal state:

For a thermal state we have: and . Using Eq. (13) the covariance matrix of the two output modes and can be written as:

(15) |

with , , and .

The values of and in this case are equal. In Fig. 2 all possible thermal states lie on the line and thus no entanglement is possible according to the criterion described in Eq. (9). In fact, in order for the modes to be entangled, and should at least have different signs [25]. Although these states do not show entanglement, there are still quantum correlations between the output modes in the form of quantum discord.

We calculated the discord of the output covariance matrix for the case of input mode in the thermal state by substituting the output covariance matrix Eq. (15) into the formula for the Gaussian discord in Eq. (12). In Fig. 3 we plot the dependence of the Gaussian discord on the absolute values of the transmission coefficients from the mode to the modes and . Notice that the discord is asymmetric against these coefficients since the measurement is performed on only one of the modes (on the mode in the Fig. 3). The discord increases monotonously with , but there can be a maximum in its dependence on , the position of which is defined by the number of photons in the input mode (Fig. 3b). At low photon numbers there is no maximum, and in that case the discord increases monotonously (Fig. 3a).

The mean discord obtained can be calculated noticing that, on average, the energy distributes equally among all the possible scattering channels and thus . If we assume that the scattered light follows a Rayleigh distribution we have [3], which can be substituted into Eq. (12) to get the average amount of discord in a pair of output modes in this configuration. As shown in Fig. 4, increases monotonically with , but decreases monotonically with . As a consequence the best conditions to observe the discord generated by multiple scattering of a thermal state of light are obtained for an intense light signal scattering over a system with a small number of channels. Therefore we suggest that light scattering from systems showing Anderson localization [1, 3] will show a significant amount of quantum discord.

## Appendix D Conclusions

The search for a universal criterion that captures all the nuances of non-classical correlations is an object of ongoing intensive discussion [38, 39]. The presented results contribute to this debate by providing an illustration of the differences between various measures of quantum correlations, such as the correlation function based on intensity fluctuations , entanglement and discord. We calculated the covariance matrix of two arbitrary output modes of the light elastically scattered by a disordered material for different states of the input mode, and analysed their correlation properties. Surprisingly, the results show that if the input is a thermal state then any two output modes will be (Gaussian) discorded, a signature of the quantum character of light. Moreover, it turns out that coherent states are the only Gaussian input that do not produce quantum correlations, as measured by any of the quantities considered here.

It is known that the propagation of light through a scattering medium is modified by quantum interference when the input state is entangled [37], but the effects of quantum discord on light propagation are still a largely unexplored subject. Quantum discord appears naturally from the multiple scattering of thermal light, even for large photon numbers. Such macroscopic effects can potentially be exploited to develop novel imaging techniques. Since the amount of expected quantum discord grows when the number of scattering channels is small, these effects will play a role especially in the case of strongly scattering materials, where the dimensionless conductance is small [3]. In particular we expect it to have an effect for systems that show Anderson localization.

Finally, when light undergoes multiple scattering, even very weak nonlinearities can have a dramatic effect [40]. Thus diffusion through a nonlinear system will produce a very rich landscape of possible output states, opening the possibility to generate multimode entanglement from classical input light.

## Appendix E Acknowledgements

We are grateful to M. Paternostro, D. Browne, and M. Williamson for insightful discussions. JA acknowledges support by EPSRC (EP/M009165/1). JB acknowledges support from the Leverhulme Trust’s Philip Leverhulme Prize. IS acknowledges support from EPSRC through the Centre of Doctoral Training in Metamaterials ().