Combined passive detection and ultrafast active imaging of cavitation events induced by short pulses of high-intensity ultrasound

Timing-synchronized passive ultrasound imaging of cavitation using eigenspace-based minimum variance beamforming and principal component analysis

BACKGROUND

Passive ultrasound imaging (PUI) allows to spatially resolve cavitation triggered during ultrasound irradiation, its application in therapeutic ultrasound has been gaining attention in recent years. The diffraction mode of the imaging transducer greatly limits the PUI axial resolution, which can be improved by transmit-receive synchronization and employment of delay sum beamforming (DSB) when transmitting short pulses, however, DSB yields poor performance in resolution and anti-interference.

PURPOSE

Inspired by adaptive beamforming and its low-complexity algorithm in active imaging field, this paper aims to develop an improved timing-synchronized PUI (TSPUI) algorithm for detection of short-pulse transmission-induced cavitation.

METHODS

The passive array data collected by timing synchronization is processed by minimum variance beamforming (MVB), whose weights are optimized by projection on the eigendecomposed signal subspace, that is, eigenspace-based MVB (EMVB), with the sum of the flight times on the transmitting and receiving paths as the delay. Applying principal component analysis (PCA) on the pre-collected MVB weight samples, a conversion matrix is constructed to allow the matrix inversion and eigendecomposition involved in weight calculation to be performed in a low dimension. The algorithm performance is confirmed by experiments, where a high-intensity focused ultrasound transducer and a linear-array transducer configured in a common parallel or vertical manner are employed for cavitation induction and cavitation imaging, and evaluated with the established indicators.

RESULTS

Reducing the eigenvalue threshold coefficient allows more sidelobes to be removed, and choosing an appropriate principal component number can reduce the time cost while guaranteeing the reconstruction quality. EMVB-PCA provides high resolution and anti-interference performance relative to DSB, with a reduction of over 60% in the point spread area and over 14 dB in the sidelobe and noise level, meanwhile, its time cost is considerably lower than EMVB, with a reduction of over 80%. Additionally, constructing the conversion matrix by simulation is feasible and valid, providing convenience for real imaging.

CONCLUSIONS

EMVB-PCA allows for high-quality TSPUI reconstruction of cavitation at a fast rate, providing an effective tool for detecting short-duration cavitation and further benefiting short-pulse therapeutic ultrasound applications.

Self-Sensing Cavitation Detection for Pulsed Cavitational Ultrasound Therapy

OBJECTIVES

Monitoring cavitation during ultrasound therapy is crucial for assessing the procedure safety and efficacy. This work aims to develop a self-sensing and low-complexity approach for robust cavitation detection in moving organs such as the heart.

METHODS

An analog-to-digital converter was connected onto one channel of the therapeutic transducer from a clinical system dedicated to cardiac therapy, allowing to record signals on a computer. Acquisition of successive echoes backscattered by the cavitation cloud on the therapeutic transducer was performed at a high repetition rate. Temporal variations of the backscattered echoes were analyzed with a Singular-Value Decomposition filter to discriminate signals associated to cavitation, based on its stochastic nature. Metrics were derived to classify the filtered backscattered echoes. Classification of raw backscattered echoes was also performed with a machine learning approach. The performances were evaluated on 155 in vitro acquisitions and 110 signals acquired in vivo during transthoracic cardiac ultrasound therapy on 3 swine.

RESULTS

Cavitation detection was achieved successfully in moving tissues with high signal to noise ratio in vitro (cSNR = 25±5) and in vivo (cSNR = 20±6) and outperformed conventional methods (cSNR = 11±6). Classification methods were validated with spectral analysis of hydrophone measurements. High accuracy was obtained using either the clutter filter-based method (accuracy of 1) or the neural network-based method (accuracy of 0.99).

CONCLUSION

Robust self-sensing cavitation detection was demonstrated to be possible with a clutter filter-based method and a machine learning approach.

SIGNIFICANCE

The self-sensing cavitation detection method enables robust, reliable and low complexity cavitation activity monitoring during ultrasound therapy.

Acoustic-pressure-driven ultrasonic activation of the mechanosensitive receptor RET and of cell proliferation in colonic tissue

Ultrasound generates both compressive and shear mechanical forces in soft tissues. However, the specific mechanisms by which these forces activate cellular processes remain unclear. Here we show that low-intensity focused ultrasound can activate the mechanosensitive RET signalling pathway. Specifically, in mouse colon tissues ex vivo and in vivo, focused ultrasound induced RET phosphorylation in colonic crypts cells, which correlated with markers of proliferation and stemness when using hours-long insonication. The activation of the RET pathway is non-thermal, is linearly related to acoustic pressure and is independent of radiation-force-induced shear strain in tissue. Our findings suggest that ultrasound could be used to regulate cell proliferation, particularly in the context of regenerative medicine, and highlight the importance of adhering to current ultrasound-safety regulations for medical imaging.

2D spatiotemporal passive cavitation imaging and evaluation during ultrasound thrombolysis based on diagnostic ultrasound platform

Acoustic cavitation plays a critical role in various biomedical applications. However, uncontrolled cavitation can lead to undesired damage to healthy tissues. Therefore, real-time monitoring and quantitative evaluation of cavitation dynamics is essential for understanding underlying mechanisms and optimizing ultrasound treatment efficiency and safety. The current research addressed the limitations of traditionally used cavitation detection methods by developing introduced an adaptive time-division multiplexing passive cavitation imaging (PCI) system integrated into a commercial diagnostic ultrasound platform. This new method combined real-time cavitation monitoring with B-mode imaging, allowing for simultaneous visualization of treatment progress and 2D quantitative evaluation of cavitation dosage within targeted area. An improved delay-and-sum (DAS) algorithm, optimized with a minimum variance (MV) beamformer, is utilized to minimize the side lobe effect and improve the axial resolution typically associated with PCI. In additional to visualize and quantitatively assess the cavitation activities generated under varied acoustic pressures and microbubble concentrations, this system was specifically applied to perform 2D cavitation evaluation for ultrasound thrombolysis mediated by different solutions, e.g., saline, nanodiamond (ND) and nitrogen-annealed nanodiamond (N-AND). This research aims to bridge the gap between laboratory-based research systems and real-time spatiotemporal cavitation evaluation demands in practical uses. Results indicate that this improved 2D cavitation monitoring and evaluation system could offer a useful tool for comprehensive evaluating cavitation-mediated effects (e.g., ultrasound thrombolysis), providing valuable insights into in-depth understanding of cavitation mechanisms and optimization of cavitation applications.

