Sensitivity of simulated transcranial ultrasound fields to acoustic medium property maps

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.

Strategies and safety simulations for ultrasonic cervical spinal cord neuromodulation

Objective. Focused ultrasound spinal cord neuromodulation has been demonstrated in small animals. However, most of the tested neuromodulatory exposures are similar in intensity and exposure duration to the reported small animal threshold for possible spinal cord damage. All efforts must be made to minimize the risk and assure the safety of potential human studies, while maximizing potential treatment efficacy. This requires an understanding of ultrasound propagation and heat deposition within the human spine.Approach. Combined acoustic and thermal modelling was used to assess the pressure and heat distributions produced by a 500 kHz source focused to the C5/C6 level via two approaches (a) the posterior acoustic window between vertebral posterior arches, and (b) the lateral intervertebral foramen from which the C6 spinal nerve exits. Pulse trains of fifty 0.1 s pulses (pulse repetition frequency: 0.33 Hz, free-field spatial peak pulse-averaged intensity: 10 W cm-2) were simulated for four subjects and for ±10 mm translational and ±10∘rotational source positioning errors.Main results.Target pressures ranged between 20%-70% of free-field spatial peak pressures with the posterior approach, and 20%-100% with the lateral approach. When the posterior source was optimally positioned, peak spine heating values were below 1 ∘C, but source mispositioning resulted in bone heating up to 4 ∘C. Heating with the lateral approach did not exceed 2 ∘C within the mispositioning range. There were substantial inter-subject differences in target pressures and peak heating values. Target pressure varied three to four-fold between subjects, depending on approach, while peak heating varied approximately two-fold between subjects. This results in a nearly ten-fold range between subjects in the target pressure achieved per degree of maximum heating.Significance. This study highlights the utility of trans-spine ultrasound simulation software and need for precise source-anatomy positioning to assure the subject-specific safety and efficacy of focused ultrasound spinal cord therapies.

Improved Transcranial Plane-Wave Imaging With Learned Speed-of-Sound Maps

Although transcranial ultrasound plane-wave imaging (PWI) has promising clinical application prospects, studies have shown that variable speed-of-sound (SoS) would seriously damage the quality of ultrasound images. The mismatch between the conventional constant velocity assumption and the actual SoS distribution leads to the general blurring of ultrasound images. The optimization scheme for reconstructing transcranial ultrasound image is often solved using iterative methods like full-waveform inversion. These iterative methods are computationally expensive and based on prior magnetic resonance imaging (MRI) or computed tomography (CT) information. In contrast, the multi-stencils fast marching (MSFM) method can produce accurate time travel maps for the skull with heterogeneous acoustic speed. In this study, we first propose a convolutional neural network (CNN) to predict SoS maps of the skull from PWI channel data. Then, use these maps to correct the travel time to reduce transcranial aberration. To validate the performance of the proposed method, numerical, phantom and intact human skull studies were conducted using a linear array transducer (L11-5v, 128 elements, pitch = 0.3 mm). Numerical simulations demonstrate that for point targets, the lateral resolution of MSFM-restored images increased by 65%, and the center position shift decreased by 89%. For the cyst targets, the eccentricity of the fitting ellipse decreased by 75%, and the center position shift decreased by 58%. In the phantom study, the lateral resolution of MSFM-restored images was increased by 49%, and the position shift was reduced by 1.72 mm. This pipeline, termed AutoSoS, thus shows the potential to correct distortions in real-time transcranial ultrasound imaging, as demonstrated by experiments on the intact human skull.

Numerical Investigation of Layered Homogeneous Skull Model for Simulations of Transcranial Focused Ultrasound

BACKGROUND AND OBJECTIVES

The influence of the intracranial pressure field must be discussed with the development of a single-element transducer for low-intensity transcranial focused ultrasound because the skull plays a significant role in blocking and dispersing ultrasound wave propagation. Ultrasound propagation is mainly affected by the structure and acoustic properties of the skull; thus, we aimed to investigate the impact of simplifying the acoustic properties of the skull on the simulation of the transcranial pressure field to present guidance for efficient skull modeling in full-wave simulations.

MATERIALS AND METHODS

We constructed a three-dimensional computational model for ultrasound transmission with the same structure but varying acoustic properties of the skull. The structural information and heterogeneous acoustic properties of the skull were acquired from computed tomography images, and we segmented the skull into three layers (3 L), including spongy and compact bones. We then assigned homogeneous acoustic properties to a single layer (1 L) or 3 L of the skull. In addition, we investigated the influence of different types of transducers and different ultrasound frequencies (1.1 MHz, 0.5 MHz, and 0.25 MHz) on the intracranial pressure field to provide a comparison of the heterogenous and homogeneous models.

RESULTS

We indicated the importance of numerical simulations in estimating the intracranial pressure field of the skull owing to beam distortions. When we simplified the skull model, both the 1 L and 3 L models showed contours of the acoustic focus comparable to those of the heterogeneous model. When we evaluated the peak pressure and volume of the acoustic focus, the 1 L model produced a better estimation of peak pressure with a difference <10%, and the 3 L model is suitable to obtain smaller errors in the volume of the acoustic focus.

CONCLUSIONS

In conclusion, we examined the possibility of simplification of skull models using 1 L and 3 L homogeneous properties in the numerical simulation for focused ultrasound. The results show that the layered homogeneous model can provide characteristics comparable to those of the acoustic focus in heterogeneous models.

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.

