Speckle-correlation imaging through a kaleidoscopic multimode fiber.

Speckle-correlation imaging techniques are widely used for noninvasive imaging through complex scattering media. While light propagation through multimode fibers and scattering media share many analogies, reconstructing images through multimode fibers from speckle correlations remains an unsolved challenge. Here, we exploit a kaleidoscopic memory effect emerging in square-core multimode fibers and demonstrate fluorescence imaging with no prior knowledge on the fiber. Experimentally, our approach simply requires to translate random speckle patterns at the input of a square-core fiber and to measure the resulting fluorescence intensity with a bucket detector. The image of the fluorescent object is then reconstructed from the autocorrelation of the measured signal by solving an inverse problem. This strategy does not require the knowledge of the fragile deterministic relation between input and output fields, which makes it promising for the development of flexible minimally invasive endoscopes.

Imaging through a square multimode fiber by scanning focused spots with the memory effect.

The existence of a shift-shift memory effect in square waveguides, whereby any translation of the input field induces translations in the output field in four symmetrical directions, has been previously observed by correlation measurements. Here we demonstrate that this memory effect is also observed in real space and can be put to use for imaging purposes. First, a focus is created at the output of a square-core multimode fiber, by wavefront shaping based on feedback from a guide-star. Then, because of the memory effect, four symmetrical spots can be scanned at the fiber output by shifting the wavefront at the fiber input. We demonstrate that this property can be exploited to perform fluorescence imaging through the multimode fiber, without requiring the measurement of a transmission matrix.

Optical memory effect in square multimode fibers.

We demonstrate experimentally the existence of a translational optical memory effect in square-core multimode fibers. We found that symmetry properties of square-core waveguides lead to speckle patterns shifting along four directions at the fiber output for any given shift direction at the input. A simple theoretical model based on a perfectly reflective square waveguide is introduced to predict and interpret this phenomenon. We report experimental results obtained with 532-nm coherent light propagating through a square-core step-index multimode fiber, demonstrating that this translational memory effect can be observed for shift distances up to typically 10 µm after propagation through several centimeters of fiber.

Optimal Control of Coherent Light Scattering for Binary Decision Problems.

Because of quantum noise fluctuations, the rate of error achievable in decision problems involving several possible configurations of a scattering system is subject to a fundamental limit known as the Helstrom bound. Here, we present a general framework to calculate and minimize this bound using coherent probe fields with tailored spatial distributions. As an example, we experimentally study a target located in between two disordered scattering media. We first show that the optimal field distribution can be directly identified using a general approach based on scattering matrix measurements. We then demonstrate that this optimal light field successfully probes the presence of the target with a number of photons that is reduced by more than 2 orders of magnitude as compared to unoptimized fields.

Polyhedral Bubble Vibrations.

Underwater bubbles are extremely good acoustic resonators, but are freely evolving and dissolving. Recently it was found that bubbles can be stabilized in frames, but the influence of the frame shape is still undocumented. Here we first explore the vibration of polyhedral bubbles with a low number of faces, shaped as the five Platonic solids. Their resonance frequency is well approximated by the formula for spherical bubbles with the same volume. Then we extend these results to shapes with a larger number of faces using fullerenes, paving the way to obtain arbitrary large resonant bubbles.

Acoustics of Cubic Bubbles: Six Coupled Oscillators.

We introduce cubic bubbles that are pinned to 3D printed millimetric frames immersed in water. Cubic bubbles are more stable over time and space than standard spherical bubbles, while still allowing large oscillations of their faces. We find that each face can be described as a harmonic oscillator coupled to the other ones. These resonators are coupled by the gas inside the cube but also by acoustic interactions in the liquid. We provide an analytical model and 3D numerical simulations predicting the resonance with very good agreement. Acoustically, cubic bubbles prove to be good monopole subwavelength emitters, with nonemissive secondary surface modes.

Photoacoustic imaging beyond the acoustic diffraction-limit with dynamic speckle illumination and sparse joint support recovery.

