# Selectively exciting quasi-normal modes in open disordered systems

Matthieu Davy, Azriel Genack

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

## To cite this version:

Matthieu Davy, Azriel Genack. Selectively exciting quasi-normal modes in open disordered systems. Nature Communications, 2018, 9 (1), pp.4714. 10.1038/s41467-018-07180-3. hal-02134651

HAL Id: hal-02134651

https://hal.science/hal-02134651

Submitted on 12 Jul 2019

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

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

<figure><img src="image_2.png" /></figure>
ARTICLE

DOI: 10.1038/s41467-018-07180-3

OPEN

# Selectively exciting quasi-normal modes in open disordered systems

Matthieu Davy<fnref n="1" /> & Azriel Z. Genack²

Transmission through disordered samples can be controlled by illuminating a sample with waveforms corresponding to the eigenchannels of the transmission matrix (TM). But can the TM be exploited to selectively excite quasi-normal modes and so control the spatial profile and dwell time inside the medium? We show in microwave and numerical studies that spectra of the TM can be analyzed into modal transmission matrices of rank unity. This makes it possible to enhance the energy within a sample by a factor equal to the number of channels. Limits to modal selectivity arise, however, from correlation in the speckle patterns of neighboring modes. In accord with an effective Hamiltonian model, the degree of modal speckle correlation grows with increasing modal spectral overlap and non-orthogonality of the modes of non-Hermitian systems. This is observed when the coupling of a sample to its surroundings increases, as in the crossover from localized to diffusive waves.

<footnote n="1">¹Institut d'Electronique et de Télécommunications de Rennes, University of Rennes 1, 35042 Rennes, France. ²Department of Physics, Queens College and Graduate Center of the City University of New York, Flushing, NY 11367, USA. Correspondence and requests for materials should be addressed to A.Z.G. (email: genack@qc.edu)</footnote>
Transmission of a monochromatic wave through a static sample is fully described by the transmission matrix (TM), t. The TM is a subset of the scattering matrix which provides the coupling of a system to its surroundings. In the last decade, there has been growing interest in measuring the TM to control the flow of waves through random systems and optical fibers¹,². Shaping the incident waveform illuminating a sample makes it possible to manipulate the net transmitted flux and its spatial intensity profile for applications in medical imaging and communications. For instance, a scattering medium can appear to be transparent or opaque when the incoming wavefront is adjusted to correspond to the first or last transmission eigenchannels³⁻⁹. The intensity can also be focused through random media at a selected point in the output by adjusting the incident wave so that all transmission channels interfere constructively at that point¹,¹⁰.

The elements $t_{ba}$ of the TM are the field transmission coefficients between the N channels leading toward and away from opposite ends of a sample, a and b, respectively. The TM was initially studied in order to explain the scaling of the conductance of wires at zero temperature¹¹,¹². Classical and quantum transport are connected by the dimensionless conductance, which is the conductance in units of the quantum of conductance, ($e^2/h$). The dimensionless conductance is equal to the average transmission, $g = \langle T \rangle$, where $\langle \cdots \rangle$ represents the average over a random ensemble. The crossover to Anderson localization occurs at $g=1$; waves are localized for $g<1$ and diffusive for $1<g<N/2$. The transmittance is the sum of all flux transmission coefficients, $|t_{ba}|^2$, which equals the sum of the N eigenvalues $\tau_n$, $T = \sum_{a,b=1}^N |t_{ba}|^2 = \text{tr}(\mathbf{t}^\dagger\mathbf{t}) = \sum_{n=1}^N \tau_n^{11,12}$.

In principle, the degree of control over transmission in diffusive samples is strong because the distribution of transmission eigenvalues is wide. This distribution is bimodal with a peak near unity containing g "open" channels and a second peak corresponding to "closed" channels with values that are exponentially small in the ratio of the sample length and the transport mean free path, $L/(6^{12-14})$. In practice, however, measurements of the TM are incomplete so that the dynamic range over which transmission can be controlled is limited⁴,¹⁵⁻¹⁷. Because the dwell time and the energy density profile inside a medium excited in a specific eigenchannel are correlated with the corresponding transmission eigenvalue, exciting transmission eigenchannels also provides a measure of control over the dwell time and the spatial distribution of energy within a random medium³,⁵⁻⁹,¹⁸. The full diversity of dwell times is given by the eigenvalues of the Wigner-Smith time-delay matrix, known as the proper delay times, which are constructed from the spectrum of the scattering matrix¹⁹⁻²².

Another approach to controlling propagation within random or structured media might be to manipulate the incident wave to preferentially excite specific quasi-normal modes²³ which have different lifetimes and spatial profiles. Modes of open systems are solutions of the wave equation over the volume of the random medium with outgoing radiation boundary conditions²⁴⁻²⁶. In resonating structures for which the complex eigenvalues and eigenvectors can be found analytically or numerically, the field for any source excitation can be reconstructed from the coherent superposition of modal contributions²⁶,²⁷. Beyond the independent contribution of each mode, the resultant field depends critically upon the interference between the fields of modes that overlap spectrally and spatially. Modal coupling plays a key role in describing the physics of photonic systems such as chaotic cavities²⁸⁻³², coupled cavities or waveguides³³,³⁴, optical resonators²⁷,³⁵, quantum plasmonic³⁶,³⁷ or disordered media³⁰,³⁸⁻⁴⁰.

In large complex systems, it is generally not possible to solve for the eigenvectors of the wave equation, but important

properties of a system and its coupling to its surroundings can be determined from the statistics of scattering spectra and their analysis into modes or energy levels. Great emphasis has been placed on the probability distributions of level spacings⁴¹⁻⁴³ and level widths²⁸,³²,⁴⁴. However, the statistics of level widths and spacings do not directly yield the statistics of scattering because the scattered wave also reflects the interference between modes and the degree to which modal speckle patterns are correlated.

Here we consider the degree of modal selectivity that can be achieved in random media by manipulating the incident waveform. We approach the problem by analyzing the spectrum of the TM into its modal components in locally 2D N-channel samples. The complex modal frequencies and amplitudes are found by decomposing the spectra of the elements of the TM into a superposition of spectral lines via Breit-Wigner theory³⁰,⁴⁵,⁴⁶

$$t_{ba}(\omega) = \sum_n \frac{t_{ba}^n}{\omega - \omega_n + i\Gamma_n/2} = \sum_n t_{ba}^n \varphi_n(\omega). \quad (1)$$

Here $\varphi_n(\omega) = (\omega - \omega_n + i\Gamma_n/2)^{-1}$ is the frequency variation of excitation of the field associated with the mode with central frequency $\omega_n$ and linewidth $\Gamma_n$, and $t_{ba}^n$ is the complex field transmission coefficient associated with the nth resonance. Each resonance is then associated with a modal transmission matrix (MTM), $t_n$, which is built upon the coefficients $t_{ba}^n$, and is the contribution of a mode of the scattering medium to the TM²³. An MTM therefore provides the incoming wavefront that maximally enhances the energy in a specific mode. However, modal selectivity becomes more challenging as the degree of modal overlap increases in non-Hermitian media.

## Results

**Modal decomposition in the effective Hamiltonian formalism.** The coupling of a system to its surroundings can be analyzed in terms of an effective Hamiltonian. The coupling is described via the $2N \times 2N$ scattering matrix $S$ expressed in terms of the $M \times M$ effective Hamiltonian $\mathbf{H}_{\text{eff}}$ as²⁸⁻³⁰,⁴⁷⁻⁴⁹