A rapid element pressure field simulation method for transcranial phase correction in focused ultrasound therapy

Transcranial focused ultrasound ablation has emerged as a promising technique for treating neurological disorders. The clinical system exclusively employed the ray tracing method to compute phase aberrations induced by the human skull, taking into account computational time constraints. However, this method compromises slightly on accuracy compared to simulation-based methods. This study evaluates a fast simulation method that simulates the time-harmonic pressure field within the region of interest for effective phase correction. Experimental validation was carried out using a 512-element, 670 kHz hemispherical transducer for fourex vivoskulls. The ray tracing method achieved a restoration ratio of 64.81% ± 4.33% of acoustic intensity normalized to hydrophone measurements. In comparison, the rapid simulation method demonstrated improved results with a restoration ratio of 73.10% ± 7.46%, albeit slightly lower than the full-wave simulation which achieved a restoration ratio of 75.87% ± 5.40%. The rapid simulation methods exhibited computational times that were less than five minutes for parallel computation with 8 threads. The incident angle was calculated, and a maximum difference of 6.8 degrees was found when the fixed position of the skull was changed. Meanwhile, the restoration ratio of acoustic intensity was validated to be above 70% for different target positions away from the geometrical focus of the transducer. The favorable balance between time consumption and correction accuracy makes this method valuable for clinical treatment applications.

Dual-Probe Transcranial Full-Waveform Inversion: A Brain Phantom Feasibility Study

OBJECTIVE

Despite being a low-cost, portable and safe medical imaging technique, transcranial ultrasound imaging is not used widely in adults because of the severe degradation and distortion of signals caused by the skull. Full-waveform inversion (FWI) has recently been found to have potential as an effective method for transcranial ultrasound tomography to obtain high-quality, subwavelength-resolution acoustic models of the brain using low-frequency ultrasound data. In this study is the first demonstration of this method in recovering a high-resolution 2-D reconstruction of a brain and skull ultrasound imaging phantom using experimentally acquired data.

METHODS

A 2:5 scale brain phantom encased within a 3-D-printed skull-mimicking layer was created to simulate a clinical transcranial imaging target. To obtain tomographic ultrasound data on the brain and skull phantom, a tomographic ultrasound acquisition system was designed and implemented using commercially available low-frequency cardiac probes. FWI reconstructions of the brain and skull phantom were performed using the acquired tomographic data and were compared with corresponding synthetic reconstructions. This comparison was used to evaluate the feasibility of the proposed imaging system when employing different transducer array configurations.

RESULTS

We demonstrate the successful FWI reconstruction of the brain phantom within the skull mimic from experimentally acquired tomographic ultrasound data. To mitigate the effects of the skull-mimicking material, a reflection-matching algorithm was applied to model the morphology of the skull layer prior to performing the inversion.

CONCLUSION

The findings of this study provide a promising step toward the clinical use of FWI for transcranial ultrasound imaging in adults.

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.

Computational modeling and minimization of unintended neuronal excitation in a LIFU stimulation

The neuromodulation effect of low-intensity focused ultrasound (LIFU) is highly target-specific. Unintended off-target neuronal excitation can be elicited when the beam focusing accuracy and resolution are limited, whereas the resulted side effect has not been evaluated quantitatively. There is also a lack of methods addressing the minimization of such side effects. Therefore, this work introduces a computational model of unintended neuronal excitation during LIFU neuromodulation, which evaluates the off-target activation area (OTAA) by integrating an ultrasound field model with the neuronal spiking model. In addition, a phased array beam focusing scheme called constrained optimal resolution beamforming (CORB) is proposed to minimize the off-target neuronal excitation area while ensuring effective stimulation in the target brain region. A lower bound of the OTAA is analytically approximated in a simplified homogeneous medium, which could guide the selection of transducer parameters such as aperture size and operating frequency. Simulations in a human head model using three transducer setups show that CORB markedly reduces the OTAA compared with two benchmark beam focusing methods. The high neuromodulation resolution demonstrates the capability of LIFU to effectively limit the side effects during neuromodulation, allowing future clinical applications such as treatment of neuropsychiatric disorders.

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 simulation study on the sensitivity of transcranial ray-tracing ultrasound modeling to skull properties

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

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.

Effects of Blood Pressure on Brain Tissue Pulsation Amplitude in a Phantom Model

OBJECTIVE

The precise mechanism and determinants of brain tissue pulsations (BTPs) are poorly understood, and the impact of blood pressure (BP) on BTPs is relatively unexplored. This study aimed to explore the relationship between BP parameters (mean arterial pressure [MAP] and pulse pressure [PP]) and BTP amplitude, using a transcranial tissue Doppler prototype.

METHODS

A phantom brain model generating arterial-induced BTPs was developed to observe BP changes in the absence of confounding variables and cerebral autoregulation feedback processes. A regression model was developed to investigate the relationship between bulk BTP amplitude and BP. The separate effects of PP and MAP were evaluated and quantified.

RESULTS

The regression model (R² = 0.978) revealed that bulk BTP amplitude measured from 27 gates significantly increased with PP but not with MAP. Every 1 mm Hg increase in PP resulted in a bulk BTP amplitude increase of 0.29 µm.

CONCLUSION

