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

Histotripsy: A Method for Mechanical Tissue Ablation with Ultrasound

Histotripsy is a relatively new therapeutic ultrasound technology to mechanically liquefy tissue into subcellular debris using high-amplitude focused ultrasound pulses. In contrast to conventional high-intensity focused ultrasound thermal therapy, histotripsy has specific clinical advantages: the capacity for real-time monitoring using ultrasound imaging, diminished heat sink effects resulting in lesions with sharp margins, effective removal of the treated tissue, a tissue-selective feature to preserve crucial structures, and immunostimulation. The technology is being evaluated in small and large animal models for treating cancer, thrombosis, hematomas, abscesses, and biofilms; enhancing tumor-specific immune response; and neurological applications. Histotripsy has been recently approved by the US Food and Drug Administration to treat liver tumors, with clinical trials undertaken for benign prostatic hyperplasia and renal tumors. This review outlines the physical principles of various types of histotripsy; presents major parameters of the technology and corresponding hardware and software, imaging methods, and bioeffects; and discusses the most promising preclinical and clinical applications.

Characterization of weakly nonlinear effects in relationship to transducer parameters in focused ultrasound therapy

BACKGROUND

Focused ultrasound therapy has been widely used for the treatment of various diseases, employing different types of transducers. The focused ultrasound pressure fields inevitably exhibit nonlinear effects, which can influence the ablation region. However, the nonlinear effects exhibit noticeable variations across different applications. The characterization of the nonlinear pressure fields of ultrasound is important for the effective implementation of focused ultrasound therapy.

PURPOSE

The traditional angular spectrum method (ASM) was extended to accurately and efficiently simulate the propagation of weakly nonlinear ultrasound in heterogeneous mediums of clinical model. The nonlinear effects were further analyzed in relationship to the transducer parameters that are different in various applications.

METHODS

The pressure fields were simulated using the extended ASM, incorporating calculations for phase aberration in the frequency domain and magnitude compensation in the spatial domain to account for heterogeneous acoustic impedance mismatch. Validation was performed by comparison to k-Wave simulation results using two simplified clinical models, an abdominal soft tissue and a transcranial skull model. The nonlinear effects were then analyzed in relation to the transducer parameters of f-number and effective source area based on the same acoustic output power. The analysis of nonlinear effects was conducted under both homogeneous medium and the clinical models.

RESULTS

The simulation results demonstrated a maximum error of 3.93% in the calculated harmonic pressure of the abdominal model, and a maximum error of 4.89% within the transcranial model when comparing the extended ASM simulation results to those obtained from k-Wave simulations. The characterization of the nonlinear effects reveals a strong correlation with the transducer parameters. Specifically, the results indicate that the nonlinear effects intensify with an increase in the effective source area and f-number, under the same acoustic output power of the transducer. However, the clinical model also showed an influence on the nonlinear effects in relation to the f-number.

CONCLUSION

The extended ASM was demonstrated as an accurate and efficient simulation tool, and the simulation results provide a reference for evaluating the intensity of nonlinear effects in various transducer designs.

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.

A simulation study on the sensitivity of transcranial ray-tracing ultrasound modeling to skull properties

In transcranial focused ultrasound therapies, such as treating essential tremor via thermal ablation in the thalamus, acoustic energy is focused through the skull using a phased-array transducer. Ray tracing is a computationally efficient method that can correct skull-induced phase aberrations via per-element phase delay calculations using patient-specific computed tomography (CT) data. However, recent studies show that variations in CT-derived Hounsfield unit may account for only 50% of the speed of sound variability in human skull specimens, potentially limiting clinical transcranial ultrasound applications. Therefore, understanding the sensitivity of treatment planning methods to material parameter variations is essential. The present work uses a ray-tracing simulation model to explore how imprecision in model inputs, arising from clinically significant uncertainties in skull properties or considerations of acoustic phenomena, affects acoustic focusing quality through the skull. We propose and validate new methods to optimize ray-tracing skull simulations for clinical treatment planning, relevant for predicting intracranial target's thermal rise, using experimental data from ex-vivo human skulls.

