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

The role of focused ultrasound for pediatric brain tumors: current insights and future implications on treatment strategies

INTRODUCTION

Focused ultrasound (FUS) is an innovative and emerging technology for the treatment of adult and pediatric brain tumors and illustrates the intersection of various specialized fields, including neurosurgery, neuro-oncology, radiation oncology, and biomedical engineering.

OBJECTIVE

The authors provide a comprehensive overview of the application and implications of FUS in treating pediatric brain tumors, with a special focus on pediatric low-grade gliomas (pLGGs) and the evolving landscape of this technology and its clinical utility.

METHODS

The fundamental principles of FUS include its ability to induce thermal ablation or enhance drug delivery through transient blood-brain barrier (BBB) disruption, emphasizing the adaptability of high-intensity focused ultrasound (HIFU) and low-intensity focused ultrasound (LIFU) applications.

RESULTS

Several ongoing clinical trials explore the potential of FUS in offering alternative therapeutic strategies for pathologies where conventional treatments fall short, specifically centrally-located benign CNS tumors and diffuse intrinsic pontine glioma (DIPG). A case illustration involving the use of HIFU for pilocytic astrocytoma is presented.

CONCLUSION

Discussions regarding future applications of FUS for the treatment of gliomas include improved drug delivery, immunomodulation, radiosensitization, and other technological advancements.

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

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

Transcranial Ultrasonic Focusing by a Phased Array Based on Micro-CT Images

In this paper, we utilize micro-computed tomography (micro-CT) to obtain micro-CT images with a resolution of 60 μm and establish a micro-CT model based on the k-wave toolbox, which can visualize the microstructures in trabecular bone, including pores and bone layers. The transcranial ultrasound phased array focusing field characteristics in the micro-CT model are investigated. The ultrasonic waves are multiply scattered in skull and time delays calculations from the transducer to the focusing point are difficult. For this reason, we adopt the pulse compression method and the linear frequency modulation Barker code to compute the time delay and implement phased array focusing in the micro-CT model. It is shown by the simulation results that ultrasonic loss is mainly caused by scattering from the microstructures of the trabecular bone. The ratio of main and side lobes of the cross-correlation calculation is improved by 5.53 dB using the pulse compression method. The focusing quality and the calculation accuracy of time delay are improved. Meanwhile, the beamwidth at the focal point and the sound pressure amplitude decrease with the increase in the signal frequency. Focusing at different depths indicates that the beamwidth broadens with the increase in the focusing depth, and beam deflection focusing maintains good consistency in the focusing effect at a distance of 9 mm from the focal point. This indicates that the phased-array method has good focusing results and focus tunability in deep cranial brain. In addition, the sound pressure at the focal point can be increased by 8.2% through amplitude regulation, thereby enhancing focusing efficiency. The preliminary experiment verification is conducted with an ex vivo skull. It is shown by the experimental results that the phased array focusing method using pulse compression to calculate the time delay can significantly improve the sound field focusing effect and is a very effective transcranial ultrasound focusing method.

Real-Time Acoustic Simulation Framework for tFUS: A Feasibility Study Using Navigation System

Transcranial focused ultrasound (tFUS), in which acoustic energy is focused on a small region in the brain through the skull, is a non-invasive therapeutic method with high spatial resolution and depth penetration. Image-guided navigation has been widely utilized to visualize the location of acoustic focus in the cranial cavity. However, this system is often inaccurate because of the significant aberrations caused by the skull. Therefore, acoustic simulations using a numerical solver have been widely adopted to compensate for this inaccuracy. Although the simulation can predict the intracranial acoustic pressure field, real-time application during tFUS treatment is almost impossible due to the high computational cost. In this study, we propose a neural network-based real-time acoustic simulation framework and test its feasibility by implementing a simulation-guided navigation (SGN) system. Real-time acoustic simulation is performed using a 3D conditional generative adversarial network (3D-cGAN) model featuring residual blocks and multiple loss functions. This network was trained by the conventional numerical acoustic simulation program (i.e., k-Wave). The SGN system is then implemented by integrating real-time acoustic simulation with a conventional image-guided navigation system. The proposed system can provide simulation results with a frame rate of 5 Hz (i.e., about 0.2 s), including all processing times. In numerical validation (3D-cGAN vs. k-Wave), the average peak intracranial pressure error was 6.8 ± 5.5%, and the average acoustic focus position error was 5.3 ± 7.7 mm. In experimental validation using a skull phantom (3D-cGAN vs. actual measurement), the average peak intracranial pressure error was 4.5%, and the average acoustic focus position error was 6.6 mm. These results demonstrate that the SGN system can predict the intracranial acoustic field according to transducer placement in real-time.

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.

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.

