Passive imaging with pulsed ultrasound insonations

Real-time passive cavitation mapping and B-mode fusion imaging via hybrid adaptive beamformer with modified diagnostic ultrasound platform

The implementation of real-time, convenient and high-resolution passive cavitation imaging (PCM) is crucial for ensuring the safety and effectiveness of ultrasound applications related to cavitation effects. However, the current B-mode ultrasound imaging system cannot achieve these functions. By developing a hybrid adaptive beamforming algorithm, the current work presented a real-time PCM and B-mode fusion imaging technique, using a modified diagnostic ultrasound platform enabling time-division multiplexing external triggering function. The proposed hybrid adaptive beamformer combined the advantages of delay-multiply-and-sum (DMAS) and minimum variance (MV) methods to effectively suppress the side lobe and tail-like artifacts, improving the resolution of PCM images. A high-pass filter was applied to selectively detect cavitation-specific signals while removing the interference from the tissue scatters. The system enabled synchronous visualization of tissue structure and cavitation activity under ultrasound exposure. Both numerical and experimental studies demonstrated that, compared with DAS, MV-DAS and DMAS methods, the proposed MV-DMAS algorithm performed better in both axial and lateral resolutions. This work represented a significant advancement in achieving high-quality real-time B-mode and PCM fusion imaging utilizing commercial medical ultrasound system, providing a powerful tool for synchronous monitoring and manipulating cavitation activity, which would enhance the safety and efficacy of cavitation-based applications.

Passive Acoustic Mapping for Convex Arrays With the Helical Wave Spectrum Method

Passive acoustic mapping (PAM) has emerged as a valuable imaging modality for monitoring the cavitation activity in focused ultrasound therapies. When it comes to imaging in the human abdomen, convex arrays are preferred due to their large acoustic window. However, existing PAM methods for convex arrays rely on the computationally expensive delay-and-sum (DAS) operation limiting the image reconstruction speed when the field-of-view (FOV) is large. In this work, we propose an efficient and frequency-selective PAM method for convex arrays. This method is based on projecting the helical wave spectrum (HWS) between cylindrical surfaces in the imaging field. Both the in silico and in vitro experiments showed that the HWS method has comparable image quality and similar acoustic cavitation source localization accuracy as the DAS-based methods. Compared to the frequency-domain and time-domain DAS methods, the time-complexity of the HWS method is reduced by one order and two orders of magnitude, respectively. A parallel implementation of the HWS method realized millisecond-level image reconstruction speed. We also show that the HWS method is inherently capable of mapping microbubble (MB) cavitation activity of different status, i.e., no cavitation, stable cavitation, or inertial cavitation. After compensating for the lens effects of the convex array, we further combined PAM formed by the HWS method and B-mode imaging as a real-time dual-mode imaging approach to map the anatomical location where MBs cavitate in a liver phantom experiment. This method may find use in applications where convex arrays are required for cavitation activity monitoring in real time.

Three-dimensional super resolution ultrasound imaging with a multi-frequency hemispherical phased array

BACKGROUND

High resolution imaging of the microvasculature plays an important role in both diagnostic and therapeutic applications in the brain. However, ultrasound pulse-echo sonography imaging the brain vasculatures has been limited to narrow acoustic windows and low frequencies due to the distortion of the skull bone, which sacrifices axial resolution since it is pulse length dependent.

PURPOSE

To overcome the detect limit, a large aperture 256-module sparse hemispherical transmit/receive array was used to visualize the acoustic emissions of ultrasound-vaporized lipid-coated decafluorobutane nanodroplets flowing through tube phantoms and within rabbit cerebral vasculature in vivo via passive acoustic mapping and super resolution techniques.

METHODS

Nanodroplets were vaporized with 55 kHz burst-mode ultrasound (burst length = 145 μs, burst repetition frequency = 9-45 Hz, peak negative acoustic pressure = 0.10-0.22 MPa), which propagates through overlying tissues well without suffering from severe distortions. The resulting emissions were received at a higher frequency (612 or 1224 kHz subarray) to improve the resulting spatial resolution during passive beamforming. Normal resolution three-dimensional images were formed using a delay, sum, and integrate beamforming algorithm, and super-resolved images were extracted via Gaussian fitting of the estimated point-spread-function to the normal resolution data.

RESULTS

With super resolution techniques, the mean lateral (axial) full-width-at-half-maximum image intensity was 16 ± 3 (32 ± 6) μm, and 7 ± 1 (15 ± 2) μm corresponding to ∼1/67 of the normal resolution at 612 and 1224 kHz, respectively. The mean positional uncertainties were ∼1/350 (lateral) and ∼1/180 (axial) of the receive wavelength in water. In addition, a temporal correlation between nanodroplet vaporization and the transmit waveform shape was observed, which may provide the opportunity to enhance the signal-to-noise ratio in future studies.

CONCLUSIONS

Here, we demonstrate the feasibility of vaporizing nanodroplets via low frequency ultrasound and simultaneously performing spatial mapping via passive beamforming at higher frequencies to improve the resulting spatial resolution of super resolution imaging techniques. This method may enable complete four-dimensional vascular mapping in organs where a hemispherical array could be positioned to surround the target, such as the brain, breast, or testicles.

Real-Time Passive Acoustic Mapping With Enhanced Spatial Resolution in Neuronavigation-Guided Focused Ultrasound for Blood-Brain Barrier Opening

OBJECTIVE

Passive acoustic mapping (PAM) provides the spatial information of acoustic energy emitted from microbubbles during focused ultrasound (FUS), which can be used for safety and efficacy monitoring of blood-brain barrier (BBB) opening. In our previous work with a neuronavigation-guided FUS system, only part of the cavitation signal could be monitored in real time due to the computational burden although full-burst analysis is required to detect transient and stochastic cavitation activity. In addition, the spatial resolution of PAM can be limited for a small-aperture receiving array transducer. For full-burst real-time PAM with enhanced resolution, we developed a parallel processing scheme for coherence-factor-based PAM (CF-PAM) and implemented it onto the neuronavigation-guided FUS system using a co-axial phased-array imaging transducer.

METHODS

Simulation and in-vitro human skull studies were conducted for the performance evaluation of the proposed method in terms of spatial resolution and processing speed. We also carried out real-time cavitation mapping during BBB opening in non-human primates (NHPs).

RESULTS

CF-PAM with the proposed processing scheme provided better resolution than that of traditional time-exposure-acoustics PAM with a higher processing speed than that of eigenspace-based robust Capon beamformer, which facilitated the full-burst PAM with the integration time of 10 ms at a rate of 2 Hz. In vivo feasibility of PAM with the co-axial imaging transducer was also demonstrated in two NHPs, showing the advantages of using real-time B-mode and full-burst PAM for accurate targeting and safe treatment monitoring.

SIGNIFICANCE

This full-burst PAM with enhanced resolution will facilitate the clinical translation of online cavitation monitoring for safe and efficient BBB opening.

Passive Cavitation Imaging Artifact Reduction Using Data-Adaptive Spatial Filtering

Passive cavitation imaging (PCI) with a clinical diagnostic array results in poor axial localization of bubble activity due to the size of the point spread function (PSF). The objective of this study was to determine if data-adaptive spatial filtering improved PCI beamforming performance relative to standard frequency-domain delay, sum, and integrate (DSI) or robust Capon beamforming (RCB). The overall goal was to improve source localization and image quality without sacrificing computation time. Spatial filtering was achieved by applying a pixel-based mask to DSI- or RCB-beamformed images. The masks were derived from DSI, RCB, or phase or amplitude coherence factors (ACFs) using both receiver operating characteristic (ROC) and precision-recall (PR) curve analyses. Spatially filtered passive cavitation images were formed from cavitation emissions based on two simulated sources densities and four source distribution patterns mimicking cavitation emissions induced by an EkoSonic catheter. Beamforming performance was assessed via binary classifier metrics. The difference in sensitivity, specificity, and area under the ROC curve (AUROC) differed by no more than 11% across all algorithms for both source densities and all source patterns. The computational time required for each of the three spatially filtered DSIs was two orders of magnitude less than that required for time-domain RCB and thus this data-adaptive spatial filtering strategy for PCI beamforming is preferable given the similar binary classification performance.

Passive acoustic mapping with absolute time-of-flight information and delay-multiply-sum beamforming

BACKGROUND

