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

Non-Invasive Blood-Brain Barrier Disruption Using Acoustic Holography With a Clinical Focused Ultrasound System

OBJECTIVE

Holographic methods can be used with phased array transducers to shape an ultrasound field. We tested a simple method to create holograms with a hemispherical 1024-element phased array transducer and explored how it could benefit ultrasound-mediated blood-brain barrier (BBB) disruption.

METHODS

With this method, individual acoustic simulations for each element of the transducer were simultaneously loaded into computer memory. Each element's phase was systematically modulated until the combined field matched a desired pattern. The method was evaluated with a 220 kHz transducer being tested clinically to enhance drug delivery via BBB disruption. The holograms were evaluated in a tissue-mimicking phantom and in vivo in experiments disrupting the BBB in rats and in a macaque. We also explored whether this approach could mitigate secondary reflections from the skull using simulations of transcranial focusing in clinical treatments of transcranial sonication for BBB disruption.

RESULTS

This approach can enlarge the focal volume in a patient-specific manner and could reduce the number of sonication targets needed to disrupt large volumes, improve the homogeneity of the disruption, and improve our ability to detect microbubble activity in tissues with low vascular density. Simulations suggest that the method could also mitigate secondary reflections during transcranial sonication.

Scalp Nerve Block Alleviates Headaches Associated With Sonication During Transcranial Magnetic Resonance-Guided Focused Ultrasound

BACKGROUND AND OBJECTIVES

In magnetic resonance-guided focused ultrasound (MRgFUS) procedures, headache is a frequent symptom and cause of treatment discontinuation. Herein, we assessed the efficacy of scalp nerve block (SNB) for alleviating headache during MRgFUS procedures.

METHODS

The effect of SNB on intraprocedural headache was examined by retrospectively comparing 2 patient cohorts at a single institution. During the study period from April 2020 to February 2022, an SNB protocol for all patients with a skull density ratio ≤0.55 was instituted on October 6, 2021. The number of patients with a skull density ratio ≤0.55 was 34 before the protocol and 36 afterward. Headache intensity was evaluated using a numerical rating scale (NRS) after each sonication. To evaluate the effect of SNB on headache intensity, multiple regression analysis was performed per patient and per sonication. In the per-patient analysis, the effect of SNB was evaluated using the maximum NRS, mean NRS, and NRS at the first ultrasound exposure that reached 52.5°C. In the per-sonication analysis, the effect of SNB was evaluated not only for the entire sonication but also for sonications classified into ≤9999 J, 10 000 to 29 999 J, and ≥30 000 J energy doses.

RESULTS

With SNB, headache alleviation was observed in the NRS after the first sonication that reached 52.5°C in each patient (β = -2.40, 95% CI -4.05 to -0.758, P = .00499), in the NRS when all sonications were evaluated (β = -0.647, 95% CI -1.19 to -0.106, P = .0201), and in the NRS when all sonications were classified into 10 000 to 29 999 J (β = -1.83, 95% CI -3.17 to -0.485, P = .00889).

CONCLUSION

SNB significantly reduced headache intensity during MRgFUS, especially that caused by sonication with a moderate-energy dose. These findings suggest that scalp nerves play a role in headache mechanisms during MRgFUS.

CMUT as a Transmitter for Microbubble-Assisted Blood-Brain Barrier Opening

Focused ultrasound (FUS) combined with microbubbles (MBs) has emerged as a promising strategy for transiently opening the blood-brain barrier (BBB) to enhance drug permeability in the brain. Current FUS systems for BBB opening use piezoelectric transducers as transmitters and receivers. While capacitive micromachined ultrasonic transducers (CMUTs) have been suggested as an FUS receiver alternative due to their broad bandwidth, their capabilities as transmitters have not been investigated. This is mainly due to the intrinsic nonlinear behavior of CMUTs, which complicates the detection of MB generated harmonic signals and their low-pressure output at FUS frequencies. Various methods have been proposed to mitigate CMUT nonlinearity; however, these approaches have primarily targeted contrast enhanced ultrasound imaging. In this study, we propose the use of polyphase modulation (PM) technique to isolate MB emissions when CMUTs are employed as transmitters for BBB opening. Our calculations for a human scale FUS system with multiple CMUT transmitters show that 10-kPa peak negative pressure (PNP) at 150-mm focal distance will be sufficient for MB excitation for BBB opening. Experimental findings indicate that this pressure level can be easily generated at 400-800 kHz using a readily available CMUT. Furthermore, more than 50-dB suppression of the fundamental harmonic signal is obtained in free field and transcranial hydrophone measurements by processing receive signals in response to phase-modulated transmit waveforms. In vitro validation of PM is also conducted using Definity MB flowing through a tube phantom. MB-filled tube phantoms show adequate nonlinear signal isolation and SNR for MB harmonic detection. Together our findings indicate that PM can effectively mitigate CMUT harmonic generation, thereby creating new opportunities for wideband transmission and receive operation for BBB opening in clinical and preclinical applications.