Dynamic mode decomposition based Doppler monitoring of de novo cavitation induced by pulsed HIFU: an in vivo feasibility study

Pulsed high-intensity focused ultrasound (pHIFU) has the capability to induce de novo cavitation bubbles, offering potential applications for enhancing drug delivery and modulating tissue microenvironments. However, imaging and monitoring these cavitation bubbles during the treatment presents a challenge due to their transient nature immediately following pHIFU pulses. A planewave bubble Doppler technique demonstrated its potential, yet this Doppler technique used conventional clutter filter that was originally designed for blood flow imaging. Our recent study introduced a new approach employing dynamic mode decomposition (DMD) to address this in an ex vivo setting. This study demonstrates the feasibility of the application of DMD for in vivo Doppler monitoring of the cavitation bubbles in porcine liver and identifies the candidate monitoring metrics for pHIFU treatment. We propose a fully automated bubble mode identification method using k-means clustering and an image contrast-based algorithm, leading to the generation of DMD-filtered bubble images and corresponding Doppler power maps after each HIFU pulse. These power Doppler maps are then correlated with the extent of tissue damage determined by histological analysis. The results indicate that DMD-enhanced power Doppler map can effectively visualize the bubble distribution with high contrast, and the Doppler power level correlates with the severity of tissue damage by cavitation. Further, the temporal characteristics of the bubble modes, specifically the decay rates derived from DMD, provide information of the bubble dissolution rate, which are correlated with tissue damage level-slower rates imply more severe tissue damage.

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.

The fundamental mechanisms of the Korotkoff sounds generation

Blood pressure measurement is the most widely performed clinical exam to predict mortality risk. The gold standard for its noninvasive assessment is the auscultatory method, which relies on listening to the so-called "Korotkoff sounds" in a stethoscope placed at the outlet of a pneumatic arm cuff. However, more than a century after their discovery, the origin of these sounds is still debated, which implies a number of clinical limitations. We imaged the Korotkoff sound generation in vivo at thousands of images per second using ultrafast ultrasound. We showed with both experience and theory that Korotkoff sounds are paradoxically not sound waves emerging from the brachial artery but rather shear vibrations conveyed in surrounding tissues by the nonlinear pulse wave propagation. When these shear vibrations reached the stethoscope, they were synchronous, correlated, and comparable in intensity with the Korotkoff sounds. Understanding this mechanism could ultimately improve blood pressure measurement and provide additional understanding of arterial mechanical properties.

Cavitation with low-energy sonication during focused ultrasound thalamotomy for a patient with tremor-dominant Parkinson's disease: a potential risk

We report a patient with tremor-dominant Parkinson's disease who had a mild cavitation bioeffect during magnetic resonance-guided focused ultrasound thalamotomy. During the aligning phase with low-energy sonication, cavitation caused mild dysarthria and paresthesia, prompting treatment cessation. At the same time, tremor and rigidity improved. MRI revealed extensive high-intensity lesions in the thalamus 1 day after the procedure followed by steroid infusion, which resulted in resolution of adverse events. Tremor and rigidity improved 1.5 years after the procedure. Although cavitation can relieve tremors and rigidity, it should be carefully monitored due to potential permanent adverse events by unpredictable and unknown behaviors.

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.

Bubble-Enhanced Mixing Induced by Standing Surface Acoustic Waves (SSAWs) in Microchannel

BAW-based micromixers usually achieve mixing enhancement with acoustic-induced bubbles, while SAW-based micromixers usually enhance mixing efficiency by varying the configuration of IDTs and microchannels. In this paper, bubble-enhanced acoustic mixing induced by standing surface acoustic waves (SSAWs) in a microchannel is proposed and experimentally demonstrated. Significant enhancement in the mixing efficiency was achieved after the bubbles were stimulated in our acoustofluidic microdevice. With an applied voltage of 5 V, 50 times amplified, the proposed mixing microdevice could achieve 90.8% mixing efficiency within 60 s at a flow rate of 240 μL/h. The bubbles were generated from acoustic cavitation assisted by the temperature increase resulting from the viscous absorption of acoustic energy. Our results also suggest that a temperature increase is harmful to microfluidic devices and temperature monitoring. Regulation is essential, especially in chemical and biological applications.

Slow-Flow Ultrasound Localization Microscopy Using Recondensation of Perfluoropentane Nanodroplets

Ultrasound localization microscopy (ULM) is an emerging, super-resolution imaging technique for detailed mapping of the microvascular structure and flow velocity via subwavelength localization and tracking of microbubbles. Because microbubbles rely on blood flow for movement throughout the vascular space, acquisition times can be long in the smallest, low-flow microvessels. In addition, detection of microbubbles in low-flow regions can be difficult because of minimal separation of microbubble signal from tissue. Nanoscale, phase-change contrast agents (PCCAs) have emerged as a switchable, intermittent or persisting contrast agent for ULM via acoustic droplet vaporization (ADV). Here, the focus is on characterizing the spatiotemporal contrast properties of less volatile perfluoropentane (PFP) PCCAs. The results indicate that at physiological temperature, nanoscale PFP PCCAs with diameters less than 100 nm disappear within microseconds after ADV with high-frequency ultrasound (16 MHz, 5- to 6-MPa peak negative pressure) and that nanoscale PFP PCCAs have an inherent deactivation mechanism via immediate recondensation after ADV. This "blinking" on-and-off contrast signal allowed separation of flow in an in vitro flow phantom, regardless of flow conditions, although with a need for some replenishment at very low flow conditions to maintain count rate. This blinking behavior allows for rapid spatial mapping in areas of low or no flow with ULM, but limits velocity tracking because there is no stable bubble formation with nanoscale PFP PCCAs.

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.

High-spatial-resolution, instantaneous passive cavitation imaging with temporal resolution in histotripsy: a simulation study

Purpose

