The histotripsy spectrum: differences and similarities in techniques and instrumentation.

Since its inception about two decades ago, histotripsy - a non-thermal mechanical tissue ablation technique - has evolved into a spectrum of methods, each with distinct potentiating physical mechanisms: intrinsic threshold histotripsy, shock-scattering histotripsy, hybrid histotripsy, and boiling histotripsy. All methods utilize short, high-amplitude pulses of focused ultrasound delivered at a low duty cycle, and all involve excitation of violent bubble activity and acoustic streaming at the focus to fractionate tissue down to the subcellular level. The main differences are in pulse duration, which spans microseconds to milliseconds, and ultrasound waveform shape and corresponding peak acoustic pressures required to achieve the desired type of bubble activity. In addition, most types of histotripsy rely on the presence of high-amplitude shocks that develop in the pressure profile at the focus due to nonlinear propagation effects. Those requirements, in turn, dictate aspects of the instrument design, both in terms of driving electronics, transducer dimensions and intensity limitations at surface, shape (primarily, the F-number) and frequency. The combination of the optimized instrumentation and the bio-effects from bubble activity and streaming on different tissues, lead to target clinical applications for each histotripsy method. Here, the differences and similarities in the physical mechanisms and resulting bioeffects of each method are reviewed and tied to optimal instrumentation and clinical applications.

Noninvasive acoustic manipulation of objects in a living body.

In certain medical applications, transmitting an ultrasound beam through the skin to manipulate a solid object within the human body would be beneficial. Such applications include, for example, controlling an ingestible camera or expelling a kidney stone. In this paper, ultrasound beams of specific shapes were designed by numerical modeling and produced using a phased array. These beams were shown to levitate and electronically steer solid objects (3-mm-diameter glass spheres), along preprogrammed paths, in a water bath, and in the urinary bladders of live pigs. Deviation from the intended path was on average <10%. No injury was found on the bladder wall or intervening tissue.

Maximizing mechanical stress in small urinary stones during burst wave lithotripsy.

Unlike shock wave lithotripsy, burst wave lithotripsy (BWL) uses tone bursts, consisting of many periods of a sinusoidal wave. In this work, an analytical theoretical approach to modeling mechanical stresses in a spherical stone was developed to assess the dependence of frequency and stone size on stress generated in the stone. The analytical model for spherical stones is compared against a finite-difference model used to calculate stress in nonspherical stones. It is shown that at low frequencies, when the wavelength is much greater than the diameter of the stone, the maximum principal stress is approximately equal to the pressure amplitude of the incident wave. With increasing frequency, when the diameter of the stone begins to exceed about half the wavelength in the surrounding liquid (the exact condition depends on the material of the stone), the maximum stress increases and can be more than six times greater than the incident pressure. These results suggest that the BWL frequency should be elevated for small stones to improve the likelihood and rate of fragmentation.

Holographic extraction of plane waves from an ultrasound beam for acoustic characterization of an absorbing layer of finite dimensions.

For the acoustic characterization of materials, a method is proposed for interpreting experiments with finite-sized transducers and test samples in terms of the idealized situation in which plane waves are transmitted through an infinite plane-parallel layer. The method uses acoustic holography, which experimentally provides complete knowledge of the wave field by recording pressure waveforms at points on a surface intersected by the acoustic beam. The measured hologram makes it possible to calculate the angular spectrum of the beam to decompose the field into a superposition of plane waves propagating in different directions. Because these waves cancel one another outside the beam, the idealized geometry of an infinite layer can be represented by a sample of finite size if its lateral dimensions exceed the width of the acoustic beam. The proposed method relies on holograms that represent the acoustic beam with and without the test sample in the transmission path. The method is described theoretically, and its capabilities are demonstrated experimentally for silicone rubber samples by measuring their frequency-dependent phase velocities and absorption coefficients in the megahertz frequency range.

A review on B/A measurement methods with a clinical perspective.