$$\mathbf{S} = \mathbf{1} - i\mathbf{V}^T \frac{1}{\omega - \mathbf{H}_{\text{eff}}} \mathbf{V}. \quad (2)$$

Here $V$ is a real $M \times 2N$ matrix describing the coupling of the $M$ modes of the closed system to the exterior via the $2N$ channels in the leads on both sides of the sample. The non-Hermitian effective Hamiltonian is

$$\mathbf{H}_{\text{eff}} = \mathbf{H}_0 - \frac{i}{2} \mathbf{V} \mathbf{V}^{\text{T}}, \quad (3)$$

where $\mathbf{H}_0$ is the Hermitian Hamiltonian of the closed system. The poles of the $S$ matrix occur at the complex eigenvalues $\tilde{\omega}_n$ of $\mathbf{H}_{\text{eff}}$, $\tilde{\omega}_n = \omega_n - i\Gamma_n/2$. Two sets of eigenvectors are associated with these eigenvalues. These are the right $|\phi_n\rangle$ and left $\langle\phi_n|$ eigenvectors, which are the transpose of one another, $\langle\phi_n| = (|\phi_n\rangle)^T$. The eigenfunctions of $\mathbf{H}_{\text{eff}}$ are bi-orthogonal and satisfy the following orthogonality condition due to the time-reversal symmetry of $\mathbf{H}_{\text{eff}}$, $\langle\phi_n^*|\phi_m\rangle = \delta_{nm}$. As a result of the imaginary part of $\mathbf{H}_{\text{eff}}$, the eigenfunctions are complex.

The modal decomposition of the TM may be expressed in terms of the complex vectors $|\mathbf{W}_{Ln}|^2$ and $|\mathbf{W}_{Rn}|^2$ which couple the eigenstates $\varphi_n$ to the scattering wavefunctions $\xi_L$ and $\xi_R$ in the left
and right leads, respectively²⁹,

$$t(\omega) = -i \sum_{n=1}^{M} \frac{\mathbf{W}_{Rn} \mathbf{W}_{Ln}^{T}}{\omega - \omega_{n} + i \frac{\Gamma_{n}}{2}}. \quad (4)$$

We identify the MTM of the nth mode at resonance as $t_n = -i\mathbf{W}_{Rn}\mathbf{W}_{Ln}^T/(\Gamma_n/2)$. The MTM for each mode is of unit rank since it is the product of the vectors $\mathbf{W}_{Rn}$ and $\mathbf{W}_{Ln}^T$. In principle, the coupling of the eigenfunctions of the closed system to the leads, and therefore vectors $\mathbf{W}_{Rn}$ and $\mathbf{W}_{Ln}$, depend on frequency. However, in the case of resonances with high quality factors $Q_n = 2\omega_n/\Gamma_n$, which are explored in the experiments described below, we can take $\mathbf{W}_{Rn}(\omega) = \mathbf{W}_{Rn}(\omega_n) = \mathbf{W}_{Rn}$ and $\mathbf{W}_{Ln}(\omega) = \mathbf{W}_{Ln}(\omega_n) = \mathbf{W}_{Ln}$. The TM is then expressed as a superposition of MTMs with modal transmission coefficients for Lorentzian lines defined at the resonance frequency. The decomposition of the TM into MTMs and the properties of the MTMs are demonstrated below in microwave experiments.

**Experimental setup.** Measurements of the TM are performed in a two-dimensional cavity containing randomly positioned disks

(see Fig. 1a and Methods for details). Spectra of the $N \times N$ TM are measured between two arrays of $N=8$ emitting and receiving antennas on the left and right sides of the cavity, respectively. Measurements are carried out in the frequency range 10.7–11.7 GHz in a scattering sample of 300 randomly distributed 6-mm-diameter aluminum disks. The sample is weakly localized and modal spectral overlap is moderate. The coupling strength $\tilde{T}_a$ of the antennas to the sample is determined using the mean value of the reflection parameter at each antenna, $\langle S_{aa} \rangle$, $\tilde{T}_a = 1 - |\langle S_{aa} \rangle|^2$ gives $\tilde{T}_a \sim 0.99$ so that the antennas are strongly coupled to the cavity.

**Decomposition into MTMs.** The modes can be found from an analysis of the spectrum of the TM as a superposition of MTMs using Eq. (1). The set $\{\omega_n, \Gamma_n\}$ is extracted via the Harmonic inversion (HI) method from the inverse Fourier transform of spectra of transmission coefficients $t_{ba}(\omega)^{45,50,51}$ (see Methods). The coefficients $t_{ba}^n$ of $t_n$ are then found from the fit of transmission coefficients in the time domain. The flux transmission coefficient between two channels, $|t_{ba}(\omega)|^2$, and the underlying modal transmission coefficients between the channels for each mode, $|t_{ba}^n|^2|\varphi_n(\omega)|^2$ are shown in Fig. 1b. The transmittance and

a

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

b

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

c

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

d

<figure><img src="image_6.png" /><figcaption>**Fig. 1** Experimental setup and decomposition of the TM into modes. **a** Experimental setup with the top plate removed to show the cavity with 300 randomly-positioned aluminum disks. The TM is measured between the 8 antennas on the left and right sides of the cavity. **b** Transmission (blue curve) between two channels, $|t_{ba}(\omega)|^2$, and its reconstruction found using HI (dashed red curve) in the [11-11.5] GHz range. The thin lines are modal strengths in transmission, $|t_{ba}^n|^2|\varphi_n(\omega)|^2$. **c** The measured transmittance, $T(\omega) = \sum_{ab} |t_{ba}(\omega)|^2$, (blue curve) and its reconstruction from the MTMs (dashed red curves), as well as the modal contributions $T_n(\omega)$ to the transmittance (thin lines). **d** The spectrum of transmission eigenvalues $\tau_i(\omega)$ and the first (blue crosses) and second (orange crosses) eigenvalues of the MTM on resonance are shown in a semilog plot. The second eigenvalue of the measured MTM is typically two orders of magnitude smaller than the first</figcaption></figure>
the contribution $T_n(\omega) = \sum_{ba} |t_{ba}^n|^2 |\phi_n(\omega)|^2$ of each mode to $T(\omega)$ are shown in Fig. 1c. The reconstructions of the transmission and of the transmittance from the modes and the measurements of these quantities are in excellent agreement.

The degree of modal overlap may be expressed as the ratio of the mode width and spacing $\delta = \delta\omega/\Delta\omega$, where $\delta\omega = <\Gamma_n>$ is the average linewidth of modes and $\Delta\omega = <\omega_{m+1} - \omega_m>$ is the typical spacing between neighboring modes. In random media, the degree to which modes overlap spectrally tracks the spatial extent of the eigenstates in the interior of the sample and so the crossover from diffusion to localization$^{44,46,52,53}$. When reflection at the interface is weak, the ensemble average of the ratio of the level width to level spacing gives the Thouless number, which is equal to the conductance $g^{44}$. Modes of the medium are generally exponentially peaked within the sample when $\delta < 1$ and extended when $\delta > 1$. Here the average linewidth is $\langle\Gamma_n\rangle \sim 9$ MHz and the degree of modal overlap is $\delta \sim 1.2$. The high average modal quality factor $Q = 3300$ justifies the assumption that the coupling vectors between the quasi-normal modes and the antennas are independent of frequency.