Physics-informed neural networks for transcranial ultrasound wave propagation

Transcranial ultrasound imaging has been playing an increasingly important role in the non-invasive treatment of brain disorders. However, the conventional mesh-based numerical wave solvers, which are an integral part of imaging algorithms, suffer from limitations such as high computational cost and discretization error in predicting the wavefield passing through the skull. In this paper, we explore the use of physics-informed neural networks (PINNs) for predicting the transcranial ultrasound wave propagation. The wave equation, two sets of time snapshots data and a boundary condition (BC) are embedded as physical constraints in the loss function during training. The proposed approach has been validated by solving the two-dimensional (2D) acoustic wave equation under three increasingly complex spatially varying velocity models. Our cases demonstrate that due to the meshless nature of PINNs, they can be flexibly applied to different wave equations and types of BCs. By adding physics constraints to the loss function, PINNs can predict wavefields far outside the training data, providing ideas for improving the generalization capability of existing deep learning methods. The proposed approach offers exciting perspectives because of the powerful framework and simple implementation. We conclude with a summary of the strengths, limitations and further research directions of this work.

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.

A Modular, Kerf-Minimizing Approach for Therapeutic Ultrasound Phased Array Construction

A novel method for fabricating a modular, kerf-minimizing histotripsy phased array was developed and tested. The method utilizes arbitrarily shaped elements, 3-D printing, water jet cutting, and a thin, 125- [Formula: see text] electrically insulating epoxy coating to maximize aperture utilization while allowing for replacement of individual transducer modules. The method was used to fabricate a 750-kHz truncated circular aperture array (165 mm ×234 mm) transducer with a focal length of 142 mm. The aperture was segmented into 260 arc-shaped modular elements, each approximately 11.5 mm ×11.5 mm, arranged in concentric rings. The resulting aperture utilization was 92%. The full-width-half-maximum (FWHM) focal zone of the array was measured to be 1.6 mm ×1.1 mm ×4.5 mm, and the FWHM electrical steering range was measured to be 38.5 mm ×33 mm 40 mm. The array was estimated to be capable of generating approximately 120-MPa peak negative pressure at the geometric focus. In addition, the array was used to ablate a 5-cm3 volume of tissue with electric focal steering.

A Comparative Feasibility Study for Transcranial Extracorporeal Shock Wave Therapy

The potential beneficial regenerative and stimulatory extracorporeal shock wave therapy (ESWT) applications to the central nervous system have garnered interest in recent years. Treatment zones for these indications are acoustically shielded by bones, which heavily impact generated sound fields. We present the results of high-resolution tissue-realistic simulations, comparing the viability of different ESWT applicators in their use for transcranial applications. The performances of electrohydraulic, electromagnetic, and piezoelectric transducers for key reflector geometries are compared. Based on density information obtained from CT imaging of the head, we utilized the non-linear wave propagation toolset Matlab k-Wave to obtain spatial therapeutic sound field geometries and waveforms. In order to understand the reliability of results on the appropriate modeling of the skull, three different bone attenuation models were compared. We find that all currently clinically ESWT applicator technologies show significant retention of peak pressures and energies past the bone barrier. Electromagnetic transducers maintain a significantly higher energy flux density compared to other technologies while low focusing strength piezoelectric applicators have the weakest transmissions. Attenuation estimates provide insights into sound field degradation and energy losses, indicating that effective transcranial therapies can readily be attained with current applicators. Furthermore, the presented approach will allow for future targeted in silico development and the design of applicators and therapy plans to ultimately improve therapeutic outcomes.

Transcranial ultrasound stimulation of the human motor cortex

