Acoustic Attenuation: Multifrequency Measurement and Relationship to CT and MR Imaging

Ultra-short time-echo based ray tracing for transcranial focused ultrasound aberration correction in human calvaria

Objective.Magnetic resonance guided transcranial focused ultrasound holds great promises for treating neurological disorders. This technique relies on skull aberration correction which requires computed tomography (CT) scans of the skull of the patients. Recently, ultra-short time-echo (UTE) magnetic resonance (MR) sequences have unleashed the MRI potential to reveal internal bone structures. In this study, we measure the efficacy of transcranial aberration correction using UTE images.Approach.We compare the efficacy of transcranial aberration correction using UTE scans to CT based correction on four skulls and two targets using a clinical device (Exablate Neuro, Insightec, Israel). We also evaluate the performance of a custom ray tracing algorithm using both UTE and CT estimates of acoustic properties and compare these against the performance of the manufacturer's proprietary aberration correction software.Main results.UTE estimated skull maps in Hounsfield units (HU) had a mean absolute error of 242 ± 20 HU (n= 4). The UTE skull maps were sufficiently accurate to improve pressure at the target (no correction: 0.44 ± 0.10, UTE correction: 0.79 ± 0.05, manufacturer CT: 0.80 ± 0.05), pressure confinement ratios (no correction: 0.45 ± 0.10, UTE correction: 0.80 ± 0.05, manufacturer CT: 0.81 ± 0.05), and targeting error (no correction: 1.06 ± 0.42 mm, UTE correction 0.30 ± 0.23 mm, manufacturer CT: 0.32 ± 0.22) (n= 8 for all values). When using CT, our ray tracing algorithm performed slightly better than UTE based correction with pressure at the target (UTE: 0.79 ± 0.05, CT: 0.84 ± 0.04), pressure confinement ratios (UTE: 0.80 ± 0.05, CT: 0.84 ± 0.04), and targeting error (UTE: 0.30 ± 0.23 mm, CT: 0.17 ± 0.15).Significance.These 3D transcranial measurements suggest that UTE sequences could replace CT scans in the case of MR guided focused ultrasound with minimal reduction in performance which will avoid ionizing radiation exposure to the patients and reduce procedure time and cost.

Enhancing transcranial ultrasound stimulation planning with MRI-derived skull masks: a comparative analysis with CT-based processing

Objective.Transcranial ultrasound stimulation (TUS) presents challenges in ultrasound wave transmission through the skull, affecting study outcomes due to aberration and attenuation. While planning strategies incorporating 3D computed tomography (CT) scans help mitigate these issues, they expose participants to radiation, which can raise ethical concerns. A solution involves generating skull masks from participants' anatomical magnetic resonance imaging (MRI). This study aims to compare ultrasound field predictions between CT-derived and MRI-derived skull masks in TUS planning.Approach.Five participants with a range of skull density ratios (SDRs: 0.31, 0.42, 0.55, 0.67, and 0.79) were selected, each having both CT and T1/T2-weighted MRI scans. Ultrasound simulations were performed using BabelBrain software with a single-element transducer (diameter = 50 mm,F# = 1) at 250, 500, and 750 kHz frequencies. CT scans were used to generate maps of the skull's acoustic properties. The MRI scans were processed using the Charm segmentation tool from the SimNIBS tool suite using default and custom settings adapted for better skull segmentation. Ultrasound was adjusted to target 30 mm below the skull's surface at 54 electroencephalogram (EEG) locations.Main Results.The custom setting in Charm significantly improved the Dice coefficient between MRI- and CT-derived masks when compared to the default setting (p< 0.001). The maximum pressure error significantly decreased in the custom setting compared to the default setting (p< 0.001). Additionally, the focus location error median across different SDRs averaged 2.32, 1.45, and 1.57 mm in default and 2.08, 1.38, and 1.44 mm in custom conditions for 250 kHz, 500 kHz, and 750 kHz respectively.Significance.MRI-derived skull masks offer satisfactory accuracy at many EEG sites, and using custom settings can further enhance this accuracy. However, significant errors at specific locations highlight the importance of carefully considering stimulation location when choosing between CT- and MRI-derived skull modeling.

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.

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.

Cavitation monitoring, treatment strategy, and acoustic simulations of focused ultrasound blood-brain barrier disruption in patients with glioblastoma

PURPOSE

We report our experience disrupting the blood-brain barrier (BBB) to improve drug delivery in glioblastoma patients receiving temozolomide chemotherapy. The goals of this retrospective analysis were to compare MRI-based measures of BBB disruption and vascular damage to the exposure levels, acoustic emissions data, and acoustic simulations. We also simulated the cavitation detectors.

METHODS

Monthly BBB disruption (BBBD) was performed using a 220 kHz hemispherical phased array focused ultrasound system (Exablate Neuro, InSightec) and Definity microbubbles (Lantheus) over 38 sessions in nine patients. Exposure levels were actively controlled via the cavitation dose obtained by monitoring subharmonic acoustic emissions. The acoustic field and sensitivity profile of the cavitation detection system were simulated. Exposure levels and cavitation metrics were compared to the level of BBBD evident in contrast-enhanced MRI and to hypointense regions in T2*-weighted MRI.

RESULTS