In histotripsy, a shock wave is transmitted, and the resulting inertial bubble cavitation that disrupts tissue is used for treatment. Therefore, it is necessary to detect when cavitation occurs and track the position of cavitation occurrence using a new passive cavitation (PC) imaging method.

Methods

An integrated PC image, which is constructed by collecting the focused signals at all times, does not provide information on when cavitation occurs and has poor spatial resolution. To solve this problem, we constructed instantaneous PC images by applying delay and sum beamforming at instantaneous time instants. By calculating instantaneous PC images at all data acquisition times, the proposed method can detect cavitation when it occurs by using the property that when signals from the cavitation are focused, their amplitude becomes large, and it can obtain a high-resolution PC image by masking out side lobes in the vicinity of cavitation.

Results

Ultrasound image simulation confirmed that the proposed method has higher resolution than conventional integrated PC imaging and showed that it can determine the position and time of cavitation occurrence as well as the signal strength.

Conclusion

Since the proposed novel PC imaging method can detect each cavitation separately when the incidence of cavitations is low, it can be used to monitor the treatment process of shock wave therapy and histotripsy, in which cavitation is an important mechanism of treatment.

Passive Cavitation Detection With a Needle Hydrophone Array

Therapeutic ultrasound and microbubble technologies seek to drive systemically administered microbubbles into oscillations that safely manipulate tissue or release drugs. Such procedures often detect the unique acoustic emissions from microbubbles with the intention of using this feedback to control the microbubble activity. However, most sensor systems reported introduce distortions to the acoustic signal. Acoustic shockwaves, a key emission from microbubbles, are largely absent in reported recording, possibly due to the sensors being too large or too narrowband, or having strong phase distortions. Here, we built a sensor array that countered such limitations with small, broadband sensors and a low-phase distorting material. We built eight needle hydrophones with polyvinylidene fluoride (PVDF, diameter: 2 mm) then fit them into a 3-D-printed scaffold in a two-layered, staggered arrangement. Using this array, we monitored microbubbles exposed to therapeutically relevant ultrasound pulses (center frequency: 0.5 MHz, peak-rarefactional pressure: 130-597 kPa, pulselength: four cycles). Our tests revealed that the hydrophones were broadband with the best having a sensitivity of -224.8 dB ± 3.2 dB re 1 V/ μ Pa from 1 to 15 MHz. The array was able to capture shockwaves generated by microbubbles. The signal-to-noise ratio (SNR) of the array was approximately two times higher than individual hydrophones. Also, the array could localize microbubbles (-3 dB lateral resolution: 2.37 mm) and determine the cavitation threshold (between 161 and 254 kPa). Thus, the array accurately monitored and localized microbubble activities, and may be an important technological step toward better feedback control methods and safer and more effective treatments.

Time-Resolved Passive Cavitation Mapping Using the Transient Angular Spectrum Approach

Passive cavitation mapping (PCM), which generates images using bubble acoustic emission signals, has been increasingly used for monitoring and guiding focused ultrasound surgery (FUS). PCM can be used as an adjunct to magnetic resonance imaging to provide crucial information on the safety and efficacy of FUS. The most widely used algorithm for PCM is delay-and-sum (DAS). One of the major limitations of DAS is its suboptimal computational efficiency. Although frequency-domain DAS can partially resolve this issue, such an algorithm is not suitable for imaging the evolution of bubble activity in real time and for cases in which cavitation events occur asynchronously. This study investigates a transient angular spectrum (AS) approach for PCM. The working principle of this approach is to backpropagate the received signal to the domain of interest and reconstruct the spatial-temporal wavefield encoded with the bubble location and collapse time. The transient AS approach is validated using an in silico model and water bath experiments. It is found that the transient AS approach yields similar results to DAS, but it is one order of magnitude faster. The results obtained by this study suggest that the transient AS approach is promising for fast and accurate PCM.

Passive Cavitation Mapping by Cavitation Source Localization From Aperture-Domain Signals-Part I: Theory and Validation Through Simulations

Passive cavitation mapping (PCM) algorithms for diagnostic ultrasound arrays based on time exposure acoustics (TEA) exhibit poor axial resolution, which is in part due to the diffraction-limited point spread function of the imaging system and poor rejection by the delay-and-sum beamformer. In this article, we adapt a method for speed of sound estimation to be utilized as a cavitation source localization (CSL) approach. This method utilizes a hyperbolic fit to the arrival times of the cavitation signals in the aperture domain, and the coefficients of the fit are related to the position of the cavitation source. Wavefronts exhibiting poor fit to the hyperbolic function are corrected to yield improved source localization. We demonstrate through simulations that this method is capable of accurate estimation of the origin of coherent spherical waves radiating from cavitation/point sources. The average localization error from simulated microbubble sources was 0.12 ± 0.12mm ( 0.15 ± 0.14λ₀ for a 1.78-MHz transmit frequency). In simulations of two simultaneous cavitation sources, the proposed technique had an average localization error of 0.2mm ( 0.23λ₀ ), whereas conventional TEA had an average localization error of 0.81mm ( 0.97λ₀ ). The reconstructed PCM-CSL image showed a significant improvement in resolution compared with the PCM-TEA approach.

Real-Time Visualization of a Focused Ultrasound Beam Using Ultrasonic Backscatter

Focused ultrasound (FUS) therapies induce therapeutic effects in localized tissues using either temperature elevations or mechanical stresses caused by an ultrasound wave. During an FUS therapy, it is crucial to continuously monitor the position of the FUS beam in order to correct for tissue motion and keep the focus within the target region. Toward the goal of achieving real-time monitoring for FUS therapies, we have developed a method for the real-time visualization of an FUS beam using ultrasonic backscatter. The intensity field of an FUS beam was reconstructed using backscatter from an FUS pulse received by an imaging array and then overlaid onto a B-mode image captured using the same imaging array. The FUS beam visualization allows one to monitor the position and extent of the FUS beam in the context of the surrounding medium. Variations in the scattering properties of the medium were corrected in the FUS beam reconstruction by normalizing based on the echogenicity of the coaligned B-mode image. On average, normalizing by echogenicity reduced the mean square error between FUS beam reconstructions in nonhomogeneous regions of a phantom and baseline homogeneous regions by 21.61. FUS beam visualizations were achieved, using a single diagnostic imaging array as both an FUS source and an imaging probe, in a tissue-mimicking phantom and a rat tumor in vivo with a frame rate of 25-30 frames/s.