It has been 40 years since the report of long-term synaptic plasticity on the rodent brain. Transcranial ultrasound stimulation (TUS) shows advantages in spatial resolution and penetration depth when compared with electrical or magnetic stimulation. The repetitive TUS (rTUS) can induce cortical excitability alteration on animals, and persistent aftereffects were observed. However, the effects of rTUS on synaptic plasticity in humans remain unelucidated. In the current study, we applied a 15-min rTUS protocol to stimulate left primary motor cortex (l-M1) in 24 male healthy participants. The single-pulsed transcranial magnetic stimulation-evoked motor evoked potential and Stop-signal task was applied to measure the rTUS aftereffects. Here, we report that conditioning the human motor cortex using rTUS may produce long-lasting and statistically significant effects on motor cortex excitability as well as motor behavior, without harmful side effects observed. These findings suggest a considerable potential of rTUS in cortical plasticity modulation and clinical intervention for impulsivity-related disorders.

"HIFU Beam:" A Simulator for Predicting Axially Symmetric Nonlinear Acoustic Fields Generated by Focused Transducers in a Layered Medium

"HIFU beam" is a freely available software tool that comprises a MATLAB toolbox combined with a user-friendly interface and binary executable compiled from FORTRAN source code (HIFU beam. (2021). Available: http://limu.msu.ru/node/3555?language=en). It is designed for simulating high-intensity focused ultrasound (HIFU) fields generated by single-element transducers and annular arrays with propagation in flat-layered media that mimic biological tissues. Numerical models incorporated in the simulator include evolution-type equations, either the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation or one-way Westervelt equation, for radially symmetric ultrasound beams in homogeneous and layered media with thermoviscous or power-law acoustic absorption. The software uses shock-capturing methods that allow for simulating strongly nonlinear acoustic fields with high-amplitude shocks. In this article, a general description of the software is given along with three representative simulation cases of ultrasound transducers and focusing conditions typical for therapeutic applications. The examples illustrate major nonlinear wave effects in HIFU fields including shock formation. Two examples simulate propagation in water, involving a single-element source (1-MHz frequency, 100-mm diameter, 90-mm radius of curvature) and a 16-element annular array (3-MHz frequency, 48-mm diameter, and 35-mm radius of curvature). The third example mimics the scenario of a HIFU treatment in a "water-muscle-kidney" layered medium using a source typical for abdominal HIFU applications (1.2-MHz frequency, 120-mm diameter, and radius of curvature). Linear, quasi-linear, and shock-wave exposure protocols are considered. It is intended that "HIFU beam" can be useful in teaching nonlinear acoustics; designing and characterizing high-power transducers; and developing exposure protocols for a wide range of therapeutic applications such as shock-based HIFU, boiling histotripsy, drug delivery, immunotherapy, and others.

Design of a Versatile Angle-Rotatable Skull-Shaped Conformal Transcranial Focused Ultrasound Transducer for Noninvasive Brain Therapy

Phase-controlled focused ultrasound transducers have opened up a new way to noninvasively treat brain diseases. However, due to the geometry and the heterogeneous medium of the human skull, the currently used hemispherical transducers cannot guarantee that the ultrasound emitted by the transducer penetrates the skull normally; consequently, the impact of the shear wave on the treatment area has to be considered. The usual approach is to turn off elements in the transducer with excessively large incident angles, but this approach reduces the efficiency of the transducer. This article presents the design of a novel transducer that can be rotated according to the different geometric shapes of the skull. The incident angles of the elements in the transducer are modified so that the effect of the shear wave on the treatment area can be ignored and the efficiency of the transducer is improved. The results of computed tomography (CT) model simulation of nine skulls verify the versatility and feasibility of the new transducer. We believe that this kind of rotatable transducer has clinical and engineering significance.

Ultrasound for the Brain: A Review of Physical and Engineering Principles, and Clinical Applications

The emergence of new ultrasound technologies has improved our understanding of the brain functions and offered new opportunities for the treatment of brain diseases. Ultrasound has become a valuable tool in preclinical animal and clinical studies as it not only provides information about the structure and function of brain tissues but can also be used as a therapy alternative for brain diseases. High-resolution cerebral flow images with high sensitivity can be acquired using novel functional ultrasound and super-resolution ultrasound imaging techniques. The noninvasive treatment of essential tremors has been clinically approved and it has been demonstrated that the ultrasound technology can revolutionize the currently existing treatment methods. Microbubble-mediated ultrasound can remotely open the blood-brain barrier enabling targeted drug delivery in the brain. More recently, ultrasound neuromodulation received a great amount of attention due to its noninvasive and deep penetration features and potential therapeutic benefits. This review provides a thorough introduction to the current state-of-the-art research on brain ultrasound and also introduces basic knowledge of brain ultrasound including the acoustic properties of the brain/skull and engineering techniques for ultrasound. Ultrasound is expected to play an increasingly important role in the diagnosis and therapy of brain diseases.

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.

Systematic examination of low-intensity ultrasound parameters on human motor cortex excitability and behavior

Low-intensity transcranial ultrasound (TUS) can non-invasively modulate human neural activity. We investigated how different fundamental sonication parameters influence the effects of TUS on the motor cortex (M1) of 16 healthy subjects by probing cortico-cortical excitability and behavior. A low-intensity 500 kHz TUS transducer was coupled to a transcranial magnetic stimulation (TMS) coil. TMS was delivered 10 ms before the end of TUS to the left M1 hotspot of the first dorsal interosseous muscle. Varying acoustic parameters (pulse repetition frequency, duty cycle, and sonication duration) on motor-evoked potential amplitude were examined. Paired-pulse measures of cortical inhibition and facilitation, and performance on a visuomotor task was also assessed. TUS safely suppressed TMS-elicited motor cortical activity, with longer sonication durations and shorter duty cycles when delivered in a blocked paradigm. TUS increased GABAA-mediated short-interval intracortical inhibition and decreased reaction time on visuomotor task but not when controlled with TUS at near-somatosensory threshold intensity.

Nonlinear ultrasound simulation in an axisymmetric coordinate system using a k-space pseudospectral method

A full-wave model for nonlinear ultrasound propagation through a heterogeneous and absorbing medium in an axisymmetric coordinate system is developed. The model equations are solved using a nonstandard or k-space pseudospectral time domain method. Spatial gradients in the axial direction are calculated using the Fourier collocation spectral method, and spatial gradients in the radial direction are calculated using discrete trigonometric transforms. Time integration is performed using a k-space corrected finite difference scheme. This scheme is exact for plane waves propagating linearly in the axial direction in a homogeneous and lossless medium and significantly reduces numerical dispersion in the more general case. The implementation of the model is described, and performance benchmarks are given for a range of grid sizes. The model is validated by comparison with several analytical solutions. This includes one-dimensional absorption and nonlinearity, the pressure field generated by plane-piston and bowl transducers, and the scattering of a plane wave by a sphere. The general utility of the model is then demonstrated by simulating nonlinear transcranial ultrasound using a simplified head model.

Focused Ultrasound Platform for Investigating Therapeutic Neuromodulation Across the Human Hippocampus

Pulsed low-intensity focused ultrasound (PLIFUS) has shown promise in inducing neuromodulation in several animal and human studies. Therefore, it is of clinical interest to develop experimental platforms to test repetitive PLIFUS as a therapeutic modality in humans with neurologic disorders. In the study described here, our aim was to develop a laboratory-built experimental device platform intended to deliver repetitive PLIFUS across the hippocampus in seizure onset zones of patients with drug-resistant temporal lobe epilepsy. The system uses neuronavigation targeting over multiple therapeutic sessions. PLIFUS (548 kHz) was emitted across multiple hippocampal targets in a human subject with temporal lobe epilepsy using a mechanically steered piezoelectric transducer. Stimulation was delivered up to 2.25 W/cm2 spatial peak temporal average intensity (free-field equivalent), with 36%-50% duty cycle, 500-ms sonications and 7-s inter-stimulation intervals lasting 140 s per target and repeated for multiple sessions. A first-in-human PLIFUS course of treatment was successfully delivered using the device platform with no adverse events.