Strategies and safety simulations for ultrasonic cervical spinal cord neuromodulation

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

MR-guided focused ultrasound in movement disorders and beyond: Lessons learned and new frontiers

The development of MR-guided focused ultrasound (MRgFUS) has provided a new therapeutic tool for neuropsychiatric disorders. In contrast to previously available neurosurgical techniques, MRgFUS allows precise impact on deep brain structures without the need for incision and yields an immediate effect. In its high-intensity modality (MRgHIFU), it produces accurate therapeutic thermoablation in previously selected brain targets. Importantly, the production of the lesion is progressive and highly controlled in real-time by both neuroimaging and clinical means. MRgHIFU ablation is already an accepted and widely used treatment for medically-refractory Parkinson's disease and essential tremor. Notably, other neurological disorders and diverse brain targets, including bilateral treatments, are currently under examination. Conversely, the low-intensity modality (MRgLIFU) shows promising prospects in neuromodulation and transient blood-brain barrier opening (BBBO). In the former circumstance, MRgLIFU could serve as a powerful clinical and research tool for non-invasively modulating brain activity and function. BBBO, on the other hand, emerges as a potentially impactful method to influence disease pathogenesis and progression by increasing brain target engagement of putative therapeutic agents. While promising, these applications remain experimental. As a recently developed technology, MRgFUS is not without challenges and questions to be addressed. Further developments and broader experience are necessary to enhance MRgFUS capabilities in both research and clinical practice, as well as to define device constraints. This clinical mini-review aims to provide an overview of the main evidence of MRgFUS application and to highlight unmet needs and future potentialities of the technique.

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.

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.

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

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

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

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

A simulation study on the sensitivity of transcranial ray-tracing ultrasound modeling to skull properties

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

Influence of cerebrospinal fluid on power absorption during transcranial magnetic resonance-guided focused ultrasound treatment

BACKGROUND

Ultrasound beam aberration correction is vital when focusing ultrasound through the skull bone in transcranial magnetic resonance-guided focused ultrasound (tcMRgFUS) applications. Current methods make transducer element phase adjustments to compensate for the variation in skull properties (shape, thickness, and acoustic properties), but do not account for variations in the internal brain anatomy.

PURPOSE

Our objective is to investigate the effect of cerebrospinal fluid (CSF) and brain anatomy on beam focusing in tcMRgFUS treatments.

METHODS

Simulations were conducted with imaging data from 20 patients previously treated with focused ultrasound for disabling tremor. The Hybrid Angular Spectrum (HAS) method was used to test the effect of including cerebral spinal fluid (CSF) and brain anatomy in determining the element phases used for aberration correction and beam focusing. Computer tomography (CT) and magnetic resonance imaging (MRI) images from patient treatments were used to construct a segmented model of each patient's head. The segmented model for treatment simulation consisted of water, skin, fat, brain, CSF, diploë, and cortical bone. Transducer element phases used for treatment simulation were determined using time reversal from the desired focus, generating a set of phases assuming a homogeneous brain in the intracranial volume, and a second set of phases assigning CSF acoustic properties to regions of CSF. In addition, for three patients, the relative effect of separately including CSF speed of sound values compared to CSF attenuation values was found.

RESULTS

We found that including CSF acoustic properties (speed of sound and attenuation) during phase planning compared to phase correction without considering CSF increased the absorbed ultrasound power density ratios at the focus over a range of 1.06 to 1.29 (mean of 17% ± 6%) for 20 patients. Separately considering the CSF speed of sound and CSF attenuation showed that the increase was due almost entirely to including the CSF speed of sound; considering only the CSF attenuation had a negligible effect.