Passive acoustic mapping (PAM) is showing increasing application potential in monitoring ultrasound therapy by spatially resolving cavitation activity. PAM with the relative time-of-flight information leads to poor axial resolution when implemented with ultrasound diagnostic transducers. Through utilizing the absolute time-of-flight information preserved by the transmit-receive synchronization and applying the common delay-sum (DS) beamforming algorithm, PAM axial resolution can be greatly improved in the short-pulse excitation scenario, as with active ultrasound imaging. However, PAM with the absolute time-of-flight information (referred as AtPAM) suffers from low imaging resolution and weak interference suppression when the DS algorithm is applied.

PURPOSE

This study aims to propose an enhanced AtPAM algorithm based on delay-multiply-sum (DMS) beamforming, to address the shortcomings of the DS-based AtPAM algorithm.

METHODS

In DMS beamforming, the element signals delayed by the absolute time delays are first processed with a signed square-root operation and then multiplied in pairs and finally summed, the resulting beamformed output is further band-pass filtered. The performances of DS- and DMS-based AtPAMs are compared by experiments, in which an ultrasound diagnostic transducer (a linear array) is employed to passively sense the wire signals generated by an unfocused ultrasound transducer and the cavitation signals generated by a focused therapeutic ultrasound transducer in a flow phantom. The AtPAM image quality is assessed by main-lobe width (MLW), intensity valley value (IVV), area of pixels (AOP), signal-to-interference ratio (SIR), and signal-to-noise ratio (SNR).

RESULTS

The single-wire experimental results show that compared to the DS algorithm, the DMS algorithm leads to an enhanced AtPAM image with a decreased transverse MLW of 0.15 mm and an improved SIR and SNR of 31.50 and 18.77 dB. For the four-wire images, the transverse (axial) IVV is decreased by 18.37 dB (13.11 dB) and the SIR (the SNR) is increased by 26.13 dB (18.47 dB) when using the DMS algorithm. The cavitation activity is better highlighted by DMS-based AtPAM, which decreases the AOP by 0.81 mm² (-10-dB level) and 4.43 mm² (-20-dB level) and increases the SIR and SNR by 20.14 and 10.48 dB respectively. The pixel distributions of AtPAM images of both wires and cavitation activity also indicate a better suppression of the DMS algorithm in sidelobe and noise.

CONCLUSIONS

The experimental results illustrate that the DMS algorithm can improve the image quality of AtPAM compared to the DS algorithm. DMS-based AtPAM is beneficial for detecting cavitation activity during short-pulse ultrasound exposure with high resolution, and further for monitoring short-pulse ultrasound therapy.

Reconstruction of thermoacoustic emission sources induced by proton irradiation using numerical time reversal

Objective.Mapping of dose delivery in proton beam therapy can potentially be performed by analyzing thermoacoustic emissions measured by ultrasound arrays. Here, a method is derived and demonstrated for spatial mapping of thermoacoustic sources using numerical time reversal, simulating re-transmission of measured emissions into the medium.Approach.Spatial distributions of thermoacoustic emission sources are shown to be approximated by the analytic-signal form of the time-reversed acoustic field, evaluated at the time of the initial proton pulse. Given calibration of the array sensitivity and knowledge of tissue properties, this approach approximately reconstructs the acoustic source amplitude, equal to the product of the time derivative of the radiation dose rate, mass density, and Grüneisen parameter. This approach was implemented using two models for acoustic fields of the array elements, one modeling elements as line sources and the other as rectangular radiators. Thermoacoustic source reconstructions employed previously reported measurements of emissions from proton energy deposition in tissue-mimicking phantoms. For a phantom incorporating a bone layer, reconstructions accounted for the higher sound speed in bone. Dependence of reconstruction quality on array aperture size and signal-to-noise ratio was consistent with previous acoustic simulation studies.Main results.Thermoacoustic source distributions were successfully reconstructed from acoustic emissions measured by a linear ultrasound array. Spatial resolution of reconstructions was significantly improved in the azimuthal (array) direction by incorporation of array element diffraction. Source localization agreed well with Monte Carlo simulations of energy deposition, and was improved by incorporating effects of inhomogeneous sound speed.Significance.The presented numerical time reversal approach reconstructs thermoacoustic sources from proton beam radiation, based on straightforward processing of acoustic emissions measured by ultrasound arrays. This approach may be useful for ranging and dosimetry of clinical proton beams, if acoustic emissions of sufficient amplitude and bandwidth can be generated by therapeutic proton sources.

Assessment of bubble activity generated by histotripsy combined with echogenic liposomes

Objective.Histotripsy is a form of focused ultrasound therapy that uses the mechanical activity of bubbles to ablate tissue. While histotripsy alone degrades the cellular content of tissue, recent studies have demonstrated it effectively disrupts the extracellular structure of pathologic conditions such as venous thrombosis when combined with a thrombolytic drug. Rather than relying on standard administration methods, associating thrombolytic drugs with an ultrasound-triggered echogenic liposome vesicle will enable targeted, systemic drug delivery. To date, histotripsy has primarily relied on nano-nuclei inherent to the medium for bubble cloud generation, and microbubbles associated with echogenic liposomes may alter the histotripsy bubble dynamics. The objective of this work was to investigate the interaction of histotripsy pulse with echogenic liposomes.Approach.Bubble clouds were generated using a focused source in anin vitromodel of venous flow. Acoustic emissions generated during the insonation were passively acquired to assess the mechanical activity of the bubble cloud. High frame rate, pulse inversion imaging was used to track the change in echogenicity of the liposomes following histotripsy exposure.Main results.For peak negative pressures less than 20 MPa, acoustic emissions indicative of stable and inertial bubble activity were observed. As the peak negative pressure of the histotripsy excitation increased, harmonics of the excitation were observed in OFP t-ELIP solutions and plasma alone. Additional observations with high frame rate imaging indicated a transition of bubble behavior as the pulse pressure transitioned to shock wave formation.Significance.These observations suggest that a complex interaction between histotripsy pulses and echogenic liposomes that may be exploited for combination treatment approaches.

A Novel Approach for the Detection of Every Significant Collapsing Bubble in Passive Cavitation Imaging

Passive cavitation image (PCI) shows the power distribution of the acoustic emissions resulting from cavitation bubble collapses. The conventional PCI convolves the emitted cavitation signals with the point spread function of an imaging system, and it suffers from a low spatial resolution and contrast due to the increased sidelobe artifacts accumulated during the temporal integral process. To overcome the problems, the present study considers a 3-D time history of instantaneous PCIs where cavitation occurs at the local maxima of the main lobes of the beamformed cavitation field surrounded by the sidelobes largely spreading out in a time-space domain. A spatial and temporal gating technique was employed to detect the local maxima so that cavitation bubbles can be identified with their collapsing strength. The proposed approach was verified by the simulation for single and multiple cavitation bubbles, proving that it accurately detects the location and strength of the collapsing bubbles. An experimental test was also carried out for the cavitation bubbles produced by a clinical extracorporeal shock wave therapeutic device, which underpins that the proposed method successfully identifies every individual cavitation bubble.

Experimental Demonstration of Trans-Skull Volumetric Passive Acoustic Mapping With the Heterogeneous Angular Spectrum Approach

Real-time, 3-D, passive acoustic mapping (PAM) of microbubble dynamics during transcranial focused ultrasound (FUS) is essential for optimal treatment outcomes. The angular spectrum approach (ASA) potentially offers a very efficient method to perform PAM, as it can reconstruct specific frequency bands pertinent to microbubble dynamics and may be extended to correct aberrations caused by the skull. Here, we experimentally assess the abilities of heterogeneous ASA (HASA) to perform trans-skull PAM. Our experimental investigations demonstrate that the 3-D PAMs of a known 1-MHz source, constructed with HASA through an ex vivo human skull segment, reduced both the localization error (from 4.7 ± 2.3 to 2.3 ± 1.6 mm) and the number, size, and energy of spurious lobes caused by aberration, with the modest additional computational expense. While further improvements in the localization errors are expected with arrays with denser elements and larger aperture, our analysis revealed that experimental constraints associated with the array pitch and aperture (here, 1.8 mm and 2.5 cm, respectively) can be ameliorated by interpolation and peak finding techniques. Beyond the array characteristics, our analysis also indicated that errors in the registration (translation and rotation of ±5 mm and ±5°, respectively) of the skull segment to the array can lead to peak localization errors of the order of a few wavelengths. Interestingly, errors in the spatially dependent speed of sound in the skull (±20%) caused only subwavelength errors in the reconstructions, suggesting that registration is the most important determinant of point source localization accuracy. Collectively, our findings show that HASA can address source localization problems through the skull efficiently and accurately under realistic conditions, thereby creating unique opportunities for imaging and controlling the microbubble dynamics in the brain.

