Alexandre AUBRY

Résumé: As light travels through a disordered medium such as biological tissues, it undergoes multiple scattering events. This phenomenon is detrimental to in-depth optical microscopy, as it causes a drastic degradation of contrast, resolution and brightness of the resulting image beyond a few scattering mean free paths. However, the information about the inner reflectivity of the sample is not lost; only scrambled. To recover this information, a matrix approach of optical imaging can be fruitful. Here, we report on a de-scanned measurement of a high-dimension reflection matrix R via low coherence interferometry. Then, we show how a set of independent focusing laws can be extracted for each medium voxel through an iterative multi-scale analysis of wave distortions contained in R. It enables an optimal and local compensation of forward multiple scattering paths and provides a three-dimensional confocal image of the sample as the latter one had become digitally transparent. The proof-of-concept experiment is performed on a human opaque cornea and an extension of the penetration depth by a factor five is demonstrated compared to the state-of-the-art.

Résumé: Matrix imaging paves the way towards a next revolution in wave physics. Based on the response matrix recorded between a set of sensors, it enables an optimized compensation of aberration phenomena and multiple scattering events that usually drastically hinder the focusing process in heterogeneous media. Although it gave rise to spectacular results in optical microscopy or seismic imaging, the success of matrix imaging has been so far relatively limited with ultrasonic waves because wave control is generally only performed with a linear array of transducers. In this paper, we extend ultrasound matrix imaging to a 3D geometry. Switching from a 1D to a 2D probe enables a much sharper estimation of the transmission matrix that links each transducer and each medium voxel. Here, we first present an experimental proof of concept on a tissue-mimicking phantom through ex-vivo tissues and then, show the potential of 3D matrix imaging for transcranial applications.

Résumé: This is the first article in a series of two dealing with a matrix approach for aberration quantification and correction in ultrasound imaging. Advanced synthetic beamforming relies on a double focusing operation at transmission and reception on each point of the medium. Ultrasound matrix imaging (UMI) consists in decoupling the location of these transmitted and received focal spots. The response between those virtual transducers form the so-called focused reflection matrix that actually contains much more information than a confocal ultrasound image. In this paper, a time-frequency analysis of this matrix is performed, which highlights the single and multiple scattering contributions as well as the impact of aberrations in the monochromatic and broadband regimes. Interestingly, this analysis enables the measurement of the incoherent input-output point spread function at any pixel of this image. A fitting process enables the quantification of the single scattering, multiple scattering and noise components in the image. From the single scattering contribution, a focusing criterion is defined, and its evolution used to quantify the amount of aberration throughout the ultrasound image. In contrast to the state-of-the-art coherence factor, this new indicator is robust to multiple scattering and electronic noise, thereby providing a contrasted map of the focusing quality at a much better transverse resolution. After a validation of the proof-of-concept based on time-domain simulations, UMI is applied to the in-vivo study of a human calf. Beyond this specific example, UMI opens a new route for speed-of-sound and scattering quantification in ultrasound imaging.

Résumé: The last decade has seen the development of a wide set of tools, such as wavefront shaping, computational or fundamental methods, that allow us to understand and control light propagation in a complex medium, such as biological tissues or multimode fibers. A vibrant and diverse community is now working in this field, which has revolutionized the prospect of diffraction-limited imaging at depth in tissues. This roadmap highlights several key aspects of this fast developing field, and some of the challenges and opportunities ahead.

Résumé: We succeeded in freeze-drying monodisperse microbubbles without degrading their performance, that is, their monodispersity in size and echogenicity. We used microfluidic technology to generate cryoprotected highly monodisperse microbubbles (coefficient of variation [CV] <5%). By using a novel retrieval technique, we were able to freeze-dry the microbubbles and resuspend them without degradation, that is, keeping their size distribution narrow (CV <6%). Acoustic characterization performed in two geometries (a centimetric cell and a millichannel) revealed that the resuspended bubbles conserved the sharpness of the backscattered resonance peak, leading to CVs ranging between 5% and 10%, depending on the geometry. As currently observed with monodisperse bubbles, the peak amplitudes are one order of magnitude higher than those of commercial ultrasound contrast agents. Our work thus solves the question of storage and transportation of highly monodisperse bubbles. This work might open pathways toward novel clinical non-invasive measurements, such as local pressure, impossible to carry out with the existing commercial ultrasound contrast agents.