Our treatment strategy evolved from using a relatively high cavitation dose goal to a lower goal and longer sonication duration and ultimately resulted in BBBD across the treatment volume with minimal petechiae. Subsonication-level feedback control of the exposure using acoustic emissions also improved consistency. Simulations of the acoustic field suggest that reflections and standing waves appear when the focus is placed near the skull, but their effects can be mitigated with aberration correction. Simulating the cavitation detectors suggest variations in the sensitivity profile across the treatment volume and between patients. A correlation was observed with the cavitation dose, BBBD and petechial hemorrhage in 8/9 patients, but substantial variability was evident. Analysis of the cavitation spectra found that most bursts did not contain wideband emissions, a signature of inertial cavitation, but biggest contribution to the cavitation dose - the metric used to control the procedure - came from bursts with wideband emissions.

CONCLUSION

Using a low subharmonic cavitation dose with a longer duration resulted in BBBD with minimal petechiae. The correlation between cavitation dose and outcomes demonstrates the benefits of feedback control based on acoustic emissions, although more work is needed to reduce variability. Acoustic simulations could improve focusing near the skull and inform our analysis of acoustic emissions. Monitoring additional frequency bands and improving the sensitivity of the cavitation detection could provide signatures of microbubble activity associated with BBB disruption that were undetected here and could improve our ability to achieve BBB disruption without vascular damage.

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.

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.

Controlled noninvasive modulation of deep brain regions in humans

Transcranial focused ultrasound provides noninvasive and reversible approaches for precise and personalized manipulations of brain circuits, with the potential to transform our understanding of brain function and treatments of brain dysfunction. However, effective applications in humans have been limited by the human head, which attenuates and distorts ultrasound severely and unpredictably. This has led to uncertain ultrasound intensities delivered into the brain. Here, we address this lingering barrier using a direct measurement approach that can be repeatedly applied to the human brain. The approach uses an ultrasonic scan of the head to measure and compensate for the attenuation of the ultrasound by all obstacles within the ultrasound path. No other imaging modality is required and the method is parameter-free and personalized to each subject. The approach accurately restores operators’ intended intensities inside ex-vivo human skulls. Moreover, the approach is critical for effective modulation of deep brain regions in humans. When applied, the approach modulates fMRI Blood Oxygen Level Dependent (BOLD) activity in disease-relevant deep brain regions. This tool unlocks the potential of emerging approaches based on low-intensity ultrasound for controlled manipulations of neural circuits in humans.

Transcranial focused ultrasound has had limited applications in humans due to the unpredictable distortions of ultrasound by the human head. Thomas Riis and colleagues report an approach which enables direct correction for the attenuation of ultrasound by the skull and hair, thus enabling controlled ultrasound therapies in humans.

The influence of bone model geometries on the determination of skull acoustic properties

In this study, we investigated the impact of various simulated skull bone geometries on the determination of skull speed of sound and acoustic attenuation values via optimization using transmitted pressure amplitudes beyond the bone. Using the hybrid angular spectrum method (HAS), we simulated ultrasound transmission through four model sets of different geometries involving sandwiched layers of diploë and cortical bone in addition to three models generated from CT images of ex-vivo human skull-bones. We characterized cost-function solution spaces for each model and, using optimization, found that when a model possessed appreciable variations in resolvable layer thickness, the predefined attenuation coefficients could be found with low error (RMSE < 0.01 Np/cm). However, we identified a spatial frequency cutoff in the models' geometry beyond which the accuracy of the property determination begins to fail, depending on the frequency of the ultrasound source. There was a large increase in error of the attenuation coefficients determined by the optimization when the variations in layer thickness were above the identified spatial frequency cutoffs, or when the lateral variations across the model were relatively low in amplitude. For our limited sample of three CT-image derived bone models, the attenuation coefficients were determined successfully. The speed of sound values were determined with low error for all models (including the CT-image derived models) that were tested (RMSE < 0.4 m/s). These results illustrate that it is possible to determine the acoustic properties of two-component models when the internal bone structure is taken into account and the structure satisfies the spatial frequency constraints discussed.

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.

Transcranial ultrasound simulations: A review

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

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

BACKGROUND

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

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.

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.

Acoustic properties across the human skull

Transcranial ultrasound is emerging as a noninvasive tool for targeted treatments of brain disorders. Transcranial ultrasound has been used for remotely mediated surgeries, transient opening of the blood-brain barrier, local drug delivery, and neuromodulation. However, all applications have been limited by the severe attenuation and phase distortion of ultrasound by the skull. Here, we characterized the dependence of the aberrations on specific anatomical segments of the skull. In particular, we measured ultrasound propagation properties throughout the perimeter of intact human skulls at 500 kHz. We found that the parietal bone provides substantially higher transmission (average pressure transmission 31 ± 7%) and smaller phase distortion (242 ± 44 degrees) than frontal (13 ± 2%, 425 ± 47 degrees) and occipital bone regions (16 ± 4%, 416 ± 35 degrees). In addition, we found that across skull regions, transmission strongly anti-correlated (R=-0.79) and phase distortion correlated (R=0.85) with skull thickness. This information guides the design, positioning, and skull correction functionality of next-generation devices for effective, safe, and reproducible transcranial focused ultrasound therapies.