On the threshold of 1/2 order subharmonic emissions in the oscillations of ultrasonically excited bubbles

The pressure threshold for 1/2 order subharmonic (SH) emissions and period doubling during the oscillations of ultrasonically excited bubbles is thought to be minimum when the bubble is sonicated with twice its resonance frequency (f_r). This estimate is based on studies that simplified or neglected the effects of thermal damping. In this work, the nonlinear dynamics of ultrasonically excited bubbles is investigated accounting for the thermal dissipation. Results are visualized using bifurcation diagrams as a function of pressure. Here we show that, and depending on the gas, the pressure threshold for 1/2 order SHs can be minimum at a frequency between 0.5f_r≤f≤0.6f_r. In this frequency range, the generation of 1/2 order SHs are due to the occurrence of 5/2 order ultra-harmonic resonance. The stability of such oscillations is size dependent. For an air bubble immersed in water, only bubbles bigger than 1 μm in diameter are able to emit non-destructive SHs in these frequency ranges.

Multifocused Ultrasound Therapy for Controlled Microvascular Permeabilization and Improved Drug Delivery

Focused ultrasound (FUS) exposure of micro-bubble (MB) contrast agents can transiently increase microvascular permeability allowing anticancer drugs to extravasate into a targeted tumor tissue. Either fixed or mechanically steered in space, most studies to date have used a single element focused transducer to deliver the ultrasound (US) energy. The goal of this study was to investigate various multi-FUS strategies implemented on a programmable US scanner (Vantage 256, Verasonics Inc.) equipped with a linear array for image guidance and a 128-element therapy transducer (HIFUPlex-06, Sonic Concepts). The multi-FUS strategies include multi-FUS with sequential excitation (multi-FUS-SE) and multi-FUS with temporal sequential excitation (multi-FUS-TSE) and were compared to single-FUS and sham treatment. This study was performed using athymic mice implanted with breast cancer cells ( N = 20 ). FUS therapy experiments were performed for 10 min after a solution containing MBs (Definity, Lantheus Medical Imaging Inc.) and near-infrared (NIR, surrogate drug) dye were injected via the tail vein. The fluorescent signal was monitored using an in vivo optical imaging system (Pearl Trilogy, LI-COR) to quantify intratumoral dye accumulation at baseline and again at 0.1, 24, and 48 h after receiving US therapy. Animals were then euthanized for ex vivo dye extraction analysis. At 48 h, fluorescent tracer accumulation within the tumor space for the multi-FUS-TSE therapy group animals was found to be 67.3%, 50.3%, and 36.2% higher when compared to sham, single-FUS, and multi-FUS-SE therapy group measures, respectively. Also, dye extraction and fluorescence measurements from excised tumor tissue found increases of 243.2%, 163.1%, and 68.1% for the multi-FUS-TSE group compared to sham, single-FUS, and multi-FUS-SE therapy group measures, respectively. In summary, experimental results revealed that for a multi-FUS sequence, increased microvascular permeability was considerably influenced by both the spatial and temporal aspects of the applied US therapy.

Nonlinear dynamics and bifurcation structure of ultrasonically excited lipid coated microbubbles

In many applications, microbubbles (MBs) are encapsulated by a lipid coating to increase their stability. However, the complex behavior of the lipid coating including buckling and rupture sophisticates the dynamics of the MBs and as a result the dynamics of the lipid coated MBs (LCMBs) are not well understood. Here, we investigate the nonlinear behavior of the LCMBs by analyzing their bifurcation structure as a function of acoustic pressure. We show that, the LC can enhance the generation of period 2 (P2), P3, higher order subharmonics (SH), superharmonics and chaos at very low excitation pressures (e.g. 1 kPa). For LCMBs sonicated by their SH resonance frequency and in line with experimental observations with increasing pressure, P2 oscillations exhibit three stages: generation at low acoustic pressures, disappearance and re-generation. Within non-destructive oscillation regimes and by pressure amplitude increase, LCMBs can also exhibit two saddle node (SN) bifurcations resulting in possible abrupt enhancement of the scattered pressure. The first SN resembles the pressure dependent resonance phenomenon in uncoated MBs and the second SN resembles the pressure dependent SH resonance. Depending on the initial surface tension of the LCMBs, the nonlinear behavior may also be suppressed for a wide range of excitation pressures.

Use of the Cross-Spectral Density Matrix for Enhanced Passive Ultrasound Imaging of Cavitation

Passive ultrasound imaging is of great interest for cavitation monitoring. Spatiotemporal monitoring of cavitation bubbles in therapeutic applications is possible using an ultrasound imaging probe to passively receive the acoustic signals from the bubbles. Fourier-domain (FD) beamformers have been proposed to process the signals received into maps of the spatial localization of cavitation activity, with reduced computing times with respect to the time-domain approach, and to take advantage of frequency selectivity for cavitation regime characterization. The approaches proposed have been mainly nonadaptive, and these have suffered from low resolution and contrast, due to the many reconstruction artifacts. Inspired by the array-processing literature and in the context of passive ultrasound imaging of cavitation, we propose here a robust estimation of the second-order statistics of data through spatial covariance matrices in the FD or cross-spectral density matrices (CSMs). The benefits of such formalism are illustrated using advanced reconstruction algorithms, such as the robust Capon beamformer, the Pisarenko class beamformer, and the multiple signal classification approach. Through both simulations and experiments in a water tank, we demonstrate that enhanced localization of cavitation activity (i.e., improved resolution and contrast with respect to nonadaptive approaches) is compatible with the rapid and frequency-selective approaches of the FD. Robust estimation of the CSM and the derived adaptive beamformers paves the way to the development of powerful passive ultrasound imaging tools.

Real-Time Passive Acoustic Mapping Using Sparse Matrix Multiplication

Passive acoustic mapping enables the spatiotemporal monitoring of cavitation with circulating microbubbles during focused ultrasound (FUS)-mediated blood-brain barrier opening. However, the computational load for processing large data sets of cavitation maps or more complex algorithms limit the visualization in real-time for treatment monitoring and adjustment. In this study, we implemented a graphical processing unit (GPU)-accelerated sparse matrix-based beamforming and time exposure acoustics in a neuronavigation-guided ultrasound system for real-time spatiotemporal monitoring of cavitation. The system performance was tested in silico through benchmarking, in vitro using nonhuman primate (NHP) and human skull specimens, and demonstrated in vivo in NHPs. We demonstrated the stability of the cavitation map for integration times longer than 62.5 [Formula: see text]. A compromise between real-time displaying and cavitation map quality obtained from beamformed RF data sets with a size of 2000 ×128 ×30 (axial [Formula: see text]) was achieved for an integration time of [Formula: see text], which required a computational time of 0.27 s (frame rate of 3.7 Hz) and could be displayed in real-time between pulses at PRF = 2 Hz. Our benchmarking tests show that the GPU sparse-matrix algorithm processed the RF data set at a computational rate of [Formula: see text]/pixel/sample, which enables adjusting the frame rate and the integration time as needed. The neuronavigation system with real-time implementation of cavitation mapping facilitated the localization of the cavitation activity and helped to identify distortions due to FUS phase aberration. The in vivo test of the method demonstrated the feasibility of GPU-accelerated sparse matrix computing in a close to a clinical condition, where focus distortions exemplify problems during treatment. These experimental conditions show the need for spatiotemporal monitoring of cavitation with real-time capability that enables the operator to correct or halt the sonication in case substantial aberrations are observed.

Imaging With Therapeutic Acoustic Wavelets-Short Pulses Enable Acoustic Localization When Time of Arrival is Combined With Delay and Sum