Résumé: Fault zones are associated with multi-scale heterogeneities of rock properties. Large scale variations may be imaged with conventional seismic reflection methods that detect offsets in geological units, and tomographic techniques that provide average seismic velocities in resolved volumes. However, characterizing elementary localized inhomogeneities of fault zones, such as cracks and fractures, constitutes a challenge for conventional techniques. Resolving these small-scale heterogeneities can provide detailed information for structural and mechanical models of fault zones. Recently, the reflection matrix approach utilizing body wave reflections in ambient noise cross-correlations was extended with the introduction of aberration corrections to handle the actual lateral velocity variations in the fault zone (Touma et al., 2021). Here this method is applied further to analyze the distribution of scatterers in the first few kilometers of the crust in the San Jacinto Fault Zone at the Sage Brush Flat (SGB) site, southeast of Anza, California. The matrix approach allows us to image not only specular reflectors but also to resolve the presence, location and reflectivity of scatterers for seismic waves starting with a simple homogeneous background velocity model of the medium. The derived three-dimensional image of the fault zone resolves lateral variations of scattering properties in the region within and around the surface fault traces, as well as differences between the Northwest (NW) and the Southeast (SE) parts of the study area. A localized intense damage zone at depth is observed in the SE section, suggesting that a geometrical complexity of the fault zone at depth induces ongoing generation of rock damage.

Résumé: Reflection seismic imaging usually suffers from a loss of resolution and contrast because of the fluctuations of the wave velocities in the Earth's crust. In the literature, phase distortion issues are generally circumvented by means of a background wave velocity model. However, it requires a prior tomography of the wave velocity distribution in the medium, which is often not possible, especially in depth. In this paper, a matrix approach of seismic imaging is developed to retrieve a 3-D image of the subsoil, despite a rough knowledge of the background wave velocity. To do so, passive noise cross-correlations between geophones of a seismic array are investigated under a matrix formalism. They form a reflection matrix that contains all the information available on the medium. A set of matrix operations can then be applied in order to extract the relevant information as a function of the problem considered. On the one hand, the background seismic wave velocity can be estimated and its fluctuations quantified by projecting the reflection matrix in a focused basis. It consists in investigating the response between virtual sources and detectors synthesized at any point in the medium. The minimization of their cross-talk can then be used as a guide star for approaching the actual wave velocity distribution. On the other hand, the detrimental effect of wave velocity fluctuations on imaging is overcome by introducing a novel mathematical object: The distortion matrix. This operator essentially connects any virtual source inside the medium with the distortion that a wavefront, emitted from that point, experiences due to heterogeneities. A time reversal analysis of the distortion matrix enables the estimation of the transmission matrix that links each real geophone at the surface and each virtual geophone in depth. Phase distortions can then be compensated for any point of the underground. Applied to passive seismic data recorded along the Clark branch of the San Jacinto fault zone (SJFZ), the present method is shown to provide an image of the fault until a depth of 4 km over the frequency range 10-20Hz with an horizontal resolution of 80 m. Strikingly, this resolution is almost one eighth below the diffraction limit imposed by the geophone array aperture. The heterogeneities of the subsoil play the role of a scattering lens and of a transverse waveguide which increase drastically the array aperture. The contrast is also optimized since most of the incoherent noise is eliminated by the iterative time reversal process. Beyond the specific case of the SJFZ, the reported approach can be applied to any scales and areas for which a reflection matrix is available at a spatial sampling satisfying the Nyquist criterion.

Résumé: © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC). In optical imaging, light propagation is affected by the inhomogeneities of the medium. Sample-induced aberrations and multiple scattering can strongly degrade the image resolution and contrast. On the basis of a dynamic correction of the incident and/or reflected wavefronts, adaptive optics has been used to compensate for those aberrations. However, it only applies to spatially invariant aberrations or to thin aberrating layers. Here, we propose a global and noninvasive approach based on the distortion matrix concept. This matrix basically connects any focusing point of the image with the distorted part of its wavefront in reflection. A singular value decomposition of the distortion matrix allows to correct for high-order aberrations and forward multiple scattering over multiple isoplanatic modes. Proof-of-concept experiments are performed through biological tissues including a turbid cornea. We demonstrate a Strehl ratio enhancement up to 2500 and recover a diffraction-limited resolution until a depth of 10 scattering mean free paths.

