Accurate simulation of transcranial ultrasound propagation for ultrasonic neuromodulation and stimulation

Model-based navigation of transcranial focused ultrasound neuromodulation in humans: Application to targeting the amygdala and thalamus

BACKGROUND

Transcranial focused ultrasound (tFUS) neuromodulation has shown promise in animals but is challenging to translate to humans because of the thicker skull that heavily scatters ultrasound waves.

OBJECTIVE

We develop and disseminate a model-based navigation (MBN) tool for acoustic dose delivery in the presence of skull aberrations that is easy to use by non-specialists.

METHODS

We pre-compute acoustic beams for thousands of virtual transducer locations on the scalp of the subject under study. We use the hybrid angular spectrum solver mSOUND, which runs in ∼4 s per solve per CPU yielding pre-computation times under 1 h for scalp meshes with up to 4000 faces and a parallelization factor of 5. We combine this pre-computed set of beam solutions with optical tracking, thus allowing real-time display of the tFUS beam as the operator freely navigates the transducer around the subject' scalp. We assess the impact of MBN versus line-of-sight targeting (LOST) positioning in simulations of 13 subjects.

RESULTS

Our navigation tool has a display refresh rate of ∼10 Hz. In our simulations, MBN increased the acoustic dose in the thalamus and amygdala by 8-67 % compared to LOST and avoided complete target misses that affected 10-20 % of LOST cases. MBN also yielded a lower variability of the deposited dose across subjects than LOST.

CONCLUSIONS

MBN may yield greater and more consistent (less variable) ultrasound dose deposition than transducer placement with line-of-sight targeting, and thus could become a helpful tool to improve the efficacy of tFUS neuromodulation.

Low-intensity pulsed ultrasound stimulation (LIPUS) modulates microglial activation following intracortical microelectrode implantation

Microglia are important players in surveillance and repair of the brain. Implanting an electrode into the cortex activates microglia, produces an inflammatory cascade, triggers the foreign body response, and opens the blood-brain barrier. These changes can impede intracortical brain-computer interfaces performance. Using two-photon imaging of implanted microelectrodes, we test the hypothesis that low-intensity pulsed ultrasound stimulation can reduce microglia-mediated neuroinflammation following the implantation of microelectrodes. In the first week of treatment, we found that low-intensity pulsed ultrasound stimulation increased microglia migration speed by 128%, enhanced microglia expansion area by 109%, and a reduction in microglial activation by 17%, indicating improved tissue healing and surveillance. Microglial coverage of the microelectrode was reduced by 50% and astrocytic scarring by 36% resulting in an increase in recording performance at chronic time. The data indicate that low-intensity pulsed ultrasound stimulation helps reduce the foreign body response around chronic intracortical microelectrodes.

A novel transcranial MR Guided focused ultrasound method based on the ultrashort echo time skull acoustic model and phase retrieval techniques

Transcranial ultrasound stimulation (TUS) has been clinically applied as a neuromodulation tool. Particularly, the phase array ultrasound can be applied in TUS to non-invasively focus on the cortex or deep brain. However, the vital phase distortion of the ultrasound induced by the skull limits its clinical application. In the current study, we aimed to develop a hybrid method that combines the ultrashort echo time (UTE) magnetic resonance imaging (MRI) sequences with the prDeep technique to achieve focusing ventral intermediate thalamic nucleus (VIM). The time-reversal (TR) approach of the UTE numerical acoustic model of the skull combined with the prDeep algorithm was used to reduce the number of iterations. The skull acoustic model simulation therapy process was establish to valid this method's prediction and focus performance, and the classical TR method were considered as the gold standard (GS). Our approach could restore 75% of the GS intensity in 25 iteration steps, with a superior the noise immunity. Our findings demonstrate that the phase aberration caused by the skull can be estimated using phase retrieval techniques to achieve a fast and accurate transcranial focus. The method has excellent adaptability and anti-noise capacity for satisfying complex and changeable scenarios.

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.

An efficient method for transcranial ultrasound focus correction based on the coupling of boundary integrals and finite elements

Transcranial focused ultrasound is a novel technique for the noninvasive treatment of brain diseases. The success of the treatment greatly depends on achieving precise and efficient intraoperative focus. However, compensating for aberrated ultrasound waves caused by the skull through numerical simulation-based phase corrections is a challenging task due to the significant computational burden involved in solving the acoustic wave equation. In this article, we propose a promising strategy using the coupling of the boundary integral equation method (BIEM) and the finite element method (FEM) to overcome the above limitation. Specifically, we adopt the BIEM to obtain the Robin-to-Dirichlet maps on the boundaries of the skull and then couple the maps to the FEM matrices via a dual interpolation technique, resulting in a computational domain including only the skull. Three simulation experiments were conducted to evaluate the effectiveness of the proposed method, including a convergence test and two skull-induced aberration corrections in 2D and 3D ultrasound. The results show that the method's convergence is guaranteed as the element size decreases, leading to a decrease in pressure error. The computation times for simulating a 500 kHz ultrasound field on a regular desktop computer were found to be 0.47 ± 0.01 s in the 2D case and 43.72 ± 1.49 s in the 3D case, provided that lower-upper decomposition (approximately 13 s in 2D and 2.5 h in 3D) was implemented in advance. We also demonstrated that more accurate transcranial focusing can be achieved by phase correction compared to the noncorrected results (with errors of 1.02 mm vs. 6.45 mm in 2D and 0.28 mm vs. 3.07 mm in 3D). The proposed strategy is valuable for enabling online ultrasound simulations during treatment, facilitating real-time adjustments and interventions.

Transcranial Phase Correction Using Pulse-Echo Ultrasound and Deep Learning: A 2-D Numerical Study

Phase aberration caused by human skulls severely degrades the quality of transcranial ultrasound images, posing a major challenge in the practical application of transcranial ultrasound techniques in adults. Aberration can be corrected if the skull profile (i.e., thickness distribution) and speed of sound (SOS) are known. However, accurately estimating the skull profile and SOS using ultrasound with a physics-based approach is challenging due to the complexity of the interaction between ultrasound and the skull. A deep learning approach is proposed herein to estimate the skull profile and SOS using ultrasound radiofrequency (RF) signals backscattered from the skull. A numerical study was performed to test the approach's feasibility. Realistic numerical skull models were constructed from computed tomography (CT) scans of five ex vivo human skulls in this numerical study. Acoustic simulations were performed on 3595 skull segments to generate array-based ultrasound backscattered signals. A deep learning model was developed and trained to estimate skull thickness and SOS from RF channel data. The trained model was shown to be highly accurate. The mean absolute error (MAE) was 0.15 mm (2% error) for thickness estimation and 13 m/s (0.5% error) for SOS estimation. The Pearson correlation coefficient between the estimated and ground-truth values was 0.99 for thickness and 0.95 for SOS. Aberration correction performed using deep-learning-estimated skull thickness and SOS values yielded significantly improved beam focusing (e.g., narrower beams) and transcranial imaging quality (e.g., improved spatial resolution and reduced artifacts) compared with no aberration correction. The results demonstrate the feasibility of the proposed approach for transcranial phase aberration correction.

Low-intensity pulsed ultrasound stimulation (LIPUS) modulates microglial activation following intracortical microelectrode implantation