Cerebrospinal Fluid Cavitation as a Mechanism of Blast-Induced Traumatic Brain Injury: A Review of Current Debates, Methods, and Findings

Cavitation has gained popularity in recent years as a potential mechanism of blast-induced traumatic brain injury (bTBI). This review presents the most prominent debates on cavitation; how bubbles can form or exist within the cerebrospinal fluid (CSF) and brain vasculature, potential mechanisms of cellular, and tissue level damage following the collapse of bubbles in response to local pressure fluctuations, and a survey of experimental and computational models used to address cavitation research questions. Due to the broad and varied nature of cavitation research, this review attempts to provide a necessary synthesis of cavitation findings relevant to bTBI, and identifies key areas where additional work is required. Fundamental questions about the viability and likelihood of CSF cavitation during blast remain, despite a variety of research regarding potential injury pathways. Much of the existing literature on bTBI evaluates cavitation based off its prima facie plausibility, while more rigorous evaluation of its likelihood becomes increasingly necessary. This review assesses the validity of some of the common assumptions in cavitation research, as well as highlighting outstanding questions that are essential in future work.

Hydrophone Spatial Averaging Correction for Acoustic Exposure Measurements From Arrays-Part I: Theory and Impact on Diagnostic Safety Indexes

This article reports underestimation of mechanical index (MI) and nonscanned thermal index for bone near focus (TIB) due to hydrophone spatial averaging effects that occur during acoustic output measurements for clinical linear and phased arrays. TIB is the appropriate version of thermal index (TI) for fetal imaging after ten weeks from the last menstrual period according to the American Institute of Ultrasound in Medicine (AIUM). Spatial averaging is particularly troublesome for highly focused beams and nonlinear, nonscanned modes such as acoustic radiation force impulse (ARFI) and pulsed Doppler. MI and variants of TI (e.g., TIB), which are displayed in real-time during imaging, are often not corrected for hydrophone spatial averaging because a standardized method for doing so does not exist for linear and phased arrays. A novel analytic inverse-filter method to correct for spatial averaging for pressure waves from linear and phased arrays is derived in this article (Part I) and experimentally validated in a companion article (Part II). A simulation was developed to estimate potential spatial-averaging errors for typical clinical ultrasound imaging systems based on the theoretical inverse filter and specifications for 124 scanner/transducer combinations from the U.S. Food and Drug Administration (FDA) 510(k) database from 2015 to 2019. Specifications included center frequency, aperture size, acoustic output parameters, hydrophone geometrical sensitive element diameter, etc. Correction for hydrophone spatial averaging using the inverse filter suggests that maximally achievable values for MI, TIB, thermal dose ( t ₄₃ ), and spatial-peak-temporal-average intensity ( [Formula: see text]) for typical clinical systems are potentially higher than uncorrected values by (means ± standard deviations) 9% ± 4% (ARFI MI), 19% ± 15% (ARFI TIB), 50% ± 41% (ARFI t ₄₃ ), 43% ± 39% (ARFI [Formula: see text]), 7% ± 5% (pulsed Doppler MI), 15% ± 11% (pulsed Doppler TIB), 42% ± 31% (pulsed Doppler t ₄₃ ), and 33% ± 27% (pulsed Doppler [Formula: see text]). These values correspond to frequencies of 3.2 ± 1.3 (ARFI) and 4.1 ± 1.4 MHz (pulsed Doppler), and the model predicts that they would increase with frequency. Inverse filtering for hydrophone spatial averaging significantly improves the accuracy of estimates of MI, TIB, t ₄₃ , and [Formula: see text] for ARFI and pulsed Doppler signals.

Tri-modality cavitation mapping in shock wave lithotripsy

Shock wave lithotripsy (SWL) has been widely used for non-invasive treatment of kidney stones. Cavitation plays an important role in stone fragmentation, yet it may also contribute to renal injury during SWL. It is therefore crucial to determine the spatiotemporal distributions of cavitation activities to maximize stone fragmentation while minimizing tissue injury. Traditional cavitation detection methods include high-speed optical imaging, active cavitation mapping (ACM), and passive cavitation mapping (PCM). While each of the three methods provides unique information about the dynamics of the bubbles, PCM has most practical applications in biological tissues. To image the dynamics of cavitation bubble collapse, we previously developed a sliding-window PCM (SW-PCM) method to identify each bubble collapse with high temporal and spatial resolution. In this work, to further validate and optimize the SW-PCM method, we have developed tri-modality cavitation imaging that includes three-dimensional high-speed optical imaging, ACM, and PCM seamlessly integrated in a single system. Using the tri-modality system, we imaged and analyzed laser-induced single cavitation bubbles in both free field and constricted space and shock wave-induced cavitation clusters. Collectively, our results have demonstrated the high reliability and spatial-temporal accuracy of the SW-PCM approach, which paves the way for the future in vivo applications on large animals and humans in SWL.

Critical instability at moving keyhole tip generates porosity in laser melting

Laser powder bed fusion is a dominant metal 3D printing technology. However, porosity defects remain a challenge for fatigue-sensitive applications. Some porosity is associated with deep and narrow vapor depressions called keyholes, which occur under high-power, low-scan speed laser melting conditions. High-speed x-ray imaging enables operando observation of the detailed formation process of pores in Ti-6Al-4V caused by a critical instability at the keyhole tip. We found that the boundary of the keyhole porosity regime in power-velocity space is sharp and smooth, varying only slightly between the bare plate and powder bed. The critical keyhole instability generates acoustic waves in the melt pool that provide additional yet vital driving force for the pores near the keyhole tip to move away from the keyhole and become trapped as defects.

Evaluation of the Uncertainty of Passive Cavitation Measurements for Blood-Brain Barrier Disruption Monitoring

Exposure to ultrasound combined with intravenous injection of microbubbles is a technique that can be used to temporarily disrupt the blood-brain barrier. Transcranial monitoring of cavitation can be done with one or more passive cavitation detectors (PCDs). However, the positioning of the PCDs relative to the cavitation site and the attenuation of these signals by the skull are two sources of error in the quantification of cavitation activity. The aim of this study was to evaluate in vitro the amplitude variation of cavitation signals that can be expected for an excised porcine skull model. The variation caused by the relative positioning of the PCD with respect to the cavitation site was quantified. A position-based correction of the signal amplitude was evaluated. Pig skull samples were used to assess variation in signal amplitude caused by bone. The overall coefficient of variation of the signals owing to these measurement biases was estimated at 30.8%.