The nonlinear parameter of ultrasound B/A has shown to be a useful diagnostic parameter, reflecting medium content, structure, and temperature. Despite its recognized values, B/A is not yet used as a diagnostic tool in the clinic due to the limitations of current measurement and imaging techniques. This review presents an extensive and comprehensive overview of the techniques developed for B/A measurement of liquid and liquid-like media (e.g., tissue), identifying the methods that are most promising from a clinical perspective. This work summarizes the progress made in the field and the typical challenges on the way to B/A estimation. Limitations and problems with the current techniques are identified, suggesting directions that may lead to further improvement. Since the basic theory of the physics behind the measurement strategies is presented, it is also suited for a reader who is new to nonlinear ultrasound.

An investigation of elastic waves producing stone fracture in burst wave lithotripsy.

Burst wave lithotripsy is a method to noninvasively fragment urinary stones by short pulses of focused ultrasound. In this study, physical mechanisms of stone fracture during burst wave lithotripsy were investigated. Photoelasticity imaging was used to visualize elastic wave propagation in model stones and compare results to numerical calculations. Epoxy and glass stone models were made into rectangular, cylindrical, or irregular geometries and exposed in a degassed water bath to focused ultrasound bursts at different frequencies. A high-speed camera was used to record images of the stone during exposure through a circular polariscope backlit by a monochromatic flash source. Imaging showed the development of periodic stresses in the stone body with a pattern dependent on frequency. These patterns were identified as guided wave modes in cylinders and plates, which formed standing waves upon reflection from the distal surfaces of the stone model, producing specific locations of stress concentration in the models. Measured phase velocities compared favorably to numerically calculated modes dependent on frequency and material. Artificial stones exposed to bursts produced cracks at positions anticipated by this mechanism. These results support guided wave generation and reflection as a mechanism of stone fracture in burst wave lithotripsy.

Modeling of photoelastic imaging of mechanical stresses in transparent solids mimicking kidney stones.

Theoretical and numerical models were developed to calculate the polariscopic integrated light intensity that forms a projection of the dynamic stress within an axisymmetric elastic object. Although the model is general, this paper addressed its application to measurements of stresses in model kidney stones from a burst wave lithotripter for stone fragmentation. The stress was calculated using linear elastic equations, and the light propagation was modeled in the instantaneous case by integrating over the volume of the stone. The numerical model was written in finite differences. The resulting images agreed well with measured images. The measured images corresponded to the maximum shear stress distribution, although other stresses were also plotted. Comparison of the modeled and observed polariscope images enabled refinement of the photoelastic constant by minimizing the error between the calculated and measured fields. These results enable quantification of the stress within the polariscope images, determination of material properties, and the modes and mechanisms of stress production within a kidney stone. Such a model may help in interpreting elastic waves in structures, such as stones, toward improving lithotripsy procedures.

Bilayer aberration-inducing gel phantom for high intensity focused ultrasound applications.

Aberrations induced by soft tissue inhomogeneities often complicate high-intensity focused ultrasound (HIFU) therapies. In this work, a bilayer phantom made from polyvinyl alcohol hydrogel and ballistic gel was built to mimic alternating layers of water-based and lipid tissues characteristic of an abdominal body wall and to reproducibly distort HIFU fields. The density, sound speed, and attenuation coefficient of each material were measured using a homogeneous gel layer. A surface with random topographical features was designed as an interface between gel layers using a 2D Fourier spectrum approach and replicating different spatial scales of tissue inhomogeneities. Distortion of the field of a 256-element 1.5 MHz HIFU array by the phantom was characterized through hydrophone measurements for linear and nonlinear beam focusing and compared to the corresponding distortion induced by an ex vivo porcine body wall of the same thickness. Both spatial shift and widening of the focal lobe were observed, as well as dramatic reduction in focal pressures caused by aberrations. The results suggest that the phantom produced levels of aberration that are similar to a real body wall and can serve as a research tool for studying HIFU effects as well as for developing algorithms for aberration correction.

Simulation of nonlinear trans-skull focusing and formation of shocks in brain using a fully populated ultrasound array with aberration correction.