Microglia are important players in surveillance and repair of the brain. Their activation mediates neuroinflammation caused by intracortical microelectrode implantation, which impedes the application of intracortical brain-computer interfaces (BCIs). While low-intensity pulsed ultrasound stimulation (LIPUS) can attenuate microglial activation, its potential to modulate the microglia-mediated neuroinflammation and enhance the bio-integration of microelectrodes remains insufficiently explored. We found that LIPUS increased microglia migration speed from 0.59±0.04 to 1.35±0.07 µm/hr on day 1 and enhanced microglia expansion area from 44.50±6.86 to 93.15±8.77 µm 2 /min on day 7, indicating improved tissue healing and surveillance. Furthermore, LIPUS reduced microglial activation by 17% on day 6, vessel-associated microglia ratio from 70.67±6.15 to 40.43±3.87% on day 7, and vessel diameter by 20% on day 28. Additionally, microglial coverage of the microelectrode was reduced by 50% in week 1, indicating better tissue-microelectrode integration. These data reveal that LIPUS helps resolve neuroinflammation around chronic intracortical microelectrodes.

Generating Patient-Specific Acoustic Simulations for Transcranial Focused Ultrasound Procedures Based on Optical Tracking Information

Optical tracking is a real-time transducer positioning method for transcranial focused ultrasound (tFUS) procedures, but the predicted focus from optical tracking typically does not incorporate subject-specific skull information. Acoustic simulations can estimate the pressure field when propagating through the cranium but rely on accurately replicating the positioning of the transducer and skull in a simulated space. Here, we develop and characterize the accuracy of a workflow that creates simulation grids based on optical tracking information in a neuronavigated phantom with and without transmission through an ex vivo skull cap. The software pipeline could replicate the geometry of the tFUS procedure within the limits of the optical tracking system (transcranial target registration error (TRE): 3.9 ± 0.7 mm). The simulated focus and the free-field focus predicted by optical tracking had low Euclidean distance errors of 0.5±0.1 and 1.2±0.4 mm for phantom and skull cap, respectively, and some skull-specific effects were captured by the simulation. However, the TRE of simulation informed by optical tracking was 4.6±0.2, which is as large or greater than the focal spot size used by many tFUS systems. By updating the position of the transducer using the original TRE offset, we reduced the simulated TRE to 1.1 ± 0.4 mm. Our study describes a software pipeline for treatment planning, evaluates its accuracy, and demonstrates an approach using MR-acoustic radiation force imaging as a method to improve dosimetry. Overall, our software pipeline helps estimate acoustic exposure, and our study highlights the need for image feedback to increase the accuracy of tFUS dosimetry.

Evaluation of synthetically generated computed tomography for use in transcranial focused ultrasound procedures

Purpose

Transcranial focused ultrasound (tFUS) is a therapeutic ultrasound method that focuses sound through the skull to a small region noninvasively and often under magnetic resonance imaging (MRI) guidance. CT imaging is used to estimate the acoustic properties that vary between individual skulls to enable effective focusing during tFUS procedures, exposing patients to potentially harmful radiation. A method to estimate acoustic parameters in the skull without the need for CT is desirable.

Approach

We synthesized CT images from routinely acquired T1-weighted MRI using a 3D patch-based conditional generative adversarial network and evaluated the performance of synthesized CT (sCT) images for treatment planning with tFUS. We compared the performance of sCT with real CT (rCT) images for tFUS planning using Kranion and simulations using the acoustic toolbox, k-Wave. Simulations were performed for 3 tFUS scenarios: (1) no aberration correction, (2) correction with phases calculated from Kranion, and (3) phase shifts calculated from time reversal.

Results

From Kranion, the skull density ratio, skull thickness, and number of active elements between rCT and sCT had Pearson's correlation coefficients of 0.94, 0.92, and 0.98, respectively. Among 20 targets, differences in simulated peak pressure between rCT and sCT were largest without phase correction (12.4 % ± 8.1 % ) and smallest with Kranion phases (7.3 % ± 6.0 % ). The distance between peak focal locations between rCT and sCT was < 1.3    mm for all simulation cases.

Conclusions

Real and synthetically generated skulls had comparable image similarity, skull measurements, and acoustic simulation metrics. Our work demonstrated similar results for 10 testing cases comparing MR-sCTs and rCTs for tFUS planning. Source code and a docker image with the trained model are available at https://github.com/han-liu/SynCT_TcMRgFUS.

A head template for computational dose modelling for transcranial focused ultrasound stimulation

Transcranial focused Ultrasound Stimulation (TUS) at low intensities is emerging as a novel non-invasive brain stimulation method with higher spatial resolution than established transcranial stimulation methods and the ability to selectively stimulate also deep brain areas. Accurate control of the focus position and strength of the TUS acoustic waves is important to enable a beneficial use of the high spatial resolution and to ensure safety. As the human skull causes strong attenuation and distortion of the waves, simulations of the transmitted waves are needed to accurately determine the TUS dose distribution inside the cranial cavity. The simulations require information of the skull morphology and its acoustic properties. Ideally, they are informed by computed tomography (CT) images of the individual head. However, suited individual imaging data is often not readily available. For this reason, we here introduce and validate a head template that can be used to estimate the average effects of the skull on the TUS acoustic wave in the population. The template was created from CT images of the heads of 29 individuals of different ages (between 20-50 years), gender and ethnicity using an iterative non-linear co-registration procedure. For validation, we compared acoustic and thermal simulations based on the template to the average of the simulation results of all 29 individual datasets. Acoustic simulations were performed for a model of a focused transducer driven at 500 kHz, placed at 24 standardized positions by means of the EEG 10-10 system. Additional simulations at 250 kHz and 750 kHz at 16 of the positions were used for further confirmation. The amount of ultrasound-induced heating at 500 kHz was estimated for the same 16 transducer positions. Our results show that the template represents the median of the acoustic pressure and temperature maps from the individuals reasonably well in most cases. This underpins the usefulness of the template for the planning and optimization of TUS interventions in studies of healthy young adults. Our results further indicate that the amount of variability between the individual simulation results depends on the position. Specifically, the simulated ultrasound-induced heating inside the skull exhibited strong interindividual variability for three posterior positions close to the midline, caused by a high variability of the local skull shape and composition. This should be taken into account when interpreting simulation results based on the template.

A phased array ultrasound system with a robotic arm for neuromodulation

BACKGROUND

Ultrasonic neuromodulation (UNMOD) provides a non-invasive brain stimulation. However, the high-resolution region-specificity of UNMOD with a single element transducer combined with a mechanical positioning system could have limits due to the intrinsic positioning error from mechanical systems.

OBJECTIVE/HYPOTHESIS

A phased array system could lead to highly selective neuromodulation with electronic control.

METHODS

A specialized phased-array system with a robotic arm is implemented for a rhesus monkey model. Various primary motor cortex areas related to tail, hand, and mouth were stimulated with a 200 μm step size. The ultrasonic parameters were ISPTA of 840 mW/cm2, pulse repetition frequency of 100 Hz, and a 5% duty factor at 600 kHz. The induced movement were recorded and analyzed.

RESULTS

Separate digits, mouth, and tongue motions were successfully induced by electronically controlling the focus. The identical body part movement could be induced when the focus was moved back to the identical primary motor cortex with electronic control. Accordingly, the reproducibility of UNMOD could be partially validated with rhesus monkey model.

CONCLUSION

A phased-array system appears to have a potential for the non-invasive and region-selective neuromodulation method.

Acoustoelectric Time-Reversal for Ultrasound Phase-Aberration Correction

