Design aspects of focal beams from high-intensity arrays

Element Position Calibration for Matrix Array Transducers with Multiple Disjoint Piezoelectric Panels

Two-dimensional ultrasound transducers enable the acquisition of fully volumetric data that have been demonstrated to provide greater diagnostic information in the clinical setting and are a critical tool for emerging ultrasound methods, such as super-resolution and functional imaging. This technology, however, is not without its limitations. Due to increased fabrication complexity, some matrix probes with disjoint piezoelectric panels may require initial calibration. In this manuscript, two methods for calibrating the element positions of the Vermon 1024-channel 8 MHz matrix transducer are detailed. This calibration is a necessary step for acquiring high resolution B-mode images while minimizing transducer-based image degradation. This calibration is also necessary for eliminating vessel-doubling artifacts in super-resolution images and increasing the overall signal-to-noise ratio (SNR) of the image. Here, we show that the shape of the point spread function (PSF) can be significantly improved and PSF-doubling artifacts can be reduced by up to 10 dB via this simple calibration procedure.

Simulation study on high spatio-temporal resolution acousto-electrophysiological neuroimaging

Objective.Acousto-electrophysiological neuroimaging (AENI) is a technique hypothesized to record electrophysiological activity of the brain with millimeter spatial and sub-millisecond temporal resolution. This improvement is obtained by tagging areas with focused ultrasound (fUS). Due to mechanical vibration with respect to the measuring electrodes, the electrical activity of the marked region will be modulated onto the ultrasonic frequency. The region's electrical activity can subsequently be retrieved via demodulation of the measured signal. In this study, the feasibility of this hypothesized technique is tested.Approach.This is done by calculating the forward electroencephalography response under quasi-static assumptions. The head is simplified as a set of concentric spheres. Two sizes are evaluated representing human and mouse brains. Moreover, feasibility is assessed for wet and dry transcranial, and for cortically placed electrodes. The activity sources are modeled by dipoles, with their current intensity profile drawn from a power-law power spectral density.Results.It is shown that mechanical vibration modulates the endogenous activity onto the ultrasonic frequency. The signal strength depends non-linearly on the alignment between dipole orientation, vibration direction and recording point. The strongest signal is measured when these three dependencies are perfectly aligned. The signal strengths are in the pV-range for a dipole moment of 5 nAm and ultrasonic pressures within Food and Drug Administration (FDA)-limits. The endogenous activity can then be accurately reconstructed via demodulation. Two interference types are investigated: vibrational and static. Depending on the vibrational interference, it is shown that millimeter resolution signal detection is possible also for deep brain regions. Subsequently, successful demodulation depends on the static interference, that at MHz-range has to be sub-picovolt.Significance.Our results show that mechanical vibration is a possible underlying mechanism of acousto-electrophyisological neuroimaging. This paper is a first step towards improved understanding of the conditions under which AENI is feasible.

Design of a 1-MHz Therapeutic Ultrasound Array for Small Volume Blood-Brain Barrier Opening at Cortical Targets in Macaques

[[gabstract]][] Focused ultrasound (FUS) can temporarily open the blood-brain barrier (BBB) and increase the delivery of chemotherapeutics, viral vectors, and other agents to the brain parenchyma. To limit FUS BBB opening to a single brain region, the transcranial acoustic focus of the ultrasound transducer must not be larger than the region targeted. In this work, we design and characterize a therapeutic array optimized for BBB opening at the frontal eye field (FEF) in macaques. We used 115 transcranial simulations in four macaques varying f-number and frequency to optimize the design for focus size, transmission, and small device footprint. The design leverages inward steering for focus tightening, a 1-MHz transmit frequency, and can focus to a simulation predicted 2.5- ± 0.3-mm lateral and 9.5- ± 1.0-mm axial full-width at half-maximum spot size at the FEF without aberration correction. The array is capable of steering axially 35 mm outward, 26 mm inward, and laterally 13 mm with 50% the geometric focus pressure. The simulated design was fabricated, and we characterized the performance of the array using hydrophone beam maps in a water tank and through an ex vivo skull cap to compare measurements with simulation predictions, achieving a 1.8-mm lateral and 9.5-mm axial spot size with a transmission of 37% (transcranial, phase corrected). The transducer produced by this design process is optimized for BBB opening at the FEF in macaques.