Increments in BP were significantly associated with increments in bulk BTP amplitude. Further work should aim to confirm the relationship between BP and BTPs in the presence of cerebral autoregulation and explore further physiological factors having an impact on BTP measurements, such as cerebral blood flow volume, tissue distensibility and intracranial pressure.

Soft Tissue Aberration Correction for Histotripsy Using Acoustic Emissions From Cavitation Cloud Nucleation and Collapse

OBJECTIVE

Phase aberration from soft tissue limits the efficacy of histotripsy, a therapeutic ultrasound technique based on acoustic cavitation. Previous work has shown that the acoustic emissions from cavitation can serve as "point sources" for aberration correction (AC). This study compared the efficacy of soft tissue AC for histotripsy using acoustic cavitation emissions (ACE) from bubble cloud nucleation and collapse.

METHODS

A 750-kHz, receive-capable histotripsy array was pulsed to generate cavitation in ex vivo porcine liver through an intervening abdominal wall. Received ACE signals were used to determine the arrival time differences to the focus and compute corrective delays. Corrections from single pulses and from the median of multiple pulses were tested.

DISCUSSION

On average, ACE AC obtained 96% ± 3% of the pressure amplitude obtained by hydrophone-based correction (compared with 71% ± 5% without AC). Both nucleation- and collapse-based corrections obtained >96% of the hydrophone-corrected pressure when using medians of ≥10 pulses. When using single-pulse corrections, nucleation obtained a range of 49%-99% of the hydrophone-corrected pressure, while collapse obtained 95%-99%.

CONCLUSION

The results suggest that (i) ACE AC can recover nearly all pressure amplitude lost owing to soft tissue aberration and that (ii) the collapse signal permits robust AC using a small number of pulses.

Physics-informed neural networks for transcranial ultrasound wave propagation

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

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.

Simultaneous multi-target ultrasound neuromodulation in freely-moving mice based on a single-element ultrasound transducer

Objective.Ultrasound neuromodulation has become an emerging method for the therapy of neurodegenerative and psychiatric diseases. The phased array ultrasonic transducer enables multi-target ultrasound neuromodulation in small animals, but the relatively large size and mass and the thick cables of the array limit the free movement of small animals. Furthermore, spatial interference may occur during multi-target ultrasound brain stimulation with multiple micro transducers.Approach.In this study, we developed a miniature power ultrasound transducer and used the virtual source time inversion method and 3D printing technology to design, optimize, and manufacture the acoustic holographic lens to construct a multi-target ultrasound neuromodulation system for free-moving mice. The feasibility of the system was verified byin vitrotranscranial ultrasound field measurements,in vivodual-target blood-brain barrier (BBB) opening experiments, andin vivodual-target ultrasound neuromodulation experiments.Main results.The developed miniature transducer had a diameter of 4.0 mm, a center frequency of 1.1 MHz, and a weight of 1.25 g. The developed miniature acoustic holographic lens had a weight of 0.019 g to generate dual-focus transcranial ultrasound. The ultrasonic field measurements' results showed that the bifocal's horizontal distance was 3.0 mm, the -6 dB focal spot width in thex-direction was 2.5 and 2.25 mm, and 2.12 and 2.24 mm in they-direction. Finally, thein vivoexperimental results showed that the system could achieve dual-target BBB opening and ultrasound neuromodulation in freely-moving mice.Significance.The ultrasonic neuromodulation system based on a miniature single-element transducer and the miniature acoustic holographic lens could achieve dual-target neuromodulation in awake small animals, which is expected to be applied to the research of non-invasive dual-target ultrasonic treatment of brain diseases in awake small animals.

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.

Three-layer model with absorption for conservative estimation of the maximum acoustic transmission coefficient through the human skull for transcranial ultrasound stimulation

Transcranial ultrasound stimulation (TUS) has been shown to be a safe and effective technique for non-invasive superficial and deep brain stimulation. Safe and efficient translation to humans requires estimating the acoustic attenuation of the human skull. Nevertheless, there are no international guidelines for estimating the impact of the skull bone. A tissue independent, arbitrary derating was developed by the U.S. Food and Drug Administration to take into account tissue absorption (0.3 dB/cm-MHz) for diagnostic ultrasound. However, for the case of transcranial ultrasound imaging, the FDA model does not take into account the insertion loss induced by the skull bone, nor the absorption by brain tissue. Therefore, the estimated absorption is overly conservative which could potentially limit TUS applications if the same guidelines were to be adopted. Here we propose a three-layer model including bone absorption to calculate the maximum pressure transmission through the human skull for frequencies ranging between 100 kHz and 1.5 MHz. The calculated pressure transmission decreases with the frequency and the thickness of the bone, with peaks for each thickness corresponding to a multiple of half the wavelength. The 95th percentile maximum transmission was calculated over the accessible surface of 20 human skulls for 12 typical diameters of the ultrasound beam on the skull surface, and varies between 40% and 78%. To facilitate the safe adjustment of the acoustic pressure for short ultrasound pulses, such as transcranial imaging or transcranial ultrasound stimulation, a table summarizes the maximum pressure transmission for each ultrasound beam diameter and each frequency.

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.

Modelling transcranial ultrasound neuromodulation: an energy-based multiscale framework