Optimizing transcranial ultrasound delivery at large incident angles by leveraging cranial leaky guided wave dispersion

We investigate the role of leaky guided waves in transcranial ultrasound transmission in temporal and parietal bones at large incidence angles. Our numerical and experimental results show that the dispersion characteristics of the fundamental leaky guided wave mode with longitudinal polarization can be leveraged to estimate the critical angle above which efficient shear mode conversion takes place, and below which major transmission drops can be expected. Simulations that employ a numerical propagator matrix and a Semi-Analytical approach establish the transcranial dispersion characteristics and transmission coefficients at different incident angles. Experimental transmission tests conducted at 500 kHz and radiation tests performed in the 200-800 kHz range confirm the numerical findings in terms of transmitted peak pressure and frequency-radiation angle spectra, based on which the connection between critical angles, dispersion and transmission is demonstrated. Our results support the identification of transcranial ultrasound strategies that leverage shear mode conversion, which is less sensitive to phase aberrations compared to normal incidence ultrasound. These findings can also enable higher transmission rates in cranial bones with low porosity by leveraging dispersion information extracted through signal processing, without requiring measurement of geometric and mechanical properties of the cranial bone.

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.

Binary acoustic metasurfaces for dynamic focusing of transcranial ultrasound

Transcranial focused ultrasound (tFUS) is a promising technique for non-invasive and spatially targeted neuromodulation and treatment of brain diseases. Acoustic lenses were designed to correct the skull-induced beam aberration, but these designs could only generate static focused ultrasound beams inside the brain. Here, we designed and 3D printed binary acoustic metasurfaces (BAMs) for skull aberration correction and dynamic ultrasound beam focusing. BAMs were designed by binarizing the phase distribution at the surface of the metasurfaces. The phase distribution was calculated based on time reversal to correct the skull-induced phase aberration. The binarization enabled the ultrasound beam to be dynamically steered along wave propagation direction by adjusting the operation frequency of the incident ultrasound wave. The designed BAMs were manufactured by 3D printing with two coding bits, a polylactic acid unit for bit “1” and a water unit for bit “0.” BAMs for single- and multi-point focusing through the human skull were designed, 3D printed, and validated numerically and experimentally. The proposed BAMs with subwavelength scale in thickness are simple to design, easy to fabric, and capable of correcting skull aberration and achieving dynamic beam steering.

Non-Invasive Thermal Therapy for Tissue Engineering and Regenerative Medicine

Owing to the development of nanotechnology and noninvasive treatment, thermal therapy in combination with external stimuli has been applied for tissue engineering and regenerative medicine (TERM), which has attracted more and more attention in recent years. In this review, the recent progress of applying a variety of non-invasive thermal therapeutic modalities for TERM, including photothermal therapy, magnetic thermotherapy, and ultrasound thermotherapy, as well as other thermal therapeutics are discussed. The parameters and conditions that need to be considered and regulated to realize a well-controlled thermal therapy for tissue regeneration are also discussed. Afterwards, the current concerns and challenges of putting thermal therapy into clinical applications are pointed out. At last, perspectives are provided for the future development directions, aiming to providing opportunities and a novel pathway for TERM.

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

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

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

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

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.

Establishing density-dependent longitudinal sound speed in the vertebral lamina

Focused ultrasound treatments of the spinal cord may be facilitated using a phased array transducer and beamforming to correct spine-induced focal aberrations. Simulations can non-invasively calculate aberration corrections using x-ray computed tomography (CT) data that are correlated to density (ρ) and longitudinal sound speed (c_L). We aimed to optimize vertebral lamina-specific c_L(ρ) functions at a physiological temperature (37 °C) to maximize time domain simulation accuracy. Odd-numbered ex vivo human thoracic vertebrae were imaged with a clinical CT-scanner (0.511 × 0.511 × 0.5 mm), then sonicated with a transducer (514 kHz) focused on the canal via the vertebral lamina. Vertebra-induced signal time shifts were extracted from pressure waveforms recorded within the canals. Measurements were repeated 5× per vertebra, with 2.5 mm vertical vertebra shifts between measurements. Linear functions relating c_L with CT-derived density were optimized. The optimized function was c_L(ρ)=0.35(ρ-ρ_w)+ c_L,w m/s, where w denotes water, giving the tested laminae a mean bulk density of 1600 ± 30 kg/m³ and a mean bulk c_L of 1670 ± 60 m/s. The optimized lamina c_L(ρ) function was accurate to λ/16 when implemented in a multi-layered ray acoustics model. This modelling accuracy will improve trans-spine ultrasound beamforming.