CONCLUSIONS

Based on HAS simulations, treatment planning phase determination using morphologically realistic CSF and brain anatomy yielded an increase of up to 29% in the ultrasound focal absorbed power density. Future work will be required to validate the CSF simulations.

A numerical investigation of passive acoustic mapping for monitoring bubble-mediated focused ultrasound treatment of the spinal cord

Focused ultrasound (FUS) combined with intravenous microbubbles (MBs) has been shown to increase drug delivery to the spinal cord in animal models. Eventual clinical translation of such a technique in the sensitive spinal cord requires robust treatment monitoring to ensure efficacy, localization, safety, and provide key intraprocedural feedback. Here, the use of passive acoustic mapping (PAM) of MB emissions with a spine-specific detector array in the context of transvertebral FUS sonications is investigated in silico. Using computed tomography-derived human vertebral geometry, transvertebral detection of MBs is evaluated over varying source locations with and without phase and amplitude corrections (PACs). The impact of prefocal cavitation is studied by simulating concurrent cavitation events in the canal and pre-laminar region. Spatially sensitive application of phase and amplitude is used to balance signal strengths emanating from different axial depths in combination with multiple dynamic ranges to elicit multisource viewing. Collectively, the results of this study encourage the use of PAM in transvertebral FUS applications with PACs to not only localize sources originating in the spinal canal but also multiple sources of innate amplitude mismatches when corrective methods are applied.

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.

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

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

A Comparative Feasibility Study for Transcranial Extracorporeal Shock Wave Therapy

The potential beneficial regenerative and stimulatory extracorporeal shock wave therapy (ESWT) applications to the central nervous system have garnered interest in recent years. Treatment zones for these indications are acoustically shielded by bones, which heavily impact generated sound fields. We present the results of high-resolution tissue-realistic simulations, comparing the viability of different ESWT applicators in their use for transcranial applications. The performances of electrohydraulic, electromagnetic, and piezoelectric transducers for key reflector geometries are compared. Based on density information obtained from CT imaging of the head, we utilized the non-linear wave propagation toolset Matlab k-Wave to obtain spatial therapeutic sound field geometries and waveforms. In order to understand the reliability of results on the appropriate modeling of the skull, three different bone attenuation models were compared. We find that all currently clinically ESWT applicator technologies show significant retention of peak pressures and energies past the bone barrier. Electromagnetic transducers maintain a significantly higher energy flux density compared to other technologies while low focusing strength piezoelectric applicators have the weakest transmissions. Attenuation estimates provide insights into sound field degradation and energy losses, indicating that effective transcranial therapies can readily be attained with current applicators. Furthermore, the presented approach will allow for future targeted in silico development and the design of applicators and therapy plans to ultimately improve therapeutic outcomes.

Two-step aberration correction: application to transcranial histotripsy

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

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.

Lesion location and lesion creation affect outcomes after focused ultrasound thalamotomy

MRI-guided focused ultrasound thalamotomy has been shown to be an effective treatment for medication refractory essential tremor. Here, we report a clinical-radiological analysis of 123 cases of MRI-guided focused ultrasound thalamotomy, and explore the relationships between treatment parameters, lesion characteristics and outcomes. All patients undergoing focused ultrasound thalamotomy by a single surgeon were included. The procedure was performed as previously described, and patients were followed for up to 1 year. MRI was performed 24 h post-treatment, and lesion locations and volumes were calculated. We retrospectively evaluated 118 essential tremor patients and five tremor-dominant Parkinson's disease patients who underwent thalamotomy. At 24 h post-procedure, tremor abated completely in the treated hand in 81 essential tremor patients. Imbalance, sensory disturbances and dysarthria were the most frequent acute adverse events. Patients with any adverse event had significantly larger lesions, while inferolateral lesion margins were associated with a higher incidence of motor-related adverse events. Twenty-three lesions were identified with irregular tails, often extending into the internal capsule; 22 of these patients experienced at least one adverse event. Treatment parameters and lesion characteristics changed with increasing surgeon experience. In later cases, treatments used higher maximum power (normalized to skull density ratio), accelerated more quickly to high power, and delivered energy over fewer sonications. Larger lesions were correlated with a rapid rise in both power delivery and temperature, while increased oedema was associated with rapid rise in temperature and the maximum power delivered. Total energy and total power did not significantly affect lesion size. A support vector regression was trained to predict lesion size and confirmed the most valuable predictors of increased lesion size as higher maximum power, rapid rise to high-power delivery, and rapid rise to high tissue temperatures. These findings may relate to a decrease in the energy efficiency of the treatment, potentially due to changes in acoustic properties of skull and tissue at higher powers and temperatures. We report the largest single surgeon series of focused ultrasound thalamotomy to date, demonstrating tremor relief and adverse events consistent with reported literature. Lesion location and volume impacted adverse events, and an irregular lesion tail was strongly associated with adverse events. High-power delivery early in the treatment course, rapid temperature rise, and maximum power were dominant predictors of lesion volume, while total power, total energy, maximum energy and maximum temperature did not improve prediction of lesion volume. These findings have critical implications for treatment planning in future patients.

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

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

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

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

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