Several animal and human studies have now established the potential of low intensity, low frequency transcranial ultrasound (TUS) for non-invasive neuromodulation. Paradoxically, the underlying mechanisms through which TUS neuromodulation operates are still unclear, and a consensus on the identification of optimal sonication parameters still remains elusive. One emerging hypothesis based on thermodynamical considerations attributes the acoustic-induced nerve activity alterations to the mechanical energy and/or entropy conversions occurring during TUS action. Here, we propose a multiscale modelling framework to examine the energy states of neuromodulation under TUS. First, macroscopic tissue-level acoustic simulations of the sonication of a whole monkey brain are conducted under different sonication protocols. For each one of them, mechanical loading conditions of the received waves in the anterior cingulate cortex region are recorded and exported into a microscopic cell-level 3D viscoelastic finite element model of neuronal axon embedded extracellular medium. Pulse-averaged elastically stored and viscously dissipated energy rate densities during axon deformation are finally computed under different sonication incident angles and are mapped against distinct combinations of sonication parameters of the TUS. The proposed multiscale framework allows for the analysis of vibrational patterns of the axons and its comparison against the spectrograms of stimulating ultrasound. The results are in agreement with literature data on neuromodulation, demonstrating the potential of this framework to identify optimised acoustic parameters in TUS neuromodulation. The proposed approach is finally discussed in the context of multiphysics energetic considerations, argued here to be a promising avenue towards a scalable framework for TUS in silico predictions. STATEMENT OF SIGNIFICANCE: Low-intensity transcranial ultrasound (TUS) is poised to become a leading neuromodulation technique for the treatment of neurological disorders. Paradoxically, how it operates at the cellular scale remains unknown, hampering progress in personalised treatment. To this end, models of the multiphysics of neurons able to upscale results to the organ scale are required. We propose here to achieve this by considering an axon submitted to an ultrasound wave extracted from a simulation at the organ scale. Doing so, information pertaining to both stored and dissipated axonal energies can be extracted for a given head/brain morphology. This two-scale multiphysics energetic approach is a promising scalable framework for in silico predictions in the context of personalised TUS treatment.

Effects of phase aberration on transabdominal focusing for a large aperture, lowf-number histotripsy transducer

Objective. Soft tissue phase aberration may be particularly severe for histotripsy due to large aperture and lowf-number transducer geometries. This study investigated how phase aberration from human abdominal tissue affects focusing of a large, strongly curved histotripsy transducer.Approach.A computational model (k-Wave) was experimentally validated withex vivoporcine abdominal tissue and used to simulate focusing a histotripsy transducer (radius: 14.2 cm,f-number: 0.62, central frequencyfc: 750 kHz) through the human abdomen. Abdominal computed tomography images from 10 human subjects were segmented to create three-dimensional acoustic property maps. Simulations were performed focusing at 3 target locations in the liver of each subject with ideal phase correction, without phase correction, and after separately matching the sound speed of water and fat to non-fat soft tissue.Main results.Experimental validation in porcine abdominal tissue showed that simulated and measured arrival time differences agreed well (average error, ∼0.10 acoustic cycles atfc). In simulations with human tissue, aberration created arrival time differences of 0.65μs (∼0.5 cycles) at the target and shifted the focus from the target by 6.8 mm (6.4 mm pre-focally along depth direction), on average. Ideal phase correction increased maximum pressure amplitude by 95%, on average. Matching the sound speed of water and fat to non-fat soft tissue decreased the average pre-focal shift by 3.6 and 0.5 mm and increased pressure amplitude by 2% and 69%, respectively.Significance.Soft tissue phase aberration of large aperture, lowf-number histotripsy transducers is substantial despite low therapeutic frequencies. Phase correction could potentially recover substantial pressure amplitude for transabdominal histotripsy. Additionally, different heterogeneity sources distinctly affect focusing quality. The water path strongly affects the focal shift, while irregular tissue boundaries (e.g. fat) dominate pressure loss.

Numerical Study of Acoustic Holograms for Deep-Brain Targeting through the Temporal Bone Window

Acoustic holograms can encode complex wavefronts to compensate the aberrations of a therapeutical ultrasound beam propagating through heterogeneous tissues such as the skull, and simultaneously, they can generate diffraction-limited acoustic images, that is, arbitrary shaped focal spots. In this work, we numerically study the performance of acoustic holograms focusing at the thalamic nuclei when the source is located at the temporal bone window. The temporal window is the thinnest area of the lateral skull and it is mainly hairless, so it is a desirable area through which to transmit ultrasonic waves to the deep brain. However, in targeting from this area the bilateral thalamic nuclei are not aligned with the elongated focal spots of conventional focused transducers, and in addition, skull aberrations can distort the focal spot. We found that by using patient-specific holographic lenses coupled to a single-element 650-kHz-frequency 65-mm-aperture source, the focal spot can be sharply adapted to the thalamic nuclei in a bilateral way while skull aberrations are mitigated. Furthermore, the performance of these holograms was studied under misalignment errors between the source and the skull, concluding that for misalignments up to 5°, acoustic images are correctly restored. This work paves the way to designing clinical applications of transcranial ultrasound such as blood-brain barrier opening for drug delivery or deep-brain neuromodulation using this low-cost and personalized technology, presenting desirable aspects for long-term treatments because the patient's head does not need to be shaved completely and skull heating is low.

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.

Mathematical Model of Ultrasound Attenuation With Skull Thickness for Transcranial-Focused Ultrasound

