Statistics and control of waves in disordered media

Fundamental concepts in the quasi-one-dimensional geometry of disordered wires and random waveguides in which ideas of scaling and the transmission matrix were first introduced are reviewed. We discuss the use of the transmission matrix to describe the scaling, fluctuations, delay time, density of states, and control of waves propagating through and within disordered systems. Microwave measurements, random matrix theory calculations, and computer simulations are employed to study the statistics of transmission and focusing in single samples and the scaling of the probability distribution of transmission and transmittance in random ensembles. Finally, we explore the disposition of the energy density of transmission eigenchannels inside random media.


Introduction
The scattering of waves by disorder is more the rule than the exception in condensed materials.Scattering is the source of electrical resistance and makes materials opaque to classical waves.Whereas, all ordered samples of a given material are alike, every disordered sample is disordered in its own way.This can be seen directly in the distinct grainy pattern of intensity of scattered laser light for each sample [1].The speckle pattern within or beyond an illuminated sample arises from the interference of waves following all possible paths to any point.Though the wave field appears to be completely random in disordered samples, it is still possible to manipulate the incident wavefront to control the net transmission and to focus the transmitted wave.This was first achieved by focusing light to a point utilizing feedback from the transmitted wave to manipulate the input wave [2].In principle, a considerable measure of control over the output wave can be exerted by tailoring the incident wave if the field transmission coefficients between all incident and outgoing channels are known [3].These coefficients are the elements of the transmission matrix (TM).In this review, we consider the use of the TM to describe the propagation and control of classical waves in random media.This will be placed in the electronic context in which the TM in the quasi-one-dimensional (quasi-1D) wire geometry was first utilized by Fisher and Lee [4] to describe electronic conductance and by Dorokhov to describe the scaling of the eigenvalues of the transmission matrix and the conductance [5,6].Because of the strong analogy between electronic and classical wave transport, it is possible to demonstrate fundamental predictions for electron transport in equivalent settings for classical waves in cases in which the predictions have not been demonstrated for electrons.In addition, consideration of the TM for classical waves has stimulated the investigation of universal aspects of transmission eigenvalues such as the dynamics and spatial profiles of the energy density of transmission eigenchannels, which cannot be explored directly in conductance measurements.
The TM is a powerful tool for studying transport in disordered media in which the random field pattern excited by a monochromatic wave is coherent in time throughout the sample [7][8][9].For classical waves, thermal fluctuations in position of the sample's constituents in static samples are much smaller than the wavelength so that the wave may be temporally coherent even on a macroscopic scale.In contrast, the electron wave is only phase-coherent in micron-sized conductors at ultralow temperature.Such samples lie between the microscopic atomic scale and macroscopic lengths.The suppression of diffusion [10] and the enhancement of fluctuations in conductance [11][12][13] due to interference in samples of intermediate or mesoscopic dimensions have been the focus of studies of electronic transport.
Anderson [10] showed that electrons in an unbounded lattice would be exponentially localized when the range of random values of the electron energy at a site is sufficiently large even when they are not trapped energetically.This leads to the absence of electron diffusion.Unlike the speckle pattern that is washed out by averaging over disorder, Anderson localization is an interference effect that survives averaging over random configurations.The suppression of average transport is a consequence of the constructive interference of partial waves that follow the same path returning to a coherence volume in the sample but in time-reversed order [14,15], as illustrated in Fig. 1.The transverse dimensions of the coherence volume is the field correlation length λ /2.Because the amplitudes and phases associated with such time-reversed partial waves are identical, the amplitude for return for paths following the loop for both senses is twice that for a single partial wave following the loop in either direction.The probability of return of a particle is therefore four times that for a single path and thus twice the probability given by the incoherent sum of probabilities for the two partial waves of each loop.This is the basis of weak localization [7] which reduces average transport.When the probability of return to a coherence volume reaches unity, the wave becomes localized.The effect of weak localization is observed directly as coherent backscattering in the domain of the reflected wavevector [16][17][18].Reflection is enhanced by a factor of two over the background level in the retroreflection direction.The angular distribution of backscattered light is the Fourier transform of the intensity distribution of the point spread function for light incident on the incident surface.The enhancement falls in an angular range of 1/k , where k is the wavenumber and the transport mean free path.Measurements of coherent backscattering therefore allow for the determination of the mean free path.
Ioffe and Regel [19] suggested that, an electron cannot be considered to be freely propagating if the scattering mean free path were only a fraction of the wavelength.This yields the localization criterion in three dimensions, k < 1.This condition is also given by a calculation of the return probability based on diffusion theory in three dimensions.Subsequently, Thouless [20,21] argued that localization in any dimension should depend only upon the dimensionless ratio of the average width and spacing of energy levels, which has come to be known as the Thouless number, δ = δ E/∆E.Electronic levels correspond to quasi-normal modes or resonances for classical waves in open systems.The width of levels, which in units of angular frequency is equal to the leakage rate of energy from the sample, is closely linked to the sensitivity of level energies to changes at the boundary since both are proportional to the ratio of the strengths of the mode at the boundary relative to the interior of the sample.When the wave is localized within the interior of a sample, the amplitude of the wave is exponentially small at the boundary and the mode is weakly coupled to the surrounding medium.Its lifetime is then long and its linewidth correspondingly narrow, so that, δ E < ∆E.On the other hand, when the wave is diffusive, modes extend throughout the medium; the wave then couples readily to its surroundings and the width of a typical level overlaps several modes, δ E > ∆E.Thus, the localization threshold occurs at δ = 1.
Thouless [20,21] showed that for diffusive samples, δ is equal to the conductance in units of the quantum of conductance δ = g ≡ G/(e 2 /h).Since the level width falls inversely with the square of sample length for diffusive samples while the spacing between levels falls only inversely with the sample length in the quasi-1D wire geometry of constant cross section, δ falls inversely with sample length.Thus there will always be a crossover from diffusive to localized wave transport as the length of a wire increases.The crossover will occur at a sample length L = ξ = N , at which δ = g = 1 [21].Here, ξ is the localization length of the sample and N is the number of channels.The channels can be the propagating transverse modes in the ideal leads attached to a conductor at a given voltage.The total number of channels is approximately 2πA/λ 2 including two electronic spin states, where A is the constant cross sectional area of the sample and λ is the electron wavelength.Abrahams et al. [22] argued that the variation of average g with the dimensions of the sample should depend only upon the value of g and the dimensionality.A localization transition is predicted to exist only in 3D in which g falls essentially exponentially with increasing sample size for g < 1 and increases for g > 1 [22].
Since the conductance exhibits large fluctuations for localized waves, the mean value of the conductance is insufficient to characterize the nature of electronic transport and the full distribution of the conductance needs to be given.Anderson et al. [23] hypothesized that the logarithm of the conductance of localized waves in one-dimensional (1D) samples follows a Gaussian distribution with a variance approaching the average of the logarithm of the conductance, ln g = −L/ξ .The scaling of the distribution of conductance for localized waves in 1D would then be completely specified by the single parameter L/ξ .Electronic conductance in disordered 1D conductors and transmission of classical waves in mesoscopic media are seen to have equivalent descriptions in the 1D Landauer relation, g = T , where T is the transmission [24,25].The Landauer relation was extended to multichannel mesoscopic systems via the TM [4].The elements t ba of the transmission matrix t are the field transmission coefficients between incoming channels it a and output channels b.The dimensionless conductance in quasi-1D samples is equal to the sum of N 2 pairs of intensity transmission coefficient, which is the transmittance, T , g = ∑ N a,b |t ba | 2 ≡ T .The dimensionless conductance or transmittance can also be expressed as the sum of the transmission eigenvalues, g = ∑ N n τ n .Here, the τ n are the eigenvalues of the matrix product tt † .The transmission eigenvalues can also be obtained via the singular value decomposition of the TM, t = UΛV † , where U and V are unitary matrices with columns which are the singular vectors yielding the output and input of the transmission eigenchannels and Λ is the diagonal matrix with elements √ τ n .Dorokhov [5,6] considered the scaling of g in samples made up of N parallel 1D disordered chains with transverse coupling to create a quasi-1D sample.He found that, instead of a single localization length ξ , there exists N auxiliary localization lengths, ξ n , which describe the scaling of the averages of each of the transmission eigenvalues vs. length via the relation, The inverse auxiliary localization lengths are are given by 1/ξ n = (2n − 1)/2ξ for integers n < N/2 so that these are equally spaced with spacing 1/ξ n+1 − 1/ξ n = 1/ξ [5,6,[26][27][28][29][30].
An intuitive Coulomb gas model [26,31] was developed to describe the statistics of transmission eigenvalues.This is based on the model originally introduced by Dyson [32] to visualize the logarithmic repulsion between eigenvalue of a large random Hamiltonian matrix [33].In the Coulomb gas model, the τ n are associated with the positions of parallel lines of charges at x n = L/ξ n = (2n − 1)L/2ξ ) for n < N/2 and image charges of the same sign at −x n [31].The line charges in the model are embedded in a compensating continuous charge background of opposite sign.The background charge tends to screen the line charges from the interaction with other line charges.For diffusive waves, the spacing between neighboring charges is small compared to the screening length.As a result, there exists a logarithmic repulsion between the transmission eigenvalues due to the Coulomb interaction.The roughly uniform distribution of x n imposed by the mutual repulsion of charges leads to the bimodal distribution of transmission eigenvalues τ with peaks in the probability distribution for τ near unity and at exponentially small values of τ, ρ(τ) = g 2τ √ (1−τ) [6,28,[34][35][36].The repulsion between charges for diffusive waves leads to universal conductance fluctuations.In quasi-1D disordered systems with time-reversal symmetry, var(T ) = 2/15 [28].Universal conductance fluctuations were first measured in mesoscopic conductors at low temperatures [11,37].Fluctuations of conductance of a gold wire at 10 mK as a function of magnetic field applied perpendicular to the wire are shown in Fig. 2. Dorokhov [6] showed that transport through disordered conductors is dominated by g eigenvalues with ξ n > L and so transmission through these channels is of order unity.The number of such open eigenchannels with τ n > 1/e, which is essentially the conductance for diffusive quasi-1D samples, is inversely proportional to the sample length in accord with Ohm's law.For localized waves, however, the spacing between charges in the Coulomb gas model is greater than the screening length so that the Coulomb interaction between charges is diminished.The smallest of the x n , x 1 , is larger than unity for localized waves so the transmission eigenvalues fall exponentially with n, τ n+1 / τ n = exp(−2L/ξ ).As a result, the conductance is dominated by the highest transmission eigenchannel τ 1 [6,26,27,[29][30][31].Fig. 2. Fluctuations of conductance of a gold wire at 10 mK as a function of magnetic field applied perpendicular to the wire.The wire is of length 310 nm and width 25 nm.(From Ref. [37]) In contrast to electronic systems, for which only the conductance can be measured, it is in principle possible to measure the elements of the TM for classical waves.Thus the transmission coefficients for the field, t ba , intensity T ba = |t ba | 2 , total transmission T a = ∑ N b T ba and transmittance T can be measured.Spectra of these variables for a single sample supporting localized and diffusive waves in different frequency ranges in which g =6.9 and g = 0.37 are shown in Fig. 3.The statistics of random ensembles of statistically equivalent samples can be obtained in microwave measurements for a succession of sample realizations produced by briefly rotating a copper tube containing randomly positioned dielectric spheres.Ensembles of random system can be studied optically by illuminating different regions of a random slab.The statistics of internal motion of colloidal samples can also be explored in measurements of the temporal correlation of transmitted light [38][39][40].In-and out-of-phase field transmission coefficients, field E ba / T ba / , transmission coefficients of intensity s ba = T ba / T ba , total transmission s a = T a / T a and transmittance s = T / T for microwave radiation propagating through a random waveguides for a diffusive sample (left column) and a localized sample (right column).For diffusive waves, δ > 1, many modes contribute to transmission at all frequencies and for all source and detector positions.In contrast, for localized waves δ < 1, distinct peaks appear when the incident radiation is on resonance with a quasi-normal mode.The resonance condition holds for all source and detector positions and therefore sharp peaks remain even when transmission is integrated over space.
The first question that was addressed in bounded samples was the scaling of conductance.This is practically difficult to realize in electronic measurements because of the limited scale of mesoscopic samples and the difficulty of producing statistically equivalent samples.The scaling of optical transmission in the slab geometry can be measured in the slab geometry by forming a wedge of random material and focusing a laser beam on the wedge as it is translated perpendicular to its vertex [41,42].The transmission of light through a slab of rutile titania particles in a polystyrene matrix is seen in Fig. 4. The lower envelopes of the transmission curves in Figs.4a,b fall inversely with sample thickness L up to a length equal to the absorption length L a indicating that the wave is diffusive.The absorption length is given by L a = √ Dτ a , where the diffusion coefficient is D = v e /3, v e is the energy velocity in the medium and τ a is the absorption time.Spikes in transmission occur at spots in the sample at which a pit has been gouged out of the otherwise smooth surface in the polishing process.Spectra of transmitted intensity in the far field at different sample thicknesses are shown in Fig. 4(c).The correlation frequencies of fluctuations in transmission at different sample thicknesses give the average linewidth of electromagnetic modes of the medium [43].In this diffusive sample, the correlation frequency scales as L −2 for L < L a [43].The field correlation function and the temporal profile of transmission form a Fourier transform pair and give the diffusion coefficient and absorption length which is consistent with measurement of the scaling of static transmission [43].Establishing that a wave is localized through scaling is a challenge in the presence of absorption of classical waves and of dephasing of the electronic wave function [42,[44][45][46][47][48].But the magnitude of fluctuations of intensity and total transmission normalized by their ensemble average values give a straightforward measure of the degree of localization which is robust against absorption [48].The magnitudes of relative fluctuations in electronic conductance [11][12][13], classical intensity and total transmission and the degrees of spatial, spectral and temporal correlation of flux increase through the localization transition [40,43,[49][50][51][52][53][54][55][56][57][58][59][60][61][62][63].Large fluctuations in transmitted intensity were observed in the critical regime of the Anderson localization transition in acoustic measurements in a slab of brazed aluminum beads [64].These results are consistent with the agreement found between the non-exponential decay of a transmitted pulse measured in this sample and self-consistent localization theory [65].Microwave measurements in quasi-1D samples of the time-dependence of transmission in more deeply localized samples [66] show dramatically reduced decay rates for transmission compared to both diffusion theory and the self-consistent localization theory.Since the slow-down of decay rates is associated with long-lived modes [67][68][69], this indicates that the self-consistent theory is valid near the localization transition but not for deeply localized waves [65,[70][71][72].The approach to localization can also be seen in the saturation of the spread of a wave with time delay on the output surface of the sample, which is independent of absorption.This was observed in acoustic measurements [64] and in pulsed laser experiments in random slabs of rutile titania particles [73].The local strength of scattering in the optical experiments was determined from measurements of the width of the coherent backscattering cone [73,74].Values of the product of the wave vector and the transport mean free path, k were larger than the value predicted by the Ioffe-Regel criterion [19] for localization in three dimensions of k < 1.The lowest value of k for these samples was 2.5.The Ioffe-Regel criterion is consistent with the Thouless criterion for localization [21,75].Localization in samples with k > 1 might occur in random samples with residual order [76] or in samples with resonances in individual scattering elements [77][78][79][80].In such samples, the density of states (DOS) may fall below the uniform medium value 8πν 2 /v 3 p , where ν is the frequency, v p is the phase velocity.Localization might then occur from values of k > 1.
Additional evidence that these samples are in the critical regime is seen in the continuous slowing decay of pulsed transmission in these samples.This shows that the exponential decay of transmission with time delay expected for diffusive samples breaks down at long times.Departures from an exponential decay of transmission for diffusive waves with values of g ∼ 2.5 have been observed in microwave transmission [68].This indicates that the distribution of mode lifetimes increases even in diffusive samples as the crossover to localization is approached [67,68].An important question is whether the time dependence of the transverse spatial profile represent the reduced spatial extent of some pre-localized modes as the localization transition is approached or localization [67,81,82].Transverse optical localization [83] has been observed in samples that are uniform in the longitudinal direction but disordered in two [84] or one [85] transverse dimension.Whether it is possible to localize electromagnetic radiation in three dimensions in samples without residual order remains an open question [76,86].
The TM has been recently measured for light [3,[87][88][89], microwave radiation [29] and acoustic waves [90,91] in static samples.The TM gives a full description of the degree to which the incident wavefront can be manipulated to control the transmitted wave [3,92].This may be exploited to focus radiation and enhance or suppress transmission for applications in imaging and telecommunications [3,29,[87][88][89][93][94][95][96][97][98][99][100][101][102][103][104][105][106][107][108][109][110][111][112].The ability to control propagation of classical waves in a random medium was first demonstrated in acoustics by focusing a pulse within and through a scattering medium [113].This is achieved by time-reversing the signal recorded by an array of transducers.The measured time-varying signal is played back in time and a pulse emerges at the location of the source.For monochromatic radiation, it is possible to focus the transmitted wave by shaping the incident waveform.Vellekoop and Mosk [2] demonstrated optical focusing through an opaque sample by adjusting the incident wavefront reflected from a spatial light modulator using a genetic algorithm with feedback from the intensity at the focus.When the TM of the scattering sample is determined, it is possible to focus at any desired channel on the output surface.This can be achieved by producing an incident waveform t * β a , which is the phase conjugate of the field one would obtain on the incident plane when a source is placed at a point β on the output plane.This would bring the fields from different incident channels a arriving at β in phase so that the intensity at the point is greatly increased by constructive interference.Vellekoop and Mosk [92] found that the background intensity in the transmitted wavefront for optimally focused radiation is also enhanced relative to the average transmitted intensity.In addition, the total transmission can be significantly enhanced to give nearly complete transmission when the incident wavefront corresponds to the highest transmission eigenchannel.
The TM gives the fullest account of transmission and conductance and enables the control of the transmitted flux.However, the TM cannot provide the intensity profiles of the transmission eigenchannels inside the sample.Recent simulations by Choi et al. [114] explored the spatial profile of the intensity of eigenchannels inside a single random sample.They have observed that the energy stored within high transmission eigenchannels is enhanced.This was confirmed by Gérardin et.al. [91] in measurements of flexural waves propagating in a disordered elastic waveguide.Using a laser source and a heterodyne interferometer, they have measured the full scattering matrix of the medium.After a proper normalization of the measured scattering ma-trix, they found that the transmission eigenvalues follow a bimodal distribution in agreement with theoretical prediction.They have also demonstrated that the intensity associated with high transmitting eigenchannels increases within the medium with a peak located near the middle of the sample.This was presented in Fig. 5.We have shown recently in microwave measurements that the integral of the intensity inside the sample for a specific eigenchannel yields the dwell time of the eigenchannel [115].The dwell time of the each of the eigenchannel is also equal to its contribution to the DOS of the sample so that the sum of the dwell times for N eigenchannels gives the DOS of the disordered material.The average intensity falls linearly inside the medium in diffusive samples, as required by the Fick's law.For a localized sample, wave interference cannot be neglected in calculating the average profile [116,117].The intensity profile for localized waves excited by radiation incident from one side of the sample is found to fall more rapidly toward the center of the open medium [70][71][72].This profile can be found by solving a generalized diffusion equation with a position dependent diffusion coefficient, which reflects the increasing renormalization of transport with depth into the sample due to wave interference [70].Recent simulations of intensity profiles inside the material for a random ensemble have suggested a universal expression for the profile of the energy density of the transmission eigenchannels [118].These profiles may be expressed in terms of the auxiliary localization lengths of the corresponding transmission eigenchannels.The structure of these profiles is given by solutions of the generalized diffusion equation [118].
In the sections of this review that follow, we will focus on recent microwave measurements and computer simulations of spectra of the TM in quasi-1D random waveguides in the crossover from diffusive to localized waves.In Section II, we discuss the statistics of transmission in single sample realizations and the statistics of the transmittance in ensembles of random samples.The statistics in singe samples with given transmittance T are found to depend upon only a single parameter, which is the eigenchannel participation number, The probability distribution of T is observed to change from a one-sided log-normal distribution to a log-normal distribution as M → 1 as the wave becomes deeply localized.In addition, the statistics of transmission are found to approach predictions for 1D.In particular, the single parameter scaling (SPS) hypothesis predicted for 1D random samples is found experimentally to hold for deeply localized waves in a quasi-1D geometry.Focusing of both static and pulsed radiation through and within opaque samples is discussed in Section III.Expressions for the spatial profile of the focused beam and the contrast between the focused and background intensities are demonstrated.In Section IV, we show that a complementary set of parameters to the transmission eigenvalues provides the dwell time and the contribution of each eigenchannel to the DOS.These parameters correspond to the integral of the energy density profile inside the sample for the corresponding transmission eigenchannel.A summary of recent results and prospects for applications are given in in section V.

Statistics of transmission in single samples and ensembles of random samples
The statistics of transmission in a single disordered sample are of particular interest since applications are typically in a specific sample.Such statistics are an essential part of the statistics of transmission over a collection of random samples.The transmitted field with a given polarization for an incident channel a in a single random speckle pattern is a complex Gaussian random variable.Since the transmission is the square of the field on the output surface, T ba = |t ba | 2 the corresponding probability distribution of intensity relative to the average in the speckle pattern,  Using the singular value decomposition of the TM, the variance of the total transmission for an incident channel a relative to the average over all possible incident channels can be expressed as [119], For diffusive samples without absorption, the bimodal distribution of the eigenvalues gives M = 3g/2.The value of g for diffusive waves is close to the number of open eigenchannels with transmission eigenvalues greater than 1/e [6,122].For deeply localized waves, the number of open eigenchannels falls to zero, but M approaches unity since the transmittance is dominated by a single transmission eigenvalue.These calculations apply to a single large TM and to subsets of samples with the same values of M −1 .We consider measurements of samples with small value of N and group together subsets of samples with similar values of M −1 .Figure 6 demonstrates that var(NT a /T ) is close to M −1 in subsets of samples with given M −1 .
The measurements of microwave transmission described in this review were carried out in collections of samples consisting of randomly positioned alumina spheres with refractive index 3.14 within a copper tube.The dielectric spheres were embedded in Styrofoam shells to produce an alumina volume fraction of 0.068.Microwave radiation was launched and detected by wire antennas connected to a vector network analyzer.The orientation of the wire antennas determines the polarization of the field that is detected.The source and detector antennas are mounted on 2D translation stages and moved to positions on a square grid over the open ends of the waveguide at which spectra of the in-and out-of-phase components of the transmitted electric field t ba are collected.Measurements are made over two frequency ranges of 10-10.24GHz and 14.7-14.94GHz.The lower frequency range is just above the first Mie resonance of the individual alumina spheres and so the scattering is strong and the wave is localized in a relatively short length of the random waveguide.The wave is diffusive in the higher frequency range.The number of allowed propagating waveguide modes increases with frequency; they are ∼30 and ∼66 in the centers of the two frequency ranges.Measurements of spectra of the TM in the two frequency ranges are made for sample lengths between 23 and 102 cm so that ensembles over a wide range of g are studied.The impact of absorption on the statistics of transmission is largely removed by Fourier transforming the field spectrum into the time domain and multiplying the time signal by exp(t/2τ a ), where t is the time delay and τ a is the absorption time [46,48].
Measured distributions of relative total transmission and intensity for subsets of samples with M −1 in the ranges 0.17 ± 0.01 and 0.995 ± 0.005 are shown in Fig. 7.The distributions in samples with the same values of M −1 but drawn from ensembles of different g are seen to coincide.This shows that the distributions P(NT a /T ) depend only upon the value of M. P(NT a /T ) are also presented in Fig. 7(b) for M −1 near unity in measurements for two ensembles of localized samples with g = 0.37 and 0.17.When τ 1 τ 2 , transmission is dominated by the first transmission eigenchannel so that, This gives NT a /T = |v 1a | 2 .The probability distribution P(NT a /T ) is therefore expected to be a negative exponential, P(NT a /T ) = exp(−NT a /T ).Measurements are in excellent agreement with this conjecture, as shown in Fig. 7(b).
The eigenchannel participation number M can be linked directly to the degree of intensity correlation within a single sample.In the limit N 1, infinite-range correlation due to fluctuations in T vanishes if computed for the TM of a single sample.The cumulant correlation function of transmitted intensity relative to its average, T /N 2 , is then M .The values of residual short-range and long-range correlation determined from the variances of relative intensity and total transmission are unity and M −1 , respectively, giving    The probability distributions of normalized total transmission P(NT a /T ) for samples with M −1 = 0.17 ± 0.01 from two random ensembles with g = 3.9 (Green open circles) and g = 0.17 (red filled circles) are seen to overlap, confirming that P(NT a /T ) depends only on the value of M. The solid curve is a calculation of P(T a / T a ) with var(T a / T a ) replaced with M −1 in the expression from Ref. [59,60].(b) P(NT a /T ) for M −1 in the range of for two ensembles with g = 0.37 and 0.17 This leads to var(N 2 T ba /T ) = 1 + 2M −1 .This is confirmed in measurements by the correspondence of the measured variances of relative intensity for the two values of M −1 of 0.17 and 0.995 of 1.38 and 3.04 with the calculated values of 1.34 and 3.00.Since the statistics of transmission in a single sample at a given frequency depend only on M, they do not reflect the nature of transport within a disordered material.This would require knowledge of the distribution of M. Previous studies have shown that the crossover to Anderson localization can be tracked in terms of the fluctuation of the intensity and total transmission over collections of samples in an ensemble [48, 58-60, 62, 63].In particular, for diffusive waves in the quasi-1D geometry, the variances of T ba and T a normalized by their ensemble averages are linked to the dimensionless conductance g via the relation, var(s ba = T ba / T ba ) = 1 + 4/3g and var(s a = T a / T a ) = 2/3g, respectively.var(s a ) is also equal to the degree of long range intensity correlation, κ ≡ δ s ba δ s b a ,where δ s ba = s ba − 1.In many circumstances, it is not possible to determine the absolute value of the ensemble average of the conductance, g, and in measurements in some geometries, such as optical measurements of transmission in a slab, the precise correspondence of transmission and conductance is not clear.Nonetheless, it is still possible to determine the degree of localization in a sample from measurements of the variances of normalized intensity or total transmission.Because of infinite-range mesoscopic intensity correlation, the fluctuations do not self-average and remain prominent even for the most spatially averaged transmission quantity, the transmittance T [11-13, 23, 40, 123-127].
Measurements of the TM for classical waves make it possible to explore long-standing predictions regarding the statistics of conductance in mesoscopic systems in the crossover to localization from the perspective of the transmission eigenvalues [23,[128][129][130][131][132][133][134][135].Unlike measurements of the electronic conductance, for which the measurement of current gives the flux of electrons directly, the transmittance of classical waves is measured on a grid in real space.In such measurements, it is difficult to measure transmission for all incident and outgoing channels.Goetschy and Stone [136] have considered the impact of the incomplete measurements of the TM in momentum space upon the statistics of the transmission eigenvalues.They find that the distribution of the density of transmission eigenvalues changes progressively from the bimodal distribution for diffusive waves to a distribution characteristic of Gaussian random matrices [137] as the ratio of measured channels N to the total number of channels N on the input and output sides, m 1 = N 1 /N and m 2 = N 2 /N, decreases.We find the impact of incomplete measurements upon the statistics of eigenvalues is low as long as N > g [30,136,138].Thus, the TM for localized waves determined in a single polarization in real space can faithfully represent the statistics of the transmission eigenvalues and the transmittance.T ) for g = 0.37 (red dots) and 0.045 (green asterisk).The solid black line is a Gaussian fit to the data.For g = 0.045, all of the data points are included, whereas for g = 0.37, only data to the left of the peak is used in the fit.(b) P(T ) for g = 0.37 in a semi-log plot exhibits an exponential tail.(From Ref. [30]) In Fig. 8, we present the probability distributions of ln T for microwave radiation transmitted through random ensembles of dielectric samples with g = 0.37 and 0.045.Since the measured distributions of the normalized transmittance s = T / T depends only upon g for localized samples, the value of g for the localized samples can be obtained by comparing the probability distribution P(s) with numerical simulations for a given value of g.P(ln T ) is seen in Fig. 8(a) to be close to a Gaussian distribution for g = 0.045.For g = 0.37, the low-transmission side of P(ln T ) is seen to be well fit by a Gaussian distribution, while the high transmission side falls sharply above the peak at ln T = −0.5.Above T = 1.1, P(T ) is seen in Fig. 8(b) to fall exponentially, in accord with predictions in [134].We show below that in addition to explaining the origin of the distributions of transmittance in the diffusive and localized limits, the Coulomb gas model can explain the anomalous distribution of transmission in the crossover between diffusive and localized transport [26,31].Understanding transport in terms of the charge model is of further interest because the charge model provides a parallel treatment of diffusive and localized waves and so forms the basis for a universal description of wave transport.The variation of the average positions of the charges x n associated with τ n for different po-sitions of the first charge x 1 in the random ensemble with g= 0.37 is plotted in Fig. 9.The repulsion between x 1 , which is associated with the highest transmission eigenvalue τ 1 , and its image at -x 1 enforces a ceiling for τ 1 of unity.The average spacing between x 1 and x 2 increases as the value of x 1 decreases since the charge at x 2 then interacts strongly with the charge at x 1 as well as with its nearby image at −x 1 .At the same time, the spacing between x 2 and x 3 and their average positions hardly change.This reflects the tendency to heal large fluctuations in charge positions for more remote charges.The source of the sharp cutoff in P(ln T ) can be seen in the context of the Coulomb gas model by examining the spacings of charges for different values of ln T , as shown in Fig. 9(b).A relatively high value of T is only achieved when the first charge is near the origin.This is an unlikely circumstance because this charge is strongly repelled by its image.P(ln T ) would be expected to fall off especially rapidly for values of T above unity since this would requires two charges along with their images to be close together near the origin as seen in Fig. 9(b).
The probability distribution of ln T for the deeply localized waves is seen to follow a Gaussian distribution.We further explore the relation between the first two moments of the distribution.The ratio of σ 2 = var(ln T ) and − ln T , R, vs. L/ξ is plotted in Fig. 10 and seen to approach unity for L ξ ∼ 24 cm.This occurs just as the weighted average of (a) Fig. 10.Approach to single-parameter scaling in multi-channel random systems.R ≡ −var(ln T )/ ln T , MT / T and M −1 T / T are plotted vs. L/ξ .The dashed line is the prediction of SPS in the limit of deeply localized samples in 1D disordered materials.(From Ref. [30]) and M −1 T / T approach unity.The latter average closely tracks R. Thus, SPS is approached in a multichannel quasi-1D random system when transmission is through a single eigenchannel and transport is essentially one dimensional [23,30,139,140].In this limit, the value of M is close to one in nearly every realization of the random sample except for sample configurations with especially small values of T, and the relationship between the statistics of transmittance, total transmission and intensity is particularly straightforward.We can therefore obtain the distributions of total transmission and intensity from the distribution of the transmittance by setting M = 1 and P(T, M) = P(T ).Because the distribution of normalized intensity for a given value of T a is a negative exponential, fluctuations of normalized intensity are closely linked to fluctuations of normalized total transmission, s n ba = n! s n a .Similarly, in the limit of M = 1, s n a = n! s n , where s = T / T .This gives the expression for the variance of s ba when transport is dominated by a single transmission eigenchannel, var(s ba ) = 2var(s a ) + 1 = 4var(s) + 3 [141].We find the values of the variances of s ba , s a and s to be 27.5, 13.2 and 6.4 in the sample of L = 102 cm, which is consistent with these relations.This further confirms that the transport in this random ensemble is via a single eigenchannel.

Controlling wave propagation in a scattering medium
Since a random speckle pattern develops on the output plane of a multiply scattering sample illuminated with monochromatic light, it appears that the information regarding incident wave is completely washed out [1].Nonetheless, because the transmitted field is linked to the incident field via the TM of the sample E b = ∑ N a t ba E a , it is possible to solve the inverse problem and to retrieve the incident waveform by recording the TM and the transmitted field [4,6,28,34].In particular, measuring the TM provides the specific waveform to create a focal spot [3,88,98,100,101].
Transmitted waves can be focused at a target point β by applying an incident waveform which is the phase conjugation of the field transmission coefficients between the target point β and the incident points a, t * β a [92].A peak in the intensity emerges at β , since the transmitted fields from different incident points a arrive at β in phase and interfere constructively while these fields are randomly phased at other points.The incident field is normalized by ∑ N a |t β a | 2 , so that the incident power is unity.The intensity at focal point β is equal to ∑ N a |t β a | 2 , which is the total transmission T β from the point β to the opposite surface.This yields a factor of N enhancement over the ensemble average transmission T ba .The background intensity is also increased because of the long range intensity correlation within the TM so that the contrast in focusing is smaller than N.The intensity at the focus and in the background can be written as, , respectively [119].The contrast for a single measurement of the TM is defined as, µ ≡ I β / I b , in which . . .b,β indicates the average over the background channels b for a given focusing point β .In the limit of N 1, we find, The open circles and squares represent measurements from transmission matrices N = 30 and 66 channels, respectively.The filled triangles give results for N × N matrices with N = 30 for points selected from a larger matrix of size N = 66.Phase conjugation is applied within the reduced matrix to achieve maximal focusing.Eq. 4 is represented by the solid red and dashed blue curves for N = 30 and 66, respectively.(From Ref. [119]) The results of measurements of the contrast in optimal focusing for diffusive samples of length L = 23 cm for incident waves in two orthogonal polarization are shown in Fig. 11 to be in excellent agreement with Eq. 4. When the number of measured points N is smaller than N and therefore the corresponding eigenvalue participation number M is smaller than M for the complete matrix, the contrast is given by Eq. 4 with the substitutions M → M and N → N .This is demonstrated in Fig. 11 with the contrast in samples with N = 66 but with the contrast computed only for N = 30 points falling on the curve for N = 30.These results may be applied to optical measurements of the transmission matrix, where the full matrix is not determined.
In contrast to a conventional focusing lens, in which the size of the focused beam is diffraction limited and so equal to λ /2 divided by the numerical aperture of the lens, the profile of the focused beam through a disordered medium depends only upon the property of the random system itself [56,94,101].The average intensity at a point b at a distance ∆r from the focal point I f oc (∆r) normalized by the peak intensity for a diffusive sample can be expressed in terms of degree of the long range intensity correlation κ and the square of the amplitude of the field correlation function, F(∆r) = |E(r)E * (r + ∆r)|/ I(r) I(r + ∆r) , The ratio between the background intensity F(∆r λ ) and the focus intensity F(∆r = 0) is equal to κ/(1 + κ), since F(∆r) → 0 for ∆r λ .In the diffusive limit, κ ∼ 1/M 1, this gives a contrast of µ = M and the focused profile is reduced to F(∆r).Since the intensity correlation function in the diffusive limit is F(∆r) according to the field factorization approximation, the profile is then expected to be the same as the intensity correlation, which has been demonstrated in optical measurements [94].The transmitted field for a single polarization is measured along a line with a spacing of 2 mm and the source antenna is translated along 49 points on a square grid with 9 mm spacing.F(∆r) (blue dots) is fit with the theoretical expression obtained from the Fourier transform of the specific intensity (dashed blue line).The black dashed line is proportional to T ba / I f oc (∆r = 0) = 1/N.Equation 5 is not valid for localized waves, but good agreement is obtained when κ is replaced by 1/(µ −1) in Eq. 5 with µ determined experimentally.(From Ref. [101]) The normalized profile of the focusing beam via phase conjugation averaged over all possible focusing spot and over all frequency points is presented in Fig. 12.The profile of the focused beam is in excellent agreement with Eq. 5 for diffusive samples.The black dashed line is 1/N, which is the ratio between the ensemble averaged intensity and the intensity at the focus.The background intensity in optimal focusing is seen to be much greater than the average transmission over a random ensemble and is enhanced by a factor of N/M.For localized waves, only a single transmission eigenchannel contributes so that the transmitted wave cannot be focused by applying phase conjugation.
It is also possible to focus a wave inside a random medium.Since the field inside a material is in general not accessible, many virtual guide stars have been created inside the sample in order to facilitate the focusing at a desired location.These include the use of feedback from florescent particles embedded in the sample [95], photoacoustic [106] and acousto-optical effects [97,104] and the use of many nonlinear phenomena, such second harmonic generation.The contrast of focusing inside a disordered sample has been explored in simulations of a scalar wave propagating through Q1D samples which are locally two dimensional.The wave equation 2 E(x, y) + k 2 0 ε(x, y)E(x, y) = 0 is discretized on a square grid and solved using the recursive Greens function method [142].Here, k 0 is the wave number in the vacuum and ε(x, y) is the spatial varying dielectric function.ε(x, y) is set to be unity in the two ideal leads and is equal to 1 ± δ ε in the disordered sample.δ ε is drawn from a rectangle distribution whose width represents the strength of the disorder.Both internal and external reflection are negligible since the sample is index matched to the leads.The product of the wave number in the leads k 0 and the grid spacing is unity and the dimension of the samples are measured in units of the grid spacing.The Green's function G(r, r ) is calculated between grid points r = (0, y) on the incident plane and r = (x , y ) at depth x within the interior of the sample with y the transverse coordinate.We considered an ensemble of diffusive samples with length L = 800 and width L t = 500.The fluctuation of dielectric function at each site is drawn randomly from a rectangular distribution from 0.8 to 1.2.This gives g = 16.
Optimal focusing at a location at (x 0 , y 0 ) inside the sample can be obtained by phase conjugating the Green's function between arrays of points in the incident plane and the target spot (x 0 , y 0 ).In Fig. 13(a), we present the intensity profile of the focused wave at the center of the transverse dimension L t /2 of the waveguide at two depths x = L/4 and x = L/2.In both cases, intensity peak emerges at the focal spot with a contrast of ∼ 37 and ∼ 17, respectively.The lower value of contrast for focusing at large depth into the scattering sample is related to the increasing mesoscopic intensity correlation towards the output surface of the sample.The focusing in the sample is given by Eq. 5 with κ replaced κ(x), the degree of long range intensity correlation at depth x.We expect that the contrast of focusing inside a sample is given in terms of the participation number of eigenvalues of the energy density at a depth x into the sample M(x) of the field matrix connecting the field at the incident plane and the focal plane at depth x.This hypothesis is confirmed in the plot of spatial variation of the focusing contrast µ(x) and 1/(1/M(x) − 1/N) in Fig. 13(b).Good agreement is seen once the depth into the sample is beyond the transport mean free path so that the field is completely randomized.
When a short pulse is incident on a scattering medium, the pulse is stretched in time as it is scrambled in space.Focusing of a pulse through a random medium was first achieved via the time reversal in acoustics [113].Recently, focusing of a pulse in transmission propagating through an opaque sample has been achieved by adjusting the incident wavefront to enhance transmission at a focal spot at a selected time delay [143][144][145][146].The temporal response of an open random medium to an incident pulse can be naturally described in terms of the quasi-normal-modes of the system [69,147,148].These modes are characterized by their central frequencies ω m and linewidths Γ m , and a Lorentzian shape in frequency, Γ m /2/(Γ m /2 + i(ω − ω m )).Once excited, each of the modes will release its energy at a rate determined by its linewidth Γ m , exp(−Γ m t).Spectra of the transmission coefficient between incident channel a and outgoing channel b can be expressed as a summation of the modes, Here, the matrix t m ba describes the spatial coupling into the m th mode via incident point a and coupling out to the output point b from the m th mode.We emphasize that, in contrast to the eigenchannels, the elements t m ba of two different modes can be strongly correlated [69].The Fourier transform of the spectra of the TM yields a time-dependent TM, which characterizes the spatial-temporal response of the random medium to an incoming pulse, t(t ) = ∑ m t m ba exp(−iω m t ) exp(−Γ m t /2), where t indicates the time delay.To have a pulse emerge in transmission at an output point β at a certain time delay t , one simply has to apply phase conjugation of the TM at time delay t , t * β a (t ), so that different pulses from incident channels a will be in phase at target point β at time delay t .
The spatial-temporal control of pulse transmission through a scattering medium via phase conjugation of the time-dependent TM is presented in Fig. 14.The intensity that would be delivered to a point at the center of the output surface of the waveguide β =(0,0) at a selected time t if the transmission matrix were phase conjugated at time t = 33 ns is presented in Figs.14b.A sharp pulse emerges at the selected time delay with intensity peaked at β =(0,0).The temporal profile of the focused pulse is the square of the field correlation function in time, For an incident Gaussian pulse, the square of the field correlation function is equal to the intensity profile of the incident pulse and is independent of delay time.
The spatial contrast of the focusing of the pulse is given by where M is the eigenchannel participation number of the TM with a size of N at time delay t .The time evolution of M and µ are shown in Fig. 15.Near the arrival time of the ballistic wave, the value of M is close to unity and the contrast is not described by Eq. 7.This is because ballistic wave is associated with the propagating waveguides modes with the highest speed and therefore the transmitted field is not randomized.Once the transmitted wave at the output is multiply scattered, a random speckle pattern develops and the measured contrast is in accord with Eq. 7.After the arrival of ballistic waves, the value of M (t ) is seen in Fig. 15 to increase rapidly before falling slowly.This reflects the distribution of lifetimes and the degree of correlation in the speckle patterns of quasi-normal modes [67,68].Just after the ballistic pulse, transmission is dominated by the shortest-lived modes [148][149][150].These modes are especially short lived and strongly transmitting because they are extended across the sample as a result of coupling between resonant centers.Sets of extended modes that are close in frequency could be expected to have similar speckle patterns in transmission, so that a number of such modes might then contribute to a single transmission channel [69].As a result, the number of independent eigenchannels of the transmission matrix contributing substantially to transmission would be relatively small at early times and M would be low.At late times, only the long-lived modes contribute appreciably to the transmission [151] and so the value of M will be close to 1. Thus for intermediate times, modes with wider distribution of lifetimes than at either early or late times contribute to transmission and these modes are less strongly correlated than at early times so that M and the contrast are peaked.

Dwell time and densities of states of the transmission eigenchannels
Spectra measurements of the TM allow us to explore the photon dwell time of the transmission eigenchannels.The delay time of transmission in 1D disordered samples is given by I(t)tdt/ I(t)dt [152].In the limit of vanishing bandwidth, it is equal to the derivative of the phase accumulated as the wave propagates through the sample with respect to the angular frequency ω, φ = dφ /dω [152].Because the delay time in transmission equals the delay time in reflection, φ is also the DOS divided by π at angular frequency ω, ρ(ω).For a multi-channel system, φ ba is still the delay time between input channel a and output channel b, but it is not  15.Time evolution of the maximal focusing contrast µ .µ is well described by Eq. 7 after the time of the ballistic arrival, t ∼ 21 ns.At early times, the signal-to-noise ratio is too low to analyze the transmission matrix.(From Ref. [146]) loner the DOS.This can be seen by noting that a phase shift of order unity occurs within the field correlation frequency, in which typically a single peak occurs in transmission, but approximately δ modes contribute to this single peak.Thus the number of modes is much greater than the phase shift.Krein, Birman, Yafaev and Schwinger [153][154][155] have shown that the DOS can be expressed in terms of the scattering matrix S of the system, ρ(ω) = −i/(2π)Tr(S † dS/dω), in which −iS † dS/dω is known as the Wigner-Smith delay-time matrix [156,157].Measurements of the DOS should therefore involve measurements of the complete scattering matrix.For a system with time-reversal symmetry, Brandbyge and Tsukada [158] demonstrated that the local DOS in an electronic system can be found from the composite derivative of the phase of the transmitted field with shift in the potential at a point inside the sample.We found that [115], Here, dθ n /dω = −i(u * n du n /dω − v * n dv n /dω) is eigenchannel density of states (EDOS), which is the contribution to the DOS of the n th eigenchannel and is the dwell time of the n th eigenchannel.dθ n /dω is thus a complementary set of parameters to τ n .Thus spectra of the TM yield the dynamics of transmission as well as static transmission for each eigenchannel.
Spectra of τ n and dθ n /dω for a diffusive sample with g = 6.9 are presented in Fig. 16.The average value of dθ n /dω is seen to decrease as the eigenchannel channel index n increases, indicating that the average dwell time becomes shorter as the value of transmission decreases.The sum of dθ n /dω represents the DOS of the sample, when the impact of incompleteness of the TM is low.We therefore expect that Eq. 8 gives the DOS in samples for which N M, as is readily achieved in samples in which the wave is localized.The DOS of the sample can also be found by identifying and counting the contribution of the each quasi-normal mode inside the sample.For an open system, the DOS can be expressed as, The central frequencies and linewidths of the modes have been found for localized samples where the modal overlap is small, by decomposing the transmitted field speckle pattern into a  16.Spectra of transmission eigenvalues τ n and the eigenchannel dwell times dθ n /dω for n = 1, 5, 15, 25 for a diffusive sample with g = 6.9.The dwell time for the eigenchannels is the contribution of the eigenchannel to the DOS.(From Ref. [115]) sum of the same set of modes [69].This allows us to directly measure the DOS of the disordered sample.Good agreement based upon these two approaches for the localized sample with g = 0.37 is seen and demonstrated in Fig. 17  Since the photon dwell time and the EDOS of an eigenchannel are proportional to the energy stored inside the sample [159], this suggests that the energy stored inside is smaller for low transmitting eigenchannels.This is consistent with the results found simulations by Choi et.al., [114] and measurements by Gérardin et.al. [91] in single samples.An expression of the profile of an eigenchannel within a random medium would be of importance, since it would afford a universal description for wave propagation that encompasses the wave inside the sample as well as the transmitted wave.This would also provide preliminary sense of the intensity inside individual eigenchannels channels for a host of application involving energy deposition inside random media.This may also be exploited to lower the lasing threshold of a diffusive random laser [160,161], in which the threshold is particularly high because the pump energy cannot penetrate deep into the sample [162].The profiles of eigenchannels over a wide range of transmission obtained in 2D simulations using the recursive Green's function method are shown in Fig. 18.The profiles are averaged over a collection of samples with the same transmission τ.The profile of the complete transmission eigenchannel with τ = 1 is found to be closely related to the probability of a wave returning to a cross-section at a depth x in an open disordered medium [70,[116][117][118].We have found that the profile of an eigenchannel with τ < 1 can be expressed as a product of the profile of the fully transmitting eigenchannel and a function governed only by the auxiliary localization length ξ , which was previously used to parameterize the corresponding transmission eigenvalues via τ = 1/ cosh 2 (L/ξ ).These results not only give a physical meaning to the set of ξ but also enable control of the energy deposition inside the opaque media.

Conclusions
In this review, we have discussed the power of the transmission matrix to describe the statistics of transmission and control of waves propagating through disordered systems.We have shown that the joint probability distribution P(T, M) determines the statistics of transmission over a random ensemble up to second order, while the statistics within a single sample of given T depends only on M. The statistics of transmission in multi-channel disordered systems approach predictions for 1D, when the transmission is dominated by a single transmission eigenchannel and so M → 1.In this limit, the SPS hypothesis for scaling in random 1D system holds for multi-channel disordered samples.The contrast of optimal focusing through and within a scattering sample is equal to M of the measured TM, provided that the number of the measured channels of the TM is much greater than g.The spatial profile of the focused wave is given in terms of the degree of long-range intensity correlation and the spatial field correlation function while the temporal profile of a focused pulse is equal to the square of the field-field correlation function in time.In addition to the transmission eigenvalues, which describe static transmission, there exists a complementary set of parameters, which gives the dwell time and the integrated intensity inside the sample for transmission eigenchannels, and the contribution of each eigenchannel to the DOS.The TM approach allows us to determine the DOS of a disordered material without finding each of the QNMs in the system, which can be difficult for diffusive samples because of the modal overlap.The ability to deposit energy deep within the sample via high-transmission eigenchannels may find applications in deep tissue imaging, depth profiling, photodynamic therapy and random lasing.

Fig. 3 .
Fig.3.Spectra of the transmission coefficients normalized to ensemble average values.In-and out-of-phase field transmission coefficients, field E ba / T ba / , transmission coefficients of intensity s ba = T ba / T ba , total transmission s a = T a / T a and transmittance s = T / T for microwave radiation propagating through a random waveguides for a diffusive sample (left column) and a localized sample (right column).For diffusive waves, δ > 1, many modes contribute to transmission at all frequencies and for all source and detector positions.In contrast, for localized waves δ < 1, distinct peaks appear when the incident radiation is on resonance with a quasi-normal mode.The resonance condition holds for all source and detector positions and therefore sharp peaks remain even when transmission is integrated over space.

Fig. 4 .
Fig. 4. Scaling of optical transmission through a wedged sample consisting of TiO 2 embedded in polystyrene.(a,b) Log-log and semi-log plots of the scaling of the transmission coefficient, T (L).The inverse variation of T (L) with L is seen in the log-log plot and the exponentially falloff of transmission beyond the absorption length L a is seen in the semilog plot.The absorption length L a = 112 ± 5µm of the sample is indicated by the dashed line in both figures.The solid line is a fit of the expression T (L) = 5αD/v sinh(αL) to the data.Here D is the diffusion coefficient, v is the speed of light in the medium and α is the inverse of L a .(c) Spectra of transmission for various sample lengths.The measured transmission in these experiments never reaches zero because the polarized intensity only vanishes at point singularities while measurements are made with a finite aperture placed before the photomultiplier tube.(From Ref. [41])

Fig. 5 .
Fig. 5. Absolute value of the wave field in the scattering medium at f = 0.36 MHz associated with (a), an incident plane wave, (b), a closed eigenchannel, (c), an open eigenchannel deduced from the measured S matrix, and (d), an open eigenchannel deduced from the normalized matrix Ŝ.The corresponding intensities averaged over the wave guide section (y axis) are shown versus depth x in lower panels (e)(h).They are all normalized by the intensity at the plane of sources (x = 0).(From Ref. [91]) is a negative exponential function, P(NT ba /T a ) = exp(−NT ba /T a )[120,121].Thus, the statistics of the transmitted flux relative to its average in a single TM are fully specified by the statistics of total transmission T a normalized by its mean within the TM, T a / T a a = NT a /T .Here, . . .a indicates the average over all incident channels a in a single sample.

Fig. 6 .
Fig.6.The variance of the normalized total transmission vs. the inverse of the eigenchannel participation number M is obtained from measurements by grouping single samples with the same value of M −1 from ensembles with a wide range of values of g.The solid line is var(NT a /T ) = M −1 .(From Ref.[119])

Fig. 7 .
Fig.7.Probability distributions of normalized total transmission and intensity in subensembles of random samples with similar values of M −1 .(a) The probability distributions of normalized total transmission P(NT a /T ) for samples with M −1 = 0.17 ± 0.01 from two random ensembles with g = 3.9 (Green open circles) and g = 0.17 (red filled circles) are seen to overlap, confirming that P(NT a /T ) depends only on the value of M. The solid curve is a calculation of P(T a / T a ) with var(T a / T a ) replaced with M −1 in the expression from Ref.[59,60].(b) P(NT a /T ) for M −1 in the range of for two ensembles with g = 0.37 and 0.17.The solid line is an exponential distribution, exp(−NT a /T ).(c,d) The normalized intensity distributions, P(N 2 T ba /T ) for the same as in (a,b).The solid lines are the calculations of the distribution based upon P(NT a /T ) in (a,b) and the exponential falloff of the distribution of the normalized intensity P(NT ba /T a ) = exp(−NT ba /T a ).(From Ref. [119]) Fig.7.Probability distributions of normalized total transmission and intensity in subensembles of random samples with similar values of M −1 .(a) The probability distributions of normalized total transmission P(NT a /T ) for samples with M −1 = 0.17 ± 0.01 from two random ensembles with g = 3.9 (Green open circles) and g = 0.17 (red filled circles) are seen to overlap, confirming that P(NT a /T ) depends only on the value of M. The solid curve is a calculation of P(T a / T a ) with var(T a / T a ) replaced with M −1 in the expression from Ref.[59,60].(b) P(NT a /T ) for M −1 in the range of for two ensembles with g = 0.37 and 0.17.The solid line is an exponential distribution, exp(−NT a /T ).(c,d) The normalized intensity distributions, P(N 2 T ba /T ) for the same as in (a,b).The solid lines are the calculations of the distribution based upon P(NT a /T ) in (a,b) and the exponential falloff of the distribution of the normalized intensity P(NT ba /T a ) = exp(−NT ba /T a ).(From Ref. [119])

Fig. 8 .
Fig.8.Probability distribution of microwave conductance.(a) P(ln T ) for g = 0.37 (red dots) and 0.045 (green asterisk).The solid black line is a Gaussian fit to the data.For g = 0.045, all of the data points are included, whereas for g = 0.37, only data to the left of the peak is used in the fit.(b) P(T ) for g = 0.37 in a semi-log plot exhibits an exponential tail.(From Ref.[30])

Fig. 9 .
Fig. 9. Coulomb gas model of transmission eigenvalues and conductance.(a) Average positions of charges and their images as the position of the first charge x 1 change in the random ensemble with g = 0.37.The dashed lines show the average positions of the charges for this ensemble.The curly brackets in this case represents the average over a subset of transmission matrices with the specified value of x 1 .(b) Average positions of charges vs. ln T in the same ensemble.Here, the curly brackets indicate the averaging is over a subset of transmission matrices with the specified value of ln T .(From Ref. [30])

66 N = 30 'Fig. 11 .
Fig. 11.Contrast in optimal focusing vs. eigenchannel participation number M. The open circles and squares represent measurements from transmission matrices N = 30 and 66 channels, respectively.The filled triangles give results for N × N matrices with N = 30 for points selected from a larger matrix of size N = 66.Phase conjugation is applied within the reduced matrix to achieve maximal focusing.Eq. 4 is represented by the solid red and dashed blue curves for N = 30 and 66, respectively.(From Ref.[119])

Fig. 12 .
Fig.12.The ensemble average of normalized intensity for focused radiation (blue circles) is compared to Eq. 5 (blue solid line) for L = 61 cm for diffusive and localized waves.The transmitted field for a single polarization is measured along a line with a spacing of 2 mm and the source antenna is translated along 49 points on a square grid with 9 mm spacing.F(∆r) (blue dots) is fit with the theoretical expression obtained from the Fourier transform of the specific intensity (dashed blue line).The black dashed line is proportional to T ba / I f oc (∆r = 0) = 1/N.Equation 5 is not valid for localized waves, but good agreement is obtained when κ is replaced by 1/(µ −1) in Eq. 5 with µ determined experimentally.(From Ref.[101])

Fig. 14 .
Fig. 14.Controlling pulse transmission through a random medium with a time dependent TM.Spectra of the TM for microwave radiation propagating through a random waveguide with a length L = 61 cm are measured from 14.7 to 15.7 GHz in 3200 steps.The TM is determined with a single polarization between pairs of 45 points on the incident and output surfaces.The time dependent transmission matrix is obtained from spectra of the transmitted field between all points a and b, t ba (ν).These spectra are multiplied by a Gaussian pulse centered in the measured spectrum at ν 0 =15.2GHz with bandwidth σ ν =150 MHz and then Fourier transformed into the time domain.This gives the time response at the detector to an incident Gaussian pulse launched by a source antenna with bandwidth σ t =1/2π/σ ν .(a) The time variation I ba (t ) of an incident pulse launched at the center of the input surface and detected at the center of the output surface in a single realization of the random sample and the average of the time-of-flight distribution I(t ) .(b) Phase conjugation is applied numerically to the same configuration as in (a) to focus at t = 33 ns at the center of the output surface.The Whittaker-Shannon sampling theorem is used to obtain high-resolution spatial intensity patterns shown in the inset of (b).(From Ref. [146]) Fig.15.Time evolution of the maximal focusing contrast µ .µ is well described by Eq. 7 after the time of the ballistic arrival, t ∼ 21 ns.At early times, the signal-to-noise ratio is too low to analyze the transmission matrix.(From Ref.[146]) Fig.16.Spectra of transmission eigenvalues τ n and the eigenchannel dwell times dθ n /dω for n = 1, 5, 15, 25 for a diffusive sample with g = 6.9.The dwell time for the eigenchannels is the contribution of the eigenchannel to the DOS.(From Ref.[115]) .

Fig. 17 .
Fig. 17.Comparison of two approaches for finding the DOS of a disordered material.(a) Contributions of the individual modes in Eq. 9 to the DOS.The integral of each mode over the angular frequency is unity.(b) DOS determined from the TM (red curve) by summing spectra of dθ n /dω and modes in (a).(From Ref. [115])

Fig. 18 .
Fig.18.Spatial profile of the energy distribution of transmission eigenchannel inside the sample.(a) Ensemble averages of the eigenchannel energy density profiles W τ (x) for eigenvalues τ = 1, 0.5, 0.1, and 0.001 for a diffusive sample with L/ξ = 0.05.W τ (x) is the energy integrated over the transverse dimension and is normalized to equal τ at the output surface.(b) Ensemble averages of energy density profiles for all transmission eigenchannels with the eigenchannel indices n from 1 to N. The linear falloff of the average of the energy density over all eigenchannels is in accord with the diffusion theory.(From Ref.[118])