From Anatomy to Functional and Molecular Biomarker Imaging and Therapy: Ultrasound Is Safe, Ultrafast, Portable, and Inexpensive

Ultrasound is the most widely used medical imaging modality worldwide. It is abundant, extremely safe, portable, and inexpensive. In this review, we consider some of the current development trends for ultrasound imaging, which build upon its current strength and the popularity it experiences among medical imaging professional users.Ultrasound has rapidly expanded beyond traditional radiology departments and cardiology practices. Computing power and data processing capabilities of commonly available electronics put ultrasound systems in a lab coat pocket or on a user's mobile phone. Taking advantage of new contributions and discoveries in ultrasound physics, signal processing algorithms, and electronics, the performance of ultrasound systems and transducers have progressed in terms of them becoming smaller, with higher imaging performance, and having lower cost. Ultrasound operates in real time, now at ultrafast speeds; kilohertz frame rates are already achieved by many systems.Ultrasound has progressed beyond anatomical imaging and monitoring blood flow in large vessels. With clinical approval of ultrasound contrast agents (gas-filled microbubbles) that are administered in the bloodstream, tissue perfusion studies are now routine. Through the use of modern ultrasound pulse sequences, individual microbubbles, with subpicogram mass, can be detected and observed in real time, many centimeters deep in the body. Ultrasound imaging has broken the wavelength barrier; by tracking positions of microbubbles within the vasculature, superresolution imaging has been made possible. Ultrasound can now trace the smallest vessels and capillaries, and obtain blood velocity data in those vessels.Molecular ultrasound imaging has now moved closer to clinic; the use of microbubbles with a specific affinity to endothelial biomarkers allows selective accumulation and retention of ultrasound contrast in the areas of ischemic injury, inflammation, or neoangiogenesis. This will aid in noninvasive molecular imaging and may provide additional help with real-time guidance of biopsy, surgery, and ablation procedures.The ultrasound field can be tightly focused inside the body, many centimeters deep, with millimeter precision, and ablate lesions by energy deposition, with thermal or mechanical bioeffects. Some of such treatments are already in clinical use, with more indications progressing through the clinical trial stage. In conjunction with intravascular microbubbles, focused ultrasound can be used for tissue-specific drug delivery; localized triggered release of sequestered drugs from particles in the bloodstream may take time to get to clinic. A combination of intravascular microbubbles with circulating drug and low-power ultrasound allows transient opening of vascular endothelial barriers, including blood-brain barrier; this approach has reached clinical trial stage. Therefore, the drugs that normally would not be getting to the target tissue in the brain will now have an opportunity to produce therapeutic efficacy.Overall, medical ultrasound is developing at a brisk rate, even in an environment where other imaging modalities are also advancing rapidly and may be considered more lucrative. With all the current advances that we discuss, and many more to come, ultrasound may help solve many problems that modern medicine is facing.

Focused Ultrasound for Noninvasive, Focal Pharmacologic Neurointervention

A long-standing goal of translational neuroscience is the ability to noninvasively deliver therapeutic agents to specific brain regions with high spatiotemporal resolution. Focused ultrasound (FUS) is an emerging technology that can noninvasively deliver energy up the order of 1 kW/cm² with millimeter and millisecond resolution to any point in the human brain with Food and Drug Administration-approved hardware. Although FUS is clinically utilized primarily for focal ablation in conditions such as essential tremor, recent breakthroughs have enabled the use of FUS for drug delivery at lower intensities (i.e., tens of watts per square centimeter) without ablation of the tissue. In this review, we present strategies for image-guided FUS-mediated pharmacologic neurointerventions. First, we discuss blood–brain barrier opening to deliver therapeutic agents of a variety of sizes to the central nervous system. We then describe the use of ultrasound-sensitive nanoparticles to noninvasively deliver small molecules to millimeter-sized structures including superficial cortical regions and deep gray matter regions within the brain without the need for blood–brain barrier opening. We also consider the safety and potential complications of these techniques, with attention to temporal acuity. Finally, we close with a discussion of different methods for mapping the ultrasound field within the brain and describe future avenues of research in ultrasound-targeted drug therapies.

Physical triggering strategies for drug delivery

Physically triggered systems hold promise for improving drug delivery by enhancing the controllability of drug accumulation and release, lowering non-specific toxicity, and facilitating clinical translation. Several external physical stimuli including ultrasound, light, electric fields and magnetic fields have been used to control drug delivery and they share some common features such as spatial targeting, spatiotemporal control, and minimal invasiveness. At the same time, they possess several distinctive features in terms of interactions with biological entities and/or the extent of stimulus response. Here, we review the key advances of such systems with a focus on discussing their physical mechanisms, the design rationales, and translational challenges.

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.

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.

Fast qualitative two-dimensional mapping of ultrasound fields with acoustic cavitation-enhanced ultrasound imaging

Characterization of ultrasound fields is a routine procedure for both diagnostic and therapeutic ultrasound. Quantitative field mapping with a calibrated hydrophone and multi-axis positioning system can be difficult and time consuming. In this study, the use of acoustic cavitation field mapping as a qualitative surrogate to acoustic pressure field mapping, albeit without acoustic pressure values is demonstrated. This technique allows for fast qualitative mapping of ultrasound fields and thereby functionality of the corresponding transducers, in a matter of seconds. In addition, this technique could be used to rapidly image in vivo acoustic cavitation fields during therapeutic ultrasound applications.

In vitro assessment of stiffness-dependent histotripsy bubble cloud activity in gel phantoms and blood clots