Transcranial-focused ultrasound (tFUS) has potential for both neuromodulation and neuroimaging. Due to the influence of head tissue, especially the skull, its attenuation is a key issue affecting precise focusing. The objective of the present study was to construct a mathematical model of ultrasound attenuation inclusive of skull thickness. First, combined with real skull phantom experiments and simulation experiments, tFUS attenuation of different head tissues was investigated. Furthermore, based on the system identification method, a mathematical model of ultrasound attenuation was constructed taking skull thickness into account. Finally, the performance of the mathematical model was tested, and its potential applications were investigated. For different head tissues, including scalp, skull, and brain tissue, the skull was found to be the biggest influencing factor for ultrasound attenuation, the attenuation caused by it being 4.70 times and 7.06 times that of attenuation caused by the brain and scalp, respectively. Consistent with the results of both the simulation and phantom experiments, the attenuation of the mathematical model increased as the skull thickness increased. The average error of the mathematical model was 1.87% in the phantom experiment. In addition, the experimental results show that the devised mathematical model is suitable for different initial pressures and different skulls with correlation coefficients higher than 0.99. Both simulation and phantom experiments validated the effectiveness of the proposed mathematical model. It can be concluded from this experiment that the proposed mathematical model can accurately calculate the tFUS attenuation and can significantly contribute to further research and application of tFUS.

Numerical Evaluation of the Effects of Transducer Displacement on Transcranial Focused Ultrasound in the Rat Brain

Focused ultrasound is a promising therapeutic technique, as it involves the focusing of an ultrasonic beam with sufficient acoustic energy into a target brain region with high precision. Low-intensity ultrasound transmission by a single-element transducer is mostly established for neuromodulation applications and blood-brain barrier disruption for drug delivery. However, transducer positioning errors can occur without fine control over the sonication, which can affect repeatability and lead to reliability problems. The objective of this study was to determine whether the target brain region would be stable under small displacement (0.5 mm) of the transducer based on numerical simulations. Computed-tomography-derived three-dimensional models of a rat head were constructed to investigate the effects of transducer displacement in the caudate putamen (CP) and thalamus (TH). Using three different frequencies (1.1, 0.69, and 0.25 MHz), the transducer was displaced by 0.5 mm in each of the following six directions: superior, interior, anterior, posterior, left, and right. The maximum value of the intracranial pressure field was calculated, and the targeting errors were determined by the full-width-at-half-maximum (FWHM) overlap between the free water space (FWHMwater) and transcranial transmission (FWHMbase). When the transducer was positioned directly above the target region, a clear distinction between the target regions was observed, resulting in 88.3%, 81.5%, and 84.5% FWHMwater for the CP and 65.6%, 76.3%, and 64.4% FWHMwater for the TH at 1.1, 0.69, and 0.25 MHz, respectively. Small transducer displacements induced both enhancement and reduction of the peak pressure and targeting errors, compared with when the transducer was displaced in water. Small transducer displacement to the left resulted in the lowest stability, with 34.8% and 55.0% targeting accuracy (FWHMwater) at 1.1 and 0.69 MHz in the TH, respectively. In addition, the maximum pressure was reduced by up to 11% by the transducer displacement. This work provides the targeting errors induced by transducer displacements through a preclinical study and recommends that attention be paid to determining the initial sonication foci in the transverse plane in the cases of small animals.

Factors Associated with Energy Efficiency of Focused Ultrasound Through the Skull: A Study of 3D-Printed Skull Phantoms and Its Comparison with Clinical Experiences

While focused ultrasound (FUS) is non-invasive, the ultrasound energy is attenuated by the skull which results in differences in energy efficiency among patients. In this study, we investigated the effect of skull variables on the energy efficiency of FUS. The thickness and density of the skull and proportion of the trabecular bone were selected as factors that could affect ultrasound energy transmittance. Sixteen 3D-printed skull models were designed and fabricated to reflect the three factors. The energy of each phantom was measured using an ultrasonic sound field energy measurement system. The thickness and proportion of trabecular bone affected the attenuation of transmitted energy. There was no difference in the density of the trabecular bone. In clinical data, the trabecular bone ratio showed a significantly greater correlation with dose/delivered energy than that of thickness and the skull density ratio. Currently, for clinical non-thermal FUS, the data are not sufficient, but we believe that the results of this study will be helpful in selecting patients and appropriate parameters for FUS treatment.

Measurement of the ultrasound attenuation and dispersion in 3D-printed photopolymer materials from 1 to 3.5 MHz

Over the past decade, the range of applications in biomedical ultrasound exploiting 3D printing has rapidly expanded. For wavefront shaping specifically, 3D printing has enabled a diverse range of new, low-cost approaches for controlling acoustic fields. These methods rely on accurate knowledge of the bulk acoustic properties of the materials; however, to date, robust knowledge of these parameters is lacking for many materials that are commonly used. In this work, the acoustic properties of eight 3D-printed photopolymer materials were characterised over a frequency range from 1 to 3.5 MHz. The properties measured were the frequency-dependent phase velocity and attenuation, group velocity, signal velocity, and mass density. The materials were fabricated using two separate techniques [PolyJet and stereolithograph (SLA)], and included Agilus30, FLXA9960, FLXA9995, Formlabs Clear, RGDA8625, RGDA8630, VeroClear, and VeroWhite. The range of measured density values across all eight materials was 1120-1180 kg  · m^-3, while the sound speed values were between 2020 to 2630 m  · s^-1, and attenuation values typically in the range 3-9 dB  ·  MHz^-1· cm^-1.