Résumé: The propagation of waves in complex media can be harnessed either by taming the incident wave-field impinging on the medium or by forcing waves along desired paths through its careful design. These two alternative strategies have given rise to fascinating concepts such as time reversal or negative refraction. Here, we show how these two processes are intimately linked through the negative reflection phenomenon. A negative reflecting mirror converts a wave of positive phase velocity into its negative counterpart and vice versa. In this article, we experimentally demonstrate this phenomenon with elastic waves in a 2D billiard and in a disordered plate by means of laser interferometry. Despite the complexity of such configurations, the negatively reflected wave field focuses back towards the initial source location, thereby mimicking a phase conjugation operation while being a fully passive process. The super-focusing capability of negative reflection is also highlighted in a monochromatic regime. The negative reflection phenomenon is not restricted to guided elastic waves since it can occur in zero-gap systems such as photonic crystals, chiral metamaterials or graphene. Negative reflection can thus become a tool of choice for the control of waves in all fields of wave physics.

Résumé: © 2018 American Physical Society. This paper provides a theoretical investigation of negative refraction and focusing of elastic guided waves in a freestanding plate with a steplike thickness change. Under certain conditions, a positive phase velocity (forward) Lamb mode can be converted into a negative phase velocity (backward) mode at such interface, giving rise to negative refraction. A semianalytical model is developed to study the influence of various parameters such as the material Poisson's coefficient, the steplike thickness, the frequency, and the incidence angle. To this end, all the Lamb and shear horizontal propagating modes and also a large number of their inhomogeneous and evanescent counterparts are taken into account. The boundary conditions applied to the stress-displacement fields at the thickness step yields an equation system. Its inversion provides the transmission and reflection coefficients between each mode at the interface. The steplike thickness and Poisson's ratio are shown to be key parameters to optimize the negative refraction process. In terms of material, Duralumin is found to be optimal as it leads to a nearly perfect conversion between forward and backward modes over broad frequency and angular ranges. An excellent focusing ability is thus predicted for a flat lens made of two symmetric thickness steps. A laser ultrasonic experiment quantitatively confirms those theoretical predictions. This study paves the way toward the optimization of elastic devices based on negative refraction, in particular for cloaking or superfocusing purposes.

Résumé: © 2017 Acoustical Society of America.Ultrasonic evaluation of coarse-grain materials generates multiple scattering at high frequency and large depth. Recent academic experiments with array probes showed the ability of a random matrix method [multiple scattering filter (MSF)] to reduce multiple scattering, hence improving detection. Here, MSF is applied to an industrial nickel-based alloy with coarse-grain structure. Two samples with average grain sizes 90 ± 60 μm and 750 ± 400 μm are inspected with wide-band 64-element arrays at central frequencies 2, 3, and 5 MHz. They contain cylindrical through-holes (1-mm radius) at various depths. The array transfer matrix is recorded and post-processed both in the flawless area and for eleven positions above each defect, which allows for a statistical analysis. MSF is compared with two conventional imaging techniques: the total focusing method (TFM) and the decomposition of the time-reversal operator (DORT). Several parameters to assess the performance of detection techniques are proposed and discussed. The results show the benefit of MSF, especially at high frequencies and for deep defects: at 5 MHz and 70 mm depth, i.e., more than three scattering mean-free paths, the detection rate for MSF ranges between 55% and 100% while it is found to be 0% both for TFM and DORT.

Résumé: © 2015 American Physical Society. We report on the passive measurement of time-dependent Green's functions in the optical frequency domain with low-coherence interferometry. Inspired by previous studies in acoustics and seismology, we show how the correlations of a broadband and incoherent wave field can directly yield the Green's functions between scatterers of a complex medium. Both the ballistic and multiple scattering components of the Green's function are retrieved. This approach opens important perspectives for optical imaging and characterization in complex scattering media.

Résumé: We report on the experimental measurement of the backscattering matrix of a weakly scattering medium in optics, composed of a few dispersed gold nanobeads. The decomposition of the time-reversal operator is applied to this matrix and we demonstrate selective and efficient focusing on individual scatterers, even through an aberrating layer. Moreover, we show that this approach provides the decomposition of the scattering pattern of a single nanoparticle. These results open important perspectives for optical imaging, characterization, and selective excitation of nanoparticles. © 2011 American Physical Society.