Transcranial magnetic resonance-guided focused ultrasound (tcMRgFUS) is gaining significant acceptance as a noninvasive treatment for motion disorders and shows promise for novel applications such as blood-brain barrier opening for tumor treatment. A typical procedure relies on CT-derived acoustic property maps to simulate the transfer of ultrasound through the skull. Accurate estimates of the acoustic attenuation in the skull are essential to accurate simulations, but there is no consensus about how attenuation should be estimated from CT images and there is interest in exploring MR as a predictor of attenuation in the skull. In this study, we measure the acoustic attenuation at 0.5, 1, and 2.25 MHz in 89 samples taken from two ex vivo human skulls. CT scans acquired with a variety of X-ray energies, reconstruction kernels, and reconstruction algorithms, and MR images acquired with ultrashort and zero echo time sequences are used to estimate the average Hounsfield unit value, MR magnitude, and T2^* value in each sample. The measurements are used to develop a model of attenuation as a function of frequency and each individual imaging parameter.

Transcranial focused ultrasound phase correction using the hybrid angular spectrum method

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

Technical Principles and Clinical Workflow of Transcranial MR-Guided Focused Ultrasound

Transcranial MR-guided focused ultrasound (MRgFUS) is a rapidly developing technology in neuroscience for manipulating brain structure and function without open surgery. The effectiveness of transcranial MRgFUS for thermoablation is well established, and the technique is actively employed worldwide for movement disorders including essential tremor. A growing number of centers are also investigating the potential of microbubble-mediated focused ultrasound-induced opening of the blood-brain barrier (BBB) for targeted drug delivery to the brain. Here, we provide a technical overview of the principles, clinical workflow, and operator considerations of transcranial MRgFUS procedures for both thermoablation and BBB opening.

Predicting Bone Marrow Damage in the Skull After Clinical Transcranial MRI-Guided Focused Ultrasound With Acoustic and Thermal Simulations

Transcranial MRI-guided focused ultrasound (TcMRgFUS) thermal ablation is a noninvasive functional neurosurgery technique. Previous reports have shown that damage in the skull bone marrow can occur at high acoustic energies. While this damage is asymptomatic, it would be desirable to avoid it. Here we examined whether acoustic and thermal simulations can predict where the thermal lesions in the marrow occurred. Post-treatment imaging was obtained at 3-15 months after 40 clinical TcMRgFUS procedures, and bone marrow lesions were observed after 16 treatments. The presence of lesions was predicted by the acoustic energy with a threshold of 18.1-21.1 kJ (maximum acoustic energy used) and 97-112 kJ (total acoustic energy applied over the whole treatment). The size of the lesions was not always predicted by the acoustic energy used during treatment alone. In contrast, the locations, sizes, and shapes of the heated regions estimated by the acoustic and thermal simulations were qualitatively similar to those of the lesions. The lesions generally appeared in areas that were predicted to have high temperatures. While more work is needed to validate the temperature estimates in and around the skull, being able to predict the locations and onset for lesions in the bone marrow could allow for better distribution of the acoustic energy over the skull. Understanding skull absorption characteristics of TcMRgFUS could also be useful in optimizing transcranial focusing.