In deep tissue photoacoustic imaging the spatial resolution is inherently limited by the acoustic wavelength. Recently, it was demonstrated that it is possible to surpass the acoustic diffraction limit by analyzing fluctuations in a set of photoacoustic images obtained under unknown speckle illumination patterns. Here, we purpose an approach to boost reconstruction fidelity and resolution, while reducing the number of acquired images by utilizing a compressed sensing computational reconstruction framework. The approach takes into account prior knowledge of the system response and sparsity of the target structure. We provide proof of principle experiments of the approach and demonstrate that improved performance is obtained when both speckle fluctuations and object priors are used. We numerically study the expected performance as a function of the measurement's signal to noise ratio and sample spatial-sparsity. The presented reconstruction framework can be applied to analyze existing photoacoustic experimental data sets containing dynamic fluctuations.

Overcoming the acoustic diffraction limit in photoacoustic imaging by the localization of flowing absorbers.

The resolution of photoacoustic imaging deep inside scattering media is limited by the acoustic diffraction limit. In this Letter, taking inspiration from super-resolution imaging techniques developed to beat the optical diffraction limit, we demonstrate that the localization of individual optical absorbers can provide super-resolution photoacoustic imaging well beyond the acoustic diffraction limit. As a proof-of-principle experiment, photoacoustic cross-sectional images of microfluidic channels were obtained with a 15 MHz linear capacitive micromachined ultrasonic transducer array, while absorbing beads were flown through the channels. The localization of individual absorbers allowed us to obtain a super-resolved cross-sectional image of the channels by reconstructing both the channel width and position with an accuracy better than λ/10. Given the discrete nature of endogenous absorbers such as red blood cells, or that of exogenous particular contrast agents, localization is a promising approach to push the current resolution limits of photoacoustic imaging.

Improving photoacoustic-guided optical focusing in scattering media by spectrally filtered detection.

We study the potential of photoacoustic guidance for light focusing through scattering samples via wavefront-shaping and iterative optimization. We experimentally demonstrate that the focusing efficiency on an extended absorber can be improved by iterative optimization of the high frequency components of the broadband photoacoustic signal detected with a spherically focused transducer. We demonstrate more than 12-fold increase in the photoacoustic signal generated by a 30 µm wire using a narrow frequency band around 60 MHz. By monitoring the speckle pattern evolution during the optimization process with a CCD camera, we experimentally confirm that such optimization leads to a smaller optical focus than what would be obtained by optimizing lower frequencies of the photoacoustic feedback.

Light focusing and two-dimensional imaging through scattering media using the photoacoustic transmission matrix with an ultrasound array.

We implement the photoacoustic transmission matrix approach on a two-dimensional photoacoustic imaging system, using a 15 MHz linear ultrasound array. Using a black leaf skeleton as a complex absorbing structure, we demonstrate that the photoacoustic transmission matrix approach allows to reveal structural features that are invisible in conventional photoacoustic images, as well as to selectively control light focusing on absorbing targets, leading to a local enhancement of the photoacoustic signal.

Accurate measurement of guided modes in a plate using a bidirectional approach.

Measuring guided wave propagation in long bones is of interest to the medical community. When an inclination exists between the probe and the tested specimen surface, a bias is introduced on the guided mode wavenumbers. The aim of this study was to generalize the bidirectional axial transmission technique initially developed for the first arriving signal. Validation tests were performed on academic materials such a bone-mimicking plate covered with either a silicon or fat-mimicking layer. For any inclination, the wavenumbers measured with the probe parallel to the waveguide surface can be obtained by averaging the wavenumbers measured in two opposite directions.

Imaging light fluence in blood vessels by combining photoacoustic fluctuation imaging and ultrasound power Doppler.

Objectives.Numerous optical biomedical imaging or therapeutic modalities suffer from unknown light fluence distribution at depths. Photoacoustic (PA) imaging, which enables imaging blood vessels at the acoustic resolution, probes the product between the fluence and effective optical absorption that depends on the size or density of blood vessels. In the case of unresolved vessels, fluence and absorption can not be decoupled using PA imaging alone without the use of inverse problems. Thus, we propose combining two modalities that are sensitive to blood vessels to directly image fluence maps within vascularized areas, including in unresolved vessels.Approach.To achieve fluence imaging, the combination of photoacoustic fluctuation (PAFI) and Ultrasound Power Doppler (UPD) images is considered. After exposing a new theoretical expression of the UPD image, we establish a fluence imaging method giving quantitative fluence in blood vessels. Fluence imaging involves resolution compensation with a PSF filter that is compared to alternative simpler corrections.Main results.This method universally applies to arbitrary hematocrit and multi-scale vessel imaging. Using a spherical sparse array, we demonstrate 3D fluence imaging within blood vessels in simulation and experiments which is not possible with PAFI alone.Significance.Overall, we show that combining PAFI and UPD has the potential for real-time light dosimetry or could enhance quantitative inverse problems in PA imaging.