Passive acoustic mapping (PAM) is an algorithm that reconstructs the location of acoustic sources using an array of receivers. This technique can monitor therapeutic ultrasound procedures to confirm the spatial distribution and amount of microbubble activity induced. Current PAM algorithms have an excellent lateral resolution but have a poor axial resolution, making it difficult to distinguish acoustic sources within the ultrasound beams. With recent studies demonstrating that short-length and low-pressure pulses-acoustic wavelets-have the therapeutic function, we hypothesized that the axial resolution could be improved with a quasi-pulse-echo approach and that the resolution improvement would depend on the wavelet's pulse length. This article describes an algorithm that resolves acoustic sources axially using time of flight and laterally using delay-and-sum beamforming, which we named axial temporal position PAM (ATP-PAM). The algorithm accommodates a rapid short pulse (RaSP) sequence that can safely deliver drugs across the blood-brain barrier. We developed our algorithm with simulations (k-wave) and in vitro experiments for one-, two-, and five-cycle pulses, comparing our resolution against that of two current PAM algorithms. We then tested ATP-PAM in vivo and evaluated whether the reconstructed acoustic sources mapped to drug delivery within the brain. In simulations and in vitro, ATP-PAM had an improved resolution for all pulse lengths tested. In vivo, experiments in mice indicated that ATP-PAM could be used to target and monitor drug delivery into the brain. With acoustic wavelets and time of flight, ATP-PAM can locate acoustic sources with a vastly improved spatial resolution.

Doppler Passive Acoustic Mapping

In therapeutic ultrasound using microbubbles, it is essential to drive the microbubbles into the correct type of activity and the correct location to produce the desired biological response. Although passive acoustic mapping (PAM) is capable of locating where microbubble activities are generated, it is well known that microbubbles rapidly move within the ultrasound beam. We propose a technique that can image microbubble movement by estimating their velocities within the focal volume. Microbubbles embedded within a wall-less channel of a tissue-mimicking material were sonicated using 1-MHz focused ultrasound. The acoustic emissions generated by the microbubbles were captured with a linear array (L7-4). PAM with robust Capon beamforming was used to localize the microbubble acoustic emissions. We spectrally analyzed the time trace of each position and isolated the higher harmonics. Microbubble velocity maps were constructed from the position-dependent Doppler shifts at different time points during sonication. Microbubbles moved primarily away from the transducer at velocities on the order of 1 m/s due to primary acoustic radiation forces, producing a time-dependent velocity distribution. We detected microbubble motion both away and toward the receiving array, revealing the influence of acoustic radiation forces and fluid motion due to the ultrasound exposure. High-speed optical images confirmed the acoustically measured microbubble velocities. Doppler PAM enables passive estimation of microbubble motion and may be useful in therapeutic applications, such as drug delivery across the blood-brain barrier, sonoporation, sonothrombolysis, and drug release.

Acoustic source localization with the angular spectrum approach in continuously stratified media

The angular spectrum approach (ASA)-a frequency domain method to calculate the acoustic field-enables highly efficient passive source localization and modeling forward propagation in homogeneous media. If the medium is continuously stratified, a first-order analytical solution may be obtained for the field at arbitrary depth. Simulations show that the proposed stratified ASA solution enables accurate source localization as compared to the uncorrected ASA (error from 1.2 ± 0.3 to 0.49 ± 0.3 wavelengths) at scalings relevant to biomedical, underwater, and atmospheric acoustic applications, and requiring milliseconds on nonspecialized hardware. The results suggest the proposed correction enables efficient and accurate localization in stratified environments.

Critical corrections to models of nonlinear power dissipation of ultrasonically excited bubbles

Current models for calculating nonlinear power dissipation during the oscillations of acoustically excited bubbles generate non-physical values for the radiation damping (D_rd) term for some frequency and pressure regions that include near resonance oscillations. Moreover, the ratio of the dissipated powers significantly deviate from the values that are calculated by the linear model at low amplitude oscillations (acoustic excitation pressure of P_A=1 kPa and expansion ratio of <≊1.01). In high amplitude oscillation regimes (P_a⩾20 kPa), the dissipated power due to D_rd deviates largely from the dissipated power as calculated by the widely accepted approach that uses the scattered power by the bubbles. We provide critical corrections to the present models. The validity of the results was examined in regimes of low amplitude oscillations and high amplitude oscillations. In the low amplitude regime, the ratio of the dissipated powers as calculated by the current and proposed model were compared with the linear model predictions. At higher amplitude oscillations, the dissipated power by radiation loss as calculated by the current and the proposed models were compared with the dissipated power calculated using the scattered power by the bubbles. We show that non-physical values are absent in the proposed model. Moreover, predictions of the proposed approach are identical to the predictions of the linear model and the dissipated power estimated using the scattered pressure by the bubble. We show that damping due to thermal effects, liquid viscosity and radiation heavily depend on the excitation pressure and that the linear model estimations are not valid even at pressures as low as 20 kPa.

A new frequency domain passive acoustic mapping method using passive Hilbert beamforming to reduce the computational complexity of fast Fourier transform

Passive acoustic mapping (PAM) is the current state-of-the-art imaging tool for monitoring cavitation activity during focused ultrasound therapy such as blood-brain barrier opening. However, PAM incurs huge computational complexity. To address this issue, frequency-domain PAM (FD-PAM) was proposed. Nevertheless, FD-PAM still requires a large number of fast Fourier transforms (FFTs) to produce the frequency components utilized for cavitation monitoring with PAM. Hence, in this paper, we proposes a frequency domain PAM method using passive Hilbert beamforming (PHB-PAM), which can significantly reduce the number of input samples for FFT by down-sampling the analytic signal of the received RF samples at each channel at a rate equal to the bandwidth of the frequency components of interest. The experimental results show that the proposed PHB-PAM provides comparable image quality to that of FD-PAM (correlation coefficient > 0.98). Additionally, the study experimentally verifies that the pre-processing block for generating the decimated analytic signal and FFT in PHB-PAM can be realized using lesser logic resources than FFT in FD-PAM when implemented in an FPGA. Especially, with 128-fold decimation, PHB-PAM reduces the amount of LUTs and DSP slices to implement the pre-processing block by 72.16% and 53.4%, respectively, compared to those of FD-PAM, which allows the 64-channel implementation of the pre-processing block in a low-cost single FPGA. Finally, a hardware-efficient architecture for the pre-processing block of PHB-PAM is described, which can be implemented by replacing the two lowpass filters of an off-the-shelf analog front-end component for ultrasound imaging with a pair of band-pass filters. If PHB-PAM is realized using such a component, it can truly minimize the computational complexity of FD-PAM.

In Vitro Thrombolytic Efficacy of Single- and Five-Cycle Histotripsy Pulses and rt-PA

Although primarily known as an ablative modality, histotripsy can increase the efficacy of lytic therapy in a retracted venous clot model. Bubble cloud oscillations are the primary mechanism of action for histotripsy, and the type of bubble activity is dependent on the pulse duration. A retracted human venous clot model was perfused with and without the thrombolytic recombinant tissue plasminogen activator (rt-PA). The clot was exposed to histotripsy pulses of single- or five-cycle duration and peak negative pressures of 0-30 MPa. Bubble activity within the clot was monitored via passive cavitation imaging. The combination of histotripsy and rt-PA was more efficacious than rt-PA alone for single- and five-cycle pulses with peak negative pressures of 25 and 20 MPa, respectively. For both excitation schemes, the detected acoustic emissions correlated with the degree of thrombolytic efficacy. These results indicate that rt-PA and single- or multicycle histotripsy pulses enhance thrombolytic therapy.

A Clinical System for Non-invasive Blood-Brain Barrier Opening Using a Neuronavigation-Guided Single-Element Focused Ultrasound Transducer

Focused ultrasound (FUS)-mediated blood-brain barrier (BBB) opening is currently being investigated in clinical trials. Here, we describe a portable clinical system with a therapeutic transducer suitable for humans, which eliminates the need for in-line magnetic resonance imaging (MRI) guidance. A neuronavigation-guided 0.25-MHz single-element FUS transducer was developed for non-invasive clinical BBB opening. Numerical simulations and experiments were performed to determine the characteristics of the FUS beam within a human skull. We also validated the feasibility of BBB opening obtained with this system in two non-human primates using U.S. Food and Drug Administration (FDA)-approved treatment parameters. Ultrasound propagation through a human skull fragment caused 44.4 ± 1% pressure attenuation at a normal incidence angle, while the focal size decreased by 3.3 ± 1.4% and 3.9 ± 1.8% along the lateral and axial dimension, respectively. Measured lateral and axial shifts were 0.5 ± 0.4 mm and 2.1 ± 1.1 mm, while simulated shifts were 0.1 ± 0.2 mm and 6.1 ± 2.4 mm, respectively. A 1.5-MHz passive cavitation detector transcranially detected cavitation signals of Definity microbubbles flowing through a vessel-mimicking phantom. T₁-weighted MRI confirmed a 153 ± 5.5 mm³ BBB opening in two non-human primates at a mechanical index of 0.4, using Definity microbubbles at the FDA-approved dose for imaging applications, without edema or hemorrhage. In conclusion, we developed a portable system for non-invasive BBB opening in humans, which can be achieved at clinically relevant ultrasound exposures without the need for in-line MRI guidance. The proposed FUS system may accelerate the adoption of non-invasive FUS-mediated therapies due to its fast application, low cost and portability.