As a bubble-based ablative therapy, the efficacy of histotripsy has been demonstrated in healthy or acutely diseased models. Chronic conditions associated with stiff tissues may require additional bubble activity prior to histotripsy liquefaction. In this study, histotripsy pulses were generated in agarose phantoms of Young's moduli ranging from 12.3 to 142 kPa, and in vitro clot models with mild and strong platelet-activated retraction. Bubble cloud emissions were tracked with passive cavitation imaging, and the threshold acoustic power associated with phantom liquefaction was extracted with receiver operator characteristic analysis. The power of histotripsy-generated emissions and the degree of liquefaction were tabulated for both clot models. For the agarose phantoms, the acoustic power associated with liquefaction increased with Young's modulus. When grouped based on agarose concentration, only two arms displayed a significant difference in the liquefaction threshold acoustic power (22.1 kPa versus 142 kPa Young's modulus). The bubble cloud dynamics tracked with passive cavitation imaging indicated no strong changes in the bubble dynamics based on the phantom stiffness. For identical histotripsy exposure, the power of acoustic emissions and degree of clot lysis did not vary based on the clot model. Overall, these results indicate that a fixed threshold acoustic power mapped with passive cavitation imaging can be utilized for predicting histotripsy liquefaction over a wide range of tissue stiffness.

Compensation of array lens effects for improved co-registration of passive acoustic mapping and B-mode images for cavitation monitoring

Passive acoustic mapping (PAM) techniques offer a simple means of spatio-temporal cavitation monitoring during therapeutic ultrasound procedures. Implementation with a conventional diagnostic ultrasound system allows natural integration of PAM with B-mode imaging. However, the refracting properties of diagnostic array lenses may introduce PAM image registration errors that could lead to inaccuracies in treatment monitoring and guidance. To address these concerns, this paper presents lens characterization of two different array designs, analytical estimation of lens-induced source mapping errors in simple media, and experimental demonstration and correction of lens effects, reducing the depth-averaged image co-registration errors to no more than 0.52 mm.

Compensation of array lens effects for improved co-registration of passive acoustic mapping and B-mode images for cavitation monitoring

Passive acoustic mapping (PAM) techniques offer a simple means of spatio-temporal cavitation monitoring during therapeutic ultrasound procedures. Implementation with a conventional diagnostic ultrasound system allows natural integration of PAM with B-mode imaging. However, the refracting properties of diagnostic array lenses may introduce PAM image registration errors that could lead to inaccuracies in treatment monitoring and guidance. To address these concerns, this paper presents lens characterization of two different array designs, analytical estimation of lens-induced source mapping errors in simple media, and experimental demonstration and correction of lens effects, reducing the depth-averaged image co-registration errors to no more than 0.52 mm.

Observation and modulation of the dissolution of histotripsy-induced bubble clouds with high-frame rate plane wave imaging

Focused ultrasound therapies are a noninvasive means to ablate tissue. Histotripsy utilizes short ultrasound pulses with sufficient tension to nucleate bubble clouds that impart lethal strain to the surrounding tissues. Tracking bubble cloud dissolution between the application of histotripsy pulses is critical to ensure treatment efficacy. In this study, plane wave B-mode imaging was employed to monitor bubble cloud motion and grayscale at frame rates up to 11.25 kHz. Minimal changes in the area or position of the bubble clouds were observed 50 ms post excitation. The bubble cloud grayscale was observed to decrease with the square root of time, indicating a diffusion-driven process. These results were qualitatively consistent with an analytic model of gas diffusion during the histotripsy process. Finally, the rate of bubble cloud dissolution was found to be dependent on the output of the imaging pulse, indicating an interaction between the bubble cloud and imaging parameters. Overall, these results highlight the utility of plane wave B-mode imaging for monitoring histotripsy bubble clouds.

SVD-Based Separation of Stable and Inertial Cavitation Signals Applied to Passive Cavitation Mapping During HIFU

Detection of inertial and stable cavitation is important for guiding high-intensity focused ultrasound (HIFU). Acoustic transducers can passively detect broadband noise from inertial cavitation and the scattering of HIFU harmonics from stable cavitation bubbles. Conventional approaches to cavitation noise diagnostics typically involve computing the Fourier transform of the time-domain noise signal, applying a custom comb filter to isolate the frequency components of interest, followed by an inverse Fourier transform. We present an alternative technique based on singular value decomposition (SVD) that efficiently separates the broadband emissions and HIFU harmonics. Spatiotemporally resolved cavitation detection was achieved using a 128-element, 5-MHz linear-array ultrasound imaging system operating in the receive mode at 15 frames/s. A 1.1-MHz transducer delivered HIFU to tissue-mimicking phantoms and excised liver tissue for a duration of 5 s. Beamformed radio frequency signals corresponding to each scan line in a frame were assembled into a matrix, and SVD was performed. Spectra of the singular vectors obtained from a tissue-mimicking gel phantom were analyzed by computing the peak ratio ( R ), defined as the ratio of the peak of its fifth-order polynomial fit and the maximum spectral peak. Singular vectors that produced an were classified as those representing stable cavitation, i.e., predominantly containing harmonics of HIFU. The projection of data onto this singular base reproduced stable cavitation signals. Similarly, singular vectors that produced an were classified as those predominantly containing broadband noise associated with inertial cavitation. These singular vectors were used to isolate the inertial cavitation signal. The R -value thresholds determined using gel data were then employed to analyze cavitation data obtained from bovine liver ex vivo. The SVD-based method faithfully reproduced the structural details in the spatiotemporal cavitation maps produced using the more cumbersome comb-filter approach with a maximum root-mean-squared error of 10%.

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.

Closed-loop cavitation control for focused ultrasound-mediated blood-brain barrier opening by long-circulating microbubbles

Focused ultrasound (FUS) exposure in the presence of microbubbles (MBs) has been successfully used in the delivery of various sizes of therapeutic molecules across the blood-brain barrier (BBB). While acoustic pressure is correlated with the BBB opening size, real-time control of BBB opening to avoid vascular and neural damage is still a challenge. This arises mainly from the variability of FUS-MB interactions due to the variations of animal-specific metabolic environment and specific experimental setup. In this study, we demonstrate a closed-loop cavitation control framework to induce BBB opening for delivering large therapeutic molecules without causing macro tissue damages. To this end, we performed in mice long-term (5 min) cavitation monitoring facilitated by using long-circulating MBs. Monitoring the long-term temporal kinetics of the MBs under varying level of FUS pressure allowed to identify in situ, animal specific activity regimes forming pressure-dependent activity bands. This enables to determine the boundaries of each activity band (i.e. steady oscillation, transition, inertial cavitation) independent from the physical and physiological dynamics of the experiment. However, such a calibration approach is time consuming and to speed up characterization of the in situ, animal specific FUS-MB dynamics, we tested a novel method called 'pre-calibration' that closely reproduces the results of long-term monitoring but with a much shorter duration. Once the activity bands are determined from the pre-calibration method, an operation band can be selected around the desired cavitation dose. To drive cavitation in the selected operation band, we developed an adaptive, closed-loop controller that updates the acoustic pressure between each sonication based on measured cavitation dose. Finally, we quantitatively assessed the safety of different activity bands and validated the proposed methods and controller framework. The proposed framework serves to optimize the FUS pressure instantly to maintain the targeted cavitation level while improving safety control.

MRI-guided transurethral insonation of silica-shell phase-shift emulsions in the prostate with an advanced navigation platform

PURPOSE

In this study, the efficacy of transurethral prostate ablation in the presence of silica-shell ultrasound-triggered phase-shift emulsions (sUPEs) doped with MR contrast was evaluated. The influence of sUPEs on MR imaging assessment of the ablation zone was also investigated.

METHODS

sUPEs were doped with a magnetic resonance (MR) contrast agent, Gd2 O3 , to assess ultrasound transition. Injections of saline (sham), saline and sUPEs alone, and saline and sUPEs with Optison microbubbles were performed under guidance of a prototype interventional MRI navigation platform in a healthy canine prostate. Treatment arms were evaluated for differences in lesion size, T1  contrast, and temperature. In addition, non-perfused areas (NPAs) on dynamic contrast-enhanced (DCE) MRI, 55°C isotherms, and areas of 240 cumulative equivalent minutes at 43°C (CEM43 ) dose or greater computed from MR thermometry were measured and correlated with ablated areas indicated by histology.

RESULTS

For treatment arms including sUPEs, the computed correlation coefficients between the histological ablation zone and the NPA, 55°C isotherm, and 240 CEM43 area ranged from 0.96-0.99, 0.98-0.99, and 0.91-0.99, respectively. In the absence of sUPEs, the computed correlation coefficients between the histological ablation zone and the NPA, 55°C isotherm, and 240 CEM43 area were 0.69, 0.54, and 0.50, respectively. Across all treatment arms, the areas of thermal tissue damage and NPAs were not significantly different (P = 0.47). Areas denoted by 55°C isotherms and 240 CEM43 dose boundaries were significantly larger than the areas of thermal damage, again for all treatment arms (P = 0.009 and 0.003, respectively). No significant differences in lesion size, T1 contrast, or temperature were observed between any of the treatment arms (P > 0.0167). Lesions exhibiting thermal fixation on histological analysis were present in six of nine insonations involving sUPE injections and one of five insonations involving saline sham injections. Significantly larger areas (P = 0.002), higher temperatures (P = 0.004), and more frequent ring patterns of restricted diffusion on ex vivo diffusion-weighted imaging (P = 0.005) were apparent in lesions with thermal fixation.

CONCLUSIONS

T1 contrast suggesting sUPE transition was not evident in sUPE treatment arms. The use of MR imaging metrics to predict prostate ablation was not diminished by the presence of sUPEs. Lesions generated in the presence of sUPEs exhibited more frequent thermal fixation, though there were no significant changes in the ablation areas when comparing arms with and without sUPEs. Thermal fixation corresponded to some qualitative imaging features.

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.

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.

Coalescence of residual histotripsy cavitation nuclei using low-gain regions of the therapy beam during electronic focal steering

Following collapse of a histotripsy cloud, residual microbubbles may persist for seconds, distributed throughout the focus. Their presence can attenuate and scatter subsequent pulses, hindering treatment speed and homogeneity. Previous studies have demonstrated use of separate low-amplitude (~1 MPa) pulses interleaved with histotripsy pulses to drive bubble coalescence (BC), significantly improving treatment speed without sacrificing homogeneity. We propose that by using electronic focal steering (EFS) to direct the therapy focus throughout specially-designed EFS sequences, it is possible to use low-gain regions of the therapy beam to accomplish BC during EFS without any additional acoustic sequence. First, to establish proof of principle for an isolated focus, a 50-foci EFS sequence was constructed with the first position isolated near the geometric focus and remaining positions distributed post-focally. EFS sequences were evaluated in tissue-mimicking phantoms with gas concentrations of 20% and 100% with respect to saturation. Results using an isolated focus demonstrated that at 20% gas concentration, 49 EFS pulses were sufficient to achieve BC in all samples for pulse repetition frequency (PRF)  ⩽  800 Hz and 84.1%  ±  3.0% of samples at 5 kHz PRF. For phantoms prepared with 100% gas concentration, BC was achieved by 49 EFS pulses in 39.2%  ±  4.7% of samples at 50 Hz PRF and 63.4%  ±  15.3% of samples at 5 kHz. To show feasibility of using the EFS-BC method to ablate a large volume quickly, a 1000-foci EFS sequence covering a volume of approximately 27 ml was tested. Results indicate that the BC effect was similarly present. A treatment rate of 27  ±  6 ml min^-1 was achieved, which is signficantly faster than standard histotripsy and ultrasound thermal ablation. This study demonstrates that histotripsy with EFS can achieve BC without employing a separate acoustic sequence which has the potential to accelerate large-volume ablation while minimizing energy deposition.

Broadband Ultrasonic Attenuation Estimation and Compensation With Passive Acoustic Mapping

Several active and passive techniques have been developed to detect, localize, and quantify cavitation activity during therapeutic ultrasound procedures. Much of the prior cavitation monitoring research has been conducted using lossless in vitro systems or small animal models in which path attenuation effects were minimal. However, the performance of these techniques may be substantially degraded by attenuation between the internal therapeutic target and the external monitoring system. As a further step toward clinical application of passive acoustic mapping (PAM), this paper presents methods for attenuation estimation and compensation based on broadband cavitation data measured with a linear ultrasound array. Soft tissue phantom experiment results are used to illustrate: 1) the impact of realistic attenuation on PAM; 2) the possibility of estimating attenuation from cavitation data; 3) cavitation source energy estimation following attenuation compensation; and 4) the impact of sound speed uncertainty on PAM-related processing. Cavitation-based estimates of attenuation were within 1.5%-6.2% of the values found from conventional through-transmission measurements. Tissue phantom attenuation reduced the PAM energy estimate by an order of magnitude, but array data compensation using the cavitation-based attenuation spectrum enabled recovery of the PAM energy estimate to within 2.9%-5.9% of the values computed in the absence of the phantom. Sound speed uncertainties were found to modestly impact attenuation-compensated PAM energies, inducing errors no larger than 28% for a 40-m/s path-averaged speed error. Together, the results indicate the potential to significantly enhance the quantitative capabilities of PAM for ensuring therapeutic safety and efficacy.

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.