Acoustoelectric imaging (AEI) is a technique that combines ultrasound (US) with radio frequency recording to detect and map local current source densities. This study demonstrates a new method called acoustoelectric time reversal (AETR), which uses AEI of a small current source to correct for phase aberrations through a skull or other US-aberrating layers with applications to brain imaging and therapy. Simulations conducted at three different US frequencies (0.5, 1.5, and 2.5 MHz) were performed through media layered with different sound speeds and geometries to induce aberrations of the US beam. Time delays of the acoustoelectric (AE) signal from a monopole within the medium were calculated for each element to enable corrections using AETR. Uncorrected aberrated beam profiles were compared with those after applying AETR corrections, which demonstrated a strong recovery (29%-100%) of lateral resolution and increases in focal pressure up to 283%. To further demonstrate the practical feasibility of AETR, we further conducted bench-top experiments using a 2.5 MHz linear US array to perform AETR through 3-D-printed aberrating objects. These experiments restored lost lateral restoration up to 100% for the different aberrators and increased focal pressure up to 230% after applying AETR corrections. Cumulatively, these results highlight AETR as a powerful tool for correcting focal aberrations in the presence of a local current source with applications to AEI, US imaging, neuromodulation, and therapy.

Full-Waveform Inversion for Breast Ultrasound Tomography Using Line-Shape Modeled Elements

OBJECTIVE

The objective of the work described here was to incorporate the spatial shapes of the transducer elements into the framework of the full-waveform inversion.

METHODS

An element is treated as its cross-section in the 2-D imaging plane, that is, a line segment. The elements are not simply modeled as a set of point sources on their surface to avoid staircasing artifacts. By use of the Fourier collocation method, an element is spatially represented as the discrete convolution between its spatial distribution and a band-limited delta function. The excitation pulses on the emitters and recorded signals on the receivers are then weighted based on the discrete convolution results. Digital and physical experiments are implemented to validate the method.

DISCUSSION

It is meaningful to model the shapes of the elements if their spatial sizes are similar to or larger than the acoustic wavelengths. It should, however, be noted that because this article focuses on 2-D imaging, the inter-plane effects are not considered.

CONCLUSION

The approach helps reduce the root mean square errors and increase the structural similarity of the reconstructed images. It also helps to improve the stability of convergence and to accelerate the convergence speed.

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.

Effect of skull porosity on ultrasound transmission and wave mode conversion at large incidence angles

BACKGROUND

Transcranial ultrasound imaging and therapy depend on the efficient transmission of acoustic energy through the skull. Multiple previous studies have concluded that a large incidence angle should be avoided during transcranial-focused ultrasound therapy to ensure transmission through the skull. Alternatively, some other studies have shown that longitudinal-to-shear wave mode conversion might improve transmission through the skull when the incidence angle is increased above the critical angle (i.e., 25° to 30°).

PURPOSE

The effect of skull porosity on the transmission of ultrasound through the skull at varying incidence angles was investigated for the first time to elucidate why transmission through the skull at large angles of incidence is decreased in some cases but improved in other cases.

METHODS

Transcranial ultrasound transmission at varying incidence angles (0°-50°) was investigated in phantoms and ex vivo skull samples with varying bone porosity (0% to 28.54% ± 3.36%) using both numerical and experimental methods. First, the elastic acoustic wave transmission through the skull was simulated using micro-computed tomography data of ex vivo skull samples. The trans-skull pressure was compared between skull segments having three levels of porosity, that is, low porosity (2.65% ± 0.03%), medium porosity (13.41% ± 0.12%), and high porosity (26.9%). Next, transmission through two 3D-printed resin skull phantoms (compact vs. porous phantoms) was experimentally measured to test the effect of porous microstructure alone on ultrasound transmission through flat plates. Finally, the effect of skull porosity on ultrasound transmission was investigated experimentally by comparing transmission through two ex vivo human skull segments having similar thicknesses but different porosities (13.78% ± 2.05% vs. 28.54% ± 3.36%).

RESULTS

Numerical simulations indicated that an increase in transmission pressure occurs at large incidence angles for skull segments having low porosities but not for those with high porosity. In experimental studies, a similar phenomenon was observed. Specifically, for the low porosity skull sample (13.78% ± 2.05%), the normalized pressure was 0.25 when the incidence angle increased to 35°. However, for the high porosity sample (28.54% ± 3.36%), the pressure was no more than 0.1 at large incidence angles.

CONCLUSIONS

These results indicate that the skull porosity has an evident effect on the transmission of ultrasound at large incidence angles. The wave mode conversion at large, oblique incidence angles could enhance the transmission of ultrasound through parts of the skull having lower porosity in the trabecular layer. However, for transcranial ultrasound therapy in the presence of highly porous trabecular bone, transmission at a normal incidence angle is preferable relative to oblique incidence angles due to the higher transmission efficiency.

Transcranial ultrasound simulations: A review

Transcranial ultrasound is more and more used for therapy and imaging of the brain. However, the skull is a highly attenuating and aberrating medium, with different structures and acoustic properties among samples and even within a sample. Thus, case-specific simulations are needed to perform transcranial focused ultrasound interventions safely. In this article, we provide a review of the different methods used to model the skull and to simulate ultrasound propagation through it.

Numerical and experimental evaluation of low-intensity transcranial focused ultrasound wave propagation using human skulls for brain neuromodulation

BACKGROUND

Low-intensity transcranial focused ultrasound (tFUS) has gained considerable attention as a promising noninvasive neuromodulatory technique for human brains. However, the complex morphology of the skull hinders scholars from precisely predicting the acoustic energy transmitted and the region of the brain impacted during the sonication. This is due to the fact that different ultrasound frequencies and skull morphology variations greatly affect wave propagation through the skull.

PURPOSE

Although the acoustic properties of human skull have been studied for tFUS applications, such as tumor ablation using a multielement phased array, there is no consensus about how to choose a single-element focused ultrasound (FUS) transducer with a suitable frequency for neuromodulation. There are interests in exploring the magnitude and dimension of tFUS beam through human parietal bone for modulating specific brain lobes. Herein, we aim to investigate the wave propagation of tFUS on human skulls to understand and address the concerns above.

METHODS

Both experimental measurements and numerical modeling were conducted to investigate the transmission efficiency and beam pattern of tFUS on five human skulls (C3 and C4 regions) using single-element FUS transducers with six different frequencies (150-1500 kHz). The degassed skull was placed in a water tank, and a calibrated hydrophone was utilized to measure acoustic pressure past it. The cranial computed tomography scan data of each skull were obtained to derive a high-resolution acoustic model (grid point spacing: 0.25 mm) in simulations. Meanwhile, we modified the power-law exponent of acoustic attenuation coefficient to validate numerical modeling and enabled it to be served as a prediction tool, based on the experimental measurements.

RESULTS

The transmission efficiency and -6 dB beamwidth were evaluated and compared for various frequencies. An exponential decrease in transmission efficiency and a logarithmic decrease of -6 dB beamwidth with an increase in ultrasound frequency were observed. It is found that a >750 kHz ultrasound leads to a relatively lower tFUS transmission efficiency (<5%), whereas a <350 kHz ultrasound contributes to a relatively broader beamwidth (>5 mm). Based on these observations, we further analyzed the dependence of tFUS wave propagation on FUS transducer aperture size.

CONCLUSIONS

We successfully studied tFUS wave propagation through human skulls at different frequencies experimentally and numerically. The findings have important implications to predict tFUS wave propagation for ultrasound neuromodulation in clinical applications, and guide researchers to develop advanced ultrasound transducers as neural interfaces.

Pulse-Echo Technique to Compensate for Laminate Membrane Transmission Loss in Phantom-Based Ultrasonic Attenuation Coefficient Measurements

OBJECTIVES

Accurately measuring the attenuation coefficient (AC) of reference phantoms is critical in clinical applications of quantitative ultrasound. Phantom AC measurement requires proper compensation of membrane transmission loss. Conventional methods require separate membrane samples to obtain membrane transmission loss. Unfortunately, separate membrane samples are often unavailable. A pulse-echo approach is proposed herein to compensate for membrane transmission loss without requiring separate membrane samples.