Multiple-Focus Patterns of Sparse Random Array Using Particle Swarm Optimization for Ultrasound Surgery

This study aims to investigate the feasibility and potential of sparse random arrays driven by the particle swarm optimization (PSO) algorithm to generate multiple-focus patterns and a large scanning range without grating lobes, which extends the scanning range of focused ultrasound in the treatment of brain tumors, opening the blood-brain barrier, and neuromodulation. Operating at 1.1 MHz, a random spherical array with 200 square elements (sparseness 58%) and a sparse random array with 660 square elements (sparseness 41%) driven by PSO are employed to simulate different focus patterns. With the same radius of curvature and diameter of transducer and element size, the scanning range of the off-axis single focus of a random 200-element array is two times that of an ordinary array using symmetric arrangement. The focal volume of multiple-focus patterns of the random array is 18 times that of the single focus. The single focus of the sparse random array with 660 elements could steer up to ±23 mm in the radial direction, without grating lobes. The maximum distance between two foci in a multiple-focus "S"-shaped deflection is approximately 25 mm. Simulation results illustrate the capability of a focused beam steered in 3-D space. Multiple-focus patterns could significantly increase the focal volume and shorten the treatment time for large target volumes. Simulation results show the feasibility and potential of the method combining PSO with a sparse random array to generate flexible focus patterns that can adapt to different needs in different tissue treatments.

Transcranial Neuromodulation Array With Imaging Aperture for Simultaneous Multifocus Stimulation in Nonhuman Primates

Even simple behaviors arise from the simultaneous activation of multiple regions in the brain. Thus, the ability to simultaneously stimulate multiple regions within a brain circuit should allow for better modulation of function. However, performing simultaneous multifocus ultrasound neuromodulation introduces challenges to transducer design. Using 3-D Fullwave simulations, we have designed an ultrasound neuromodulation array for nonhuman primates that: 1) can simultaneously focus on multiple targets and 2) include an imaging aperture for additional functional imaging. This design is based on a spherical array, with 128 15-mm elements distributed in a spherical helix pattern. It is shown that clustering the elements tightly around the 65-mm imaging aperture located at the top of the array improves targeting at shallow depths, near the skull surface. Spherical arrays have good focusing capabilities through the skull at the center of the array, but focusing on off-center locations is more challenging due to the natural geometric configuration and the angle of incidence with the skull. In order to mitigate this, the 64 elements closest to the aperture were rotated toward and focusing on a shallow target, and the 64 elements farthest from the aperture were rotated toward and focusing on a deeper target. Data illustrated that this array produced focusing on the somatosensory cortex with a gain of 4.38 and to the thalamus with a gain of 3.82. To improve upon this, the array placement was optimized based on phase aberration simulations, allowing for the elements with the largest impact on the gain at each focal point to be found. This optimization resulted in an array design that can focus on the somatosensory cortex with a gain of 5.19 and the thalamus with a gain of 4.45. Simulations were also performed to evaluate the ability of the array to focus on 28 additional brain regions, showing that off-center target regions can be stimulated, but those closer to the skull will require corrective steps to deliver the same amount of energy to those locations. This simulation and design process can be adapted to an individual monkey or human skull morphologies and specific target locations within individuals by using orientable 3-D printing of the transducer case and by electronic phase aberration correction.

A Prototype Therapy System for Boiling Histotripsy in Abdominal Targets Based on a 256-Element Spiral Array

Boiling histotripsy (BH) uses millisecond-long ultrasound (US) pulses with high-amplitude shocks to mechanically fractionate tissue with potential for real-time lesion monitoring by US imaging. For BH treatments of abdominal organs, a high-power multielement phased array system capable of electronic focus steering and aberration correction for body wall inhomogeneities is needed. In this work, a preclinical BH system was built comprising a custom 256-element 1.5-MHz phased array (Imasonic, Besançon, France) with a central opening for mounting an imaging probe. The array was electronically matched to a Verasonics research US system with a 1.2-kW external power source. Driving electronics and software of the system were modified to provide a pulse average acoustic power of 2.2 kW sustained for 10 ms with a 1-2-Hz repetition rate for delivering BH exposures. System performance was characterized by hydrophone measurements in water combined with nonlinear wave simulations based on the Westervelt equation. Fully developed shocks of 100-MPa amplitude are formed at the focus at 275-W acoustic power. Electronic steering capabilities of the array were evaluated for shock-producing conditions to determine power compensation strategies that equalize BH exposures at multiple focal locations across the planned treatment volume. The system was used to produce continuous volumetric BH lesions in ex vivo bovine liver with 1-mm focus spacing, 10-ms pulselength, five pulses/focus, and 1% duty cycle.

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/cm² 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.