The Effects of the Structural and Acoustic Parameters of the Skull Model on Transcranial Focused Ultrasound

Transcranial focused ultrasound (tFUS) has great potential in brain imaging and therapy. However, the structural and acoustic differences of the skull will cause a large number of technical problems in the application of tFUS, such as low focus energy, focal shift, and defocusing. To have a comprehensive understanding of the skull effect on tFUS, this study investigated the effects of the structural parameters (thickness, radius of curvature, and distance from the transducer) and acoustic parameters (density, acoustic speed, and absorption coefficient) of the skull model on tFUS based on acrylic plates and two simulation methods (self-programming and COMSOL). For structural parameters, our research shows that as the three factors increase the unit distance, the attenuation caused from large to small is the thickness (0.357 dB/mm), the distance to transducer (0.048 dB/mm), and the radius of curvature (0.027 dB/mm). For acoustic parameters, the attenuation caused by density (0.024 dB/30 kg/m³) and acoustic speed (0.021 dB/30 m/s) are basically the same. Additionally, as the absorption coefficient increases, the focus acoustic pressure decays exponentially. The thickness of the structural parameters and the absorption coefficient of the acoustic parameters are the most important factors leading to the attenuation of tFUS. The experimental and simulation trends are highly consistent. This work contributes to the comprehensive and quantitative understanding of how the skull influences tFUS, which further enhances the application of tFUS in neuromodulation research and treatment.

Simulation study on the effects of cancellous bone structure in the skull on ultrasonic wave propagation

The transcranial Doppler method (TCD) enables the measurement of cerebral blood flow velocity and detection of emboli by applying an ultrasound probe to the temporal bone window, or the orbital or greater occipital foramina. TCD is widely used for evaluation of cerebral vasospasm after subarachnoid hemorrhage, early detection of patients with arterial stenosis, and the assessment of brain death. However, measurements often become difficult in older women. Among various factors contributing to this problem, we focused on the effect of the diploe in the skull bone on the penetration of ultrasound into the brain. In particular, the effect of the cancellous bone structure in the diploe was investigated. Using a 2D digital bone model, wave propagation through the skull bone was investigated using the finite-difference time-domain (FDTD) method. We fabricated digital bone models with similar structure but different BV/TV (bone volume/total volume) values in the diploe. At a BV/TV of approximately 50–60% (similar to that of older women), the minimum ultrasound amplitude was observed as a result of scattering and multiple reflections in the cancellous diploe. These results suggest that structural changes such as osteoporosis may be one factor hampering TCD measurements.

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.

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.

The impact of CT image parameters and skull heterogeneity modeling on the accuracy of transcranial focused ultrasound simulations

Objective. Low-intensity transcranial ultrasound stimulation (TUS) is a promising non-invasive brain stimulation (NIBS) technique. TUS can reach deeper areas and target smaller regions in the brain than other NIBS techniques, but its application in humans is hampered by the lack of a straightforward and reliable procedure to predict the induced ultrasound exposure. Here, we examined how skull modeling affects computer simulations of TUS.Approach. We characterized the ultrasonic beam after transmission through a sheep skull with a hydrophone and performed computed tomography (CT) image-based simulations of the experimental setup. To study the skull model's impact, we varied: CT acquisition parameters (tube voltage, dose, filter sharpness), image interpolation, segmentation parameters, acoustic property maps (speed-of-sound, density, attenuation), and transducer-position mismatches. We compared the impact of modeling parameter changes on model predictions and on measurement agreement. Spatial-peak intensity and location, total power, and the Gamma metric (a measure for distribution differences) were used as quantitative criteria. Modeling-based sensitivity analysis was also performed for two human head models.Main results. Sheep skull attenuation assignment and transducer positioning had the most important impact on spatial peak intensity (overestimation up to 300%, respectively 30%), followed by filter sharpness and tube voltage (up to 20%), requiring calibration of the mapping functions. Positioning and skull-heterogeneity-structure strongly affected the intensity distribution (gamma tolerances exceeded in>80%, respectively>150%, of the focus-volume in water), necessitating image-based personalized modeling. Simulation results in human models consistently demonstrate a high sensitivity to the skull-heterogeneity model, attenuation tuning, and transducer shifts, the magnitude of which depends on the underlying skull structure complexity.Significance. Our study reveals the importance of properly modeling the skull-heterogeneity and its structure and of accurately reproducing the transducer position. The results raise red flags when translating modeling approaches among clinical sites without proper standardization and/or recalibration of the imaging and modeling parameters.

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.

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.

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

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

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.

Computationally Efficient Transcranial Ultrasonic Focusing: Taking Advantage of the High Correlation Length of the Human Skull