An additional check is placed on the accuracy of the modal decomposition of the TM when the fit to transmission is carried out simultaneously at several points in the sample: the rank of MTMs found in the fits of transmission must be close to unit rank, as predicted by Eq. (4). This additional check is not possible when a single spectrum of a transmission coefficient is decomposed into modes with use of HI, as has been done in studies of chaotic cavities$^{45}$. Spectra of the transmission eigenvalues of the measured TM, $\tau_i(\omega)$, are compared in Fig. 1d to the transmission eigenvalues of modes found from the diagonalization of $t_n t_n^\dagger$. The second modal eigenvalue is typically smaller than the first by a factor of $10^{-2}$. The dominance of the first modal eigenvalue supports the predicted decomposition of the TM into MTMs of unit rank. This is further confirmed in measurements with smaller modal overlap (Supplementary Note 1) in which the ratios of the second and first modal eigenvalue are substantially smaller. The ratio is still smaller in simulations for samples with modal overlap comparable to that in experiments indicating that the quality of the modal decomposition is degraded by noise in the measurements.

We find that when the first and second eigenvalues of the MTM are close in value, the results are likely to be spurious. The modal analysis is limited here to samples with moderate modal overlap. For higher modal overlap, spurious resonances may appear in the modal analysis due to the contributions of modes with large linewidth which cannot be resolved and to modes that lie outside the frequency range but still contribute since their linewidth is broad.

**Modal selectivity.** The strength of excitation of an individual mode is maximized when the incoming wave on the left or right excites the sample with the optimal modal patterns $W_{Ln}^*$ and $W_{Rn}^*$, respectively. These are the complex conjugates or the time-reversal of modal speckle patterns at the sample boundaries. Using Eq. (4), the vector of the transmitted field for an excitation of the sample from the left with the normalized optimal wave-form, $v_n = W_{Ln}^*/W_{Ln}$, can be expressed as

$$ \mathbf{E}_{\text{max}}(\omega) = \mathbf{t}(\omega)v_n = -i \frac{\mathbf{W}_{Ln}^2}{\omega - \omega_n + i \frac{\Gamma_n}{2}} \frac{\mathbf{W}_{Rn}}{\mathbf{W}_{Ln}} -i \sum_{m \neq n} \frac{\mathbf{W}_{Ln}^T \mathbf{W}_{Ln}^* \mathbf{W}_{Rn}}{(\omega - \omega_m + i \frac{\Gamma_m}{2}) \mathbf{W}_{Ln}} \quad (5) $$

The first term in Eq. (5) gives the contribution of the $n^{th}$ mode to transmission for maximal coupling and the sum in the second

term gives the contributions of other modes. Apart from the Lorentzian function, the energy in the mode to which the field is maximally coupled is equal to $||W_{Ln}||^2 ||W_{Rn}||^2$. This can be compared to the average energy for a normalized random incoming waveform $v_{rand}$ which is $\langle |W_{Ln}^T v_{rand}|^2 \rangle ||W_{Ln}||^2 = ||W_{Ln}||^2 ||W_{Rn}||^2/N$. The energy in the mode for maximal coupling in an $N$-channel system is therefore enhanced by a factor $N$ using the optimal modal pattern. This property is a consequence of the unit-rank of the MTMs. At the same time, the contribution to transmission of the selected mode vanishes for any incoming vector orthogonal to the optimal modal pattern $W_{Ln}$. Residual transmission is due to the contributions of neighboring modes.

Excitation of specific quasi-normal modes differs from excitation of transmission eigenchannels. The eigenchannels and eigenvalues of the TM can be found via a singular value decomposition in which the TM at a single frequency is expressed as the product of three $N \times N$ matrices, $t(\omega) = U\Lambda V^\dagger$. Here $\Lambda$ is a diagonal matrix whose elements are the singular values $\sqrt{\tau_i}$, and $V$ and $U$ are unitary matrices and correspond to the waveforms of the transmission eigenchannel on the input and output of the sample, respectively. In contrast to modes, which have a Lorentzian spectrum, the eigenchannels are defined at a specific frequency; a new set of transmission eigenvalues and eigenchannels must be computed at each frequency$^{5,23,54,55}$. However, the spectral characteristics of the channels can be obtained by decomposing the transmission eigenchannels into modes$^{23,54}$. When a single mode dominates transmission, the MTM for this mode is close to the first eigenchannel. At the resonance, the first transmission eigenvalue is, $\tau_1(\omega_n) = ||W_{Ln}||^2 ||W_{Rn}||^2 / (\Gamma_n/2)$. When several resonances overlap, however, the first transmission eigenchannel is a combination of modal contributions of several modes.

**Experimental demonstration of modal selectivity.** To explore the degree of modal selectivity for different incident waveforms, the field coefficient within the medium $e_a(x, y, \omega)$ is measured using a wire antenna inserted through subwavelength holes (see Methods). The contributions of modes inside the medium, $e_a^n(x, y)$, are then obtained from a fit of the coefficients $e_a(x, y, \omega) = \sum_n e_a^n(x, y)\phi_n(\omega)$ using the set of resonances $\{\omega_n, \Gamma_n\}$ obtained from the modal expansion of the TM. The spatial energy distribution for each mode is reconstructed from the contributions of modal field patterns due to the eight incoming channels summed to give the optimum incident modal pattern.

For isolated modes, strong modal discrimination is readily accomplished by tuning to resonance. The strength of excitation is enhanced over the average of random excitation by a factor of $N$ by adjusting the incident wavefront to the optimal modal pattern. In Fig. 2, we consider two weakly overlapping modes at $f_1 = 11.434$ GHz and $f_2 = 11.461$ GHz with linewidths of $\Gamma_1/(2\pi) = 3.75$ MHz and $\Gamma_2/(2\pi) = 4.84$ MHz. This gives a degree of modal overlap between the modes of $\delta_{12} = [\Gamma_n + \Gamma_{n+1}/2] / (\omega_{n+1} - \omega_n)$ of 0.15. The two modes are spatially distinct and peaked at different points within the sample, as seen in Fig. 2a. For maximal coupling to the first and second modes, the transmission is seen in Fig. 2c, d to be enhanced by a factor of close to $N=8$ at the resonance of the two modes in comparison transmission for a random incoming wavefront shown in Fig. 2b. The energy density inside the medium at resonance then closely matches the spatial distribution of the mode, as seen in the insets of Fig. 2c, d. In contrast, for vanishing coupling to the first mode using the third singular vector of the MTM, the first mode does not contribute to transmission and the energy density is due to the contribution of
<figure><img src="image_7.png" /><figcaption>**Fig. 2** Selectivity in modal excitation of weakly overlapping modes. **a** Spatial profiles of the two modes whose resonant frequencies are $f_1 = 11.434$ GHz and $f_2 = 11.461$ GHz are indicated with arrows in **b**. **b-e** Transmission and strengths of modes in transmission (dashed lines) for **b** a random incoming wavefront, **c**, **d** maximal coupling to the first mode (**c**) and second mode (**d**) indicated with arrows. **e** Vanishing coupling to the first mode. The insets in **b-e** are spatial intensity profiles determined from measurements. **f** Modal energy density profiles inside the sample averaged over the cross-section for the first mode is shown for maximal coupling (blue curve) and vanishing coupling (red curve).</figcaption></figure>

weakly overlapping modes, as seen in Fig. 2e. The energy density is concentrated at the beginning of the sample and falls rapidly into sample. Thus, specific modes can be selected using optimal modal incident wave patterns when the modes overlap weakly.

