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

A comparative study of experimental and simulated ultrasound beam propagation through cranial bones

Objective.Transcranial ultrasound is used in a variety of treatments, including neuromodulation, opening the blood-brain barrier, and high intensity focused ultrasound therapies. To ensure safety and efficacy of these treatments, numerical simulations of the ultrasound field within the brain are used for treatment planning and evaluation. This study investigates the accuracy of numerical modelling of the propagation of focused ultrasound through cranial bones.Approach.Holograms of acoustic fields after propagation through four human skull specimens were measured for frequencies ranging from 270 kHz to 1 MHz, using both quasi-continuous and pulsed modes. The open-source k-Wave toolbox was employed for simulations, using an equivalent-source hologram and a uniform bowl source with parameters that best matched the measured free-field pressure distribution.Main results.The average absolute error in k-Wave simulations with sound speed and density derived from CT scans compared to measurements was 15% for the spatial-peak acoustic pressure amplitude, 2.7 mm for the position of the focus, and 35% for the focal volume. Optimised uniform bowl sources achieved calculation accuracy comparable to that of the hologram sources.Significance.This method is demonstrated as a suitable tool for prediction of focal position, size and overall distribution of transcranial ultrasound fields. The accuracy of the shape and position of the focal region demonstrate the suitability of the sound speed and density mapping used here. However, large errors in pressure amplitude and transmission loss in some individual cases show that alternative methods for mapping individual skull attenuation are needed and the possibility of considerable errors in pressure amplitude should be taken into account when planning focused ultrasound studies or interventions in the human brain, and appropriate safety margins should be used.

Sonogenetics is a novel antiarrhythmic mechanism

Arrhythmia of the heart is a dangerous and potentially fatal condition. The current widely used treatment is the implantable cardioverter defibrillator (ICD), but it is invasive and affects the patient's quality of life. The sonogenetic mechanism proposed here focuses ultrasound on a cardiac tissue, controls endogenous stretch-activated Piezo1 ion channels on the focal region's cardiomyocyte sarcolemma, and restores normal heart rhythm. In contrast to anchoring the implanted ICD lead at a fixed position in the myocardium, the size and position of the ultrasound focal region can be selected dynamically by adjusting the signals of every piezoelectric chip on the ultrasonic phased array, and it allows novel and efficient defibrillations. Based on the developed interdisciplinary electro-mechanical model of sonogenetic treatment, our analysis shows that the proposed ultrasound intensity and frequency will be safe and painless for humans and well below the limits established by the U.S. Food and Drug Administration.

Treatment envelope of transcranial histotripsy: challenges and strategies to maximize the treatment location profile

A 750 kHz, 360-element ultrasound array has been built for transcranial histotripsy applications. This study aims to evaluate its performance to determine whether this array is adequate for treating a wide range of brain locations through a human skull. Treatment location profiles in 2 excised human skulls were experimentally characterized based on passive cavitation mapping. Full-wave acoustic simulations were performed in 8 human skulls to analyze the ultrasound propagation at shallow targets in skulls with different properties. Results showed that histotripsy successfully generated cavitation from deep to shallow targets within 5 mm from the skull surface in the skull with high SDR and small thickness, whereas in the skull with low SDR and large thickness, the treatment envelope was limited up to 16 mm from the skull surface. Simulation results demonstrated that the treatment envelope was highly dependent on the skull acoustic properties. Pre-focal pressure hotspots were observed in both simulation and experiments when targeting near the skull. For each skull, the acoustic pressure loss increases significantly for shallow targets compared to central targets due to high attenuation, large incident angles, and pre-focal pressure hotspots. Strategies including array design optimization, pose optimization, and amplitude correction, are proposed to broaden the treatment envelope. This study identifies the capabilities and limitations of the 360-element transcranial histotripsy array and suggests strategies for designing the next-generation transcranial histotripsy array to expand the treatment location profile for a future clinical trial.

Machine Learning-Empowered Real-Time Acoustic Trapping: An Enabling Technique for Increasing MRI-Guided Microbubble Accumulation