The phase correction necessary for transcranial ultrasound therapy requires numerical simulation to noninvasively assess the phase shift induced by the skull bone. Ideally, the numerical simulations need to be fast enough for clinical implementation in a brain therapy protocol and to provide accurate estimation of the phase shift to optimize the refocusing through the skull. In this article, we experimentally performed transcranial ultrasound focusing at 900 kHz on N = 5 human skulls. To reduce the computation time, we propose here to perform the numerical simulation at 450 kHz and use the corresponding phase shifts experimentally at 900 kHz. We demonstrate that a 450-kHz simulation restores 94.2% of the pressure when compared with a simulation performed at 900 kHz and 85.0% of the gold standard pressure obtained by an invasive time reversal procedure based on the signal recorded by a hydrophone placed at the target. From a 900- to 450-kHz simulation, the grid size is divided by 8, and the computation time is divided by 10.

Echo-Focusing in Transcranial Focused Ultrasound Thalamotomy for Essential Tremor: A Feasibility Study

BACKGROUND

Transcranial magnetic resonance-guided focused ultrasound (TcMRgFUS) systems currently employ computed tomography (CT)-based aberration corrections, which may provide suboptimal trans-skull focusing.

OBJECTIVES

The objective of this study was to evaluate a contrast agent microbubble imaging-based transcranial focusing method, echo-focusing (EF), during TcMRgFUS for essential tremor.

METHODS

A clinical trial of TcMRgFUS thalamotomy using EF for the treatment of essential tremor was conducted (NCT03935581; funded by InSightec [Tirat Carmel, Israel]). Patients (n = 12) were injected with Definity (Lantheus Medical Imaging, North Billerica, MA) microbubbles, and EF was performed using a research feature add-on to a commercial TcMRgFUS system (ExAblate Neuro, InSightec). Subablative thermal sonications carried out using (1) EF and (2) CT-based aberration corrections were compared via magnetic resonance thermometry, and the optimal focusing method for each patient was employed for TcMRgFUS thalamotomy.

RESULTS

EF aberration corrections provided increased sonication efficiency, decreased focal size, and equivalent targeting accuracy relative to CT-based focusing. EF aberration corrections were employed successfully for lesion formation in all 12 patients, 3 of whom had previously undergone unsuccessful TcMRgFUS thalamotomy via CT-based focusing. There were no adverse events related directly to the EF procedure.

CONCLUSIONS

EF is feasible and appears safe during TcMRgFUS thalamotomy for essential tremor and improves on the trans-skull focal quality provided by existing CT-based focusing methods. © 2020 International Parkinson and Movement Disorder Society.

Accurate anatomical head segmentations: a data set for biomedical simulations

Detailed computational anatomical models of the entire head are needed for accurate in silico modeling in a variety of transcranial stimulation applications. Models from different subjects help to understand and account for population variability. To this end, we have developed a new library of head models of 20 individuals, segmented from co-aligned multi-modal medical image data. The acquired image modalities allow to accurately model tissues with different material properties, such as electrical conductivity or spatially varying acoustic properties. The usefulness of the models is illustrated for two example applications.

Experimental Validation of k-Wave: Nonlinear Wave Propagation in Layered, Absorbing Fluid Media

Models of ultrasound propagation in biologically relevant media have applications in planning and verification of ultrasound therapies and computational dosimetry. To be effective, the models must be able to accurately predict both the spatial distribution and amplitude of the acoustic pressure. This requires that the models are validated in absolute terms, which for arbitrarily heterogeneous media should be performed by comparison with measurements of the acoustic field. In this article, simulations performed using the open-source k-Wave acoustics toolbox, with a measurement-based source definition, were quantitatively validated against measurements of acoustic pressure in water and layered absorbing fluid media. In water, the measured and simulated spatial-peak pressures agreed to within 3% under linear conditions and 7% under nonlinear conditions. After propagation through a planar or wedge-shaped glycerol-filled phantom, the difference in spatial-peak pressure was 8.5% and 10.7%, respectively. These differences are within or close to the expected uncertainty of the acoustic pressure measurement. The -6 dB width and length of the focus agreed to within 4% in all cases, and the focal positions were within 0.7 mm for the planar phantom and 1.2 mm for the wedge-shaped phantom. These results demonstrate that when the acoustic medium properties and geometry are well known, accurate quantitative predictions of the acoustic field can be made using k-Wave.

Elementwise approach for simulating transcranial MRI-guided focused ultrasound thermal ablation

This work explored an elementwise approach to model transcranial MRI-guided focused ultrasound (TcMRgFUS) thermal ablation, a noninvasive approach to neurosurgery. Each element of the phased array transducer was simulated individually and could be simultaneously loaded into computer memory, allowing for rapid (~2.5 s) calculation of the pressure field for different phase offsets used for beam steering and aberration correction. We simulated the pressure distribution for 431 sonications in 32 patients, applied the phase and magnitude values used during treatment, and estimated the resulting temperature rise. We systematically varied the relationship between CT (computerized tomography)-derived skull density and the acoustic attenuation and sound speed to obtain the best agreement between the predictions and MR temperature imaging (MRTI). The optimization was validated with simulations of 396 sonications from 40 additional treatments. After optimization, the predicted and measured heating agreed well (R 2: 0.74 patients 1-32; 0.71 patients 33-72). The dimensions and obliquity of the heating in the simulated temperature maps were correlated with the MRTI (R 2: 0.62, 0.74, respectively), but the measured heating was more spatially diffuse. The energy needed to achieve ablation varied by an order of magnitude (3.3-36.1 kJ). While this elementwise approach required more computation time up front (the combined simulation matrices were approximately 4.6 times higher than a single large simulation), it could be performed in parallel on a computing cluster. It allows for rapid calculation of the three-dimensional heating at the focus for different phase and magnitude values on the array. We also show how this approach can be used to optimize the relationship between CT-derived skull density and acoustic properties. While the relationships found here need further validation in a larger patient population, these results demonstrate the promise of this approach to model TcMRgFUS.