Ultrasound Localization Microscopy and Super-Resolution: A State of the Art

Because it drives the compromise between resolution and penetration, the diffraction limit has long represented an unreachable summit to conquer in ultrasound imaging. Within a few years after the introduction of optical localization microscopy, we proposed its acoustic alter ego that exploits the micrometric localization of microbubble contrast agents to reconstruct the finest vessels in the body in-depth. Various groups now working on the subject are optimizing the localization precision, microbubble separation, acquisition time, tracking, and velocimetry to improve the capacity of ultrasound localization microscopy (ULM) to detect and distinguish vessels much smaller than the wavelength. It has since been used in vivo in the brain, the kidney, and tumors. In the clinic, ULM is bound to improve drastically our vision of the microvasculature, which could revolutionize the diagnosis of cancer, arteriosclerosis, stroke, and diabetes.

Inside/outside the brain binary cavitation localization based on the lowpass filter effect of the skull on the harmonic content: a proof of concept study

Cavitation activity induced by ultrasound may occur during high intensity focused ultrasound (HIFU) treatment, due to bubble nucleation under high peak negative pressure, and during blood-brain-barrier (BBB) disruption, due to injected ultrasound contrast agents (UCAs). Such microbubble activity has to be monitored to assess the safety and efficiency of ultrasonic brain treatments. In this study, we aim at assessing whether cavitation occurs within cerebral tissue by binary discriminating cavitation activity originating from the inside or the outside of the skull. The results were obtained from both in vitro experiments mimicking BBB opening, by using UCA flow, and in vitro thermal necrosis in calf brain samples. The sonication was applied using a 1 MHz focused transducer and the acoustic response of the microbubbles was recorded with a wideband passive cavitation detector. The spectral content of the recorded signal was used to localize microbubble activity. Since the skull acts as a low pass filter, the ratio of high harmonics to low harmonics is lower for cavitation events located inside the skull compared to events outside the skull. Experiments showed that the ratio of the 5/2 ultraharmonic to the 1/2 subharmonic for binary localization cavitation activity achieves 100% sensitivity and specificity for both monkey and human skulls. The harmonic ratio of the fourth to the second harmonic provided 100% sensitivity and 96% and 46% specificity on a non-human primate for thermal necrosis and BBB opening, respectively. Nonetheless, the harmonic ratio remains promising for human applications, as the experiments showed 100% sensitivity and 100% specificity for both thermal necrosis and BBB opening through the human skull. The study requires further validation on a larger number of skull samples.

