Numerical investigation of the energy distribution of Low-intensity transcranial focused ultrasound neuromodulation for hippocampus

Focused ultrasound heating in brain tissue/skull phantoms with 1 MHz single-element transducer

PURPOSE

The study aims to provide insights on the practicality of using single-element transducers for transcranial Focused Ultrasound (tFUS) thermal applications.

METHODS

FUS sonications were performed through skull phantoms embedding agar-based tissue mimicking gels using a 1 MHz single-element spherically focused transducer. The skull phantoms were 3D printed with Acrylonitrile Butadiene Styrene (ABS) and Resin thermoplastics having the exact skull bone geometry of a healthy volunteer. The temperature field distribution during and after heating was monitored in a 3 T Magnetic Resonance Imaging (MRI) scanner using MR thermometry. The effect of the skull's thickness on intracranial heating was investigated.

RESULTS

A single FUS sonication at focal acoustic intensities close to 1580 W/cm2 for 60 s in free field heated up the agar phantom to ablative temperatures reaching about 90 °C (baseline of 37 °C). The ABS skull strongly blocked the ultrasonic waves, resulting in zero temperature increase within the phantom. Considerable heating was achieved through the Resin skull, but it remained at hyperthermia levels. Conversely, tFUS through a 1 mm Resin skull showed enhanced ultrasonic penetration and heating, with the focal temperature reaching 70 °C.

CONCLUSIONS

The ABS skull demonstrated poorer performance in terms of tFUS compared to the Resin skull owing to its higher ultrasonic attenuation and porosity. The thin Resin phantom of 1 mm thickness provided an efficient acoustic window for delivering tFUS and heating up deep phantom areas. The results of such studies could be particularly useful for accelerating the establishment of a wider range of tFUS applications.

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

BACKGROUND AND OBJECTIVES

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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

An efficient method for transcranial ultrasound focus correction based on the coupling of boundary integrals and finite elements

Transcranial focused ultrasound is a novel technique for the noninvasive treatment of brain diseases. The success of the treatment greatly depends on achieving precise and efficient intraoperative focus. However, compensating for aberrated ultrasound waves caused by the skull through numerical simulation-based phase corrections is a challenging task due to the significant computational burden involved in solving the acoustic wave equation. In this article, we propose a promising strategy using the coupling of the boundary integral equation method (BIEM) and the finite element method (FEM) to overcome the above limitation. Specifically, we adopt the BIEM to obtain the Robin-to-Dirichlet maps on the boundaries of the skull and then couple the maps to the FEM matrices via a dual interpolation technique, resulting in a computational domain including only the skull. Three simulation experiments were conducted to evaluate the effectiveness of the proposed method, including a convergence test and two skull-induced aberration corrections in 2D and 3D ultrasound. The results show that the method's convergence is guaranteed as the element size decreases, leading to a decrease in pressure error. The computation times for simulating a 500 kHz ultrasound field on a regular desktop computer were found to be 0.47 ± 0.01 s in the 2D case and 43.72 ± 1.49 s in the 3D case, provided that lower-upper decomposition (approximately 13 s in 2D and 2.5 h in 3D) was implemented in advance. We also demonstrated that more accurate transcranial focusing can be achieved by phase correction compared to the noncorrected results (with errors of 1.02 mm vs. 6.45 mm in 2D and 0.28 mm vs. 3.07 mm in 3D). The proposed strategy is valuable for enabling online ultrasound simulations during treatment, facilitating real-time adjustments and interventions.

Altered hippocampal functional connectivity after the rupture of anterior communicating artery aneurysm

Background and purpose

Aneurysmal subarachnoid hemorrhage (SAH) predisposes hippocampal injury, a major cause of follow-up cognitive impairment. Our previous study has revealed an abnormal resting-state brain network in patients after the rupture of anterior communicating artery (ACoA) aneurysm. However, the functional connectivity (FC) characteristics of the hippocampus and its relationship with cognitive performance in these patients remain unknown.

Methods

This study ultimately included 26 patients and 19 age- and sex-matched controls who completed quality control for resting-state functional magnetic resonance imaging (fMRI). The mean time series for each side of the hippocampus was extracted from individuals and then a seed-to-voxel analysis was performed. We compared the difference in FC strength between the two groups and subsequently analyzed the correlations between abnormal FC and their cognitive performance.

Results

The results of bilateral hippocampus-based FC analysis were largely consistent. Compared with the healthy controls, patients after the rupture of ACoA aneurysm exhibited significantly decreased FC between the hippocampus and other brain structures within the Papez circuit, including bilateral anterior and middle cingulate cortex (MCC), bilateral medial superior frontal gyrus, and left inferior temporal gyrus (ITG). Instead, increased FC between the hippocampus and bilateral insula was observed. Correlation analyses showed that more subjective memory complaints or lower total cognitive scores were associated with decreased connectivity in the hippocampus and several brain regions such as left anterior cingulate cortex (ACC) and frontotemporal cortex.

Conclusion

These results extend our previous findings and suggest that patients with ruptured ACoA aneurysm exist hypoconnectivity between the hippocampus and multiple brain regions within the Papez circuit. Deactivation of the Papez circuit may be a crucial neural mechanism related to cognitive deficits in patients after the rupture of ACoA aneurysm.