Improving visibility in photoacoustic imaging using dynamic speckle illumination.

In high-frequency photoacoustic imaging with uniform illumination, homogeneous photoabsorbing structures may be invisible because of their large size or limited-view issues. Here we show that, by exploiting dynamic speckle illumination, it is possible to reveal features that are normally invisible with a photoacoustic system comprised of a 20 MHz linear ultrasound array. We demonstrate imaging of a ∅5 mm absorbing cylinder and a 30 µm black thread arranged in a complex shape. The hidden structures are directly retrieved from photoacoustic images recorded for different random speckle illuminations of the phantoms by assessing the variation in the value of each pixel over the illumination patterns.

Simulations of ultrasound propagation in random arrangements of elliptic scatterers: occurrence of two longitudinal waves.

Ultrasound propagation in clusters of elliptic (two-dimensional) or ellipsoidal (three-dimensional) scatterers randomly distributed in a fluid is investigated numerically. The essential motivation for the present work is to gain a better understanding of ultrasound propagation in trabecular bone. Bone microstructure exhibits structural anisotropy and multiple wave scattering. Some phenomena remain partially unexplained, such as the propagation of two longitudinal waves. The objective of this study was to shed more light on the occurrence of these two waves, using finite-difference simulations on a model medium simpler than bone. Slabs of anisotropic, scattering media were randomly generated. The coherent wave was obtained through spatial and ensemble-averaging of the transmitted wavefields. When varying relevant medium parameters, four of them appeared to play a significant role for the observation of two waves: (i) the solid fraction, (ii) the direction of propagation relatively to the scatterers orientation, (iii) the ability of scatterers to support shear waves, and (iv) a continuity of the solid matrix along the propagation. These observations are consistent with the hypothesis that fast waves are guided by the locally plate/bar-like solid matrix. If confirmed, this interpretation could significantly help developing approaches for a better understanding of trabecular bone micro-architecture using ultrasound.

Full-visibility 3D imaging of oxygenation and blood flow by simultaneous multispectral photoacoustic fluctuation imaging (MS-PAFI) and ultrasound Doppler.

We present a method and setup that provide complementary three-dimensional (3D) images of blood oxygenation (via quantitative photoacoustic imaging) and blood flow dynamics (via ultrasound Doppler). The proposed approach is label-free and exploits blood-induced fluctuations, and is implemented on a sparse array with only 256 elements, driven with a commercially available ultrasound electronics. We first implement 3D photoacoustic fluctuation imaging (PAFI) to image chicken embryo, and obtain full-visibility images of the vascular morphology. We obtain simultaneously 3D ultrasound power Doppler with a comparable image quality. We then introduce multispectral photoacoustic fluctuation imaging (MS-PAFI), and demonstrate that it can provide quantitative measurements of the absorbed optical energy density with full visibility and enhanced contrast, as compared to conventional delay-and-sum multispectral photoacoustic imaging. We finally showcase the synergy and complementarity between MS-PAFI, which provides 3D quantitative oxygenation (SO[Formula: see text]) imaging, and 3D ultrasound Doppler, which provides quantitative information on blood flow dynamics. MS-PAFI represents a promising alternative to model-based inversions with the advantage of resolving all the visibility artefacts without prior and regularization, by use of a straightforward processing scheme.

Acousto-optical coherence tomography with a digital holographic detection scheme.

Acousto-optical coherence tomography (AOCT) consists in using random phase jumps on ultrasound and light to achieve a millimeter resolution when imaging thick scattering media. We combined this technique with heterodyne off-axis digital holography. Two-dimensional images of absorbing objects embedded in scattering phantoms are obtained with a good signal-to-noise ratio. We study the impact of the phase modulation characteristics on the amplitude of the acousto-optic signal and on the contrast and apparent size of the absorbing inclusion.

Attenuation, scattering, and absorption of ultrasound in the skull bone.