In Fig. 2f, we show the average over the cross-section of the modal strength inside the medium for maximal and vanishing coupling to the first mode at $f_1 = 11.434$ GHz. The contribution of the mode is maximum for the optimal incident modal pattern and should vanish for the orthogonal waveforms. For vanishing coupling, the intensity is seen to fall exponentially within the sample and transmission is more than two orders of magnitude below that for maximal coupling, in agreement with the ratio between the first and second modal eigenvalues.

We next consider selectivity in a case of two strongly overlapping modes (Fig. 3). The modes at $f_1 = 11.763$ GHz and $f_2 = 11.773$ GHz with linewidths $\Gamma_1/(2\pi) = 12.1$ MHz and $\Gamma_2/(2\pi) = 16.6$ MHz

have the modal overlap factor $\delta_{12} = 1.35$. The spatial profiles of the two modes are seen in Fig. 3a to be more extended than the modes discussed previously with $\delta_{12} = 0.15$ and to be very similar. In addition to enhancing the contribution of the maximally excited mode, maximal coupling is seen in Fig. 3c, d to enhance the contribution of the neighboring mode in comparison to a wavefront that has not been optimized. The modal transmission associated with the second mode is enhanced by a factor 4 for maximal coupling to the first mode.

The distributions of energy density for a random wavefront and for the first transmission eigenchannel at a frequency midway between the two resonances $\omega_0 = (\omega_n + \omega_{n+1})/2$, are seen in Fig. 4a, b to be primarily mixtures of the modal spatial patterns of the two neighboring modes shown in Fig. 3$^{54}$. Nevertheless, it is possible to preferentially excite a single mode by adjusting the incident wave to match the pattern of one of the nearby resonant
<figure><img src="image_8.png" /></figure>

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

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

<figure><img src="image_11.png" /><figcaption>**Fig. 3** Selectivity in modal excitation of strongly overlapping modes. **a** Spatial profile of the energy density of two modes at $f_1$ = 11.763 GHz and $f_2$ = 11.773 GHz. **b–d** Transmission (blue curve) and strengths of modes (dashed curves) for **b** a random incoming wavefront, and **c, d** for the wavefront which maximally couple to the modes peaked at the frequencies indicated with arrows</figcaption></figure>

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

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

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

<figure><img src="image_15.png" /><figcaption>**Fig. 4** Selectivity between resonances of overlapping modes. **a–d** Spatial intensity distributions at frequency $(f_1 + f_2)/2$ for **a** a random incoming wavefront, **b** the first transmission eigenchannel, and **c, d** maximal coupling to the first (**c**) and second (**d**) modes</figcaption></figure>

modes. In Fig. 4c, d, the energy density at the frequency between the modal resonances is shown for maximal coupling to one or the other of the modes. In each case, the energy density matches the spatial distribution of the selected mode shown in Fig. 3a. The degree of modal selectivity between two modes achieved by maximizing the input for one of the modes is reduced as a result of the hybridization and spectral broadening of the modes of the

closed system when the sample is coupled to its surroundings. We will see in the theoretical analysis and measurements below that the similarity in modal patterns in the interior of the sample seen in Fig. 3a is a consequence of the bi-orthogonality of the eigenfunctions and the correlation between them.

The degree of modal selectivity can be further enhanced at the expense of the net excitation of the mode by exciting with a
waveform that is orthogonal to the neighboring mode. In the case of two overlapping modes, this is achieved by illuminating the sample with an incident wavefront $aW_{L1}^* + bW_{L2}$ with coefficients $a$ and $b$ satisfying the condition, $aW_{L2}^T W_{L1}^* + bW_{L2}^T W_{L2}^* = 0$. The contribution of the second mode therefore vanishes so that the desired modes is perfectly selected relative to the second mode. Because of the correlation of modal speckle pattern, this also suppresses the excitation of the selected mode.

### Mixing of eigenfunctions

**Nonorthogonality matrix.** Equation (4) shows that the contribution to the output speckle pattern of the mth mode for the incident waveform that couples maximally to the nth mode at resonance, $\omega = \omega_n$, is equal to

$$C_{mn} = \frac{W_{Lm}^{\dagger} W_{Ln}}{\omega_n - \tilde{\omega}_m^*} \qquad (6)$$

The matrix $C$ involves the degree of correlation between the modal patterns at the input. It can be related to the Bell-Steinberger nonorthogonality matrix $U$, which gives the correlation over the volume between the eigenfunctions $\phi_n$ of the system$^{32,56}$. The elements of $U$ given by the scalar product, $U_{mn} = \phi_m^\dagger \phi_n$ can also be expressed in terms of the vectors $W_n$, in the absence of losses that are not due to the coupling of the antennas to the system.

$$U_{mn} = \phi_m^\dagger \phi_n = i \frac{W_m^\dagger W_n}{\tilde{\omega}_n - \tilde{\omega}_m^*} \qquad (7)$$

Here, the vector $W_n$ with 2N elements is the concatenation of the two vectors $W_{Ln}$ and $W_{Rn}$. In Hermitian systems, $U$ is the identity matrix since the modes are orthogonal over the volume. However for non-Hermitian systems, the diagonal elements of $U$, which are the inverses of the phase rigidities of the eigenfunctions, $U_{nn} = 1/\rho_n$, increase with modal overlap$^{32,57,58}$. The phase rigidity, $\rho_n \equiv \langle\phi_n^2\rangle/\langle|\phi_n|\rangle^2$, can be expressed in terms of the degree of complexness of the eigenfunctions of $H_{eff}$ $q_n^2 = \langle\mathrm{Im}(\phi_n)^2\rangle/\langle\mathrm{Re}(\phi_n)^2\rangle$, as $\rho_n = (1 - q_n^2)/(1 + q_n^2)^{39}$. For traveling waves, the real and imaginary parts of the eigenstates are the same on average so that $q_n = 1$ and $\rho_n = 0$. In contrast, for isolated resonances, the eigenfunctions coincide with the real eigenfunctions of the closed system and $\rho_n \to 1$. In random media, the degree of complexness and the phase rigidity track the transition from the diffusive to localized regime$^{59}$.

The completeness of the eigenfunctions implies the sum rule $\Sigma_m U_{nm}^2 = 1^{60}$. The positive diagonal elements of $U$ matrix and the negative off-diagonal elements are enhanced as the degree of modal overlap $\delta$ increases. In the weak coupling regime, the diagonal elements are of order $1 + \delta^{248}$, while the magnitude of off-diagonal elements increase as $-\delta^{58}$.

Equation (7) shows that the non-orthogonality of eigenfunctions over the volume yields a non-vanishing degree of correlation between modal speckle patterns at the interface. In the case of maximal coupling to a mode, the off-diagonal elements are also responsible of the non-vanishing of the contribution of the neighboring modes as shown by the similarity of Eqs. (6) and (7). However, two differences can be observed. First, $C_{nm}$ involves $W_{Lm}^\dagger W_{Ln}$ instead of $W_m^\dagger W_n = W_{Lm}^\dagger W_{Ln} + W_{Rm}^\dagger W_{Rn}$. Since $W_{Ln}$ and $W_{Rn}$ are statistically independent random variables, we may apply the central limit theorem and approximate $W_m^\dagger W_n \sim 2W_{Lm}^\dagger W_{Ln}$ for $N \gg 1$.