Collective nonlinear behavior of interacting polydisperse microbubble clusters

Acoustically excited microbubbles (MBs) have shown to exhibit rich dynamics, enabling them to be employed in various applications ranging from chemistry to medicine. Exploiting the full potential of MBs for applications requires a good understanding of their complex dynamics. Improved understanding of MB oscillations can lead to further enhancement in optimizing their efficacy in many applications and also invent new ones. Oscillating MBs have been shown to generate secondary pressure waves that modify the dynamics of the MBs in their proximity. A modified Keller-Miksis equation is used to account for inter-bubble interactions. The oscillatory dynamics of each MB within clusters was computed by numerically solving the resulting system of coupled nonlinear second order differential equations in potential fluid flow. Frequency response analysis and bifurcation diagrams were employed to track the dynamics of interacting MBs. We start with investigating the effect of inter-bubble interactions for cases of three and four MBs over a wide range of acoustic and geometric parameters. Emergent collective behavior was observed which are dominated by the dynamics of the largest MB within the cluster. The emergent dynamics of smaller MBs within clusters can be characterized by constructive and destructive inter-bubble interactions. In constructive interactions, the radial oscillations of smaller MBs matched those of the largest MB and their oscillations are amplified. In destructive interactions, the oscillations of smaller bubbles are suppressed so that their oscillations match those of the largest MB. Furthermore, a special case of constructive interactions is presented where dominant MB (largest) can force smaller MBs into period doubling and subharmonic oscillations. The collective behavior is further investigated in large MB cluster and it is shown that largest MBs, even in small numbers can force smaller ones into period doubling and subharmonic oscillations.

Delay multiply and sum beamforming method applied to enhance linear-array passive acoustic mapping of ultrasound cavitation

PURPOSE

Passive acoustic mapping (PAM) has been proposed as a means of monitoring ultrasound therapy, particularly nonthermal cavitation-mediated applications. In PAM, the most common beamforming algorithm is a delay, sum, and integrate (DSAI) approach. However, using DSAI leads to low-quality images for the case where a narrow-aperture receiving array such as a standard B-mode linear array is used. This study aims to propose an enhanced linear-array PAM algorithm based on delay, multiply, sum, and integrate (DMSAI).

METHODS

In the proposed algorithm, before summation, the delayed signals are combinatorially coupled and multiplied, which means that the beamformed output of the proposed algorithm is the spatial coherence of received acoustic emissions. We tested the performance of the proposed DMSAI using both simulated and experimental data and compared it with DSAI. The reconstructed cavitation images were evaluated quantitatively by using source location errors between the two algorithms, full width at half maximum (FWHM), size of point spread function (A₅₀ area), signal-to-noise ratio (SNR), and computational time.

RESULTS

The results of simulations and experiments for single cavitation source show that, by introducing DMSAI, the FWHM and the A₅₀ area are reduced and the SNR is improved compared with those obtained by DSAI. The simulation results for two symmetric or nonsymmetric cavitation sources and multiple cavitation sources show that DMSAI can significantly reduce the A₅₀ area and improve the SNR, therefore improving the detectability of multiple cavitation sources.

CONCLUSIONS

The results indicate that the proposed DMSAI algorithm outperforms the conventionally used DSAI algorithm. This work may have the potential of providing an appropriate method for ultrasound therapy monitoring.

A simple method to analyze the super-harmonic and ultra-harmonic behavior of the acoustically excited bubble oscillator

The bubble oscillator is a highly nonlinear system, which makes it difficult to generate a comprehensive understanding of its oscillatory behavior. One method used to investigate such complex dynamical systems is the bifurcation analysis. Numerous investigations have employed the method of bifurcation diagrams to study the effect of different control parameters on the bubble behavior. These studies, however, focused mainly on investigating the subharmonic (SH) and chaotic oscillations of the bubbles. Super-harmonic (SuH) and ultra-harmonic (UH) bubble oscillations remain under-investigated. One reason is that the conventional method used for generating bifurcation diagrams cannot reliably identify features that are responsible for the identification of SuH and UH oscillations. Additionally, the conventional method cannot distinguish between the UHs and SHs. We introduce a simple procedure to address this shortcoming. In this method, the maxima of the bubble oscillatory response were selected and plotted alongside the traditional bifurcation points for the corresponding control parameter. Results show that depending on the control parameters the conventional method or the method of maxima may miss intricate details of the oscillations. In order to have a comprehensive knowledge on the rich dynamics of the system, the two methods should be employed side by side. Through plotting the two bifurcation structures in tandem, the oscillatory behavior of the bubble was analyzed with more detail, and stable SuH and UH bubble oscillations were investigated. Based on this new analysis, the conditions for the generation and amplification of UH and SuH regimes are discussed.

Comparison study of passive acoustic mapping and high-speed photography for monitoring in situ cavitation bubbles

The spatiotemporal accuracy of passive acoustic mapping (PAM) for monitoring in situ cavitation bubbles has not been assessed directly via optical means. Here, the cavitation bubbles are monitored from two image sequences obtained simultaneously with PAM and high-speed photography (HSP). The temporal accuracy of PAM for detecting cavitation nucleation and the spatial resolution for cavitation localization are compared with those measured from HSP. The results show that PAM has a temporal accuracy of 20 μs. Mean differences in the spatial locations of PAM and HSP are as small as 10.0 and 30.5 μm along the lateral and axial directions, respectively.

Passive cavitation mapping using dual apodization with cross-correlation in ultrasound therapy monitoring

Recently, passive acoustic mapping (PAM) has been successfully applied for dynamic monitoring of ultrasound therapy by beamforming acoustic emissions of cavitation activity during ultrasound exposure. The most widely used PAM algorithm in the literature is time exposure acoustics (TEA), which is a standard delay, sum, and integrate algorithm. However, it results in large point spread function (PSF) and serious imaging artifacts for the case where a narrow-aperture receiving array such as a standard B-mode linear array is used, therefore degrading the quality of cavitation image. To address these challenges, in this paper, we proposed a novel PAM algorithm namely dual apodization with cross-correlation (DAX)-based TEA, in which DAX was originally used as a reconstruction algorithm in medical ultrasound imaging. In the proposed algorithm, two sets of signals were beamformed by two receive apodization functions with alternating elements enabled, and the cross-correlation coefficient of the two signals served as a weighting factor that would be multiplied to the sum of the two signals. The performance of the proposed algorithm was tested on simulated channel data obtained using a multi-bubble model, and experiments were also performed in an in vitro vessel phantom with flowing microbubbles as cavitation nuclei. The reconstructed cavitation images were evaluated quantitatively using established quality metrics including full width at half maximum (FWHM), A_-6dB area, and signal-to-noise ratio (SNR). The results suggested that the proposed algorithm significantly outperformed the conventionally used TEA algorithm. This work may have the potential of providing a useful tool for highly accurate localization of cavitation activity during ultrasound therapy.

Advances in acoustic monitoring and control of focused ultrasound-mediated increases in blood-brain barrier permeability

Transcranial focused ultrasound (FUS) combined with intravenously circulating microbubbles can transiently and selectively increase blood-brain barrier permeability to enable targeted drug delivery to the central nervous system, and is a technique that has the potential to revolutionize the way neurological diseases are managed in medical practice. Clinical testing of this approach is currently underway in patients with brain tumors, early Alzheimer's disease, and amyotrophic lateral sclerosis. A major challenge that needs to be addressed in order for widespread clinical adoption of FUS-mediated blood-brain barrier permeabilization to occur is the development of systems and methods for real-time treatment monitoring and control, to ensure that safe and effective acoustic exposure levels are maintained throughout the procedures. This review gives a basic overview of the oscillation dynamics, acoustic emissions, and biological effects associated with ultrasound-stimulated microbubbles in vivo, and provides a summary of recent advances in acoustic-based strategies for detecting, controlling, and mapping microbubble activity in the brain. Further development of next-generation clinical FUS brain devices tailored towards microbubble-mediated applications is warranted and required for translation of this potentially disruptive technology into routine clinical practice.