A random phased-array for MR-guided transcranial ultrasound neuromodulation in non-human primates

Transcranial focused ultrasound (FUS) is a non-invasive technique for therapy and study of brain neural activation. Here we report on the design and characterization of a new MR-guided FUS transducer for neuromodulation in non-human primates at 650 kHz. The array is randomized with 128 elements 6.6 mm in diameter, radius of curvature 7.2 cm, opening diameter 10.3 cm (focal ratio 0.7), and 46% coverage. Simulations were used to optimize transducer geometry with respect to focus size, grating lobes, and directivity. Focus size and grating lobes during electronic steering were quantified using hydrophone measurements in water and a three-axis stage. A novel combination of optical tracking and acoustic mapping enabled measurement of the 3D pressure distribution in the cortical region of an ex vivo skull to within ~3.5 mm of the surface, and allowed accurate modelling of the experiment via non-homogeneous 3D acoustic simulations. The data demonstrates acoustic focusing beyond the skull bone, with the focus slightly broadened and shifted proximal to the skull. The fabricated design is capable of targeting regions within the S1 sensorimotor cortex of macaques.

In situ bone tissue engineering via ultrasound-mediated gene delivery to endogenous progenitor cells in mini-pigs

More than 2 million bone-grafting procedures are performed each year using autografts or allografts. However, both options carry disadvantages, and there remains a clear medical need for the development of new therapies for massive bone loss and fracture nonunions. We hypothesized that localized ultrasound-mediated, microbubble-enhanced therapeutic gene delivery to endogenous stem cells would induce efficient bone regeneration and fracture repair. To test this hypothesis, we surgically created a critical-sized bone fracture in the tibiae of Yucatán mini-pigs, a clinically relevant large animal model. A collagen scaffold was implanted in the fracture to facilitate recruitment of endogenous mesenchymal stem/progenitor cells (MSCs) into the fracture site. Two weeks later, transcutaneous ultrasound-mediated reporter gene delivery successfully transfected 40% of cells at the fracture site, and flow cytometry showed that 80% of the transfected cells expressed MSC markers. Human bone morphogenetic protein-6 (BMP-6) plasmid DNA was delivered using ultrasound in the same animal model, leading to transient expression and secretion of BMP-6 localized to the fracture area. Micro-computed tomography and biomechanical analyses showed that ultrasound-mediated BMP-6 gene delivery led to complete radiographic and functional fracture healing in all animals 6 weeks after treatment, whereas nonunion was evident in control animals. Collectively, these findings demonstrate that ultrasound-mediated gene delivery to endogenous mesenchymal progenitor cells can effectively treat nonhealing bone fractures in large animals, thereby addressing a major orthopedic unmet need and offering new possibilities for clinical translation.

Development of a spherically focused phased array transducer for ultrasonic image-guided hyperthermia

A 1.5 MHz prolate spheroidal therapeutic array with 128 circular elements was designed to accommodate standard imaging arrays for ultrasonic image-guided hyperthermia. The implementation of this dual-array system integrates real-time therapeutic and imaging functions with a single ultrasound system (Vantage 256, Verasonics). To facilitate applications involving small animal imaging and therapy the array was designed to have a beam depth of field smaller than 3.5 mm and to electronically steer over distances greater than 1 cm in both the axial and lateral directions. In order to achieve the required f number of 0.69, 1-3 piezocomposite modules were mated within the transducer housing. The performance of the prototype array was experimentally evaluated with excellent agreement with numerical simulation. A focal volume (2.70 mm (axial)  ×  0.65 mm (transverse)  ×  0.35 mm (transverse)) defined by the  -6 dB focal intensity was obtained to address the dimensions needed for small animal therapy. An electronic beam steering range defined by the  -3 dB focal peak intensity (17 mm (axial)  ×  14 mm (transverse)  ×  12 mm (transverse)) and  -8 dB lateral grating lobes (24 mm (axial)  ×  18 mm (transverse)  ×  16 mm (transverse)) was achieved. The combined testing of imaging and therapeutic functions confirmed well-controlled local heating generation and imaging in a tissue mimicking phantom. This dual-array implementation offers a practical means to achieve hyperthermia and ablation in small animal models and can be incorporated within protocols for ultrasound-mediated drug delivery.