Second, the denominator of Eq. (6) depends on $\omega_n - \tilde{\omega}_m^*$ instead of $\tilde{\omega}_n - \tilde{\omega}_m^*$. However, in the case of strongly overlapping resonances, the spacing between the central frequencies $\omega_n - \tilde{\omega}_m^*$ is much smaller than the linewidths $\Gamma_n$ and $\Gamma_m$ so that $\omega_n - \tilde{\omega}_m^* \sim -i\Gamma_m/2$ and $\tilde{\omega}_n - \tilde{\omega}_m^* \sim -i(\Gamma_n + \Gamma_m)/2$. For resonances with similar linewidths, we therefore obtain $C_{nm} \sim U_{nm}$. The non-orthogonality of eigenfunctions therefore yields a non-vanishing degree of correlation between modal speckle patterns.

We define the modal selectivity for maximal modal coupling as the ratio of the strength of the selected mode in transmission over the incoherent sum of strengths of all modes

$$S_{\text{mode}} = \frac{T_n(\omega_n)}{\sum_{m=1}^{M} T_m(\omega_n)} \sim \frac{|C_{nm}|^2}{\sum_m |C_{nm}|^2}. \qquad (8)$$

A diagonal matrix $C$ would correspond to perfect modal selectivity with $S_{mode} = 1$. However, for non-Hermitian systems, modal selectivity falls below unity due to off-diagonal elements of $C_{nm}$. Thus the bi-orthogonality of the eigenfunctions of a non-Hermitian system reduces the degree of modal selectivity for maximal coupling to a mode. This is illustrated in the analytical analysis of a two-level non-Hermitian effective Hamiltonian model.

**Two-level effective Hamiltonian.** The modes of the system are expressed in the basis of the two modes of the closed cavity$^{39,40}$

$$\mathbf{H}_{\text{eff}} = \begin{pmatrix} \omega_1 & 0 \\ 0 & \omega_2 \end{pmatrix} - \frac{i}{2} \begin{pmatrix} \Gamma_{11} & \Gamma_{12} \\ \Gamma_{12} & \Gamma_{22} \end{pmatrix} \qquad (9)$$