Résumé: The propagation of surface plasmons along rough metal surfaces is investigated with transformation optics. The roughness is modeled on a nanometer scale either by partly embedding a cylinder of metal into the surface (convex rough surface) or by excavating a cylindrical cavity from it (concave rough surface). These two structures can be treated analytically by means of conformal transformation. The interaction of surface plasmons with the singularities of these structures is shown to induce extreme field enhancements. These modes dominate the surface-enhanced Raman-scattering response and enhancement factors of the order of 107 are predicted. Interestingly, concave rough surfaces are shown to be the best candidates for surface-enhanced Raman scattering due to a stronger field enhancement and a lower sensitivity to the incident light polarization. Our analytical approach also points out the influence of the contact angle between the asperities and the metal surface on the bandwidth and the efficiency of the light-harvesting process. © 2011 American Physical Society.

Résumé: The distribution of singular values of the propagation operator in a random medium is investigated, in a backscattering configuration. Experiments are carried out with pulsed ultrasonic waves around 3 MHz, using an array of 64 programmable transducers placed in front of a random scattering medium. The impulse responses between each pair of transducers are measured and form the response matrix. The evolution of its singular values with time and frequency is computed by means of a short-time Fourier analysis. The mean distribution of singular values exhibits a very different behaviour in the single and multiple scattering regimes. The results are compared with random matrix theory. Once the experimental matrix coefficients are renormalized, experimental results and theoretical predictions are found to be in a very good agreement. Two kinds of random media have been investigated: a highly scattering medium in which multiple scattering predominates and a weakly scattering medium. In both cases, residual correlations that may exist between matrix elements are shown to be a key parameter. Finally, the possibility of detecting a target embedded in a random scattering medium based on the statistical properties of the strongest singular value is discussed. © 2010 Taylor & Francis.

Résumé: The singular values distribution of the propagation operator in a random medium is investigated in a backscattering configuration. Experiments are carried out with pulsed ultrasonic waves around 3 MHz, using an array of transducers. Coherent backscattering and field correlations are taken into account. Interestingly, the distribution of singular values shows a dramatically different behavior in the single and multiple-scattering regimes. Based on a matrix separation of single and multiple-scattered waves, an experimental illustration of imaging through a highly scattering slab is presented. © 2009 The American Physical Society.

Résumé: As classical imaging fails with diffusive media, one way to image a multiple-scattering medium is to achieve local measurements of the dynamic transport properties of a wave undergoing diffusion. This paper presents a method to obtain local measurements of the diffusion constant D in a multiple-scattering medium. The experimental setup consists in an array of programmable transducers placed in front of the multiple-scattering medium to be imaged. By achieving Gaussian beamforming both at emission and reception, an array of virtual sources and receivers located in the near field is constructed. The time evolution of the incoherent component of the intensity backscattered on this virtual array is shown to represent directly the growth of the diffusive halo as Dt. A matrix treatment is proposed to separate the incoherent intensity from the coherent backscattering peak. Once the incoherent contribution is isolated, a local measurement of the diffusion constant is possible. The technique is applied to image the long-scale variations of D in a random-scattering sample made of two parts with a different concentration of cylindrical scatterers. This experimental result is obtained with ultrasonic waves around 3 MHz. It illustrates the possibility of imaging diffusive media from local measurements of the diffusion constant, based on coherent Gaussian beamforming and a matrix "antisymmetrization," which creates a virtual antireciprocity. © 2007 The American Physical Society.

Résumé: The DORT method (French acronym for decomposition of the time reversal operator) is an active remote sensing technique using an array of antennas for the detection and localization of scatterers. This method is based on the singular value decomposition of the interelement response matrix. In this paper an analytical expression of the singular vectors due to the reflection from a large rigid cylinder is provided. Depending on the array aperture, two asymptotic regimes are described. It is shown that the singular vectors correspond to Hermite-Gaussian modes for large apertures and Legendre polynomials for small ones. Using perturbation theory, the corresponding singular values are deduced. Theoretical predictions are in good agreement with experimental results. © 2006 Acoustical Society of America.