Visualization of the Intensity Field of a Focused Ultrasound Source In Situ

In an increasing number of applications of focused ultrasound (FUS) therapy, such as opening of the blood-brain barrier or collapsing microbubbles in a tumor, elevation of tissue temperature is not involved. In these cases, real-time visualization of the field distribution of the FUS source would allow localization of the FUS beam within the targeted tissue and allow repositioning of the FUS beam during tissue motion. In this paper, in order to visualize the FUS beam in situ, a 6-MHz single-element transducer ( f /2) was used as the FUS source and aligned perpendicular to a linear array which passively received scattered ultrasound from the sample. An image of the reconstructed intensity field pattern of the FUS source using bistatic beamforming was then superimposed on a registered B-mode image of the sample acquired using the same linear array. The superimposed image is used to provide anatomical context of the FUS beam in the sample being treated. The intensity field pattern reconstructed from a homogeneous scattering phantom was compared with the field characteristics of the FUS source characterized by the wire technique. The beamwidth estimates at the FUS focus using the in situ reconstruction technique and the wire technique were 1.5 and 1.2 mm, respectively. The depth-of-field estimates for the in situ reconstruction technique and the wire technique were 11.8 and 16.8 mm, respectively. The FUS beams were also visualized in a two-layer phantom and a chicken breast. The novel reconstruction technique was able to accurately visualize the field of an FUS source in the context of the interrogated medium.

2D and 3D real-time passive cavitation imaging of pulsed cavitation ultrasound therapy in moving tissues

Pulsed cavitation ultrasound therapy (PCUT) is an effective non-invasive therapeutic approach in various medical indications that relies on the mechanical effects generated by cavitation bubbles. Even though limited by the poor contrast, conventional ultrasound B-Mode imaging has been widely used for the guidance and monitoring of the therapeutic procedure, allowing the visualization of the cavitation bubble cloud. However, the visualization of the bubble cloud is often limited in deep organs such as the liver and the heart and remains moreover completely subjective for the operator. Our goal is to develop a new imaging mode to better identify the cavitation cloud. Active and passive cavitation imaging methods have been developed but none of them has been able to locate the cavitation bubble created by PCUT in real-time and in moving organs. In this paper we propose a passive ultrasound imaging approach combined with a spatiotemporal singular value decomposition filter to detect and map the bubble cloud with high sensitivity and high contrast. In moving applications at a maximal motion speed of 10 mm s^-1, the contrast-to-noise ratio for passive cavitation imaging is up to 10 times higher than for active cavitation imaging, with a temporal resolution of about 100 ms. The mapping of the bubble cloud can be overlaid in real-time to the conventional B-Mode, which permits to locate the cavitation phenomena in relation to the anatomic image. Finally, we extend the technique to volumetric imaging and show its feasibility on moving phantoms.

Weighting the Passive Acoustic Mapping Technique With the Phase Coherence Factor for Passive Ultrasound Imaging of Ultrasound-Induced Cavitation

Ultrasound (US) cavitation is currently being explored for low-invasive therapy techniques applied to a wide panel of pathologies. Because of the random behavior of cavitation, a real-time spatial monitoring system may be required. For this purpose, the US passive imaging techniques have been recently investigated. In particular, the passive acoustic mapping (PAM) beamforming method enables the reconstruction of cavitation activity maps by beamforming acoustic signals passively recorded by an array transducer. In this paper, an optimized version of PAM, PAM weighted with a phase coherence factor (PAM-PCF), is considered. A general validation process is developed including simulations on a point source and experiments on a wire. Furthermore, using a focused regulated US-induced cavitation generator, reproducible cavitation experiments are conducted in water and in agar gel. The spatial behavior of a bubble cavitation cloud is determined using the PAM-PCF beamforming method to localize the focal cavitation point in two perpendicular imaging planes.

Passive Acoustic Mapping Using Data-Adaptive Beamforming Based on Higher Order Statistics

Sources of nonlinear acoustic emissions, particularly those associated with cavitation activity, play a key role in the safety and efficacy of current and emerging therapeutic ultrasound applications, such as oncological drug delivery, blood-brain barrier opening, and histotripsy. Passive acoustic mapping (PAM) is the first technique to enable real-time and non-invasive imaging of cavitation activity during therapeutic ultrasound exposure, through the recording and passive beamforming of broadband acoustic emissions using an array of ultrasound detectors. Initial limitations in PAM spatial resolution led to the adoption of optimal data-adaptive beamforming algorithms, such as the robust capon beamformer (RCB), that provide improved interference suppression and calibration error mitigation compared to non-adaptive beamformers. However, such approaches are restricted by the assumption that the recorded signals have a Gaussian distribution. To overcome this limitation and further improve the source resolvability of PAM, we propose a new beamforming approach termed robust beamforming by linear programming (RLPB). Along with the variance, this optimization-based method uses higher-order-statistics of the recorded signals, making no prior assumption on the statistical distribution of the acoustic signals. The RLPB is found via numerical simulations to improve resolvability over time exposure acoustics and RCB. In vitro experimentation yielded improved resolvability with respect to the source-to-array distance on the order of 22% axially and 13% transversely relative to RCB, whilst successfully accounting for array calibration errors. The improved resolution and decreased dependence on accurate calibration of RLPB is expected to facilitate the clinical translation of PAM for diagnostic, including super-resolution, and therapeutic ultrasound applications.

Characterization of cavitation-radiated acoustic power using diffraction correction

A method is developed for compensating absolute pressure measurements made by a calibrated passive cavitation detector (PCD) to estimate the average acoustic power radiated from a region of interest (ROI) defined to encompass all cavitating bubbles. A diffraction correction factor for conversion of PCD-measured pressures to cavitation-radiated acoustic power per unit area or volume is derived as a simple analytic expression, accounting for position- and frequency-dependent PCD sensitivity. This approach can be applied to measurements made by any PCD without precise knowledge of the number, spatial, or temporal distribution of cavitating bubbles. The diffraction correction factor is validated in simulation for a wide range of ROI dimensions and frequencies. The correction factor is also applied to emission measurements obtained during in vitro ultrasound-enhanced sonophoresis experiments, allowing comparison of stable cavitation levels between therapeutic configurations with different source center frequencies. Results incorporating sonication at both 0.41 and 2.0 MHz indicate that increases in skin permeability correlate strongly with the acoustic power of subharmonic emissions radiated per unit skin area.

Closed Loop Spatial and Temporal Control of Cavitation Activity with Passive Acoustic Mapping

Ultrasonically actuated microbubble oscillations hold great promise for minimally invasive therapeutic interventions. While several preclinical studies have demonstrated the potential of this technology, real-time methods to control the amplitude and type of microbubble oscillations (stable vs inertial acoustic cavitation) and ensure that cavitation occurs within the targeted region are needed for their successful translation to the clinic. In this paper, we propose a real-time nonlinear state controller that uses specific frequency bands of the microbubble acoustic emissions (harmonic, ultra-harmonic, etc.) to control cavitation activity (observer states). To attain both spatial and temporal control of cavitation activity with high signal to noise ratio, we implement a controller using fast frequency-selective passive acoustic mapping (PAM) based on the angular spectrum approach. The controller includes safety states based on the recorded broadband signal level and is able to reduce sensing inaccuracies with the inclusion of multiple frequency bands. In its simplest implementation the controller uses the peak intensity of the passive acoustic maps, reconstructed using the 3rd harmonic (4.896 × 0.019 MHz) of the excitation frequency. Our results show that the proposed real-time nonlinear state controller based on PAM is able to reach the targeted level of observer state (harmonic emissions) in less than 6 seconds and remain within 10 % of tolerance for the duration of the experiment (45 seconds). Similar response was observed using the acoustic emissions from single element passive cavitation detection, albeit with higher susceptibility to background noise and lack of spatial information. Importantly, the proposed PAM-based controller was able to control cavitation activity with spatial selectivity when cavitation existed simultaneously in multiple regions. The robustness of the controller is demonstrated using a range of controller parameters, multiple observer states concurrently (harmonic, ultra-harmonic, and broadband), noise levels (°6 to 12 dB SNR), and bubble concentrations (0.3 to 180 × 103 bubbles per microliter). More research in this direction under preclinical and clinical conditions is warranted.