Multi-element high-intensity focused ultrasound phased arrays in the shape of hemispheres are currently used in clinics for thermal lesioning in deep brain structures. Certain side effects of overheating non-targeted tissues and skull bones have been revealed. Here, an approach is developed to mitigate these effects. A specific design of a fully populated 256-element 1-MHz array shaped as a spherical segment (F-number, F# = 1) and filled by randomly distributed equal-area polygonal elements is proposed. Capability of the array to generate high-amplitude shock fronts at the focus is tested in simulations by combining three numerical algorithms for linear and nonlinear field modeling and aberration correction. The algorithms are based on the combination of the Rayleigh integral, a linear pseudo-spectral time domain Kelvin-Voigt model, and nonlinear Westervelt model to account for the effects of inhomogeneities, aberrations, reflections, absorption, nonlinearity, and shear waves in the skull. It is shown that the proposed array can generate nonlinear waveforms with shock amplitudes >60 MPa at the focus deep inside the brain without exceeding the existing technical limitation on the intensity of 40 W/cm2 at the array elements. Such shock amplitudes are sufficient for mechanical ablation of brain tissues using the boiling histotripsy approach and implementation of other shock-based therapies.

Ultrasound-based cell sorting with microbubbles: A feasibility study.

The isolation and sorting of cells is an important process in research and hospital labs. Most large research and commercial labs incorporate fluorescently or magnetically labeled antibodies adherent to cell surface antigens for cell identification and separation. In this paper, a process is described that merges biochemical labeling with ultrasound-based separation. Instead of lasers and fluorophore tags, or magnets and magnetic particle tags, the technique uses ultrasound and microbubble tags. Streptavidin-labeled microbubbles were mixed with a human acute lymphoblastic leukemia cell line, CCL 119, conjugated with biotinylated anti-CD7 antibodies. Tagged cells were forced under ultrasound, and their displacement and velocity quantified. Differential displacement in a flow stream was quantified against erythrocytes, which showed almost no displacement under ultrasound. A model for the acoustic radiation force on the conjugated pairs compares favorably with observations. This technology may improve on current time-consuming and costly purification procedures.

Modeling of the acoustic radiation force in elastography.

Elastography is a non-invasive imaging technique that can assess in vivo tissue stiffness. In shear wave elastography imaging, the acoustic radiation force (ARF) produced by focused ultrasound generates a local force that produces shear waves. The authors compare three existing formulations for the ARF: its full expression in the second-order approximation and two simplified formulations using a quasi-plane wave and an attenuated plane wave approximation. Analytical expressions for the ARF are derived for the special cases of a concave spherical source and a quasi-Gaussian beam. They provide expressions for the resulting ARF and show discrepancies between the different formulations. For strongly divergent or highly focused beams the ARF expressed by the second-order approximation significantly differs from both simplified formulations. However, despite those differences the second-order and quasi-plane wave approximations create identical shear displacements in tissue. To compute the ARF and the displacements produced by a conventional ultrasound probe, the three formulations were incorporated into the k-Wave simulation package. The second-order and quasi-plane wave approximations give different forces but nearly identical displacements while the plane wave approximation significantly differs. It is concluded that to properly take into account the ultrasound field structure, the second-order or quasi-plane wave approximations should be preferably used.

Shock formation and nonlinear saturation effects in the ultrasound field of a diagnostic curvilinear probe.

Newer imaging and therapeutic ultrasound technologies may benefit from in situ pressure levels higher than conventional diagnostic ultrasound. One example is the recently developed use of ultrasonic radiation force to move kidney stones and residual fragments out of the urinary collecting system. A commercial diagnostic 2.3 MHz C5-2 array probe has been used to deliver the acoustic pushing pulses. The probe is a curvilinear array comprising 128 elements equally spaced along a convex cylindrical surface. The effectiveness of the treatment can be increased by using higher transducer output to provide a stronger pushing force; however nonlinear acoustic saturation can be a limiting factor. In this work nonlinear propagation effects were analyzed for the C5-2 transducer using a combined measurement and modeling approach. Simulations were based on the three-dimensional Westervelt equation with the boundary condition set to match low power measurements of the acoustic pressure field. Nonlinear focal waveforms simulated for different numbers of operating elements of the array at several output power levels were compared to fiber-optic hydrophone measurements and were found to be in good agreement. It was shown that saturation effects do limit the acoustic pressure in the focal region of a diagnostic imaging probe.

Fragmentation of urinary calculi in vitro by burst wave lithotripsy.

PURPOSE: We developed a new method of lithotripsy that uses short, broadly focused bursts of ultrasound rather than shock waves to fragment stones. We investigated the characteristics of stone comminution by burst wave lithotripsy in vitro. MATERIALS AND METHODS: Artificial and natural stones (mean ± SD size 8.2 ± 3.0 mm, range 5 to 15) were treated with ultrasound bursts using a focused transducer in a water bath. Stones were exposed to bursts with focal pressure amplitude of 6.5 MPa or less at a 200 Hz burst repetition rate until completely fragmented. Ultrasound frequencies of 170, 285 and 800 kHz were applied using 3 transducers, respectively. Time to fragmentation for each stone type was recorded and fragment size distribution was measured by sieving. RESULTS: Stones exposed to ultrasound bursts were fragmented at focal pressure amplitudes of 2.8 MPa or greater at 170 kHz. Fractures appeared along the stone surface, resulting in fragments that separated at the surface nearest to the transducer until the stone was disintegrated. All natural and artificial stones were fragmented at the highest focal pressure of 6.5 MPa with a mean treatment duration of 36 seconds for uric acid stones to 14.7 minutes for cystine stones. At a frequency of 170 kHz the largest artificial stone fragments were less than 4 mm. Exposure at 285 and 800 kHz produced only fragments less than 2 mm and less than 1 mm, respectively. CONCLUSIONS: Stone comminution with burst wave lithotripsy is feasible as a potential noninvasive treatment method for nephrolithiasis. Adjusting the fundamental ultrasound frequency allows for stone fragment size to be controlled.

Acoustic holography as a metrological tool for characterizing medical ultrasound sources and fields.

Acoustic holography is a powerful technique for characterizing ultrasound sources and the fields they radiate, with the ability to quantify source vibrations and reduce the number of required measurements. These capabilities are increasingly appealing for meeting measurement standards in medical ultrasound; however, associated uncertainties have not been investigated systematically. Here errors associated with holographic representations of a linear, continuous-wave ultrasound field are studied. To facilitate the analysis, error metrics are defined explicitly, and a detailed description of a holography formulation based on the Rayleigh integral is provided. Errors are evaluated both for simulations of a typical therapeutic ultrasound source and for physical experiments with three different ultrasound sources. Simulated experiments explore sampling errors introduced by the use of a finite number of measurements, geometric uncertainties in the actual positions of acquired measurements, and uncertainties in the properties of the propagation medium. Results demonstrate the theoretical feasibility of keeping errors less than about 1%. Typical errors in physical experiments were somewhat larger, on the order of a few percent; comparison with simulations provides specific guidelines for improving the experimental implementation to reduce these errors. Overall, results suggest that holography can be implemented successfully as a metrological tool with small, quantifiable errors.

Dynamic Mode Decomposition for Transient Cavitation Bubbles Imaging in Pulsed High-Intensity Focused Ultrasound Therapy.

Pulsed high-intensity focused ultrasound (pHIFU) can induce sparse de novo inertial cavitation without the introduction of exogenous contrast agents, promoting mild mechanical disruption in targeted tissue. Because the bubbles are small and rapidly dissolve after each HIFU pulse, mapping transient bubbles and obtaining real-time quantitative metrics correlated with tissue damage are challenging. Prior work introduced Bubble Doppler, an ultrafast power Doppler imaging method as a sensitive means to map cavitation bubbles. The main limitation of that method was its reliance on conventional wall filters used in Doppler imaging and its optimization for imaging blood flow rather than transient scatterers. This study explores Bubble Doppler enhancement using dynamic mode decomposition (DMD) of a matrix created from a Doppler ensemble for mapping and extracting the characteristics of transient cavitation bubbles. DMD was first tested in silico with a numerical dataset mimicking the spatiotemporal characteristics of backscattered signal from tissue and bubbles. The performance of DMD filter was compared to other widely used Doppler wall filter-singular value decomposition (SVD) and infinite impulse response (IIR) high-pass filter. DMD was then applied to an ex vivo tissue dataset where each HIFU pulse was immediately followed by a plane wave Doppler ensemble. In silico DMD outperformed SVD and IIR high-pass filter and ex vivo provided physically interpretable images of the modes associated with bubbles and their corresponding temporal decay rates. These DMD modes can be trackable over the duration of pHIFU treatment using k-means clustering method, resulting in quantitative indicators of treatment progression.

Dynamic Mode Decomposition for Transient Cavitation Bubbles Imaging in Pulsed High Intensity Focused Ultrasound Therapy.

Pulsed high-intensity focused ultrasound (pHIFU) can induce sparse de novo inertial cavitation without the introduction of exogenous contrast agents, promoting mild mechanical disruption in targeted tissue. Because the bubbles are small and rapidly dissolve after each HIFU pulse, mapping transient bubbles and obtaining real-time quantitative metrics correlated to tissue damage are challenging. Prior work introduced Bubble Doppler, an ultrafast power Doppler imaging method as a sensitive means to map cavitation bubbles. The main limitation of that method was its reliance on conventional wall filters used in Doppler imaging and optimized for imaging blood flow rather than transient scatterers. This study explores Bubble Doppler enhancement using dynamic mode decomposition (DMD) of a matrix created from a Doppler ensemble for mapping and extracting the characteristics of transient cavitation bubbles. DMD was first tested in silico with a numerical dataset mimicking the spatiotemporal characteristics of backscattered signal from tissue and bubbles. The performance of DMD filter was compared to other widely used Doppler wall filters - singular value decomposition (SVD) and infinite impulse response (IIR) highpass filter. DMD was then applied to an ex vivo tissue dataset where each HIFU pulse was immediately followed by a plane wave Doppler ensemble. In silico DMD outperformed SVD and IIR high pass filter and ex vivo provided physically interpretable images of the modes associated with bubbles and their corresponding temporal decay rates. These DMD modes can be trackable over the duration of pHIFU treatment using k-means clustering method, resulting in quantitative indicators of treatment progression.

Treatment Planning and Aberration Correction Algorithm for HIFU Ablation of Renal Tumors.

High-intensity focused ultrasound (HIFU) applications for thermal or mechanical ablation of renal tumors often encounter challenges due to significant beam aberration and refraction caused by oblique beam incidence, inhomogeneous tissue layers, and presence of gas and bones within the beam. These losses can be significantly mitigated through sonication geometry planning, patient positioning, and aberration correction using multielement phased arrays. Here, a sonication planning algorithm is introduced, which uses the simulations to select the optimal transducer position and evaluate the effect of aberrations and acoustic field quality at the target region after aberration correction. Optimization of transducer positioning is implemented using a graphical user interface (GUI) to visualize a segmented 3-D computed tomography (CT)-based acoustic model of the body and to select sonication geometry through a combination of manual and automated approaches. An HIFU array (1.5 MHz, 256 elements) and three renal cell carcinoma (RCC) cases with different tumor locations and patient body habitus were considered. After array positioning, the correction of aberrations was performed using a combination of backpropagation from the focus with an ordinary least squares (OLS) optimization of phases at the array elements. The forward propagation was simulated using a combination of the Rayleigh integral and k-space pseudospectral method (k-Wave toolbox). After correction, simulated HIFU fields showed tight focusing and up to threefold higher maximum pressure within the target region. The addition of OLS optimization to the aberration correction method yielded up to 30% higher maximum pressure compared to the conventional backpropagation and up to 250% higher maximum pressure compared to the ray-tracing method, particularly in strongly distorted cases.

Elastic Properties of Aging Human Hematoma Model In Vitro and Its Susceptibility to Histotripsy Liquefaction.

OBJECTIVE: Tissue susceptibility to histotripsy disintegration has been reported to depend on its elastic properties. This work was aimed at investigation of histotripsy efficiency for liquefaction of human hematomas, depending on their stiffness and degree of retraction over time (0-10 d). METHODS: As an in vitro hematoma model, anticoagulated human blood samples (200 mL) were recalcified at different temperatures. In one set of samples, the shear modulus was measured by shear wave elastography during blood clotting at 10℃, 22℃ and 37℃, and then daily during further aging. The ultrastructure of the samples was analyzed daily with scanning electron microscopy (SEM). Another set of blood samples (50-200 mL) were recalcified at 37℃ for density and retraction measurements over aging and exposed to histotripsy at varying time points. Boiling histotripsy (2.5 ms pulses) and hybrid histotripsy (0.2 ms pulses) exposures (2 MHz, 1% dc, P+/P-/As = 182/-27/207 MPa in situ) were used to produce either individual cigar-shaped or volumetric (0.8-3 mL) lesions in samples incubated for 3 h, 5 d and 10 d. The obtained lesions were sized, then the lysate aspirated under B-mode guidance was analyzed ultrastructurally and diluted in distilled water for sizing of residual fragments. RESULTS: It was found that clotting time decreased from 113 to 25 min with the increase in blood temperature from 10℃ to 37℃. The shear modulus increased to 0.53 ± 0.17 kPa during clotting and remained constant within 8 d of incubation at 2℃. Sample volumes decreased by 57% because of retraction within 10 d. SEM revealed significant echinocytosis but unchanged ultrastructure of the fibrin meshwork. Liquefaction rate and lesion dimensions produced with the same histotripsy protocols correlated with the increase in the degree of retraction and were lower in retracted samples versus freshly clotted samples. More than 80% of residual fibrin fragments after histotripsy treatment were shorter than 150 µm; the maximum length was 208 µm, allowing for unobstructed aspiration of the lysate with most clinically used needles. CONCLUSION: The results indicate that hematoma susceptibility to histotripsy liquefaction is not entirely determined by its stiffness, and correlates with the retraction degree.

Radiation force of an arbitrary acoustic beam on an elastic sphere in a fluid.

A theoretical approach is developed to calculate the radiation force of an arbitrary acoustic beam on an elastic sphere in a liquid or gas medium. First, the incident beam is described as a sum of plane waves by employing conventional angular spectrum decomposition. Then, the classical solution for the scattering of a plane wave from an elastic sphere is applied for each plane-wave component of the incident field. The net scattered field is expressed as a superposition of the scattered fields from all angular spectrum components of the incident beam. With this formulation, the incident and scattered waves are superposed in the far field to derive expressions for components of the radiation stress tensor. These expressions are then integrated over a spherical surface to analytically describe the radiation force on an elastic sphere. Limiting cases for particular types of incident beams are presented and are shown to agree with known results. Finally, the analytical expressions are used to calculate radiation forces associated with two specific focusing transducers.

Phase holograms for the three-dimensional patterning of unconstrained microparticles.

Acoustic radiation forces can remotely manipulate particles. Forces from a standing wave field align microscale particles along the nodal or anti-nodal locations of the field to form three-dimensional (3D) patterns. These patterns can be used to form 3D microstructures for tissue engineering applications. However, standing wave generation requires more than one transducer or a reflector, which is challenging to implement in vivo. Here, a method is developed and validated to manipulate microspheres using a travelling wave from a single transducer. Diffraction theory and an iterative angular spectrum approach are employed to design phase holograms to shape the acoustic field. The field replicates a standing wave and aligns polyethylene microspheres in water, which are analogous to cells in vivo, at pressure nodes. Using Gor'kov potential to calculate the radiation forces on the microspheres, axial forces are minimized, and transverse forces are maximized to create stable particle patterns. Pressure fields from the phase holograms and resulting particle aggregation patterns match predictions with a feature similarity index > 0.92, where 1 is a perfect match. The resulting radiation forces are comparable to those produced from a standing wave, which suggests opportunities for in vivo implementation of cell patterning toward tissue engineering applications.