Frequency considerations for deep ablation with high-intensity focused ultrasound: A simulation study

PURPOSE

The objective of this study was to explore frequency considerations for large-volume, deep thermal ablations with focused ultrasound. Though focal patterns, focal steering rate, and the size of focal clusters have all been explored in this context, frequency studies have generally explored shallower depths and hyperthermia applications. This study examines both treatment efficiency and near-field heating rate as functions of frequency and depth.

METHODS

Flat, 150 mm transducer arrays were simulated to operate at frequencies of 250, 500, 750, 1000, 1250, and 1500 kHz. Each array had λ2 interelement spacing yielding arrays of 2000-70 000 piston-shaped elements arranged in concentric rings. Depths of 50, 100, and 150 mm were explored, with attenuation (α) values of 2.5-10 (Np/m)/MHz. Ultrasound propagation was simulated with the Rayleigh-Sommerfeld integral over a volume of homogeneous simulated tissue. Absorbed power density was determined from the acoustic pressure which, in turn, was modeled with the Pennes bioheat transfer equation. Using this knowledge of temperature over time, thermal dose function of Sapareto and Dewey was used to model the resulting bioeffect of each simulated sonication. Initially, single foci at each depth, frequency, and α were examined with either fixed peak temperatures or fixed powers. Based on the size of the resulting, single foci lesions, larger compound sonications were designed with foci packed together in multiple layers and rings. For each depth, focal patterns were chosen to produce a similar total ablated volume for each frequency. These compound sonications were performed with a fixed peak temperature at each focus. The resulting energy efficiency (volume ablated per acoustic energy applied), near-field heating rate (temperature increase in the anterior third of the simulation space per unit volume ablated), and near- and far-field margins were assessed.

RESULTS

Lesions of comparable volume were created with different frequencies at different depths. The results reflect the interconnected nature of frequency as it effects focal size (decreasing with frequency), peak pressure (generally increasing with frequency), and attenuation (also increasing with frequency). The ablation efficiency was the highest for α = 5 (Np/m)/MHz at a frequency of 750 kHz at each depth. For α = 10 (Np/m)/MHz, efficiency was the highest at 750 kHz for a depth of 50 mm, and 500 kHz at depths of 100 and 150 mm. At all sonication depths, near-field heating was minimized with lower frequencies of 250 and 500 kHz.

CONCLUSIONS

Large-volume ablations are most efficient at frequencies of 500-750 kHz at depths of 100-150 mm. When one considers that near-field heat accumulation tends to be the rate limiting factor in large-volume ablations like uterine fibroid surgery, the results show that frequencies as low as 500 kHz are favored for their ability to reduce heating in the near-field.

Flexible integration of high-imaging-resolution and high-power arrays for ultrasound-induced thermal strain imaging (US-TSI)

Ultrasound-induced thermal strain imaging (USTSI) for carotid artery plaque detection requires both high imaging resolution (<100 μm) and sufficient US-induced heating to elevate the tissue temperature (~1°C to 3°C within 1 to 3 cardiac cycles) to produce a noticeable change in sound speed in the targeted tissues. Because the optimization of both imaging and heating in a monolithic array design is particularly expensive and inflexible, a new integrated approach is presented which utilizes independent ultrasound arrays to meet the requirements for this particular application. This work demonstrates a new approach in dual-array construction. A 3-D printed manifold was built to support both a high-resolution 20 MHz commercial imaging array and 6 custom heating elements operating in the 3.5 to 4 MHz range. For the application of US-TSI in carotid plaque characterization, the tissue target site is 20 to 30 mm deep, with a typical target volume of 2 mm (elevation) × 8 mm (azimuthal) × 5 mm (depth). The custom heating array performance was fully characterized for two design variants (flat and spherical apertures), and can easily deliver 30 W of total acoustic power to produce intensities greater than 15 W/cm(2) in the tissue target region.