PURPOSE: Measured values of ultrasound attenuation in bone represent a combination of different loss mechanisms. As a wave is transmitted from a fluid into bone, reflections occur at the interface. In the bone, mode conversion occurs between longitudinal and shear modes and the mechanical wave is scattered by its complex internal microstructure. Finally, part of the wave energy is absorbed by the bone and converted into heat. Due to the complexity of the wave propagation and the difficulty in performing measurements that are capable of separating the various loss mechanisms, there are currently no estimates of the absorption in bone. The aim of this paper is, thus, to quantify the attenuation, scattering, and thermal absorption in bone. METHODS: An attenuating model of wave propagation in bone is established and used to develop a three-dimensional finite difference time domain numerical algorithm. Hydrophone and optical heterodyne interferometer measurements of the acoustic field as well as a x-ray microtomography of the bone sample are used to drive the simulations and to measure the attenuation. The acoustic measurements are performed concurrently with an infrared camera that can measure the temperature elevation during insonication. A link between the temperature and the absorption via a three-dimensional thermal simulation is then used to quantify the absorption coefficients for longitudinal and shear waves in cortical bone. RESULTS: We demonstrate that only a small part of the attenuation is due to absorption in bone and that the majority of the attenuation is due to reflection, scattering, and mode conversion. In the nine samples of a human used for the study, the absorption time constant for cortical bone was determined to be 1.04 µs ± 28%. This corresponds to a longitudinal absorption of 2.7 dB/cm and a shear absorption of 5.4 dB/cm. The experimentally measured attenuation across the approximately 8 mm thick samples was 13.3 ± 0.97 dB/cm. CONCLUSIONS: This first measurement of ultrasound absorption in bone can be used to estimate the amount of heat deposition based on knowledge of the acoustic field.

Compensating for visibility artefacts in photoacoustic imaging with a deep learning approach providing prediction uncertainties.

Conventional photoacoustic imaging may suffer from the limited view and bandwidth of ultrasound transducers. A deep learning approach is proposed to handle these problems and is demonstrated both in simulations and in experiments on a multi-scale model of leaf skeleton. We employed an experimental approach to build the training and the test sets using photographs of the samples as ground truth images. Reconstructions produced by the neural network show a greatly improved image quality as compared to conventional approaches. In addition, this work aimed at quantifying the reliability of the neural network predictions. To achieve this, the dropout Monte-Carlo procedure is applied to estimate a pixel-wise degree of confidence on each predicted picture. Last, we address the possibility to use transfer learning with simulated data in order to drastically limit the size of the experimental dataset.

Super-resolution photoacoustic and ultrasound imaging with sparse arrays.

It has previously been demonstrated that model-based reconstruction methods relying on a priori knowledge of the imaging point spread function (PSF) coupled to sparsity priors on the object to image can provide super-resolution in photoacoustic (PA) or in ultrasound (US) imaging. Here, we experimentally show that such reconstruction also leads to super-resolution in both PA and US imaging with arrays having much less elements than used conventionally (sparse arrays). As a proof of concept, we obtained super-resolution PA and US cross-sectional images of microfluidic channels with only 8 elements of a 128-elements linear array using a reconstruction approach based on a linear propagation forward model and assuming sparsity of the imaged structure. Although the microchannels appear indistinguishable in the conventional delay-and-sum images obtained with all the 128 transducer elements, the applied sparsity-constrained model-based reconstruction provides super-resolution with down to only 8 elements. We also report simulation results showing that the minimal number of transducer elements required to obtain a correct reconstruction is fundamentally limited by the signal-to-noise ratio. The proposed method can be straigthforwardly applied to any transducer geometry, including 2D sparse arrays for 3D super-resolution PA and US imaging.

Single-shot hybrid photoacoustic-fluorescent microendoscopy through a multimode fiber with wavefront shaping.

We present a minimally-invasive endoscope based on a multimode fiber that combines photoacoustic and fluorescence sensing. From the measurement of a transmission matrix during a prior calibration step, a focused spot is produced and raster-scanned over a sample at the distal tip of the fiber by use of a fast spatial light modulator. An ultra-sensitive fiber-optic ultrasound sensor for photoacoustic detection placed next to the fiber is combined with a photodetector to obtain both fluorescence and photoacoustic images with a distal imaging tip no larger than 250 µm. The high signal-to-noise ratio provided by wavefront shaping based focusing and the ultra-sensitive ultrasound sensor enables imaging with a single laser shot per pixel, demonstrating fast two-dimensional hybrid in vitro imaging of red blood cells and fluorescent beads.