Design and Implementation of a Transmit/Receive Ultrasound Phased Array for Brain Applications

Focused ultrasound phased array systems have attracted increased attention for brain therapy applications. However, such systems currently lack a direct and real-time method to intraoperatively monitor ultrasound pressure distribution for securing treatment. This study proposes a dual-mode ultrasound phased array system design to support transmit/receive operations for concurrent ultrasound exposure and backscattered focal beam reconstruction through a spherically focused ultrasound array. A 256-channel ultrasound transmission system was used to transmit focused ultrasonic energy (full 256 channels), with an extended implementation of multiple-channel receiving function (up to 64 channels) using the same 256-channel ultrasound array. A coherent backscatter-received beam formation algorithm was implemented to map the point spread function (PSF) and focal beam distribution under a free-field/transcranial environment setup, with the backscattering generated from a strong scatterer (a point reflector or a microbubble-perfused tube) or a weakly scattered tissue-mimicking graphite phantom. Our results showed that PSF and focal beam can be successfully reconstructed and visualized in free-field conditions and can also be transcranially reconstructed following skull-induced aberration correction. In vivo experiments were conducted to demonstrate its capability to preoperatively and semiquantitatively map a focal beam to guide blood-brain barrier opening. The proposed system may have potential for real-time guidance of ultrasound brain intervention, and may facilitate the design of a dual-mode ultrasound phased array for brain therapeutic applications.

Pulse inversion enhances the passive mapping of microbubble-based ultrasound therapy

Therapeutic ultrasound combined with preformed circulating microbubbles has enabled non-invasive and targeted drug delivery into the brain, tumors, and blood clots. Monitoring the microbubble activity is essential for the success of such therapies; however, skull and tissues limit our ability to detect low acoustic signals. Here, we show that by emitting consecutive therapeutic pulses of inverse polarity, the sensitivity in the detection of weak bubble acoustic signals during blood-brain barrier opening is enhanced compared to therapeutic pulses of the same polarity. Synchronous passive mapping of the cavitation activity was conducted using delay-and-sum beamforming with absolute time delays, which offers superior spatial resolution compared to the existing asynchronous passive imaging techniques. Sonication with pulse inversion allowed filter-free suppression of the tissue signals by up to 8 dB in a tissue-mimicking phantom and by 7 dB in vivo, compared to exposure without pulse inversion, enabling enhanced passive mapping of microbubble activity. Both therapeutic schemes resulted in similar free-field microbubble activation in vitro and efficient blood-brain barrier opening in vivo.

Frequency-sum beamforming for passive cavitation imaging

Beamforming includes a variety of spatial filtering techniques that may be used for determining sound source locations from near-field sensor array recordings. For this scenario, beamforming resolution depends on the acoustic frequency, array geometry, and target location. Random scattering in the medium between the source and the array may degrade beamforming resolution with higher frequencies being more susceptible to degradation. The performance of frequency-sum (FS) beamforming for reducing such sensitivity to mild scattering while increasing resolution is reported here. FS beamforming was used with a data-dependent [minimum variance (MV)] or data-independent (delay-and-sum, DAS) weight vector to produce higher frequency information from lower frequency signal components via a quadratic product of complex signal amplitudes. The current findings and comparisons are based on simulations and passive cavitation imaging experiments using 3 MHz and 6 MHz emissions recorded by a 128-element linear array. FS beamforming results are compared to conventional DAS and MV beamforming using four metrics: point spread function (PSF) size, axial and lateral contrast, and computation time. FS beamforming produces a smaller PSF than conventional DAS beamforming with less computation time than MV beamforming in free space and mild scattering environments. However, it may fail when multiple unknown sound sources are present.

A dual-mode hemispherical sparse array for 3D passive acoustic mapping and skull localization within a clinical MRI guided focused ultrasound device

Previous work has demonstrated that passive acoustic imaging may be used alongside MRI for monitoring of focused ultrasound therapy. However, past implementations have generally made use of either linear arrays originally designed for diagnostic imaging or custom narrowband arrays specific to in-house therapeutic transducer designs, neither of which is fully compatible with clinical MR-guided focused ultrasound (MRgFUS) devices. Here we have designed an array which is suitable for use within an FDA-approved MR-guided transcranial focused ultrasound device, within the bore of a 3 Tesla clinical MRI scanner. The array is constructed from 5  ×  0.4 mm piezoceramic disc elements arranged in pseudorandom fashion on a low-profile laser-cut acrylic frame designed to fit between the therapeutic elements of a 230 kHz InSightec ExAblate 4000 transducer. By exploiting thickness and radial resonance modes of the piezo discs the array is capable of both B-mode imaging at 5 MHz for skull localization, as well as passive reception at the second harmonic of the therapy array for detection of cavitation and 3D passive acoustic imaging. In active mode, the array was able to perform B-mode imaging of a human skull, showing the outer skull surface with good qualitative agreement with MR imaging. Extension to 3D showed the array was able to locate the skull within  ±2 mm/2° of reference points derived from MRI, which could potentially allow registration of a patient to the therapy system without the expense of real-time MRI. In passive mode, the array was able to resolve a point source in 3D within a  ±10 mm region about each axis from the focus, detect cavitation (SNR ~ 12 dB) at burst lengths from 10 cycles to continuous wave, and produce 3D acoustic maps in a flow phantom. Finally, the array was used to detect and map cavitation associated with microbubble activity in the brain in nonhuman primates.

Power cavitation-guided blood-brain barrier opening with focused ultrasound and microbubbles

Image-guided monitoring of microbubble-based focused ultrasound (FUS) therapies relies on the accurate localization of FUS-stimulated microbubble activity (i.e. acoustic cavitation). Passive cavitation imaging with ultrasound arrays can achieve this, but with insufficient spatial resolution. In this study, we address this limitation and perform high-resolution monitoring of acoustic cavitation-mediated blood-brain barrier (BBB) opening with a new technique called power cavitation imaging. By synchronizing the FUS transmit and passive receive acquisition, high-resolution passive cavitation imaging was achieved by using delay and sum beamforming with absolute time delays. Since the axial image resolution is now dependent on the duration of the received acoustic cavitation emission, short pulses of FUS were used to limit its duration. Image sets were acquired at high-frame rates for calculation of power cavitation images analogous to power Doppler imaging. Power cavitation imaging displays the mean intensity of acoustic cavitation over time and was correlated with areas of acoustic cavitation-induced BBB opening. Power cavitation-guided BBB opening with FUS could constitute a standalone system that may not require MRI guidance during the procedure. The same technique can be used for other acoustic cavitation-based FUS therapies, for both safety and guidance.

Passive acoustic mapping of cavitation using eigenspace-based robust Capon beamformer in ultrasound therapy

Pulse-echo imaging technique can only play a role when high intensity focused ultrasound (HIFU) is turned off due to the interference between the primary HIFU signal and the transmission pulse. Passive acoustic mapping (PAM) has been proposed as a tool for true real-time monitoring of HIFU therapy. However, the most-used PAM algorithm based on time exposure acoustic (TEA) limits the quality of cavitation image. Recently, robust Capon beamformer (RCB) has been used in PAM to provide improved resolution and reduced artifacts over TEA-based PAM, but the presented results have not been satisfactory. In the present study, we applied an eigenspace-based RCB (EISRCB) method to further improve the PAM image quality. The optimal weighting vector of the proposed method was found by projecting the RCB weighting vector onto the desired vector subspace constructed from the eigenstructure of the covariance matrix. The performance of the proposed PAM was validated by both simulations and in vitro histotripsy experiments. The results suggested that the proposed PAM significantly outperformed the conventionally used TEA and RCB-based PAM. The comparison results between pulse-echo images of the residual bubbles and cavitation images showed the potential of our proposed PAM in accurate localization of cavitation activity during HIFU therapy.

Post Hoc Analysis of Passive Cavitation Imaging for Classification of Histotripsy-Induced Liquefaction in Vitro