Acoustic trap, using ultrasound interference to ensnare bioparticles, has emerged as a versatile tool for life sciences due to its non-invasive nature. Bolstered by magnetic resonance imaging's advances in sensing acoustic interference and tracking drug carriers (e.g., microbubble), acoustic trap holds promise for increasing MRI-guided microbubbles (MBs) accumulation in target microvessels, improving drug carrier concentration. However, accurate trap generation remains challenging due to complex ultrasound propagation in tissues. Moreover, the MBs' short lifetime demands high computation efficiency for trap position adjustments based on real-time MRI-guided carrier monitoring. To this end, we propose a machine learning-based model to modulate the transducer array. Our model delivers accurate prediction of both time-of-flight (ToF) and pressure amplitude, achieving low average prediction errors for ToF (-0.45 µs to 0.67 µs, with only a few isolated outliers) and amplitude (-0.34% to 1.75%). Compared with the existing methods, our model enables rapid prediction (<10 ms), achieving a four-order of magnitude improvement in computational efficiency. Validation results based on different transducer sizes and penetration depths support the model's adaptability and potential for future ultrasound treatments.

tFUSFormer: Physics-Guided Super-Resolution Transformer for Simulation of Transcranial Focused Ultrasound Propagation in Brain Stimulation

Transcranial focused ultrasound (tFUS) has emerged as a new mode of non-invasive brain stimulation (NIBS), with its exquisite spatial precision and capacity to reach the deep regions of the brain. The placement of the acoustic focus onto the desired part of the brain is critical for successful tFUS procedures; however, acoustic wave propagation is severely affected by the skull, distorting the focal location/shape and the pressure level. High-resolution (HR) numerical simulation allows for monitoring of acoustic pressure within the skull but with a considerable computational burden. To address this challenge, we employed a 4× super-resolution (SR) Swin Transformer method to improve the precision of estimating tFUS acoustic pressure field, targeting operator-defined brain areas. The training datasets were obtained through numerical simulations at both ultra-low (2.0 [Formula: see text]) and high (0.5 [Formula: see text]) resolutions, conducted on in vivo CT images of 12 human skulls. Our multivariable datasets, which incorporate physical properties of the acoustic pressure field, wave velocity, and skull CT images, were utilized to train three-dimensional SR models. We found that our method yielded 87.99 ± 4.28% accuracy in terms of focal volume conformity under foreseen skull data, and accuracy of 82.32 ± 5.83% for unforeseen skulls, respectively. Moreover, a significant improvement of 99.4% in computational efficiency compared to the traditional 0.5 [Formula: see text] HR numerical simulation was shown. The presented technique, when adopted in guiding the placement of the FUS transducer to engage specific brain targets, holds great potential in enhancing the safety and effectiveness of tFUS therapy.

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.

Sustained reduction of essential tremor with low-power non-thermal transcranial focused ultrasound stimulations in humans

BACKGROUND

Transcranial ultrasound stimulation (TUS) is a non-invasive brain stimulation technique; when skull aberrations are compensated for, this technique allows, with millimetric accuracy, circumvention of the invasive surgical procedure associated with deep brain stimulation (DBS) and the limited spatial specificity of transcranial magnetic stimulation.

OBJECTIVE

/hypothesis: We hypothesize that MR-guided low-power TUS can induce a sustained decrease of tremor power in patients suffering from medically refractive essential tremor.

METHODS

The dominant hand only was targeted, and two anatomical sites were sonicated in this exploratory study: the ventral intermediate nucleus of the thalamus (VIM) and the dentato-rubro-thalamic tract (DRT). Patients (N = 9) were equipped with MR-compatible accelerometers attached to their hands to monitor their tremor in real-time during TUS.

RESULTS

VIM neurostimulations followed by a low-duty cycle (5 %) DRT stimulation induced a substantial decrease in the tremor power in four patients, with a minimum of 89.9 % reduction when compared with the baseline power a few minutes after the DRT stimulation. The only patient stimulated in the VIM only and with a low duty cycle (5 %) also experienced a sustained reduction of the tremor (up to 93.4 %). Four patients (N = 4) did not respond. The temperature at target was 37.2 ± 1.4 °C compared to 36.8 ± 1.4 °C for a 3 cm away control point.

CONCLUSIONS

MR-guided low power TUS can induce a substantial and sustained decrease of tremor power. Follow-up studies need to be conducted to reproduce the effect and better to understand the variability of the response amongst patients. MR thermometry during neurostimulations showed no significant thermal rise, supporting a mechanical effect.

New semi-analytical method for fast transcranial ultrasonic field simulation