METHODS

The proposed method consists of the following steps. First, the insertion loss, caused by phantom attenuation and membrane transmission loss, is measured. Second, the membrane reflection coefficient is measured. Third, the unknown acoustic parameters of the membrane and phantom material are estimated by fitting theoretical reflection coefficient to the measured one. Finally, the fitted parameters are used to estimate membrane transmission loss and phantom AC. The proposed method was validated through k-Wave simulations and phantom experiments. Experimental AC measurements were repeated on 5 distinct phantoms by 2 operators to assess the repeatability and reproducibility of the proposed method. Five transducers were used to cover a broad bandwidth (0.7-16 MHz).

RESULTS

The acquired AC in the simulations had a maximum error of 0.06 dB/cm-MHz for simulated phantom AC values ranging from 0.5 to 1 dB/cm-MHz. The acquired AC in the experiments had a maximum error of 0.045 dB/cm-MHz for phantom AC values ranging from 0.28 to 1.48 dB/cm-MHz. Good repeatability and cross-operator reproducibility were observed with a mean coefficient of variation below 0.054.

CONCLUSION

The proposed method simplifies phantom AC measurement while providing satisfactory accuracy and precision.

S.M.A.R.T. F.U.S: Surrogate Model of Attenuation and Refraction in Transcranial Focused Ultrasound