The parameters $\Gamma_{nm}$ are given by, $\Gamma_{nm} = \sum_{c=1}^{2N} V_n^c V_m^c$, where the vectors $V_n$ represent the coupling of the closed system to the leads (see Eq. (1))$^{39}$. The parameter $\gamma = \Gamma_{12}$ is the coupling parameter between the resonances. Diagonalizing $\mathbf{H}_{\text{eff}}$ gives the eigenvalues $\tilde{\omega}_{1,2} = \frac{1}{2}(\omega_2 + \omega_1) \mp \frac{1}{2}\sqrt{\epsilon^2 - \gamma^2} - \frac{i}{4}(\Gamma_{11} + \Gamma_{22})$, with $\epsilon = (\omega_2 - \omega_1) - \frac{i}{2}(\Gamma_{11} - \Gamma_{22})$. The eigenvectors $|\phi_n\rangle$ of the effective Hamiltonian $\mathbf{H}_{\text{eff}}$ can be written in the basis $\{|\psi_n\rangle$ of the unperturbed eigenvectors of the Hamiltonian of the closed system $H^{48}$

$$|\phi_1\rangle = \frac{1}{\sqrt{1-f^2}} \begin{pmatrix} 1 \\ -if \end{pmatrix}, \quad |\phi_2\rangle = \frac{1}{\sqrt{1-f^2}} \begin{pmatrix} if \\ 1 \end{pmatrix}. \qquad (10)$$

The mixing of the two eigenstates depends on the single parameter, $f = \gamma/(\epsilon + \sqrt{\epsilon^2 - \gamma^2})$. The degree of complexness of the eigenfunctions is the same for both modes, $q_1^2 = q_2^2 = f^2$. In the present case of the two-level Hamiltonian, $q_1^2$ increases from 0 for isolated modes ($\gamma \ll \epsilon$) to unity for $f=1$, which is the case of an exceptional point at which the eigenvalues $\tilde{\omega}_1$ and $\tilde{\omega}_2$ coalesce$^{61}$.

The two-level Hamiltonian model is illustrated in Fig. 5 by the eigenfunctions of two hybridized modes at $f_1 = 11.260$ GHz and $f_2 = 11.266$ GHz, which are isolated from other modes but overlap strongly with a degree of overlap between the modes of $\delta_{12} = 5.5$. The eigenfunctions are normalized following the bi-orthogonality condition, $\int^d r q_n^2(r) = 1$. There is a strong similarity between $\mathrm{Re}(\phi_1)$ and $\mathrm{Im}(\phi_2)$ and between $\mathrm{Im}(\phi_1)$ and $\mathrm{Re}(\phi_2)$, as anticipated in Eq. (10). The eigenfunctions give $q_1^2 = 0.86$ and $q_2^2 = 0.4$. The two values are not equal because of the weak overlap with other modes. We observe in Fig. 5e, $f$ that when the incident wave is maximally coupled to the first or second mode, the transmission spectra and the contribution of
<figure><img src="image_16.png" /><figcaption>**Fig. 5** Modal mixing. **a-d** The square of the real and imaginary parts, Re($\phi_n$)$^2$ (**a**, **c**) and Im($\phi_n$)$^2$ (**b**, **d**), of two strongly overlapping modes with resonance at $f_1 = 11.260$ GHz and $f_2 = 11.266$ GHz with linewidths $\Gamma_1/(2\pi) = 21.7$ MHz and $\Gamma_2/(2\pi) = 12.3$ MHz, respectively. **e, f** Transmission spectra (blue line) for maximal coupling to the first (**e**) and the second (**f**) mode. The dashed lines are the modal strengths. The small value of transmission in comparison to modal strengths in **e**, **f** is a result of strong destructive interference between the highly correlated speckle patterns of strongly overlapping modes</figcaption></figure>

the two modes are very similar. The correlation between the incident waveforms $W_{L1}$ and $W_{L2}$, $|W_{L1}^\dagger W_{L2}|/(||W_{L1}|| ||W_{L2}||)$ is 0.98. The modal mixing of two strongly overlapping modes and its impact on the modal selectivity are further confirmed in finite-element simulations in Supplementary Note 2.

The decrease of $S_{\text{mode}}$ with increasing degree of modal mixing is demonstrated analytically in Supplementary Note 3 within the framework of the two-level Hamiltonian model. By expressing the average degree of correlation between vectors $W_{L1}$ and $W_{L2}$ as a function of $f$ in the limit $N \gg 1$ and $f \ll 1$, we find using Eq. (4) that $S_{\text{mode}}$ only depends upon the modal overlap and $f$ with

$$S_{\text{mode}} \sim \left[ 1 + \frac{4f^2 \Gamma_1^2}{4\Delta_{12}^2 + \Gamma_2^2} \right]^{-1}, \quad \text{where } \Delta_{12} = \omega_2 - \omega_1.$$

**Average modal selectivity.** In order to investigate the average modal selectivity in a large number of samples, we carry out simulations utilizing the recursive Green's function method$^{62}$ in random quasi-1D samples connected to leads supporting $N$ channels to find $t(\omega)$ (see Methods). Four ensembles with modal overlap $\delta$ equal to 0.08, 0.11, 0.64, and 1.13 are studied. The number of channels is $N=10$ for $\delta = 0.08$, $N=16$ for $\delta = 0.11$ and 0.64, and $N=33$ for diffusive samples with $\delta > 1$. The HI method is applied to more than one hundred samples for each ensemble giving more than 4,000 modes.

The modal selectivity for maximal coupling $S_{\text{mode}}$ is shown as a function of the modal overlap $\delta$ in Fig. 6, and compared to $S_{\text{ran}}$ computed for a random incident wavefront. As expected, maximal coupling enhances modal selectivity, but $S_{\text{mode}}$ falls below the value expected in the case of isolated modes of unity,

<figure><img src="image_17.png" /><figcaption>**Fig. 6** Average modal selectivity. The average modal selectivity for maximal coupling $S_{\text{mode}}$ (blue circles) and for a random wavefront $S_{\text{ran}}$ (red crosses) is plotted for four ensembles as a function of the modal overlap $\delta$</figcaption></figure>

even for localized waves. For $\delta = 0.08$, modes overlap and interfere giving $S_{\text{mode}} = 0.92$. Modal selectivity is seen to decrease with increasing $\delta$ as a result of greater spectral overlap with a larger number of modes and consequent increased correlation between their MTMs.

## Discussion

We have considered the inverse problem of characterizing and controlling the modes within the sample on the basis of the properties of waves scattered from the sample. We have seen that the road to the control of modes and of wave properties related to modes runs through the MTM which can be obtained from
spectra of the TM. We have demonstrated that a single mode can be excited within the sample even in the case of a moderate modal overlap. The incident vector of the MTM of unit rank couples maximally to the mode with an enhancement by a factor N over excitation by a random wavefront. However, as modal overlap increases as a result of greater coupling of the modes to the environment through the boundaries of the sample, the bi-orthogonality of the eigenfunctions leads to increasing modal correlation so that the degree of modal control is reduced.

This investigation of selective excitation of quasi-normal modes in disordered system has been carried out in systems with resonances with high quality factors. The decomposition of the TM into MTMs is more challenging for diffusive waves when the degree of modal overlap is high. However, as Alpeggiani et al. have demonstrated analytically, the scattering matrix can be reconstructed from the far-field properties of the eigenmodes for any degree of overlap²⁶. The modal coefficients depend on all other contributing modes through a coupling matrix, but MTMs are still of unit rank. Future studies will be dedicated to exploring the characteristics of modes and the limits of selectivity in systems with strong modal overlap. This will advance a more comprehensive understanding of the relationship between eigenchannels, time-delay eigenstates and quasi-normal modes of open system.

Selecting specific modes in samples in which the wave is localized would make it possible to deliver energy to specific regions of a sample. In the case of moderate modal overlap in open random systems, the modes extend over the entire sample, even in absorbing samples, and so if a single mode or a small number of modes is selected, it would be possible to deliver energy to the center of the sample. In contrast, when, many modes overlap, the average profile of energy density within the sample is determined by the diffusion equation and is concentrated within an absorption depth of the sample $L_a = (Dt_a)^{1/2}$, where D is the diffusion coefficient and $\tau_a$ is the absorption time, and particularly the absorption associated with exciting gain in the medium⁶³. As a result, light emitted from excited samples with weakly overlapping modes, in which energy penetrates more deeply into a random medium, is longer lived than in samples in which stronger modal overlap. There is therefore greater opportunity for the emitted photons to stimulate emission before escaping the sample. The lasing threshold will consequently be lowered and narrow-line emission will be observed⁶⁴. In general the smaller the degree of modal overlap, the more it is possible to deposit energy into a random absorbing medium. Finding the MTM and then pumping from the front of the sample⁶⁴,⁶⁵ could thereby lower the threshold of random lasers via coherent feedback.

Modal decomposition of the TM has been demonstrated with microwave radiation but is in principle possible in optics. Measurement of spectrally resolved TM has indeed been reported recently⁶⁶. Modal selectivity could then be utilized to enhance light-matter interactions in photonic materials⁶⁷, solar cells⁶⁸, or biomedical optics⁶⁹. The use of MTMs is not restricted to random media but can also be applied to structured media such as optical microcavities³¹, or photonic crystals⁷⁰.

## Methods

**Experimental setup.** The aluminum cavity has length $L = 500$ mm, width $W = 268$ mm, and height $H = 8$ mm and only supports a single waveguide mode in the vertical dimension over the frequency range of the measurements. The antennas are waveguide to coax adapters designed for the Ku band (12–18 GHz). Measurements of individual elements of the TM between two antenna arrays are made with use of electro-mechanical switches and a vector network analyzer. The channels of the switches that are turned off are matched to a 50 Ω load so that the boundary conditions of the system do not change when different emitting source

and receiving antennas are used. The spacing between two antennas on the left and right sides of the sample is metallic as seen in Fig. 1a.

The field inside the waveguide is detected by inserting an antenna sequentially into a square grid of holes which are 4 mm in diameter and spaced by 8 mm on a side. The transmission coefficient $e_a(x, y, \omega)$ is measured between each source antenna and a wire antenna inserted 0.5 mm below the bottom of the 6-mm thick aluminum cover. The penetration depth of the antenna is small enough that it does not distort the field profile in the waveguide. The spatial energy distribution $I(x, y)$ for any incoming vector $v$ is then reconstructed from the coherent superposition of the fields arising from each source antenna, $I(x, y) = |\sum_a e_a(x, y)v_a|^2$.

Once the decomposition into modes of the field inside the sample, $e_a(x, y, \omega) = \sum_n e_n^m(x, y)\phi_n(\omega)$, has been obtained, the modal spatial profile is found from an average on incoming channels $I_n(x, y) = \langle |e_n^m(x, y)|^2 \rangle$. We also present in Fig. 2f the modal energy density profile for incoming vector $v$ computed from $W_n(x, y) = |\sum_a e_a^n(x, y)v_a|^2$. For maximal modal coupling $v = W_{Ln}||W_{Ln}$| and for vanishing coupling $v$ is orthogonal to $W_{Ln}$.

**Harmonic inversion.** The modal analysis of the TM is performed using the HI method to obtain the complex modal frequencies $\omega_m - i\Gamma_m$ within a range $\omega_{min} < \omega_m < \omega_{max}$. Following the algorithm described in ref. ⁵⁰, we choose a Fourier-type Krylov basis to extract the resonances. The HI method consists of solving a generalized eigenvalue problem applied to two matrices $U^{(0)}$ and $U^{(1)}$ of dimension $J \times J$ with $J = \tilde{N}'dt(\omega_{max} - \omega_{min})/(4\pi)$ which are created from a time signal of length $N'$ and time step $dt$. The generalized eigenvalue problem is $U^{(1)}B_n = u_nU^{(0)}B_n$, where $u_n$ and $B_n$ are the eigenvalues and the eigenvectors, respectively. The M non-zero eigenvalues yield the M unknown complex frequencies; the associated complex amplitudes can be computed from the eigenvectors.

However, we may miss a few modes by applying HI to a single spectrum of the TM since the modal speckle patterns are random vectors. The strength of a particular mode could thus vanish for a single spectrum when the incoming and outgoing channels correspond to two nodes of this mode. We therefore include several spectra from the TM in the generalized eigenvalue problem. We first perform an inverse Fourier transform of nine field transmission spectra randomly chosen from the $N^2$ spectra of the TM. The time window of the inverse Fourier transform is taken to be the time domain for which the signals are above the noise level. We then create 9 matrices $U_i^{(0)}$ and $U_i^{(1)}$ from each time signals $s_i(t)$. Those 9 matrices of dimensions $J \times J$ are then concatenated into two matrices of dimension $3J \times 3J$, $U_T^{(0)}$ and $U_T^{(1)}$. Because we are seeking for the same set of resonance time signals $s_i(t)$, we solve the generalized eigenvalue problem $U_T^{(1)}B_n = u_nU_T^{(0)}B_n$. We then extract the resonances $\tilde{\omega}_n$ from the significant eigenvalues $u_n$. However, this does not directly provide the associated modal transmission coefficients $t_{ba}^n$ which are the elements of the MTMs. These are then obtained from a simple inverse problem by fitting in the time domain each transmission coefficient of the TM using Eq. (4) of the main text: $t_{ba}(\omega) = \sum_n \frac{t_{ba}^n}{\omega - \omega_n + i\Gamma_n/2}$.

**Recursive Green's function simulations.** The Green's functions between points at the input and output surfaces of a waveguide are obtained by solving the two-dimensional wave equation $\nabla^2\psi(x, y, \omega) + k_n^2\epsilon(x, y)\psi(x, y, \omega) = 0$ using the recursive Green's function method. The random dielectric permittivity $\epsilon(x, y)$ is drawn from a rectangular distribution centered on unity. The field transmission coefficients between each of the incoming modes $a$ and outgoing modes $b$ at frequency $\omega$, $t_{ba}(\omega)$, are then calculated using the projection of the Green's function onto the modes of the empty waveguide. Spectra of the TM is then obtained by computing the TM for each frequency over a frequency range which is much larger than the typical linewidth of the resonances.

## Data availability

The authors declare that all data that support the findings of this study are available from Matthieu Davy at matthieu.davy@univ-rennes1.fr upon reasonable request.

Received: 28 June 2018 Accepted: 17 October 2018

Published online: 09 November 2018

## References

1. Mosk, A. P., Lagendijk, A., Lerosey, G. & Fink, M. Controlling waves in space and time for imaging and focusing in complex media. *Nat. Photon.* **6**, 283–292 (2012).

2. Rotter, S. & Gigan, S. Light fields in complex media: mesoscopic scattering meets wave control. *Rev. Mod. Phys.* **89**, 015005 (2017).

3. Choi, W., Mosk, A. P., Park, Q. H. & Choi, W. Transmission eigenchannels in a disordered medium. *Phys. Rev. B* **83**, 134207 (2011).
4. Shi, Z. & Genack, A. Z. Transmission eigenvalues and the bare conductance in the crossover to anderson localization. *Phys. Rev. Lett.* **108**, 043901 (2012).

5. Liew, S. F., Popoff, S. M., Mosk, A. P., Vos, W. L. & Cao, H. Transmission channels for light in absorbing random media: from diffusive to ballistic-like transport. *Phys. Rev. B* **89**, 224202 (2014).

6. Gérardin, B., Laurent, J., Derode, A., Prada, C. & Aubry, A. Full transmission and reflection of waves propagating through a maze of disorder. *Phys. Rev. Lett.* **113**, 173901 (2014).

7. DavyM., & ShiZ. & ParkJ. & TianC. & GenackA. Z. Universal structure of transmission eigenchannels inside opaque media. *Nat. Commun.* **6**, 6893 (2015).

8. Yamilov, A., Petrenko, S., Sarma, R. & Cao, H. Shape dependence of transmission, reflection, and absorption eigenvalue densities in disordered waveguides with dissipation. *Phys. Rev. B* **93**, 100201 (2016).

9. Sarma, R., Yamilov, A. G., Petrenko, S., Bromberg, Y. & Cao, H. Control of energy density inside a disordered medium by coupling to open or closed channels. *Phys. Rev. Lett.* **117**, 086803 (2016).

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

11. Fisher, D. S. & Lee, P. A. Relation between conductivity and transmission matrix. *Phys. Rev. B* **23**, 6851-6854 (1981).

12. Dorokhov, O. N. On the coexistence of localized and extended electronic states in the metallic phase. *Solid State Commun.* **51**, 381-384 (1984).

13. Pretre, J. B. P. A. Ma. A. B. Maximal fluctuations—a new phenomenon in disordered system. *Phys. A* **168**, 7 (1990).

14. Imry, Y. Active transmission channels and universal conductance fluctuations. *Europhys. Lett.* **1**, 249-256 (1986).

15. Kim, M. et al. Maximal energy transport through disordered media with the implementation of transmission eigenchannels. *Nat. Photon.* **6**, 583-587 (2012).

16. Goetschy, A. & Stone, A. D. Filtering random matrices: the effect of incomplete channel control in multiple scattering. *Phys. Rev. Lett.* **111**, 063901 (2013).

17. Popoff, S. M., Goetschy, A., Liew, S. F., Stone, A. D. & Cao, H. Coherent control of total transmission of light through disordered media. *Phys. Rev. Lett.* **112**, 133903 (2014).

18. Davy, M., Shi, Z., Wang, J., Cheng, X. & Genack, A. Z. Transmission eigenchannels and the densities of states of random media. *Phys. Rev. Lett.* **114**, 033901 (2015).

19. Wigner, E. P. Lower limit for the energy derivative of the scattering phase shift. *Phys. Rev.* **98**, 145-147 (1955).

20. Smith, F. T. Lifetime matrix in collision theory. *Phys. Rev.* **118**, 349-356 (1960).

21. Rotter, S., Ambichl, P. & Libisch, F. Generating particlelike scattering states in wave transport. *Phys. Rev. Lett.* **106**, 120602 (2011).

22. Gérardin, B. et al. Particlelike wave packets in complex scattering systems. *Phys. Rev. B* **94**, 014209 (2016).

23. Shi, Z. & Genack, A. Z. Dynamic and spectral properties of transmission eigenchannels in random media. *Phys. Rev. B* **92**, 184202 (2015).

24. Ching, E. S. C. et al. Quasinormal-mode expansion for waves in open systems. *Rev. Mod. Phys.* **70**, 1545-1554 (1998).

25. Leung, P. T., Liu, S. Y. & Young, K. Completeness and orthogonality of quasinormal modes in leaky optical cavities. *Phys. Rev. A* **49**, 3057-3067 (1994).

26. Alpeggiani, F., Parappurath, N., Verhagen, E. & Kuipers, L. Quasinormal-mode expansion of the scattering matrix. *Phys. Rev. X* **7**, 021035 (2017).

27. Lalanne, P., Yan, W., Vynck, K., Sauvan, C. & Hugonin, J.-P. Light interaction with photonic and plasmonic resonances. *Laser Photon. Rev.* **12**, 170013 (2018).

28. Fyodorov, Y. V. & Sommers, H.-J. Statistics of resonance poles, phase shifts and time delays in quantum chaotic scattering: random matrix approach for systems with broken time-reversal invariance. *J. Math. Phys.* **38**, 1918-1981 (1997).

29. Rotter, I. Dynamics of quantum systems. *Phys. Rev. E* **64**, 036213 (2001).

30. Rotter, I. A non-Hermitian Hamilton operator and the physics of open quantum systems. *J. Phys. A* **42**, 153001 (2009).

31. Cao, H. & Wiersig, J. Dielectric microcavities: model systems for wave chaos and non-Hermitian physics. *Rev. Mod. Phys.* **87**, 61-111 (2015).

32. Sokolov, V. V. & Zelevinsky, V. G. Dynamics and statistics of unstable quantum states. *Nucl. Phys. A* **504**, 562-588 (1989).

33. Huang, W.-P. Coupled-mode theory for optical waveguides: an overview. *J. Opt. Soc. Am. A* **11**, 963-983 (1994).

34. Suh, W., Wang, Z. & Fan, S. Temporal coupled-mode theory and the presence of non-orthogonal modes in lossless multimode cavities. *IEEE J. Quantum Electron.* **40**, 1511-1518 (2004).

35. Vial, B. & Hao, Y. A coupling model for quasi-normal modes of photonic resonators. *J. Opt.* **18**, 115004 (2016).

36. Yang, J., Perrin, M. & Lalanne, P. Analytical formalism for the interaction of two-level quantum systems with metal nanoresonators. *Phys. Rev. X* **5**, 021008 (2015).

37. H. Lourenço-Martins, P. Das, L. H. Tizei, R. Weil and M. Kociak, Self-hybridization within non-Hermitian localized plasmonic systems, *Nat. Phys.* **14** 360–364 (2018).

38. Leung, P. T., Suen, W. M., Sun, C. P. & Young, K. Waves in open systems via a biorthogonal basis. *Phys. Rev. E* **57**, 6101–6104 (1998).

39. Savin, D. V., Legrand, O. & Mortessagne, F. Inhomogeneous losses and complexness of wave functions in chaotic cavities. *Europhys. Lett.* **76**, 774 (2006).

40. Xeridat, O., Poli, C., Legrand, O., Mortessagne, F. & Sebbah, P. Quasimodes of a chaotic elastic cavity with increasing local losses. *Phys. Rev. E* **80**, 035201 (2009).

41. Wigner, E. P. On the statistical distribution of the widths and spacings of nuclear resonance levels. *Math. Proc. Camb. Philos. Soc.* **47**, 790-798 (1951).

42. M. L. Mehta, Random Matrices, 3rd ed. (Academic Press, New York, 2004).

43. F. J. Dyson. A brownian-motion model for the eigenvalues of a random matrix. *J. Math. Phys.* **3**, 1191 (1962).

44. Thouless, D. J. Maximum metallic resistance in thin wires. *Phys. Rev. Lett.* **39**, 1167 (1977).

45. Kuhl, U., Höhmann, R., Main, J. & Stöckmann, H. J. Resonance widths in open microwave cavities studied by harmonic inversion. *Phys. Rev. Lett.* **100**, 254101 (2008).

46. Wang, J. & Genack, A. Z. Transport through modes in random media. *Nature* **471**, 345-348 (2011).

47. Brouwer, P. W. Wave function statistics in open chaotic billiards. *Phys. Rev. E* **68**, 046205 (2003).

48. Poli, C., Savin, D. V., Legrand, O. & Mortessagne, F. Statistics of resonance states in open chaotic systems: a perturbative approach. *Phys. Rev. E* **80**, 046203 (2009).

49. Kumar, S. et al. Distribution of scattering matrix elements in quantum chaotic scattering. *Phys. Rev. Lett.* **111**, 030403 (2013).

50. Mandelshtam, V. A. & Taylor, H. S. Harmonic inversion of time signals and its applications. *J. Chem. Phys.* **107**, 6756-6769 (1997).

51. Main, J. Use of harmonic inversion techniques in semiclassical quantization and analysis of quantum spectra. *Phys. Rep.* **316**, 233-338 (1999).

52. Abrahams, E., Anderson, P. W., Licciardello, D. C. & Ramakrishnan, T. V. Scaling theory of localization: absence of quantum diffusion in two dimensions. *Phys. Rev. Lett.* **42**, 673 (1979).

53. Mirlin, A. D. Statistics of energy levels and eigenfunctions in disordered systems. *Phys. Rep.* **326**, 259-382 (2000).

54. Peña, A., Girschik, A., Libisch, F., Rotter, S. and Chabanov, A., The single-channel regime of transport through random media, *Nat. Commun.* **5**, 3488 (2014).

55. Hsu, C. W., Goetschy, A., Bromberg, Y., Stone, A. D. & Cao, H. Broadband coherent enhancement of transmission and absorption in disordered media. *Phys. Rev. Lett.* **115**, 223901 (2015).

56. Bell J. and Steinberger J. In:Proc. Oxford International Conference on *Elementary Particles* (Oxford, UK, 1965).

57. Schomerus, H., Frahm, K. M., Patra, M. & Beenakker, C. W. J. Quantum limit of the laser line width in chaotic cavities and statistics of residues of scattering matrix poles. *Phys. A* **278**, 469-496 (2000).

58. Fyodorov, Y. V. & Savin, D. V. Statistics of resonance width shifts as a signature of eigenfunction nonorthogonality. *Phys. Rev. Lett.* **108**, 184101 (2012).

59. Vanneste, C. & Sebbah, P. Complexity of two-dimensional quasimodes at the transition from weak scattering to Anderson localization. *Phys. Rev. A* **79**, 041802 (2009).

60. Chalker, J. T. & Mehlig, B. Eigenvector statistics in non-hermitian random matrix ensembles. *Phys. Rev. Lett.* **81**, 3367-3370 (1998).

61. Eleuch, H. & Rotter, I. Nearby states in non-Hermitian quantum systems I: two states. *Eur. Phys. J. D.* **69**, 229 (2015).

62. Baranger, H. U., DiVincenzo, D. P., Jalabert, R. A. & Stone, A. D. Classical and quantum ballistic-transport anomalies in microjunctions. *Phys. Rev. B* **44**, 10637 (1991).

63. Genack, A. & Drake, J. Scattering for super-radiation. *Nature* **368**, 400–401 (1994).

64. Cao, H. et al. Random laser action in semiconductor powder. *Phys. Rev. Lett.* **82**, 2278-2281 (1999).

65. Bachelard, N., Gigan, S., Noblin, X. & Sebbah, P. Adaptive pumping for spectral control of random lasers. *Nat. Phys.* **10**, 426 (2014).
66. Andreoli, D. et al. Deterministic control of broadband light through a multiply scattering medium via the multispectral transmission matrix. *Sci. Rep.* 5, 10347 (2015).

67. Riboli, F. et al. Engineering of light confinement in strongly scattering disordered media. *Nat. Mater.* **13**, 720 (2014).

68. Polman, A. & Atwater, H. A. Photonic design principles for ultrahigh-efficiency photovoltaics. *Nat. Mater.* **11**, 174 (2012).

69. Yu, H. et al. Recent advances in wavefront shaping techniques for biomedical applications. *Curr. Appl. Phys.* **15**, 632–641 (2015).

70. Akahane, Y., Asano, T., Song, B.-S. & Noda, S. High-Q photonic nanocavity in a two-dimensional photonic crystal. *Nature* **425**, 944 (2003).

## Acknowledgements

This publication was supported by the European Union through the European Regional Development Fund (ERDF), by the French region of Brittany and Rennes Métropole through the CPER Project SOPHIE/STIC & Ondes, and by the National Science Foundation under grant DMR-/BSF: 1609218. We would also like to acknowledge Cécile Leconte for her help in automating the scan and Zhou Shi, Yan Fyodorov, Ulrich Kühl, Olivier Legrand, Fabrice Mortessagne, Jing Wang, and Eli Ashoush for stimulating discussions.

## Author contributions

M.D. carried out the experiment, numerical simulations, and theoretical analysis. A.Z.G. and M.D. conceived the project, discussed the results, and wrote the manuscript.

## Additional information

**Supplementary Information** accompanies this paper at https://doi.org/10.1038/s41467-018-07180-3.

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

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons

Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018