Objective.To optimize and ensure the safety of ultrasound brain therapy, personalized transcranial ultrasound simulations are very useful. They allow to predict the pressure field, depending on the patient skull and probe position. Most transcranial ultrasound simulations are based on numerical methods which have a long computation time and a high memory usage. The goal of this study is to develop a new semi-analytical field computation method that combines realism and computation speed.Approach.Instead of the classic ray tracing, the ultrasonic paths are computed by time of flight minimization. Then the pressure field is computed using the pencil method. This method requires a smooth and homogeneous skull model. The simulation algorithm, so-called SplineBeam, was numerically validated, by comparison with existing solvers, and experimentally validated by comparison with hydrophone measured pressure fields through anex vivohuman skull.Main results.SplineBeam simulated pressure fields were close to the experimentally measured ones, with a focus position difference of the order of the positioning error and a maximum pressure difference lower than 6.02%. In addition, for those configurations, SplineBeam computation time was lower than another simulation software, k-Wave's, by two orders of magnitude, thanks to its capacity to compute the field only at the focal spot.Significance.These results show the potential of this new method to compute fast and realistic transcranial pressure fields. The combination of this two assets makes it a promising tool for real time transcranial pressure field prediction during ultrasound brain therapy interventions.

Finite difference-embedded UNet for solving transcranial ultrasound frequency-domain wavefield

Transcranial ultrasound imaging assumes a growing significance in the detection and monitoring of intracranial lesions and cerebral blood flow. Accurate solution of partial differential equation (PDE) is one of the prerequisites for obtaining transcranial ultrasound wavefields. Grid-based numerical solvers such as finite difference (FD) and finite element methods have limitations including high computational costs and discretization errors. Purely data-driven methods have relatively high demands on training datasets. The fact that physics-informed neural network can only target the same model limits its application. In addition, compared to time-domain approaches, frequency-domain solutions offer advantages of reducing computational complexity and enabling stable and accurate inversions. Therefore, we introduce a framework called FD-embedded UNet (FEUNet) for solving frequency-domain transcranial ultrasound wavefields. The PDE error is calculated using the optimal 9-point FD operator, and it is integrated with the data-driven error to jointly guide the network iterations. We showcase the effectiveness of this approach through experiments involving idealized skull and brain models. FEUNet demonstrates versatility in handling various input scenarios and excels in enhancing prediction accuracy, especially with limited datasets and noisy information. Finally, we provide an overview of the advantages, limitations, and potential avenues for future research in this study.

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.

Shedding light on ultrasound in action: Optical and optoacoustic monitoring of ultrasound brain interventions

Monitoring brain responses to ultrasonic interventions is becoming an important pillar of a growing number of applications employing acoustic waves to actuate and cure the brain. Optical interrogation of living tissues provides a unique means for retrieving functional and molecular information related to brain activity and disease-specific biomarkers. The hybrid optoacoustic imaging methods have further enabled deep-tissue imaging with optical contrast at high spatial and temporal resolution. The marriage between light and sound thus brings together the highly complementary advantages of both modalities toward high precision interrogation, stimulation, and therapy of the brain with strong impact in the fields of ultrasound neuromodulation, gene and drug delivery, or noninvasive treatments of neurological and neurodegenerative disorders. In this review, we elaborate on current advances in optical and optoacoustic monitoring of ultrasound interventions. We describe the main principles and mechanisms underlying each method before diving into the corresponding biomedical applications. We identify areas of improvement as well as promising approaches with clinical translation potential.

Improving Sonication Efficiency in Transcranial MR-Guided Focused Ultrasound Treatment: A Patient-Data Simulation Study

The potential improvement in sonication efficiency achieved by tilting the focused ultrasound (FUS) transducer of the transcranial MR-guided FUS system is presented. A total of 56 cases of patient treatment data were used. The relative position of the clinical FUS transducer to the patient's head was reconstructed, and region-specific skull density and porosity were calculated based on the patient's CT volume image. The total transmission coefficient of acoustic waves emitted from each channel was calculated. Then, the total energy penetrating the human skull-which represents the sonication efficiency-was estimated. As a result, improved sonication efficiency was by titling the FUS transducer to a more appropriate angle achieved in all 56 treatment cases. This simulation result suggests the potential improvement in transcranial-focused ultrasound treatment by simply adjusting the transducer angle.

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.

Aberration correction in abdominal histotripsy

In transabdominal histotripsy, ultrasound pulses are focused on the body to noninvasively destroy soft tissues via cavitation. However, the ability to focus is limited by phase aberration, or decorrelation of the ultrasound pulses due to spatial variation in the speed of sound throughout heterogeneous tissue. Phase aberration shifts, broadens, and weakens the focus, thereby reducing the safety and efficacy of histotripsy therapy. This paper reviews and discusses aberration effects in histotripsy and in related therapeutic ultrasound techniques (e.g., high intensity focused ultrasound), with an emphasis on aberration by soft tissues. Methods for aberration correction are reviewed and can be classified into two groups: model-based methods, which use segmented images of the tissue as input to an acoustic propagation model to predict and compensate phase differences, and signal-based methods, which use a receive-capable therapy array to detect phase differences by sensing acoustic signals backpropagating from the focus. The relative advantages and disadvantages of both groups of methods are discussed. Importantly, model-based methods can correct focal shift, while signal-based methods can restore substantial focal pressure, suggesting that both methods should be combined in a 2-step approach. Aberration correction will be critical to improving histotripsy treatments and expanding the histotripsy treatment envelope to enable non-invasive, non-thermal histotripsy therapy for more patients.

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