Histotripsy utilizes focused ultrasound to generate bubble clouds for transcutaneous tissue liquefaction. Bubble activity maps are under development to provide image guidance and monitor treatment progress. The aim of this paper was to investigate the feasibility of using plane wave B-mode and passive cavitation images to be used as binary classifiers of histotripsy-induced liquefaction. Prostate tissue phantoms were exposed to histotripsy pulses over a range of pulse durations (5- ) and peak negative pressures (12-23 MPa). Acoustic emissions were recorded during the insonation and beamformed to form passive cavitation images. Plane wave B-mode images were acquired following the insonation to detect the hyperechoic bubble cloud. Phantom samples were sectioned and stained to delineate the liquefaction zone. Correlation between passive cavitation and plane wave B-mode images and the liquefaction zone was assessed using receiver operating characteristic (ROC) curve analysis. Liquefaction of the phantom was observed for all the insonation conditions. The area under the ROC (0.94 versus 0.82), accuracy (0.90 versus 0.83), and sensitivity (0.81 versus 0.49) was greater for passive cavitation images relative to B-mode images ( ) along the azimuth of the liquefaction zone. The specificity was greater than 0.9 for both imaging modalities. These results demonstrate a stronger correlation between histotripsy-induced liquefaction and passive cavitation imaging compared with the plane wave B-mode imaging, albeit with limited passive cavitation image range resolution.

Real-time monitoring of controllable cavitation erosion in a vessel phantom with passive acoustic mapping

Cavitation erosion in blood vessel plays an important role in ultrasound thrombolysis, drug delivery, and other clinical applications. The controllable superficial vessel erosion based on ultrasonic standing wave (USW) has been used to effectively prevent vessel ruptures and haemorrhages, and optical method is used to observe the experiments. But optical method can only work in transparent media. Compared with standard B-mode imaging, passive acoustic mapping (PAM) can monitor erosion in real time and has better sensitivity of cavitation detection. However, the conventionally used PAM has limitations in imaging resolution and artifacts. In this study, a unique PAM method that combined the robust Capon beamformer (RCB) with the sign coherence factor (SCF) was proposed to monitor the superficial vessel erosion in real time. The performance of the proposed method was validated by simulations. In vitro experiments showed that the lateral (axial) resolution of the proposed PAM was 2.31±0.51 (3.19±0.38) times higher than time exposure acoustics (TEA)-based PAM and 1.73±0.38 (1.76±0.48) times higher than RCB-based PAM, and the cavitation-to-artifact ratio (CAR) of the proposed PAM could be improved by 22.5±3.2dB and 7.1±1.2dB compared with TEA and RCB-based PAM. These results showed that the proposed PAM can precisely monitor the superficial vessel erosion and the erosion shift after USW modulation. This work may have the potential of developing a useful tool for precise spatial control and real-time monitoring of the superficial vessel erosion.

Passive Microlesion Detection and Mapping for Treatment of Hypertrophic Cardiomyopathy

Intermittent high intensity ultrasound pulses with circulating contrast agent microbubbles can induce scattered microlesions of potential value for myocardial reduction therapy. This paper presents an in vitro setup imitating the treatment for monitoring development. A preclinical imaging system with a single element transducer, synchronization and receive-only imaging transducer array has been implemented on a research platform. Contrast agent microbubbles pumped in a dialysis tubing setup were exposed to high intensity focused ultrasound at 1.0/3.5 MHz center frequencies. Polystyrene spheres were employed as linear scatterers compared to contrast agents for system transfer function equalization. A cavitation mapping technique was employed to spatially locate and depict microbubble activity during treatment. For high acoustic pressure amplitudes a 5 dB difference between contrast agent and solid spheres was observed and spatially mapped. The in-plane resolution was 4.5 mm for axial and 1.5 mm laterally. In the future, this cavitation detection scheme will be applied to monitor in vivo microlesioning in real-time.

Combined passive acoustic mapping and magnetic resonance thermometry for monitoring phase-shift nanoemulsion enhanced focused ultrasound therapy

Focused ultrasound (FUS) has the potential to enable precise, image-guided noninvasive surgery for the treatment of cancer in which tumors are identified and destroyed in a single integrated procedure. However, success of the method in highly vascular organs has been limited due to heat losses to perfusion, requiring development of techniques to locally enhance energy absorption and heating. In addition, FUS procedures are conventionally monitored using MRI, which provides excellent anatomical images and can map temperature, but is not capable of capturing the full gamut of available data such as the acoustic emissions generated during this inherently acoustically-driven procedure. Here, we employed phase-shift nanoemulsions (PSNE) embedded in tissue phantoms to promote cavitation and hence temperature rise induced by FUS. In addition, we incorporated passive acoustic mapping (PAM) alongside simultaneous MR thermometry in order to visualize both acoustic emissions and temperature rise, within the bore of a full scale clinical MRI scanner. Focal cavitation of PSNE could be resolved using PAM and resulted in accelerated heating and increased the maximum elevated temperature measured via MR thermometry compared to experiments without nanoemulsions. Over time, the simultaneously acquired acoustic and temperature maps show translation of the focus of activity towards the FUS transducer, and the magnitude of the increase in cavitation and focal shift both increased with nanoemulsion concentration. PAM results were well correlated with MRI thermometry and demonstrated greater sensitivity, with the ability to detect cavitation before enhanced heating was observed. The results suggest that PSNE could be beneficial for enhancement of thermal focused ultrasound therapies and that PAM could be a critical tool for monitoring this process.

Visualizing the movement of the contact between vocal folds during vibration by using array-based transmission ultrasonic glottography

For the purpose of noninvasively visualizing the dynamics of the contact between vibrating vocal fold medial surfaces, an ultrasonic imaging method which is referred to as array-based transmission ultrasonic glottography is proposed. An array of ultrasound transducers is used to detect the ultrasound wave transmitted from one side of the vocal folds to the other side through the small-sized contact between the vocal folds. A passive acoustic mapping method is employed to visualize and locate the contact. The results of the investigation using tissue-mimicking phantoms indicate that it is feasible to use the proposed method to visualize and locate the contact between soft tissues. Furthermore, the proposed method was used for investigating the movement of the contact between the vibrating vocal folds of excised canine larynges. The results indicate that the vertical movement of the contact can be visualized as a vertical movement of a high-intensity stripe in a series of images obtained by using the proposed method. Moreover, a visualization and analysis method, which is referred to as array-based ultrasonic kymography, is presented. The velocity of the vertical movement of the contact, which is estimated from the array-based ultrasonic kymogram, could reach 0.8 m/s during the vocal fold vibration.

Passive Acoustic Mapping with the Angular Spectrum Method

In the present proof of principle study, we evaluated the homogenous angular spectrum method for passive acoustic mapping (AS-PAM) of microbubble oscillations using simulated and experimental data. In the simulated data we assessed the ability of AS-PAM to form 3D maps of a single and multiple point sources. Then, in the two dimensional limit, we compared the 2D maps from AS-PAM with alternative frequency and time domain passive acoustic mapping (FD- and TD-PAM) approaches. Finally, we assessed the ability of AS-PAM to visualize microbubble activity in vivo with data obtained during 8 different experiments of FUS-induced blood-brain barrier disruption in 3 nonhuman primates, using a clinical MR-guided FUS system. Our in silico results demonstrate AS-PAM can be used to perform 3D passive acoustic mapping. 2D AS-PAM as compared to FD- PAM and TD-PAM is 10 and 200 times faster respectively and has similar sensitivity, resolution, and localization accuracy, even when the noise was 10-fold higher than the signal. In-vivo, the AS-PAM reconstructions of emissions at frequency bands pertinent to the different types of microbubble oscillations were also found to be more sensitive than TD-PAM. AS-PAM of harmonic-only components predicted safe blood-brain barrier disruption, whereas AS-PAM of broadband emissions correctly identified MR-evident tissue damage. The disparity (3.2 mm) in the location of the cavitation activity between the three methods was within their resolution limits. These data clearly demonstrate that AS-PAM is a sensitive and fast approach for PAM, thus providing a clinically relevant method to guide therapeutic ultrasound procedures.

Quantitative Frequency-Domain Passive Cavitation Imaging

Passive cavitation detection has been an instrumental technique for measuring cavitation dynamics, elucidating concomitant bioeffects, and guiding ultrasound therapies. Recently, techniques have been developed to create images of cavitation activity to provide investigators with a more complete set of information. These techniques use arrays to record and subsequently beamform received cavitation emissions, rather than processing emissions received on a single-element transducer. In this paper, the methods for performing frequency-domain delay, sum, and integrate passive imaging are outlined. The method can be applied to any passively acquired acoustic scattering or emissions, including cavitation emissions. To compare data across different systems, techniques for normalizing Fourier transformed data and converting the data to the acoustic energy received by the array are described. A discussion of hardware requirements and alternative imaging approaches is additionally outlined. Examples are provided in MATLAB.