Ultrasound Axicon: Systematic Approach to Optimize Focusing Resolution through Human Skull Bone

The use of axicon lenses is useful in many high-resolution-focused ultrasound applications, such as mapping, detection, and have recently been extended to ultrasonic brain therapies. However, in order to achieve high spatial resolution with an axicon lens, it is necessary to adjust the separation, called stand-off (δ), between a conventional transducer and the lens attached to it. Comprehensive ultrasound simulations, using the open-source k-Wave toolbox, were performed for an axicon lens attached to a piezo-disc type transducer with a radius of 14 mm, and a frequency of about 0.5 MHz, that is within the range of optimal frequencies for transcranial transmission. The materials properties were measured, and the lens geometry was modelled. Hydrophone measurements were performed through a human skull phantom. We obtained an initial easygoing design model for the lens angle and optimal stand-off using relatively simple formulas. The skull is not an obstacle for focusing of ultrasound with optimized axicon lenses that achieve an identical resolution to spherical transducers, but with the advantage that the focusing distance is shortened. An adequate stand-off improves the lateral resolution of the acoustic beam by approximately 50%. The approach proposed provides an effective way of designing polydimethylsiloxane (PDMS)-based axicon lenses equipped transducers.

Method to optimize the placement of a single-element transducer for transcranial focused ultrasound

BACKGROUND AND OBJECTIVE

Transcranial focused ultrasound (tFUS) is a promising neuromodulation technique because of its non-invasiveness and high spatial resolution (within millimeter scale). However, the presence of the skull can lead to disrupting and shifting the acoustic focus in the brain. In this study, we propose a computationally efficient way to determine the optimal position of a single-element focused ultrasound transducer which can effectively deliver acoustic energy to the brain target. We hypothesized that the placement of a single element transducer with the lowest average reflection coefficient would be the optimal position.

METHODS

The reflection coefficient is defined by the ratio of the amplitude of the reflected wave to the incident wave. To calculate the reflection coefficient, we assumed ultrasound waves as straight lines (beam lines). At each beam line, the reflection coefficient was calculated from the incidence angle at the skull interface (outer/inner skull surfaces). The average reflection coefficient (ARC) was calculated at each possible placement of the transducer using a custom-built software. For comparison purposes, acoustic simulations (k-Wave MATLAB toolbox) which numerically solved the linear wave equation were performed with the same transducer positions used in the ARC calculation. In addition, the experimental validation of our proposed method was also performed by measuring acoustic wave propagation through the calvaria skull phantom in water. The accuracy of our method was defined as the distance between the two optimal transducer placements which were determined from the acoustic simulations and from the ARC method.

RESULT

Simulated acoustic pressure distribution corresponding to each ARC showed an inverse relationship with peak acoustic pressures produced in the brain. In comparison to the acoustic simulations, the accuracy of our method was 5.07 ± 4.27 mm when targeting the cortical region in the brain. The computing time of ARC calculations were 0.08% of the time required for acoustic pressure simulations.

CONCLUSION

We calculated the ARC to find the optimal position of the tFUS transducer used in the present study. The optimal placement of the transducer was found when the ARC was the lowest. Our numerical and experimental results showed that the proposed ARC method can effectively be used to find the optimal position of a single-element tFUS transducer for targeting the cortex region of the brain in a computationally inexpensive way.

Safety of transcranial focused ultrasound stimulation: A systematic review of the state of knowledge from both human and animal studies

BACKGROUND

Low-intensity transcranial focused ultrasound stimulation (TFUS) holds great promise as a highly focal technique for transcranial stimulation even for deep brain areas. Yet, knowledge about the safety of this novel technique is still limited.

OBJECTIVE

To systematically review safety related aspects of TFUS. The review covers the mechanisms-of-action by which TFUS may cause adverse effects and the available data on the possible occurrence of such effects in animal and human studies.

METHODS

Initial screening used key term searches in PubMed and bioRxiv, and a review of the literature lists of relevant papers. We included only studies where safety assessment was performed, and this results in 33 studies, both in humans and animals.

RESULTS

Adverse effects of TFUS were very rare. At high stimulation intensity and/or rate, TFUS may cause haemorrhage, cell death or damage, and unintentional blood-brain barrier (BBB) opening. TFUS may also unintentionally affect long-term neural activity and behaviour. A variety of methods was used mainly in rodents to evaluate these adverse effects, including tissue staining, magnetic resonance imaging, temperature measurements and monitoring of neural activity and behaviour. In 30 studies, adverse effects were absent, even though at least one Food and Drug Administration (FDA) safety index was frequently exceeded. Two studies reported microhaemorrhages after long or relatively intense stimulation above safety limits. Another study reported BBB opening and neuronal damage in a control condition, which intentionally and substantially exceeded the safety limits.

CONCLUSION

Most studies point towards a favourable safety profile of TFUS. Further investigations are warranted to establish a solid safety framework for the therapeutic window of TFUS to reliably avoid adverse effects while ensuring neural effectiveness. The comparability across studies should be improved by a more standardized reporting of TFUS parameters.