BabelBrain: An Open-Source Application for Prospective Modeling of Transcranial Focused Ultrasound for Neuromodulation Applications

BabelBrain is an open-source standalone graphic-user-interface application designed for studies of neuromodulation using transcranial-focused ultrasound (FUS). It calculates the transmitted acoustic field in the brain tissue, taking into account the distortion effects caused by the skull barrier. The simulation is prepared using scans from magnetic resonance imaging (MRI) and, if available, computed tomography (CT) and zero-echo time MRI scans. It also calculates the thermal effects based on a given ultrasound regime, such as the total duration of exposure, the duty cycle, and acoustic intensity. The tool is designed to work in tandem with neuronavigation and visualization software, such as 3-DSlicer. It uses image processing to prepare domains for ultrasound simulation and uses the BabelViscoFDTD library for transcranial modeling calculations. BabelBrain supports multiple GPU backends, including Metal, OpenCL, and CUDA, and works on all major operating systems including Linux, macOS, and Windows. This tool is particularly optimized for Apple ARM64 systems, which are common in brain imaging research. The article presents the modeling pipeline used in BabelBrain and a numerical study where different methods of acoustic properties mapping were tested to select the best method that can reproduce the transcranial pressure transmission efficiency reported in the literature.

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.

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.

Transcranial ultrasound neuromodulation of the thalamic visual pathway in a large animal model and the dose-response relationship with MR-ARFI

Neuromodulation of deep brain structures via transcranial ultrasound stimulation (TUS) is a promising, but still elusive approach to non-invasive treatment of brain disorders. The purpose of this study was to confirm that MR-guided TUS of the lateral geniculate nucleus (LGN) can modulate visual evoked potentials (VEPs) in the intact large animal; and to study the impact on cortical brain oscillations. The LGN on one side was identified with T2-weighted MRI in sheep (all male, n = 9). MR acoustic radiation force imaging (MR-ARFI) was used to confirm localization of the targeted area in the brain. Electroencephalographic (EEG) signals were recorded, and the visual evoked potential (VEP) peak-to-peak amplitude (N70 and P100) was calculated for each trial. Time-frequency spectral analysis was performed to elucidate the effect of TUS on cortical brain dynamics. The VEP peak-to-peak amplitude was reversibly suppressed relative to baseline during TUS. Dynamic spectral analysis demonstrated a change in cortical oscillations when TUS is paired with visual sensory input. Sonication-associated microscopic displacements, as measured by MR-ARFI, correlated with the TUS-mediated suppression of visual evoked activity. TUS non-invasively delivered to LGN can neuromodulate visual activity and oscillatory dynamics in large mammalian brains.

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.

Towards controlled drug delivery in brain tumors with microbubble-enhanced focused ultrasound

Brain tumors are particularly challenging malignancies, due to their location in a structurally and functionally distinct part of the human body - the central nervous system (CNS). The CNS is separated and protected by a unique system of brain and blood vessel cells which together prevent most bloodborne therapeutics from entering the brain tumor microenvironment (TME). Recently, great strides have been made through microbubble (MB) ultrasound contrast agents in conjunction with ultrasound energy to locally increase the permeability of brain vessels and modulate the brain TME. As we elaborate in this review, this physical method can effectively deliver a wide range of anticancer agents, including chemotherapeutics, antibodies, and nanoparticle drug conjugates across a range of preclinical brain tumors, including high grade glioma (glioblastoma), diffuse intrinsic pontine gliomas, and brain metastasis. Moreover, recent evidence suggests that this technology can promote the effective delivery of novel immunotherapeutic agents, including immune check-point inhibitors and chimeric antigen receptor T cells, among others. With early clinical studies demonstrating safety, and several Phase I/II trials testing the preclinical findings underway, this technology is making firm steps towards shaping the future treatments of primary and metastatic brain cancer. By elaborating on its key components, including ultrasound systems and MB technology, along with methods for closed-loop spatial and temporal control of MB activity, we highlight how this technology can be tuned to enable new, personalized treatment strategies for primary brain malignancies and brain metastases.