Three-dimensional transcranial microbubble imaging for guiding volumetric ultrasound-mediated blood-brain barrier opening

Focused ultrasound (FUS)-mediated blood-brain barrier (BBB) opening recently entered clinical testing for targeted drug delivery to the brain. Sources of variability exist in the current procedures, motivating the development of real-time monitoring and control techniques to improve treatment safety and efficacy. Here we used three-dimensional (3D) transcranial microbubble imaging to calibrate FUS exposure levels for volumetric BBB opening.

Methods: Using a sparse hemispherical transmit/receive ultrasound phased array, pulsed ultrasound was focused transcranially into the thalamus of rabbits during microbubble infusion and multi-channel 3D beamforming was performed online with receiver signals captured at the subharmonic frequency. Pressures were increased pulse-by-pulse until subharmonic activity was detected on acoustic imaging (p_sub), and tissue volumes surrounding the calibration point were exposed at 50-100%p_sub via rapid electronic beam steering.

Results: Spatially-coherent subharmonic microbubble activity was successfully reconstructed transcranially in vivo during calibration sonications. Multi-point exposures induced volumetric regions of elevated BBB permeability assessed via contrast-enhanced magnetic resonance imaging (MRI). At exposure levels ≥75%p_sub, MRI and histological examination occasionally revealed tissue damage, whereas sonications at 50%p_sub were performed safely. Substantial intra-grid variability of FUS-induced bioeffects was observed via MRI, prompting future development of multi-point calibration schemes for improved treatment consistency. Receiver array sparsity and sensor configuration had substantial impacts on subharmonic detection sensitivity, and are factors that should be considered when designing next-generation clinical FUS brain therapy systems.

Conclusion: Our findings suggest that 3D subharmonic imaging can be used to calibrate exposure levels for safe FUS-induced volumetric BBB opening, and should be explored further as a method for cavitation-mediated treatment guidance.

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.

Using the cavitation collapse time to indicate the extent of histotripsy-induced tissue fractionation

Histotripsy is an ultrasonic tissue ablation method based on acoustic cavitation. It has been shown that cavitation dynamics change depending on the mechanical properties of the host medium. During histotripsy treatment, the target-tissue is gradually fractionated and eventually liquefied to acellular homogenate. In this study, the change in the collapse time (t _col) of the cavitation bubble cloud over the course of histotripsy treatment is investigated as an indicator for progression of the tissue fractionation process throughout treatment. A 500 kHz histotripsy transducer is used to generate single-location lesions within tissue-mimicking agar phantoms of varying stiffness levels as well as ex vivo bovine liver samples. Cavitation collapse signals are acquired with broadband hydrophones, and cavitation is imaged optically using a high-speed camera in transparent tissue-mimicking phantoms. The high-speed-camera-acquired measurements of t _col validate the acoustic hydrophone measurements. Increases in t _col are observed both with decreasing phantom stiffness and throughout histotripsy treatment with increasing number of pulses applied. The increasing trend of t _col throughout the histotripsy treatment correlates well with the progression of lesion formation generated in tissue-mimicking phantoms (R ²  =  0.87). Finally, the increasing trend of t _col over the histotripsy treatment is validated in ex vivo bovine liver.