Low-intensity focused ultrasound (LIFU) is an increasingly applied method for achieving non-invasive brain stimulation. However, transmission of ultrasound through the human skull can substantially affect focal point characteristics of LIFU, including dramatic attenuation in intensity and refraction of focal point location. These effects depend on a high-dimensional parameter space, making these effects difficult to estimate from previous work. Instead, focal point properties of LIFU experiments are often estimated using numerical simulation of LIFU sonication through skull. However, this procedure presents many entry barriers to even computationally savvy investigators and often requires expensive computational hardware, impeding LIFU research. We present a novel MATLAB toolbox (data: doi:10.5068/D1QD60; Matlab Scripts: https://doi.org/10.5281/zenodo.5811122) for rapidly estimating beam properties of LIFU transmitted through bone. Users provide specific values for frequency of LIFU, bone thickness, angle at which LIFU is applied, depth of the LIFU focal point, and diameter of the transducer used and receive an estimation of the degree of refraction/attenuation expected for the given parameters.

Comparison between MR and CT imaging used to correct for skull-induced phase aberrations during transcranial focused ultrasound

Transcranial focused ultrasound with the InSightec Exablate system uses thermal ablation for the treatment of movement and mood disorders and blood brain barrier disruption for tumor therapy. The system uses computed tomography (CT) images to calculate phase corrections that account for aberrations caused by the human skull. This work investigates whether magnetic resonance (MR) images can be used as an alternative to CT images to calculate phase corrections. Phase corrections were calculated using the gold standard hydrophone method and the standard of care InSightec ray tracing method. MR binary image mask, MR-simulated-CT (MRsimCT), and CT images of three ex vivo human skulls were supplied as inputs to the InSightec ray tracing method. The degassed ex vivo human skulls were sonicated with a 670 kHz hemispherical phased array transducer (InSightec Exablate 4000). 3D raster scans of the beam profiles were acquired using a hydrophone mounted on a 3-axis positioner system. Focal spots were evaluated using six metrics: pressure at the target, peak pressure, intensity at the target, peak intensity, positioning error, and focal spot volume. Targets at the geometric focus and 5 mm lateral to the geometric focus were investigated. There was no statistical difference between any of the metrics at either target using either MRsimCT or CT for phase aberration correction. As opposed to the MRsimCT, the use of CT images for aberration correction requires registration to the treatment day MR images; CT misregistration within a range of ± 2 degrees of rotation error along three dimensions was shown to reduce focal spot intensity by up to 9.4%. MRsimCT images used for phase aberration correction for the skull produce similar results as CT-based correction, while avoiding both CT to MR registration errors and unnecessary patient exposure to ionizing radiation.

Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models

Computational models of acoustic wave propagation are frequently used in transcranial ultrasound therapy, for example, to calculate the intracranial pressure field or to calculate phase delays to correct for skull distortions. To allow intercomparison between the different modeling tools and techniques used by the community, an international working group was convened to formulate a set of numerical benchmarks. Here, these benchmarks are presented, along with intercomparison results. Nine different benchmarks of increasing geometric complexity are defined. These include a single-layer planar bone immersed in water, a multi-layer bone, and a whole skull. Two transducer configurations are considered (a focused bowl and a plane piston operating at 500 kHz), giving a total of 18 permutations of the benchmarks. Eleven different modeling tools are used to compute the benchmark results. The models span a wide range of numerical techniques, including the finite-difference time-domain method, angular spectrum method, pseudospectral method, boundary-element method, and spectral-element method. Good agreement is found between the models, particularly for the position, size, and magnitude of the acoustic focus within the skull. When comparing results for each model with every other model in a cross-comparison, the median values for each benchmark for the difference in focal pressure and position are less than 10% and 1 mm, respectively. The benchmark definitions, model results, and intercomparison codes are freely available to facilitate further comparisons.

Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models

Computational models of acoustic wave propagation are frequently used in transcranial ultrasound therapy, for example, to calculate the intracranial pressure field or to calculate phase delays to correct for skull distortions. To allow intercomparison between the different modeling tools and techniques used by the community, an international working group was convened to formulate a set of numerical benchmarks. Here, these benchmarks are presented, along with intercomparison results. Nine different benchmarks of increasing geometric complexity are defined. These include a single-layer planar bone immersed in water, a multi-layer bone, and a whole skull. Two transducer configurations are considered (a focused bowl and a plane piston operating at 500 kHz), giving a total of 18 permutations of the benchmarks. Eleven different modeling tools are used to compute the benchmark results. The models span a wide range of numerical techniques, including the finite-difference time-domain method, angular spectrum method, pseudospectral method, boundary-element method, and spectral-element method. Good agreement is found between the models, particularly for the position, size, and magnitude of the acoustic focus within the skull. When comparing results for each model with every other model in a cross-comparison, the median values for each benchmark for the difference in focal pressure and position are less than 10% and 1 mm, respectively. The benchmark definitions, model results, and intercomparison codes are freely available to facilitate further comparisons.

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.

Two-step aberration correction: application to transcranial histotripsy

Objective: Phase aberration correction is essential in transcranial histotripsy to compensate for focal distortion caused by the heterogeneity of the intact skull bone. This paper improves the 2-step aberration correction (AC) method that has been previously presented and develops an AC workflow that fits in the clinical environment, in which the computed tomography (CT)-based analytical approach was first implemented, followed by a cavitation-based approach using the shockwaves from the acoustic cavitation emission (ACE).Approach:A 700 kHz, 360-element hemispherical transducer array capable of transmit-and-receive on all channels was used to transcranially generate histotripsy-induced cavitation and acquire ACE shockwaves. For CT-AC, two ray-tracing models were investigated: a forward ray-tracing model (transducer-to-focus) in the open-source software Kranion, and an in-house backward ray-tracing model (focus-to-transducer) accounting for refraction and the sound speed variation in skulls. Co-registration was achieved by aligning the skull CT data to the skull surface map reconstructed using the acoustic pulse-echo method. For ACE-AC, the ACE signals from the collapses of generated bubbles were aligned by cross-correlation to estimate the corresponding time delays.Main results:The performance of the 2-step method was tested with 3 excised human calvariums placed at 2 different locations in the transducer array. Results showed that the 2-step AC achieved 90 ± 7% peak focal pressure compared to the gold standard hydrophone correction. It also reduced the focal shift from 0.84 to 0.30 mm and the focal volume from 10.6 to 2.0 mm3on average compared to the no AC cases.Significance:The 2-step AC yielded better refocusing compared to either CT-AC or ACE-AC alone and can be implemented in real-time for transcranial histotripsy brain therapy.

Ultrasound stimulation of the motor cortex during tonic muscle contraction

Transcranial ultrasound stimulation (tUS) shows potential as a noninvasive brain stimulation (NIBS) technique, offering increased spatial precision compared to other NIBS techniques. However, its reported effects on primary motor cortex (M1) are limited. We aimed to better understand tUS effects in human M1 by performing tUS of the hand area of M1 (M1hand) during tonic muscle contraction of the index finger. Stimulation during muscle contraction was chosen because of the transcranial magnetic stimulation-induced phenomenon known as cortical silent period (cSP), in which transcranial magnetic stimulation (TMS) of M1hand involuntarily suppresses voluntary motor activity. Since cSP is widely considered an inhibitory phenomenon, it presents an ideal parallel for tUS, which has often been proposed to preferentially influence inhibitory interneurons. Recording electromyography (EMG) of the first dorsal interosseous (FDI) muscle, we investigated effects on muscle activity both during and after tUS. We found no change in FDI EMG activity concurrent with tUS stimulation. Using single-pulse TMS, we found no difference in M1 excitability before versus after sparsely repetitive tUS exposure. Using acoustic simulations in models made from structural MRI of the participants that matched the experimental setups, we estimated in-brain pressures and generated an estimate of cumulative tUS exposure experienced by M1hand for each subject. We were unable to find any correlation between cumulative M1hand exposure and M1 excitability change. We also present data that suggest a TMS-induced MEP always preceded a near-threshold cSP.

Increasing the transmission efficiency of transcranial ultrasound using a dual-mode conversion technique based on Lamb waves

Transcranial focused ultrasound (FUS) is a noninvasive treatment for brain tumors and neuromodulation. Based on normal incidence, conventional FUS techniques use a focused or an array of ultrasonic transducers to overcome the attenuation and absorption of ultrasound in the skull; however, this remains the main limitation of using FUS. A dual-mode conversion technique based on Lamb waves is proposed to achieve high transmission efficiency. This concept was validated using the finite element analysis (FEA) and experiments based on changes in the incident angle. Aluminum, plexiglass, and a human skull were used as materials with different attenuations. The transmission loss was calculated for each material, and the results were compared with the reflectance function of the Lamb waves. Oblique incidence based on dual-mode conversion exhibited a better transmission efficiency than that of a normal incidence for all of the specimens. The total transmission losses for the materials were 13.7, 15.46, and 3.91 dB less than those associated with the normal incidence. A wedge transducer was designed and fabricated to implement the proposed method. The results demonstrated the potential applicability of the dual-mode conversion technique for the human skull.

A Spatial Multitarget Ultrasound Neuromodulation System Using High-Powered 2-D Array Transducer

Transcranial focused ultrasound (tFUS) is increasingly used in experimental neuroscience due to its neuromodulatory effectiveness in animal studies. However, achieving multitarget tFUS in small animals is typically limited by transducer size, energy transfer efficiency, and brain volume. The objective of this work was to construct an ultrasound system for multitarget neuromodulation in small animals. First, a miniaturized high-powered 2-D array transducer was developed. The phase delay of each array element was calculated based on the multifocal time-reversal method, generating multiple foci simultaneously in a 3-D field. The effects of the axial focal length, interfocus spacing (lateral distance between the two focal centers), and the number of foci on the focal properties of the pressure field were examined through numerical simulations. In-vitro ultrasonic measurements and transcranial simulations on a rat skull were conducted. The minimum interfocus spacing separating two -6-dB foci and the peak full-width at half-maximum were positively correlated with axial focal length; the relative relationship between the interfocus spacing and pressure field properties was similar for each axial focal length. The maximum acoustic pressure and spatial average intensity at focus in deionized water were 2.21 MPa and 133 W/cm², respectively. The simulated and experimental results were compared, demonstrating agreement in both peak position and focus shape. The ultrasound system can provide a neuroscientific platform for evaluating the feasibility of multitarget ultrasound stimulation treatment protocols, thus improving the understanding of functional neuroanatomy.

Multitarget Transcranial Ultrasound Therapy in Small Animals Based on Phase-Only Acoustic Holographic Lens

Transcranial ultrasound therapy has become a noninvasive method for treating neurological and psychiatric disorders, and studies have further demonstrated that multitarget transcranial ultrasound therapy is a better solution. At present, multitarget transcranial ultrasound therapy in small animals can only be achieved by the multitransducer or phased array. However, multiple transducers may cause spatial interference, and the phased array system is complicated, expensive, and especially unsuitable for small animals. This study is the first to design and fabricate a miniature acoustic holography lens for multitarget transcranial ultrasound therapy in rats. The acoustic holographic lens, working at a frequency of 1.0 MHz, with a size of 10.08 mm ×10.08 mm and a pixel resolution of 0.72 mm, was designed, optimized, and fabricated. The dual-focus transcranial ultrasound generated based on the lens was measured; the full-width at half-maximum (FWHM) of the focal spots in the y -direction was 2.15 and 2.27 mm and in the z -direction was 2.3 and 2.36 mm. The focal length was 5.4 mm, and the distance between the two focuses was 5.6 mm, close to the desired values of 5.4 and 6.0 mm. Finally, the multiple-target blood-brain barrier opening in rats' bilateral secondary visual cortex (mediolateral area, V2ML) was demonstrated using the transcranial ultrasound therapy system based on the lens. These results demonstrate the good performance of the multitarget transcranial ultrasound therapy system for small animals, including high spatial resolution, small size, and low cost.

Transcranial Theranostic Ultrasound for Pre-Planning and Blood-Brain Barrier Opening: A Feasibility Study Using an Imaging Phased Array In Vitro and In Vivo

Focused ultrasound (FUS) for blood-brain barrier (BBB) opening is a safe, reversible and non-invasive strategy for targeted drug delivery to the brain, however extensive pre-planning strategies are necessary for successful FUS-mediated BBB opening through the structurally complex primate skull.

OBJECTIVE

This work aims to demonstrate a pre-planning pipeline consisting of transcranial simulations and in vitro experimentation used to inform synchronous BBB opening and power cavitation imaging (PCI) with a single theranostic ultrasound (TUS) phased array.

METHODS

Acoustic wave propagation simulation readouts of pressure attenuation and focal shift through clinical-CT and micro-CT-based primate skull models were compared, while the latter were used to determine the impact of beam steering angle on focal shift and pressure attenuation. In vitro experimentation with a channel phantom enabled characterization of skull-induced receive focal shift (RFS), while in vivo BBB opening and PCI using in silico and in vitro pre-planning information was conducted using a custom Verasonics/MATLAB script.

RESULTS

Simulations confirmed steering angle dependent transcranial focal shift and pressure attenuation, while in vitro experiments revealed minimal (0.30-1.50 mm) skull-induced RFS. In vivo rodent experiments with overlaid primate skull fragments demonstrated successful TUS-mediated BBB opening and spatially correlated power cavitation images (PCI) with regions of BBB opening on T1-weighted magnetic resonance images (MRI).

CONCLUSION

Herein, we demonstrate the feasibility for TUS-mediated BBB opening in vivo using in silico and in vitro pre-planning information.

SIGNIFICANCE

TUS as an ultrasound-guided modality for BBB opening could be a promising alternative to current FUS-mediated BBB opening configurations in the clinic.

Prognostic value of brain tissues' volumes in patients with essential tremor treated with MRgFUS thalamotomy

MRgFUS Vim thalamotomy is a novel, effective, minimally invasive therapeutic option for patients with essential tremor (ET). Among the selection criteria, some parameters related to the patient's anatomy, such as the skull density ratio (SDR), are well recognized. The role of brain tissue interposed between the target and the ultrasound transducers has never been explored. Therefore, the purpose of our study was to evaluate the correlation and the possible predictive value between brain tissue volumes (grey matter - GM, white matter - WM, and cerebrospinal fluid - CSF) and several treatment-related variables (periprocedural parameters, MRI imaging findings, and the clinical outcome). We analysed data from thirty ET patients previously submitted to MRgFUS thalamotomy. Pre-treatment images were automatically segmented in sopra-tentorial (ST) WM, GM, and CSF using SPM 12. The most significant findings were a positive correlation of the ST-GM with the Accumulated Thermal Dose (ATD) (p < 0,001) and a negative correlation of the ATD temperature with ST-CSF and ST-TIV (p < 0,001). Ultrasound propagation speed is lower in fluids than brain tissues. Also, WM has an attenuation rate of 1.5 higher than the GM. Therefore, the difference in the ATD may be explained by the different acoustic properties of normal brain tissues interposed between the transducers and the VIM.

"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.

Simulation and performance study of circular ultrasonic array for tubes' internal inspection

A circular array (CA) transducer is developed and validated to implement internal inspection of the small diameter tubes, one of whose potential application targets is the heat exchanger tubes in Modular High Temperature Gas-cooled Reactor. A geometric model of the transducer and the delay law are proposed according to the inspection setup. Using k-Wave, an open source acoustic simulation toolbox, the beam profiles of various aperture sizes are analyzed to select a suitable aperture. The beams with different focal depths are compared to verify the delay law and to choose a preferable focus. Reflector responses are simulated to predict the transducer's detection sensitivity. To validate the simulation results, a prototype transducer composed of 64 elements is fabricated. Its electrical and acoustical performances are acquired with a pulser-receiver. A 4-element aperture with the focal depth of 1.0 mm is selected to perform internal inspection experiment according to the simulation results. The eight longitudinal grooves with widths of 0.2, 0.3, 0.4, and 0.5 mm on inner and outer tube walls are applied as typical reflectors to test circumferential detection sensitivity. The other five grooves with depths from 0.1 mm to 0.5 mm are chosen to verify radial detection sensitivity. The experiment results demonstrate the detection sensitivity of the CA transducer reaches at least 0.2 mm in circumferential direction and 0.1 mm in radial direction. The CA is proven capable of inspecting small diameter tubes.

Effect of Skull Porous Trabecular Structure on Transcranial Ultrasound Imaging in the Presence of Elastic Wave Mode Conversion at Varying Incidence Angle

With the advancement of aberration correction techniques, transcranial ultrasound imaging has exhibited great potential in applications such as imaging neurological function and guiding therapeutic ultrasound. However, the feasibility of transcranial imaging varies among individuals because of the differences in skull acoustic properties. To better understand the fundamental mechanisms underlying the variation in imaging performance, the effect of the structure of the porous trabecular bone on transcranial imaging performance (i.e., target localization errors and resolution) was investigated for the first time through the use of elastic wave simulations and experiments. Simulation studies using high-resolution computed tomography data from ex vivo skull samples revealed that imaging at large incidence angles reduced the target localization error for skulls having low porosity; however, as skull porosity increased, large angles of incidence resulted in degradation of resolution and increased target localization errors. Experimental results indicate that imaging at normal incidence introduced a localization error of 1.85 ± 0.10 mm, while imaging at a large incidence angle (40°) resulted in an increased localization error of 6.54 ± 1.33 mm and caused a single point target to no longer appear as a single, coherent target in the resulting image, which is consistent with simulation results. This first investigation of the effects of skull microstructure on transcranial ultrasound imaging indicates that imaging performance is highly dependent on the porosity of the skull, particularly at non-normal angles of incidence.

A Generalized Split-Step Angular Spectrum Method for Efficient Simulation of Wave Propagation in Heterogeneous Media

Angular spectrum (AS) methods enable efficient calculation of wave propagation from one plane to another inside homogeneous media. For wave propagation in heterogeneous media such as biological tissues, AS methods cannot be applied directly. Split-stepping techniques decompose the heterogeneous domain into homogeneous and perturbation parts, and provide a solution for forward wave propagation by propagating the incident wave in both frequency-space and frequency-wavenumber domains. Recently, a split-step hybrid angular spectrum (HAS) method was proposed for plane wave propagation of focused ultrasound beams. In this study, we extend these methods to enable simulation of acoustic pressure field for an arbitrary source distribution, by decomposing the source and reflection spectra into orthogonal propagation direction components, propagating each component separately, and summing all components to get the total field. We show that our method can efficiently simulate the pressure field of arbitrary sources in heterogeneous media. The accuracy of the method was analyzed comparing the resultant pressure field with pseudospectral time domain (PSTD) solution for breast tomography and hemispherical transcranial-focused ultrasound simulation models. Eighty times acceleration was achieved for a 3-D breast simulation model compared to PSTD solution with 0.005 normalized root mean-squared difference (NRMSD) between two solutions. For the hemispherical phased array, aberrations due to skull were accurately calculated in a single simulation run as evidenced by the resultant-focused ultrasound beam simulations, which had 0.001 NRMSD with 40 times acceleration factor compared to the PSTD method.

Discrete Approximations of Acoustic Source Distributions

The modeling of source a distributions of finite spatial extent in ultrasound and medical imaging applications is a problem of longstanding interest. In time-domain methods, such as the finite-difference time-domain or pseudospectral approaches, one requirement is the representation of such distributions over a grid, normally Cartesian. Various artifacts, including staircasing errors, can arise. In this short contribution, the problem of the representation of distribution over a grid is framed as an optimization problem in the Fourier domain over a preselected set of grid points, thus maintaining control over computational cost and allowing the fine-tuning of the optimization to the wavenumber range of interest for a particular numerical method. Numerical results are presented in the important special case of the spherical cap or bowl source.

Comparison Between Ray-Tracing and Full-Wave Simulation for Transcranial Ultrasound Focusing on a Clinical System Using the Transfer Matrix Formalism

Only one high-intensity focused ultrasound device has been clinically approved for transcranial brain surgery at the time of writing. The device operates within 650 and 720 kHz and corrects the phase distortions induced by the skull of each patient using a multielement phased array. Phase correction is estimated adaptively using a proprietary algorithm based on computed-tomography (CT) images of the patient's skull. In this article, we assess the performance of the phase correction computed by the clinical device and compare it to: 1) the correction obtained with a previously validated full-wave simulation algorithm using an open-source pseudo-spectral toolbox and 2) a hydrophone-based correction performed invasively to measure the aberrations induced by the skull at 650 kHz. For the full-wave simulation, three different mappings between CT Hounsfield units and the longitudinal speed of sound inside the skull were tested. All methods are compared with the exact same setup due to transfer matrices acquired with the clinical system for N = 5 skulls and T = 2 different targets for each skull. We show that the clinical ray-tracing software and the full-wave simulation restore, respectively, 84% ± 5% and 86% ± 5% of the pressure obtained with hydrophone-based correction for targets located in central brain regions. On the second target (off-center), we also report that the performance of both algorithms degrades when the average incident angles of the acoustic beam at the skull surface increase. When incident angles are higher than 20°, the restored pressure drops below 75% of the pressure restored with hydrophone-based correction.

Safety of focused ultrasound neuromodulation in humans with temporal lobe epilepsy

OBJECTIVE

Transcranial Focused Ultrasound (tFUS) is a promising new potential neuromodulation tool. However, the safety of tFUS neuromodulation has not yet been assessed adequately. Patients with refractory temporal lobe epilepsy electing to undergo an anterior temporal lobe resection present a unique opportunity to evaluate the safety and efficacy of tFUS neuromodulation. Histological changes in tissue after tFUS can be examined after surgical resection, while further potential safety concerns can be assessed using neuropsychological testing.

METHODS

Neuropsychological functions were assessed in eight patients before and after focused ultrasound sonication of the temporal lobe at intensities up to 5760 mW/cm². Using the BrainSonix Pulsar 1002, tFUS was delivered under MR guidance, using the Siemens Magnetom 3T Prisma scanner. Neuropsychological changes were assessed using various batteries. Histological changes were assessed using hematoxylin and eosin staining, among others.

RESULTS

With respect to safety, the histological analysis did not reveal any detectable damage to the tissue, except for one subject for whom the histology findings were inconclusive. In addition, neuropsychological testing did not show any statistically significant changes in any test, except for a slight decrease in performance on one of the tests after tFUS.

SIGNIFICANCE

This study supports the hypothesis that low-intensity Transcranial Focused Ultrasound (tFUS) used for neuromodulation of brain circuits at intensities up to 5760 mW/cm² may be safe for use in human research. However, due to methodological limitations in this study and inconclusive findings, more work is warranted to establish the safety. Future directions include greater number of sonications as well as longer exposure at higher intensity levels to further assess the safety of tFUS for modulation of neuronal circuits.

mSOUND: An Open Source Toolbox for Modeling Acoustic Wave Propagation in Heterogeneous Media

mSOUND is an open-source toolbox written in MATLAB. This toolbox is intended for modeling linear/ nonlinear acoustic wave propagation in media (primarily biological tissues) with arbitrary heterogeneities, in which, the speed of sound, density, attenuation coefficient, power-law exponent, and nonlinear coefficient are all spatially varying functions. The computational model is an iterative one-way model based on a mixed domain method. In this article, a general guideline is given along with three representative examples to illustrate how to set up simulations using mSOUND. The first example uses the transient mixed-domain method (TMDM) forward projection to compute the transient acoustic field for a given source defined on a plane. The second example uses the frequency-specific mixed-domain method (FSMDM) forward projection to rapidly obtain the pressure distribution directly at the frequencies of interest, assuming linear or weakly nonlinear wave propagation. The third example demonstrates how to use TMDM backward projection to reconstruct the initial acoustic pressure field to facilitate photoacoustic tomography (PAT). mSOUND (https://m-sound.github.io/mSOUND/home) is designed to be complementary to existing ultrasound modeling toolboxes and is expected to be useful for a wide range of applications in medical ultrasound including treatment planning, PAT, transducer design, and characterization.

Transcranial focused ultrasound phase correction using the hybrid angular spectrum method

The InSightec Exablate system is the standard of care used for transcranial focused ultrasound ablation treatments in the United States. The system calculates phase corrections that account for aberrations caused by the human skull. This work investigates whether skull aberration correction can be improved by comparing the standard of care InSightec ray tracing method with the hybrid angular spectrum (HAS) method and the gold standard hydrophone method. Three degassed ex vivo human skulls were sonicated with a 670 kHz hemispherical phased array transducer (InSightec Exablate 4000). Phase corrections were calculated using four different methods (straight ray tracing, InSightec ray tracing, HAS, and hydrophone) and were used to drive the transducer. 3D raster scans of the beam profiles were acquired using a hydrophone mounted on a 3-axis positioner system. Focal spots were evaluated using six metrics: pressure at the target, peak pressure, intensity at the target, peak intensity, positioning error, and focal spot volume. For three skulls, the InSightec ray tracing method achieved 52 ± 21% normalized target intensity (normalized to hydrophone), 76 ± 17% normalized peak intensity, and 0.72 ± 0.47 mm positioning error. The HAS method achieved 74 ± 9% normalized target intensity, 81 ± 9% normalized peak intensity, and 0.35 ± 0.09 mm positioning error. The InSightec-to-HAS improvement in focal spot targeting provides promise in improving treatment outcomes. These improvements to skull aberration correction are also highly relevant for the applications of focused ultrasound neuromodulation and blood brain barrier opening, which are currently being translated for human use.

Suppressing Grating Lobes for Transcranial Focused Ultrasound System by Frequency-Modulated Excitation

Transcranial focused ultrasound is a novel noninvasive therapeutic modality for glioblastoma and other disorders of the brain. However, because the phase aberrations caused by the skull need to be corrected with computed tomography (CT) images, the transcranial transducer is tightly fixed on the patient's head to avoid any variation in the relative position, and the focus shifting relies mainly on the capacity for electronic beam steering. Due to the presence of grating lobes and the rapid degradation of the focus quality with increasing focus-shifting distance, transcranial focus-shifting sonication may damage healthy brain tissue unintentionally. To reduce the risks associated with transcranial focused ultrasound therapy, linear frequency-modulated (FM) excitation is proposed. The k-space corrected pseudospectral time domain (PSTD) and acoustic holography approach based on the Rayleigh integral are combined to calculate the distribution of the deposited acoustic power. The corresponding simulation was performed with axial/lateral focus shifting at different distances. The distributions of the deposited acoustic power show that linear FM excitation can effectively suppress undesired prefocal grating lobes without compromising focus quality.

Transcranial focused ultrasound stimulation with high spatial resolution

BACKGROUND

Low-intensity transcranial focused ultrasound stimulation is a promising candidate for noninvasive brain stimulation and accurate targeting of brain circuits because of its focusing capability and long penetration depth. However, achieving a sufficiently high spatial resolution to target small animal sub-regions is still challenging, especially in the axial direction.

OBJECTIVE

To achieve high axial resolution, we designed a dual-crossed transducer system that achieved high spatial resolution in the axial direction without complex microfabrication, beamforming circuitry, and signal processing.

METHODS

High axial resolution was achieved by crossing two ultrasound beams of commercially available piezoelectric curved transducers at the focal length of each transducer. After implementation of the fixture for the dual-crossed transducer system, three sets of in vivo animal experiments were conducted to demonstrate high target specificity of ultrasound neuromodulation using the dual-crossed transducer system (n = 38).

RESULTS

The full-width at half maximum (FWHM) focal volume of our dual-crossed transducer system was under 0.52 μm³. We report a focal diameter in both lateral and axial directions of 1 mm. To demonstrate successful in vivo brain stimulation of wild-type mice, we observed the movement of the forepaws. In addition, we targeted the habenula and verified the high spatial specificity of our dual-crossed transducer system.

CONCLUSIONS

Our results demonstrate the ability of the dual-crossed transducer system to target highly specific regions of mice brains using ultrasound stimulation. The proposed system is a valuable tool to study the complex neurological circuitry of the brain noninvasively.

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.

Design, Development, and Operation of a Low-Intensity Focused Ultrasound Pulsation (LIFUP) System for Clinical Use

Noninvasive low-intensity focused ultrasound pulsation (LIFUP) neuromodulation provides a unique approach to both investigating and treating the brain. This work describes a well-calibrated, simple-to-use ultrasound stimulation system for neuromodulation studies. It provides a single-element 650-kHz transducer design and a straightforward control mechanism, with extensive calibration and internal electronic monitoring to prevent unwanted over or under treatment. One goal of this approach is to relieve researchers of many of the details associated with developing their own exposure equipment. A unique transducer positioning system and distinctive MRI fiducial targets simplify alignment and targeting. The system design, control software, calibration, and alignment techniques are described in detail. Examples of transducer targeting and neurostimulation using the system are provided.

Modeling continuous source distributions in wave-based virtual acoustics

All acoustic sources are of finite spatial extent. In volumetric wave-based simulation approaches (including, e.g., the finite difference time domain method among many others), a direct approach is to represent such continuous source distributions in terms of a collection of point-like sources at grid locations. Such a representation requires interpolation over the grid and leads to common staircasing effects, particularly under rotation or translation of the distribution. In this article, a different representation is shown, based on a spherical harmonic representation of a given distribution. The source itself is decoupled from any particular arrangement of grid points, and is compactly represented as a series of filter responses used to drive a canonical set of source terms, each activating a given spherical harmonic directivity pattern. Such filter responses are derived for a variety of commonly encountered distributions. Simulation results are presented, illustrating various features of such a representation, including convergence, behaviour under rotation, the extension to the time varying case, and differences in computational cost relative to standard grid-based source representations.

Transcranial magnetic stimulation, deep brain stimulation, and other forms of neuromodulation for substance use disorders: Review of modalities and implications for treatment

Given the high prevalence of individuals diagnosed with substance use disorder, along with the elevated rate of relapse following treatment initiation, investigating novel approaches and new modalities for substance use disorder treatment is of vital importance. One such approach involves neuromodulation which has been used therapeutically for neurological and psychiatric disorders and has demonstrated positive preliminary findings for the treatment of substance use disorder. The following article provides a review of several forms of neuromodulation which warrant consideration as potential treatments for substance use disorder. PubMed, PsycINFO, Ovid MEDLINE, and Web of Science were used to identify published articles and clinicaltrials.gov was used to identify currently ongoing or planned studies. Search criteria for Brain Stimulation included the following terminology: transcranial direct current stimulation, transcranial magnetic stimulation, theta burst stimulation, deep brain stimulation, vagus nerve stimulation, trigeminal nerve stimulation, percutaneous nerve field stimulation, auricular nerve stimulation, and low intensity focused ultrasound. Search criteria for Addiction included the following terminology: addiction, substance use disorder, substance-related disorder, cocaine, methamphetamine, amphetamine, alcohol, nicotine, tobacco, smoking, marijuana, cannabis, heroin, opiates, opioids, and hallucinogens. Results revealed that there are currently several forms of neuromodulation, both invasive and non-invasive, which are being investigated for the treatment of substance use disorder. Preliminary findings have demonstrated the potential of these various neuromodulation techniques in improving substance treatment outcomes by reducing those risk factors (e.g. substance craving) associated with relapse. Specifically, transcranial magnetic stimulation has shown the most promise with several well-designed studies supporting the potential for reducing substance craving. Deep brain stimulation has also shown promise, though lacks well-controlled clinical trials to support its efficacy. Transcranial direct current stimulation has also demonstrated promising results though consistently designed, randomized trials are also needed. There are several other forms of neuromodulation which have not yet been investigated clinically but warrant further investigation given their mechanisms and potential efficacy based on findings from other studied indications. In summary, given promising findings in reducing substance use and craving, neuromodulation may provide a non-pharmacological option as a potential treatment and/or treatment augmentation for substance use disorder. Further research investigating neuromodulation, both alone and in combination with already established substance use disorder treatment (e.g. medication treatment), warrants consideration.

An efficient and accurate parallel hybrid acoustic signal correction method for transcranial ultrasound

Phased-control focused ultrasound transducers provide a new and noninvasive treatment method for brain disease. However, improving the accuracy of phase correction and reducing the calculation time during treatment have always been contradictory constraints. In this paper, a hybrid acoustic signal correction (HASC) method combined with k-Wave stage and holography stage was introduced for phase correction and simulation of transcranial focused ultrasound. The k-Wave stage is mainly used to calculate the sound field in a heterogeneous medium (skull), which divides the sound field calculation process into paths that can be calculated in parallel, and the transcranial correction phase can also be obtained during the calculation. The holography stage is sufficient to simulate the acoustic field in the homogenous intracranial medium after ultrasound transmitting through the skull. The agreement of the k-space corrected pseudospectral time domain method and HASC method was assessed by statistical methods: linear regression between the two methods provided a slope of 0.9735, intercept of 0.0078, and R ² of 0.9982. The Bland-Altman method provided a bias of 0.0015 and 95% limits of agreement 0.065 apart. We demonstrated that the difference in sound intensity at the focal point corrected by HASC and time reversal phase correction method was 0.2% and 0.5% in the results of simulation and experiment, respectively. Not only that, the phase calculation time by the HASC phase correction method can be reduced to 11 min on a multi GPU array, which has clinical potential for ultrasound treatment of brain therapy.

Impact of Aperture, Depth, and Acoustic Clutter on the Performance of Coherent Multi-Transducer Ultrasound Imaging

Transducers with a larger aperture size are desirable in ultrasound imaging to improve resolution and image quality. A coherent multi-transducer ultrasound imaging system (CoMTUS) enables an extended effective aperture through the coherent combination of multiple transducers. In this study, the discontinuous extended aperture created by CoMTUS and its performance for deep imaging and through layered media are investigated by both simulations and experiments. Typical image quality metrics—resolution, contrast and contrast-to-noise ratio—are evaluated and compared with a standard single probe imaging system. Results suggest that the image performance of CoMTUS depends on the relative spatial location of the arrays. The resulting effective aperture significantly improves resolution, while the separation between the arrays may degrade contrast. For a limited gap in the effective aperture (less than a few centimetres), CoMTUS provides benefits to image quality compared to the standard single probe imaging system. Overall, CoMTUS shows higher sensitivity and reduced loss of resolution with imaging depth. In general, CoMTUS imaging performance was unaffected when imaging through a layered medium with different speed of sound values and resolution improved up to 80% at large imaging depths.

Low-intensity pulsed ultrasound ameliorates depression-like behaviors in a rat model of chronic unpredictable stress

Introduction

There is an unmet need for better nonpharmaceutical treatments for depression. Low‐intensity pulsed ultrasound (LIPUS) is a novel type of neuromodulation that could be helpful for depressed patients.

Objective

The goal of this study was to investigate the feasibility and potential mechanisms of LIPUS in the treatment of depression.

Methods

Chronic unpredictable stress (CUS) was used to generate rats with depression‐like features that were treated with four weeks of LIPUS stimulation of the ventromedial prefrontal cortex. Depression‐like behaviors were assessed with the sucrose preference, forced swim, and open field tests. BDNF/mTORC1 signaling was examined by Western blot to investigate this potential molecular mechanism. The safety of LIPUS was evaluated using hematoxylin‐eosin and Nissl staining.

Results

Four weeks of LIPUS stimulation significantly increased sucrose preference and reduced forced swim immobility time in CUS rats. LIPUS also partially reversed the molecular effects of CUS that included decreased levels of BDNF, phosphorylated tyrosine receptor kinase B (TrkB), extracellular signal‐regulated kinase (ERK), mammalian target of rapamycin complex 1 (mTORC1), and S6 kinase (S6K). Moreover, histological staining revealed no gross tissue damage.

Conclusions

Chronic LIPUS stimulation can effectively and safely improve depression‐like behaviors in CUS rats. The underlying mechanisms may be related to enhancement of BDNF/ERK/mTORC1 signaling pathways in the prefrontal cortex (PFC). LIPUS is a promising noninvasive neuromodulation tool that merits further study as a potential treatment for depression.