mSOUND: An Open Source Toolbox for Modeling Acoustic Wave Propagation in Heterogeneous Media

Acoustic-holography-patterned primary hepatocytes possess liver functions

Acoustic holography (AH), a promising approach for cell patterning, emerges as a powerful tool for constructing novel invitro 3D models that mimic organs and cancers features. However, understanding changes in cell function post-AH remains limited. Furthermore, replicating complex physiological and pathological processes solely with cell lines proves challenging. Here, we employed acoustical holographic lattice to assemble primary hepatocytes directly isolated from mice into a cell cluster matrix to construct a liver-shaped tissue sample. For the first time, we evaluated the liver functions of AH-patterned primary hepatocytes. The patterned model exhibited large numbers of self-assembled spheroids and superior multifarious core hepatocyte functions compared to cells in 2D and traditional 3D culture models. AH offers a robust protocol for long-term in vitro culture of primary cells, underscoring its potential for future applications in disease pathogenesis research, drug testing, and organ replacement therapy.

Low-intensity focused ultrasound to the posterior insula reduces temporal summation of pain

BACKGROUND

The insula and dorsal anterior cingulate cortex (dACC) are core brain regions involved in pain processing and central sensitization, a shared mechanism across various chronic pain conditions. Methods to modulate these regions may serve to reduce central sensitization, though it is unclear which target may be most efficacious for different measures of central sensitization.

OBJECTIVE/HYPOTHESIS

Investigate the effect of low-intensity focused ultrasound (LIFU) to the anterior insula (AI), posterior insula (PI), or dACC on conditioned pain modulation (CPM) and temporal summation of pain (TSP).

METHODS

N = 16 volunteers underwent TSP and CPM pain tasks pre/post a 10 min LIFU intervention to either the AI, PI, dACC or Sham stimulation. Pain ratings were collected pre/post LIFU.

RESULTS

Only LIFU to the PI significantly attenuated pain ratings during the TSP protocol. No effects were found for the CPM task for any of the LIFU targets. LIFU pressure modulated group means but did not affect overall group differences.

CONCLUSIONS

LIFU to the PI reduced temporal summation of pain. This may, in part, be due to dosing (pressure) of LIFU. Inhibition of the PI with LIFU may be a future potential therapy in chronic pain populations demonstrating central sensitization. The minimal effective dose of LIFU for efficacious neuromodulation will help to translate LIFU for therapeutic options.

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.

Enabling strain imaging in realistic Eulerian ultrasound simulation methods

Cardiovascular strain imaging is continually improving due to ongoing advances in ultrasound acquisition and data processing techniques. The phantoms used for validation of new methods are often burdensome to make and lack flexibility to vary mechanical and acoustic properties. Simulations of US imaging provide an alternative with the required flexibility and ground truth strain data. However, the current Lagrangian US strain imaging models cannot simulate heterogeneous speed of sound distributions and higher-order scattering, which limits the realism of the simulations. More realistic Eulerian modelling techniques exist but have so far not been used for strain imaging. In this research, a novel sampling scheme was developed based on a band-limited interpolation of the medium, which enables accurate strain simulation in Eulerian methods. The scheme was validated in k-Wave using various numerical phantoms and by a comparison with Field II. The method allows for simulations with a large range in strain values and was accurate with errors smaller than -60 dB. Furthermore, an excellent agreement with the Fourier theory of US scattering was found. The ability to perform simulations with heterogeneous speed of sound distributions was demonstrated using a pulsating artery model. The developed sampling scheme contributes to more realistic strain imaging simulations, in which the effect of heterogenous acoustic properties can be taken into account.

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.

Nonlinear Acoustic Holography With Adaptive Sampling

Accurate and efficient numerical simulation of highly nonlinear ultrasound propagation is essential for a wide range of therapeutic and physical ultrasound applications. However, due to large domain sizes and the generation of higher harmonics, such simulations are computationally challenging, particularly in 3-D problems with shock waves. Current numerical methods are based on computationally inefficient uniform meshes that resolve the highest harmonics across the entire spatial domain. To address this challenge, we present an adaptive numerical algorithm for computationally efficient nonlinear acoustic holography. At each propagation step, the algorithm monitors the harmonic content of the acoustic signal and adjusts its discretization parameters accordingly. This enables efficient local resolution of higher harmonics in areas of high nonlinearity while avoiding unnecessary resolution elsewhere. Furthermore, the algorithm actively adapts to the signal's nonlinearity level, eliminating the need for prior reference simulations or information about the spatial distribution of the harmonic content of the acoustic field. The proposed algorithm incorporates an upsampling process in the frequency domain to accommodate the generation of higher harmonics in forward propagation and a downsampling process when higher harmonics are decimated in backward propagation. The efficiency of the algorithm was evaluated for highly nonlinear 3-D problems, demonstrating a significant reduction in computational cost with a nearly 50-fold speedup over a uniform mesh implementation. Our findings enable a more rapid and efficient approach to modeling nonlinear high-intensity focused ultrasound (HIFU) wave propagation.

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.

Corrections to "mSOUND: An Open Source Toolbox for Modeling Acoustic Wave Propagation in Heterogeneous Media"

In the above article [1]>, the authors regret that there were some mistakes pertaining to (1)-(3) and (5).