Modeling nonlinear ultrasound propagation in heterogeneous media with power law absorption using a k-space pseudospectral method

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.

Transcranial Blood-Brain Barrier Opening in Alzheimer's Disease Patients Using a Portable Focused Ultrasound System with Real-Time 2-D Cavitation Mapping

Background

Focused ultrasound (FUS) in combination with microbubbles has recently shown great promise in facilitating blood-brain barrier (BBB) opening for drug delivery and immunotherapy in Alzheimer's disease (AD). However, it is currently limited to systems integrated within the MRI suites or requiring post-surgical implants, thus restricting its widespread clinical adoption. In this pilot study, we investigate the clinical safety and feasibility of a portable, non-invasive neuronavigation-guided FUS (NgFUS) system with integrated real-time 2-D microbubble cavitation mapping.

Methods

A phase 1 clinical study with mild to moderate AD patients (N=6) underwent a single session of microbubble-mediated NgFUS to induce transient BBB opening (BBBO). Microbubble activity under FUS was monitored with real-time 2-D cavitation maps and dosing to ensure the efficacy and safety of the NgFUS treatment. Post-operative MRI was used for BBB opening and closure confirmation as well as safety assessment. Changes in AD biomarker levels in both blood serum and extracellular vesicles (EVs) were evaluated, while changes in amyloid-beta (Aβ) load in the brain were assessed through 18F-Florbetapir PET.

Results

BBBO was achieved in 5 out of 6 subjects with an average volume of 983±626 mm3 following FUS at the right frontal lobe both in white and gray matter regions. The outpatient treatment was completed within 34.8±10.7 min. Cavitation dose significantly correlated with the BBBO volume (R2>0.9, N=4), demonstrating the portable NgFUS system's capability of predicting opening volumes. The cavitation maps co-localized closely with the BBBO location, representing the first report of real-time transcranial 2-D cavitation mapping in the human brain. Larger opening volumes correlated with increased levels of AD biomarkers, including Aβ42 (R2=0.74), Tau (R2=0.95), and P-Tau181 (R2=0.86), assayed in serum-derived EVs sampled 3 days after FUS (N=5). From PET scans, subjects showed a lower Aβ load increase in the treated frontal lobe region compared to the contralateral region. Reduction in asymmetry standardized uptake value ratios (SUVR) correlated with the cavitation dose (R2>0.9, N=3). Clinical changes in the mini-mental state examination over 6 months were within the expected range of cognitive decline with no additional changes observed as a result of FUS.

Conclusion

We showed the safety and feasibility of this cost-effective and time-efficient portable NgFUS treatment for BBBO in AD patients with the first demonstration of real-time 2-D cavitation mapping. The cavitation dose correlated with BBBO volume, a slowed increase in pathology, and serum detection of AD proteins. Our study highlights the potential for accessible FUS treatment in AD, with or without drug delivery.

Nonlinear effects of dual-frequency focused ultrasound on the on-demand regulation of acoustic droplet vaporization

Dual-frequency ultrasound has been widely employed to enhance and regulate acoustic droplet vaporization (ADV) but the role of ultrasonic nonlinear effects on it remains unclear. The main objective of this study is to investigate the influence of nonlinear effects on the control of ADV nucleation under different dual-frequency focused ultrasound conditions. ADV nucleation of PFC nanodroplets activated by nonlinear dual-frequency ultrasound was modeled and parametric studies were conducted to investigate the influence of dual-frequency ultrasound frequency and acoustic power on the degree of nonlinearity (DoN), nucleation rates and dimensions of the nucleation region in a wide parameter range. The results showed that the ultrasonic nonlinearity caused a significant decrease in peak negative pressure due to waveform distortion, which leads to a lower nucleation rate in the nonlinear model compared to that in the linear model. Furthermore, the distributions of nucleation regions were also affected by the interaction between waves of different frequencies and cloud-like spatial distributions were produced, which could be modulated by the dual-frequency ultrasound parameters and have great potentials in the spatial regulation of the ADV and customized treatment protocols in clinical applications. In addition, represented by 1.5 MHz + 3 MHz, such a dual-frequency combination of fundamental and second harmonic could effectively enhance ultrasonic nonlinear effects with relatively lower peak negative pressure and higher DoN. Therefore, nonlinear effect of the dual-frequency ultrasound plays an important role in the ADV regulation, which should be considered in the numerical model and practical applications.

Treatment Planning and Aberration Correction Algorithm for HIFU Ablation of Renal Tumors

High-intensity focused ultrasound (HIFU) applications for thermal or mechanical ablation of renal tumors often encounter challenges due to significant beam aberration and refraction caused by oblique beam incidence, inhomogeneous tissue layers, and presence of gas and bones within the beam. These losses can be significantly mitigated through sonication geometry planning, patient positioning, and aberration correction using multielement phased arrays. Here, a sonication planning algorithm is introduced, which uses the simulations to select the optimal transducer position and evaluate the effect of aberrations and acoustic field quality at the target region after aberration correction. Optimization of transducer positioning is implemented using a graphical user interface (GUI) to visualize a segmented 3-D computed tomography (CT)-based acoustic model of the body and to select sonication geometry through a combination of manual and automated approaches. An HIFU array (1.5 MHz, 256 elements) and three renal cell carcinoma (RCC) cases with different tumor locations and patient body habitus were considered. After array positioning, the correction of aberrations was performed using a combination of backpropagation from the focus with an ordinary least squares (OLS) optimization of phases at the array elements. The forward propagation was simulated using a combination of the Rayleigh integral and k-space pseudospectral method (k-Wave toolbox). After correction, simulated HIFU fields showed tight focusing and up to threefold higher maximum pressure within the target region. The addition of OLS optimization to the aberration correction method yielded up to 30% higher maximum pressure compared to the conventional backpropagation and up to 250% higher maximum pressure compared to the ray-tracing method, particularly in strongly distorted cases.

On the imaging depth limit of photoacoustic tomography in the visible and first near-infrared windows

It is well known that photoacoustic tomography (PAT) can circumvent the photon scattering problem in optical imaging and achieve high-contrast and high-resolution imaging at centimeter depths. However, after two decades of development, the long-standing question of the imaging depth limit of PAT in biological tissues remains unclear. Here we propose a numerical framework for evaluating the imaging depth limit of PAT in the visible and the first near-infrared windows. The established framework simulates the physical process of PAT and consists of seven modules, including tissue modelling, photon transportation, photon to ultrasound conversion, sound field propagation, signal reception, image reconstruction, and imaging depth evaluation. The framework can simulate the imaging depth limits in general tissues, such as the human breast, the human abdomen-liver tissues, and the rodent whole body and provide accurate evaluation results. The study elucidates the fundamental imaging depth limit of PAT in biological tissues and can provide useful guidance for practical experiments.

Evolution of acoustic nonlinearity in outdoor blast propagation from firearms: On the persistence of nonlinear behavior

Acoustic events exceeding a certain threshold of intensity cannot benefit from a linearization of the governing wave equation, posing an additional burden on the numerical modelling. Weak shock theory associates nonlinearity with the generation of high frequency harmonics that compensate for atmospheric attenuation. Overlooking the persistence of this phenomenon at large distances can lead to mispredictions in gun detection procedures, noise abatement protocols, and auditory risk assessment. The state-of-the-art mostly addresses aircraft jet noise, a stationary and largely random type of signal. The extension of such conclusions to muzzle blasts requires caution in considering their peculiar impulsive and broadband nature. A methodology based on the time and frequency analysis of an experimental dataset of eight calibres intends to find quantitative metrics linked to acoustic nonlinearity in outdoor muzzle blast propagation. Propagating three waveforms (SCAR-L 7.62 mm, Browning 9 mm, and Howitzer 105 mm) up to 300 [m] with the in-house numerical solver based on the nonlinear progressive wave equation, demonstrates that the propagation does not downgrade to truly linear.

An all-ultrasound cranial imaging method to establish the relationship between cranial FUS incidence angle and transcranial attenuation in non-human primates in 3D

Focused ultrasound (FUS) is a non-invasive and non-ionizing technique which deploys ultrasound waves to induce bio-effects. When paired with acoustically active particles such as microbubbles (MBs), it can open the blood brain barrier (BBB) to facilitate drug delivery otherwise inhibited due to the presence of BBB. One of the parameters that affects the FUS beam propagation is the beam incidence angle on the skull. Prior work by our group has shown that, as incidence angles deviate from 90°, FUS focal pressures attenuate and result in a smaller BBB opening volume. The incidence angles calculated in our prior studies were in 2D and used skull information from CT. The study presented herein develops methods to calculate incidence angle in 3D in non-human primate (NHP) skull fragments using harmonic ultrasound imaging without using ionizing radiation. Our results show that ultrasound harmonic imaging is capable of accurately depicting features such as sutures and eye-sockets of the skull. Furthermore, we were able to reproduce previously reported relationships between the incidence angle and FUS beam attenuation. We also show feasibility of performing ultrasound harmonic imaging in in-vivo non-human primates. The all-ultrasound method presented herein combined with our neuronavigation system stands to increase more widespread adoption of FUS and render it accessible by eliminating the need for CT cranial mapping.

Investigation of a method to estimate the average speed of sound using phase variances of element signals for ultrasound compound imaging

PURPOSE

In the receive beamforming of an ultrasonography system, a B-mode image is reconstructed by assuming an average speed of sound (SoS) as a constant value. In our previous studies, we proposed a method for estimating the average SoS based on the coherence factor (CF) and the reciprocal of phase variances of element signals in delay-and-sum (DAS) beamforming. In this paper, we investigate the accuracy of estimation of the average SoS for compound imaging.

METHODS

For this purpose, two numerical simulations were performed with k-Wave software. Also, the estimation methods based on the CF and the reciprocal were applied to in vivo data from the common carotid artery, and B-mode images were reconstructed using the estimated average SoS.

RESULTS

In the first numerical simulation using an inhomogeneous phantom, the relationship between the accuracy and the transmission angles for the estimation was investigated, and the root mean squared errors (RMSEs) of estimates obtained based on the CF and the reciprocal of the phase variance were 1.25 ± 0.09, and 0.765 ± 0.17% at the transmission sequence of steering angles of (- 10°, - 5°, 0°, 5°, 10°), respectively. In the second numerical simulation using a cyst phantom, lateral resolutions were improved by reconstructing the image using the estimates obtained using the proposed strategy (reciprocal). By the proposed strategy, improvement of the continuity of the lumen-intima interface in the lateral direction was observed in the in vivo experiment.

CONCLUSION

Consequently, the results indicated that the proposed strategy was beneficial for estimation of the average SoS and image reconstruction.

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

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

Hermite-scan imaging for differentiating glioblastoma from normal brain: Simulations and ex vivo studies for applications in intra-operative tumor identificationa)

Hermite-scan (H-scan) imaging is a tissue characterization technique based on the analysis of raw ultrasound radio frequency (RF) echoes. It matches the RF echoes to Gaussian-weighted Hermite polynomials of various orders to extract information related to scatterer diameter. It provides a color map of large and small scatterers in the red and blue H-scan image channels, respectively. H-scan has been previously reported for characterizing breast, pancreatic, and thyroid tumors. The present work evaluated H-scan imaging to differentiate glioblastoma tumors from normal brain tissue ex vivo. First, we conducted 2-D numerical simulations using the k-wave toolbox to assess the performance of parameters derived from H-scan images of acoustic scatterers (15-150 μm diameters) and concentrations (0.2%-1% w/v). We found that the parameter intensity-weighted percentage of red (IWPR) was sensitive to changes in scatterer diameters independent of concentration. Next, we assessed the feasibility of using the IWPR parameter for differentiating glioblastoma and normal brain tissues (n = 11 samples per group). The IWPR parameter estimates for normal tissue (44.1% ± 1.4%) were significantly different (p < 0.0001) from those for glioblastoma (36.2% ± 0.65%). These findings advance the development of H-scan imaging for potential use in differentiating glioblastoma tumors from normal brain tissue during resection surgery.

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.

A Forward Model Incorporating Elevation-Focused Transducer Properties for 3-D Full-Waveform Inversion in Ultrasound Computed Tomography

Ultrasound computed tomography (USCT) is an emerging medical imaging modality that holds great promise for improving human health. Full-waveform inversion (FWI)-based image reconstruction methods account for the relevant wave physics to produce high spatial resolution images of the acoustic properties of the breast tissues. A practical USCT design employs a circular ring-array comprised of elevation-focused ultrasonic transducers, and volumetric imaging is achieved by translating the ring-array orthogonally to the imaging plane. In commonly deployed slice-by-slice (SBS) reconstruction approaches, the 3-D volume is reconstructed by stacking together 2-D images reconstructed for each position of the ring-array. A limitation of the SBS reconstruction approach is that it does not account for 3-D wave propagation physics and the focusing properties of the transducers, which can result in significant image artifacts and inaccuracies. To perform 3-D image reconstruction when elevation-focused transducers are employed, a numerical description of the focusing properties of the transducers should be included in the forward model. To address this, a 3-D computational model of an elevation-focused transducer is developed to enable 3-D FWI-based reconstruction methods to be deployed in ring-array-based USCT. The focusing is achieved by applying a spatially varying temporal delay to the ultrasound pulse (emitter mode) and recorded signal (receiver mode). The proposed numerical transducer model is quantitatively validated and employed in computer simulation studies that demonstrate its use in image reconstruction for ring-array USCT.

Real-Time Passive Acoustic Mapping With Enhanced Spatial Resolution in Neuronavigation-Guided Focused Ultrasound for Blood-Brain Barrier Opening

OBJECTIVE

Passive acoustic mapping (PAM) provides the spatial information of acoustic energy emitted from microbubbles during focused ultrasound (FUS), which can be used for safety and efficacy monitoring of blood-brain barrier (BBB) opening. In our previous work with a neuronavigation-guided FUS system, only part of the cavitation signal could be monitored in real time due to the computational burden although full-burst analysis is required to detect transient and stochastic cavitation activity. In addition, the spatial resolution of PAM can be limited for a small-aperture receiving array transducer. For full-burst real-time PAM with enhanced resolution, we developed a parallel processing scheme for coherence-factor-based PAM (CF-PAM) and implemented it onto the neuronavigation-guided FUS system using a co-axial phased-array imaging transducer.

METHODS

Simulation and in-vitro human skull studies were conducted for the performance evaluation of the proposed method in terms of spatial resolution and processing speed. We also carried out real-time cavitation mapping during BBB opening in non-human primates (NHPs).

RESULTS

CF-PAM with the proposed processing scheme provided better resolution than that of traditional time-exposure-acoustics PAM with a higher processing speed than that of eigenspace-based robust Capon beamformer, which facilitated the full-burst PAM with the integration time of 10 ms at a rate of 2 Hz. In vivo feasibility of PAM with the co-axial imaging transducer was also demonstrated in two NHPs, showing the advantages of using real-time B-mode and full-burst PAM for accurate targeting and safe treatment monitoring.

SIGNIFICANCE

This full-burst PAM with enhanced resolution will facilitate the clinical translation of online cavitation monitoring for safe and efficient BBB opening.

Acoustoelectric Time-Reversal for Ultrasound Phase-Aberration Correction

Acoustoelectric imaging (AEI) is a technique that combines ultrasound (US) with radio frequency recording to detect and map local current source densities. This study demonstrates a new method called acoustoelectric time reversal (AETR), which uses AEI of a small current source to correct for phase aberrations through a skull or other US-aberrating layers with applications to brain imaging and therapy. Simulations conducted at three different US frequencies (0.5, 1.5, and 2.5 MHz) were performed through media layered with different sound speeds and geometries to induce aberrations of the US beam. Time delays of the acoustoelectric (AE) signal from a monopole within the medium were calculated for each element to enable corrections using AETR. Uncorrected aberrated beam profiles were compared with those after applying AETR corrections, which demonstrated a strong recovery (29%-100%) of lateral resolution and increases in focal pressure up to 283%. To further demonstrate the practical feasibility of AETR, we further conducted bench-top experiments using a 2.5 MHz linear US array to perform AETR through 3-D-printed aberrating objects. These experiments restored lost lateral restoration up to 100% for the different aberrators and increased focal pressure up to 230% after applying AETR corrections. Cumulatively, these results highlight AETR as a powerful tool for correcting focal aberrations in the presence of a local current source with applications to AEI, US imaging, neuromodulation, and therapy.

Aberration correction in diagnostic ultrasound: A review of the prior field and current directions

Medical ultrasound images are reconstructed with simplifying assumptions on wave propagation, with one of the most prominent assumptions being that the imaging medium is composed of a constant sound speed. When the assumption of a constant sound speed are violated, which is true in most in vivoor clinical imaging scenarios, distortion of the transmitted and received ultrasound wavefronts appear and degrade the image quality. This distortion is known as aberration, and the techniques used to correct for the distortion are known as aberration correction techniques. Several models have been proposed to understand and correct for aberration. In this review paper, aberration and aberration correction are explored from the early models and correction techniques, including the near-field phase screen model and its associated correction techniques such as nearest-neighbor cross-correlation, to more recent models and correction techniques that incorporate spatially varying aberration and diffractive effects, such as models and techniques that rely on the estimation of the sound speed distribution in the imaging medium. In addition to historical models, future directions of ultrasound aberration correction are proposed.

Frequency-domain full-waveform inversion-based musculoskeletal ultrasound computed tomography

Recently, full-waveform inversion (FWI) has become a promising tool for ultrasound computed tomography (USCT). However, as a computationally intensive technique, FWI suffers from computational burden, especially in conventional time-domain full-waveform inversion (TDFWI). On the contrary, frequency-domain full-waveform inversion (FDFWI) provides a relatively high computational efficiency as the propagation of discrete frequencies is much cheaper than full time-domain modeling. FDFWI has already been applied in soft tissue imaging, such as breast, but for the musculoskeletal model with high impedance contrast between hard and soft tissues, there is still a lack of an effective source estimation method. In this paper, a water-referenced data calibration method is proposed to address the source estimation challenge in the presence of bones, which achieves consistency between the measured and simulated data before the FDFWI procedure. To avoid the cycle-skipping local minimum effect and facilitate the algorithm convergence, a starting frequency criterion for musculoskeletal FDFWI is further proposed. The feasibility of the proposed method is demonstrated by numerical studies on retrieving the anatomies of the leg models and different musculoskeletal lesions. The study extends the advanced FDFWI method to the musculoskeletal system and provides an alternative solution for musculoskeletal USCT imaging with high computational efficiency.

Numerical modeling of Mach-stem formation in high-amplitude time-reversal focusing

In acoustics, time-reversal processing is commonly used to exploit multiple scatterings in reverberant environments to focus sound to a specific location. Recently, the nonlinear characteristics of time-reversal focusing at amplitudes as high as 200 dB have been reported [Patchett and Anderson, J. Acoust. Soc. Am. 151(6), 3603-3614 (2022)]. These studies were experimental in nature and suggested that converging waves nonlinearly interact in the focusing of waves, leading to nonlinear amplification. This study investigates the nonlinear interactions and subsequent characteristics from a model-based approach. Utilizing both finite difference and finite-element models, it is shown that nonlinear interactions between high-amplitude waves lead to free-space Mach-wave coalescence of the converging waves. The number of waves used in both models represents a small piece of the full aperture of converging waves experimentally. Limiting the number of waves limits the number of Mach-stem formations and reduces the nonlinear growth of the focus amplitudes when compared to experiment. However, limiting the number of waves allows the identification of individual Mach waves. Mach wave coalescence leading to Mach-stem formation appears to be the mechanism behind nonlinear amplification of peak focus amplitudes observed in high-amplitude time-reversal focusing.

Interframe Echo Intensity Variation of Subregions and Whole Plaque in Two-Dimensional Carotid Ultrasonography: Simulations and In Vivo Observations

OBJECTIVES

The risk of cardiovascular disease is associated with the echo intensity of carotid plaques in ultrasound images and their cardiac cycle-induced intensity variations. In this study, we aimed to 1) explore the underlying origin of echo intensity variations by using simulations and 2) evaluate the association between the two-dimensional (2D) spatial distribution of these echo intensity variations and plaque vulnerability.

METHODS

First, we analyzed how out-of-plane motion and compression of simulated scattering spheres of different sizes affect the ultrasound echo intensity. Next, we propose a method to analyze the features of the 2D spatial distribution of interframe plaque echo intensity in carotid ultrasound image sequences and explore their associations with plaque vulnerability in experimental data.

RESULTS

The simulations showed that the magnitude of echo intensity changes was similar for both the out-of-plane motion and compression, but for scattering objects smaller than 1 mm radius, the out-of-plane motion dominated. In experimental data, maps of the 2D spatial distribution of the echo intensity variations had a low correlation with standard B-mode echo intensity distribution, indicating complementary information on plaque tissue composition. In addition, we found the existence of ∼1 mm diameter subregions with pronounced echo intensity variations associated with plaque vulnerability.

CONCLUSIONS

The results indicate that out-of-plane motion contributes to intra-plaque regions of high echo intensity variation. The 2D echo intensity variation maps may provide complementary information for assessing plaque composition and vulnerability. Further studies are needed to verify this method's role in identifying vulnerable plaques and predicting cardiovascular disease risk.

An aberration correction approach for single and dual aperture ultrasound imaging of the abdomen

Abdominal ultrasound image quality is hampered by phase aberration, that is mainly caused by the large speed-of-sound (SoS) differences between fat and muscle tissue in the abdominal wall. The mismatch between the assumed and actual SoS distribution introduces general blurring of the ultrasound images, and acoustic refraction can lead to geometric distortion of the imaged features. Large aperture imaging or dual-transducer imaging can improve abdominal imaging at deep locations by providing increased contrast and resolution. However, aberration effects for large aperture imaging can be even more severe, which limits its full potential. In this study, a model-based aberration correction method for arbitrary acquisition schemes is introduced for delay-and-sum (DAS) beamforming and its performance was analyzed for both single- and dual-transducer ultrasound imaging. The method employs aberration corrected wavefront arrival times, using manually assigned local SoS values. Two wavefront models were compared. The first model is based on a straight ray approximation, and the second model on the Eikonal equation, which is solved by a multi-stencils fast marching method. Their accuracy for abdominal imaging was evaluated in acoustic simulations and phantom experiments involving tissue-mimicking and porcine material with large SoS contrast (∼100 m/s). The lateral resolution was improved by up to 90% in simulations and up to 65% in experiments compared to standard DAS, in which the use of Eikonal beamforming generally outperformed straight ray beamforming. Moreover, geometric distortions were mitigated in multi-aperture imaging, leading to a reduction in position error of around 80%. A study on the sensitivity of the aberration correction to shape and SoS of aberrating layers was performed, showing that even with imperfect segmentations or SoS values, aberration correction still outperforms standard DAS.

The effect of source backing materials and excitation pulse durations on laser-generated ultrasound waveforms

In this article, it is shown experimentally that a planar laser-generated ultrasound source with a hard reflective backing will generate higher acoustic pressures than a comparable source with an acoustically matched backing when the stress confinement condition is not met. Furthermore, while the source with an acoustically matched backing will have a broader bandwidth when the laser pulse is short enough to ensure stress confinement, the bandwidths of both source types will converge as the laser pulse duration increases beyond stress confinement. The explanation of the results is supported by numerical simulations.

Physics-informed neural networks for transcranial ultrasound wave propagation

Transcranial ultrasound imaging has been playing an increasingly important role in the non-invasive treatment of brain disorders. However, the conventional mesh-based numerical wave solvers, which are an integral part of imaging algorithms, suffer from limitations such as high computational cost and discretization error in predicting the wavefield passing through the skull. In this paper, we explore the use of physics-informed neural networks (PINNs) for predicting the transcranial ultrasound wave propagation. The wave equation, two sets of time snapshots data and a boundary condition (BC) are embedded as physical constraints in the loss function during training. The proposed approach has been validated by solving the two-dimensional (2D) acoustic wave equation under three increasingly complex spatially varying velocity models. Our cases demonstrate that due to the meshless nature of PINNs, they can be flexibly applied to different wave equations and types of BCs. By adding physics constraints to the loss function, PINNs can predict wavefields far outside the training data, providing ideas for improving the generalization capability of existing deep learning methods. The proposed approach offers exciting perspectives because of the powerful framework and simple implementation. We conclude with a summary of the strengths, limitations and further research directions of this work.

Multifocal skull-compensated transcranial focused ultrasound system for neuromodulation applications based on acoustic holography

Transcranial focused ultrasound stimulation is a promising therapeutic modality for human brain disorders because of its noninvasiveness, long penetration depth, and versatile spatial control capability through beamforming and beam steering. However, the skull presents a major hurdle for successful applications of ultrasound stimulation. Specifically, skull-induced focal aberration limits the capability for accurate and versatile targeting of brain subregions. In addition, there lacks a fully functional preclinical neuromodulation system suitable to conduct behavioral studies. Here, we report a miniature ultrasound system for neuromodulation applications that is capable of highly accurate multiregion targeting based on acoustic holography. Our work includes the design and implementation of an acoustic lens for targeting brain regions with compensation for skull aberration through time-reversal recording and a phase conjugation mirror. Moreover, we utilize MEMS and 3D-printing technology to implement a 0.75-g lightweight neuromodulation system and present in vivo characterization of the packaged system in freely moving mice. This preclinical system is capable of accurately targeting the desired individual or multitude of brain regions, which will enable versatile and explorative behavior studies using ultrasound neuromodulation to facilitate widespread clinical adoption.

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.

Acoustic Radiation Forces at the Crossroads of Ultrasound Exposimetry, HIFU, and Elastography

The applications of the acoustic radiation force (ARF) continue to multiply and extend from elastography into high-intensity focused ultrasound (HIFU), diagnostic imaging, lithotripsy, sonochemistry, levitation, and microsonics yet fundamental principles remain shrouded in mystery. A well-known and popular equation often used for calculating ARF in elastography is in conflict with the equation commonly employed for calculating ARF for determining acoustic power in radiation force balances (RFBs). Controversies have sparked debate for over a century concerning the physical mechanisms underlying ARFs. For over four decades, the science of ultrasound exposimetry has steadily improved and has provided clues in terms of accumulated data about the characteristics of transmitted ultrasound fields. Concurrently, the availability and capability of predicting these fields have improved significantly. The author draws on these sources to re-examine the physical principles behind ARFs. Conflicts are shown to stem from idealized configurations and simplistic assumptions. By more fully accounting for the pulse shape and spectrum, the effect of frequency power law attenuation, diffraction, and nonlinearity, more accurate equations are developed for ARFs for practical applications which are more consistent with exposimetry observations. Simulations compare well to corrected 1.5 MHz RFB data. While some questions await resolution, the approach presented here settles several longstanding conflicts and provides a new broadband framework for future ARF work.

Fast iterative reconstruction for photoacoustic tomography using learned physical model: Theoretical validation

Iterative reconstruction has demonstrated superior performance in medical imaging under compressed, sparse, and limited-view sensing scenarios. However, iterative reconstruction algorithms are slow to converge and rely heavily on hand-crafted parameters to achieve good performance. Many iterations are usually required to reconstruct a high-quality image, which is computationally expensive due to repeated evaluations of the physical model. While learned iterative reconstruction approaches such as model-based learning (MBLr) can reduce the number of iterations through convolutional neural networks, it still requires repeated evaluations of the physical models at each iteration. Therefore, the goal of this study is to develop a Fast Iterative Reconstruction (FIRe) algorithm that incorporates a learned physical model into the learned iterative reconstruction scheme to further reduce the reconstruction time while maintaining robust reconstruction performance. We also propose an efficient training scheme for FIRe, which releases the enormous memory footprint required by learned iterative reconstruction methods through the concept of recursive training. The results of our proposed method demonstrate comparable reconstruction performance to learned iterative reconstruction methods with a 9x reduction in computation time and a 620x reduction in computation time compared to variational reconstruction.

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.

Nonlinear Harmonic Distortion of Complementary Golay Codes

Recent advances in electronics miniaturization have led to the development of low-power, low-cost, point-of-care ultrasound scanners. Low-cost systems employing simple bi-level pulse generation devices need only utilize binary phase modulated coded excitations to significantly improve sensitivity; however the performance of complementary codes in the presence of nonlinear harmonic distortion has not been thoroughly investigated. Through simulation, it was found that nonlinear propagation media with little attenuative properties can significantly deteriorate the Peak Sidelobe Level (PSL) performance of complementary Golay coded pulse compression, resulting in PSL levels of -62 dB using nonlinear acoustics theory contrasted with -198 dB in the linear case. Simulations of 96 complementary pairs revealed that some pairs are more robust to sidelobe degradation from nonlinear harmonic distortion than others, up to a maximum PSL difference of 17 dB between the best and worst performing codes. It is recommended that users consider the effects of nonlinear harmonic distortion when implementing binary phase modulated complementary Golay coded excitations.

Pulse-Echo Technique to Compensate for Laminate Membrane Transmission Loss in Phantom-Based Ultrasonic Attenuation Coefficient Measurements

OBJECTIVES

Accurately measuring the attenuation coefficient (AC) of reference phantoms is critical in clinical applications of quantitative ultrasound. Phantom AC measurement requires proper compensation of membrane transmission loss. Conventional methods require separate membrane samples to obtain membrane transmission loss. Unfortunately, separate membrane samples are often unavailable. A pulse-echo approach is proposed herein to compensate for membrane transmission loss without requiring separate membrane samples.

METHODS

The proposed method consists of the following steps. First, the insertion loss, caused by phantom attenuation and membrane transmission loss, is measured. Second, the membrane reflection coefficient is measured. Third, the unknown acoustic parameters of the membrane and phantom material are estimated by fitting theoretical reflection coefficient to the measured one. Finally, the fitted parameters are used to estimate membrane transmission loss and phantom AC. The proposed method was validated through k-Wave simulations and phantom experiments. Experimental AC measurements were repeated on 5 distinct phantoms by 2 operators to assess the repeatability and reproducibility of the proposed method. Five transducers were used to cover a broad bandwidth (0.7-16 MHz).

RESULTS

The acquired AC in the simulations had a maximum error of 0.06 dB/cm-MHz for simulated phantom AC values ranging from 0.5 to 1 dB/cm-MHz. The acquired AC in the experiments had a maximum error of 0.045 dB/cm-MHz for phantom AC values ranging from 0.28 to 1.48 dB/cm-MHz. Good repeatability and cross-operator reproducibility were observed with a mean coefficient of variation below 0.054.

CONCLUSION

The proposed method simplifies phantom AC measurement while providing satisfactory accuracy and precision.

Experiments and simulations demonstrating the rapid ultrasonic rewarming of frozen tissue cryovials

The development of methods to safely rewarm large cryopreserved biological samples remains a barrier to the widespread adoption of cryopreservation. Here, experiments and simulations were performed to demonstrate that ultrasound can increase rewarming rates relative to thermal conduction alone. An ultrasonic rewarming setup based on a custom 444 kHz tubular piezoelectric transducer was designed, characterized, and tested with 2 ml cryovials filled with frozen ground beef. Rewarming rates were characterized in the -20 °C to 5 °C range. Thermal conduction-based rewarming was compared to thermal conduction plus ultrasonic rewarming, demonstrating a tenfold increase in rewarming rate when ultrasound was applied. The maximum recorded rewarming rate with ultrasound was 57° C/min, approximately 2.5 times faster than with thermal conduction alone. Coupled acoustic and thermal simulations were developed and showed good agreement with the heating rates demonstrated experimentally and were also used to demonstrate spatial heating distributions with small (<3° C) temperature differentials throughout the sample when the sample was below 0° C. The experiments and simulations demonstrate the potential for ultrasonic cryovial rewarming with a possible application to large volume rewarming, as faster rewarming rates may improve the viability of cryopreserved tissues and reduce the time needed for cells to regain normal function.

3-D Protoacoustic Imaging Through a Planar Ultrasound Array: A Simulation Workflow

Bragg peak range uncertainties are a persistent constraint in proton therapy. Pulsed proton beams generate protoacoustic emissions proportional to absorbed proton energy, thereby encoding dosimetry information in a detectable acoustic wave. Here, we seek to derive and model 3D protoacoustic imaging with an ultrasound array and examine the frequency characteristics of protoacoustic emissions. A formalism is presented through which protoacoustic signals can be characterized considering transducer bandwidth as well as pulse duration of the incident beam. We have also collected an experimental proton beam intensity signal from a Mevion S250 clinical machine to analyze our formalism. We also show that proton-acoustic image reconstruction is possible even when the noise amplitude is larger than the signal amplitude on individual transducers. We find that a 4μ s Gaussian proton pulse can generate a signal in the range of MHz as long as the spatial heating function has sufficiently high temperature gradients.

Evaluation of the effects of focused ultrasound stimulation on the central nervous system through a multiscale simulation approach

The lack of sensory feedback represents one of the main drawbacks of commercial upper limb prosthesis. Transcranial Focused Ultrasound Stimulation (tFUS) seems to be a valid non-invasive technique for restoring sensory feedback allowing to deliver acoustic energy to cortical sensory areas with high spatial resolution and depth penetration. This paper aims at studying in simulation the use of tFUS on cortical sensory areas to evaluate its effects in terms of latency ad firing rate of the cells response, for understanding if these parameters influence the safety and the efficacy of the stimulation. In this paper, in order to study the propagation of the ultrasound wave from the transducer to the cortical cells, a multiscale approach was implemented by building a macroscopic model, which estimates the pressure profile in a simplified 2D human head geometry, and coupling it with the SONIC microscale model, that describes the electrical behaviour of a cortical neuron. The influence of the stimulation parameters and of the skull thickness on the latency and the firing rate are evaluated and the obtained behaviour is linked to the sensory response obtained on human subjects. Results have shown that slight changes in the transducer position should not affect the efficacy of the stimulation; however, high skull thickness leads to lower cells activation. These results will be useful for evaluating safety and effectiveness of tFUS for sensory feedback in closed-loop prosthetic systems.

Effects of acoustic nonlinearity on communication performance in soft tissues

Acoustic communication has been gaining traction as an alternative communication method in nontraditional media, such as underwater or through tissue. Acoustic propagation is known to be a nonlinear phenomenon; nonlinear propagation of acoustic waves in soft tissues at biomedical frequencies and intensities has been widely demonstrated. However, the effects of acoustic nonlinearity on communication performance in biological tissues have not yet been examined. In this work, nonlinear propagation of a communication signal in soft tissues is analyzed. The relationship between communication parameters (signal amplitude, bandwidth, and center frequency) and nonlinear distortion of the communication signal propagating in soft tissues with different acoustic properties is investigated. Simulated experiments revealed that, unlike linear channels, bit error rates increase as signal amplitude and bandwidth increase. Linear and decision feedback equalizers fail to address the increased error rates. When tissue properties and transmission parameters can be estimated, receivers based on maximum likelihood sequence estimation approach the performance of an ideal receiver in an ideal additive white Gaussian noise channel.

Classical and Learned MR to Pseudo-CT Mappings for Accurate Transcranial Ultrasound Simulation

Model-based treatment planning for transcranial ultrasound therapy typically involves mapping the acoustic properties of the skull from an X-ray computed tomography (CT) image of the head. Here, three methods for generating pseudo-CT (pCT) images from magnetic resonance (MR) images were compared as an alternative to CT. A convolutional neural network (U-Net) was trained on paired MR-CT images to generate pCT T images from either T1-weighted or zero-echo time (ZTE) MR images (denoted tCT and zCT, respectively). A direct mapping from ZTE to pCT was also implemented (denoted cCT). When comparing the pCT and ground-truth CT images for the test set, the mean absolute error was 133, 83, and 145 Hounsfield units (HU) across the whole head, and 398, 222, and 336 HU within the skull for the tCT, zCT, and cCT images, respectively. Ultrasound simulations were also performed using the generated pCT images and compared to simulations based on CT. An annular array transducer was used targeting the visual or motor cortex. The mean differences in the simulated focal pressure, focal position, and focal volume were 9.9%, 1.5 mm, and 15.1% for simulations based on the tCT images; 5.7%, 0.6 mm, and 5.7% for the zCT; and 6.7%, 0.9 mm, and 12.1% for the cCT. The improved results for images mapped from ZTE highlight the advantage of using imaging sequences, which improves the contrast of the skull bone. Overall, these results demonstrate that acoustic simulations based on MR images can give comparable accuracy to those based on CT.

Effects of Beamforming Techniques on Quality of Ultrasound Computed Tomography Images

BACKGROUND

In body tissues, tumors generally have different speeds of sound (SOS) than normal tissues. In this respect, ultrasound computed tomography (UCT) can generate a cross-sectional SOS map as an innovative ultrasound imaging method. This technique can produce images with a resolution of millimeters and a high signal-to-noise ratio.

OBJECTIVE

This study aimed to improve UCT image quality without increasing breast cancer screening and diagnosis time.

MATERIAL AND METHODS

In this analytical study, a ring-shaped UCT breast imaging system was simulated using the K-wave toolbox of MATLAB. The system has a 20 cm diameter and 256 ultrasonic piezoelectrics placed in the ring's circumference. Different beamforming techniques imaged two designed phantoms (i.e., resolution and contrast), and the resolution and contrast to noise ratio (CNR) were calculated.

RESULTS

The results of resolution phantom imaging without any beamforming showed that only bars with the value of 0.125 and 0.167 lp/mm were distinguishable, and the 0.1 bars were not recognizable in the imaging. In addition, increasing the number of transmitters led to no noticeable change in resolution for 0.125 and 0.167 lp/mm bars. In all beamforming techniques for imaging the contrast phantom, the CNR parameter up to an object with a diameter of 8 mm increases with increasing diameter without any change.

CONCLUSION

The beamforming technique using three simultaneous transmitters improved the resolution by about 1 mm compared to the normal strategy. In addition to high-contrast images, beamforming with 9 simultaneous transmitters led to a preferable technique.

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

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

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

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

Effects of phase aberration on transabdominal focusing for a large aperture, lowf-number histotripsy transducer

Objective. Soft tissue phase aberration may be particularly severe for histotripsy due to large aperture and lowf-number transducer geometries. This study investigated how phase aberration from human abdominal tissue affects focusing of a large, strongly curved histotripsy transducer.Approach.A computational model (k-Wave) was experimentally validated withex vivoporcine abdominal tissue and used to simulate focusing a histotripsy transducer (radius: 14.2 cm,f-number: 0.62, central frequencyfc: 750 kHz) through the human abdomen. Abdominal computed tomography images from 10 human subjects were segmented to create three-dimensional acoustic property maps. Simulations were performed focusing at 3 target locations in the liver of each subject with ideal phase correction, without phase correction, and after separately matching the sound speed of water and fat to non-fat soft tissue.Main results.Experimental validation in porcine abdominal tissue showed that simulated and measured arrival time differences agreed well (average error, ∼0.10 acoustic cycles atfc). In simulations with human tissue, aberration created arrival time differences of 0.65μs (∼0.5 cycles) at the target and shifted the focus from the target by 6.8 mm (6.4 mm pre-focally along depth direction), on average. Ideal phase correction increased maximum pressure amplitude by 95%, on average. Matching the sound speed of water and fat to non-fat soft tissue decreased the average pre-focal shift by 3.6 and 0.5 mm and increased pressure amplitude by 2% and 69%, respectively.Significance.Soft tissue phase aberration of large aperture, lowf-number histotripsy transducers is substantial despite low therapeutic frequencies. Phase correction could potentially recover substantial pressure amplitude for transabdominal histotripsy. Additionally, different heterogeneity sources distinctly affect focusing quality. The water path strongly affects the focal shift, while irregular tissue boundaries (e.g. fat) dominate pressure loss.

Time-of-flight completion in ultrasound computed tomography based on the singular value threshold algorithm

Ultrasound computed tomography (USCT) has been developed for breast tumor screening. The sound-speed modal of USCT can provide quantitative sound-speed values to help tumor diagnosis. Time-of-flight (TOF) is the critical input in sound-speed reconstruction. However, we found that the missing data problem in the detected TOF causes artifacts on the reconstructed sound-speed images, which may affect the tumor identification. In this study, to address the missing TOF data problem, we first adopted the singular value threshold (SVT) algorithm to complete the TOF matrix. The threshold value in SVT is difficult to determine, so we proposed a selection strategy, that is, to enumerate the threshold values as the multiples of the maximum singular value of the incomplete matrix and then evaluate the image quality to select the proper threshold value. In the numerical breast phantom experiment, the artifacts are eliminated, and the accuracy is higher than the accuracy of the compared methods. In the in vivo experiment, we reconstructed the sound-speed image of the breast of a volunteer with invasive breast cancer, and the SVT algorithm improved the image sharpness. The completion of DTOF based on SVT gives better accuracy than the compared methods, but too large a threshold value decreases the accuracy. In the future, the selection method of the threshold value needs further research, and more USCT cases should be included in the experiments.

Multiscale approach for tFUS neurocomputational modelling

Among the non-invasive methods employed for brain stimulation, trans cranial Focused Ultrasound Stimulation (tFUS) is the technique with the best penetration into the tissues and spatial resolution. The development of computational models of US propagation in brain tissue can be useful for estimating the behaviour of neural cells subjected to mechanical stimulus due to US. This paper aims at studying the neural cell response of a cortical Regular Spiking point neuron model, for different values of stimulus Duty Cycle (DC). The main goal is to use a multiscale approach to couple the results obtained from a macroscale simulation on wave propagation in tissue, with neuron model described by Hodgkin-Huxley equations to study latency and firing rate of the RS model. The obtained results showed that latency and firing rate have slight variations along the propagation direction of the US beam, in the focal region under the skull model, for different stimulus DC.

Development of a Numerical Model of High-Intensity Focused Ultrasound Treatment in Mobile and Elastic Organs: Application to a Beating Heart

High-intensity focused ultrasound (HIFU) is a promising method used to treat cardiac arrhythmias, as it can induce lesions at a distance throughout myocardium thickness. Numerical modeling is commonly used for ultrasound probe development and optimization of HIFU treatment strategies. This study was aimed at describing a numerical method to simulate HIFU thermal ablation in elastic and mobile heart models. The ultrasound pressure field is computed on a 3-D orthonormal grid using the Rayleigh integral method, and the attenuation is calculated step by step between cells. The temperature distribution is obtained by resolution of the bioheat transfer equation on a 3-D non-orthogonally structured curvilinear grid using the finite-volume method. The simulation method is applied on two regions of the heart (atrioventricular node and ventricular apex) to compare the thermal effects of HIFU ablation depending on deformation, motion type and amplitude. The atrioventricular node requires longer sonication than the ventricular apex to reach the same lesion volume. Motion considerably influences treatment duration, lesion shape and distribution in cardiac HIFU treatment. These results emphasize the importance of considering local motion and deformation in numerical studies to define efficient and accurate treatment strategies.

SIMUS: An open-source simulator for medical ultrasound imaging. Part II: Comparison with four simulators

BACKGROUND AND OBJECTIVE

Computational ultrasound imaging has become a well-established methodology in the ultrasound community. In the accompanying paper (part I), we described a new ultrasound simulator (SIMUS) for MATLAB, which belongs to the Matlab UltraSound Toolbox (MUST). SIMUS can generate pressure fields and radiofrequency RF signals for simulations in medical ultrasound imaging. It works in a harmonic domain and uses far-field and paraxial linear equations.

METHODS

In this article (part II), we illustrate how SIMUS compares with other ultrasound simulators (Field II, k-Wave, FOCUS, and Verasonics) for a homogeneous medium. We designed different transmit sequences (focused, planar, and diverging wavefronts) and calculated the corresponding 2-D and 3-D (with elevation focusing) RMS pressure fields.

RESULTS

SIMUS produced pressure fields similar to those of Field II, FOCUS, and k-Wave. The acoustic fields provided by the Verasonics simulator were significantly different from those of SIMUS and k-Wave, although the overall appearance remained consistent.

CONCLUSION

Our simulations tend to demonstrate that SIMUS is reliable and can be used for realistic medical ultrasound simulations.

Ultrasonic backscatter coefficient estimation in nonlinear regime using an in situ calibration target

Tissue characterization based on the backscatter coefficient (BSC) can be degraded by acoustic nonlinearity. Often, this degradation is due to the method used for obtaining a reference spectrum, i.e., using a planar reference in water compared to a reference phantom approach resulted in more degradation. We hypothesize that an in situ calibration approach can improve BSC estimates in the nonlinear regime compared to using the reference phantom approach. The in situ calibration target provides a reference within the medium being interrogated and, therefore, nonlinear effects would already be contained in the in situ reference signal. Simulations and experiments in phantoms and in vivo were performed. A 2 mm diameter titanium bead was embedded in the interrogated media. An L9-4/38 probe (BK Ultrasound, Peabody, MA) and an analysis bandwidth from 4.5 to 7.4 MHz were used in experiments. Radiofrequency data from the sample, bead, and reference phantoms were acquired at a quasi-linear baseline power level and at further increments of output power. Better agreement between the BSC obtained at low power compared to high power was observed for the in situ calibration compared to the reference phantom approach.

A reduced aperture allows for transcranial focus localization at lower pressure

Localizing the focus during transcranial focused ultrasound procedures is important to ensure accurate targeting of specific brain regions and interpretation of results. Magnetic resonance acoustic radiation force imaging uses the displacement induced by the ultrasound focus in the brain to localize the beam, but the high pressure required to displace brain tissue may cause damage or confounds during subsequent neuromodulatory experiments. Here, reduced apertures were applied to a phased array transducer to generate comparable displacement to the full aperture but with 20% lower free field pressure.

Theoretical Framework to Predict Generalized Contrast-to-Noise Ratios of Photoacoustic Images With Applications to Computer Vision

The successful integration of computer vision, robotic actuation, and photoacoustic imaging to find and follow targets of interest during surgical and interventional procedures requires accurate photoacoustic target detectability. This detectability has traditionally been assessed with image quality metrics, such as contrast, contrast-to-noise ratio, and signal-to-noise ratio (SNR). However, predicting target tracking performance expectations when using these traditional metrics is difficult due to unbounded values and sensitivity to image manipulation techniques like thresholding. The generalized contrast-to-noise ratio (gCNR) is a recently introduced alternative target detectability metric, with previous work dedicated to empirical demonstrations of applicability to photoacoustic images. In this article, we present theoretical approaches to model and predict the gCNR of photoacoustic images with an associated theoretical framework to analyze relationships between imaging system parameters and computer vision task performance. Our theoretical gCNR predictions are validated with histogram-based gCNR measurements from simulated, experimental phantom, ex vivo, and in vivo datasets. The mean absolute errors between predicted and measured gCNR values ranged from 3.2 ×10^-3 to 2.3 ×10^-2 for each dataset, with channel SNRs ranging -40 to 40 dB and laser energies ranging 0.07 [Formula: see text] to 68 mJ. Relationships among gCNR, laser energy, target and background image parameters, target segmentation, and threshold levels were also investigated. Results provide a promising foundation to enable predictions of photoacoustic gCNR and visual servoing segmentation accuracy. The efficiency of precursory surgical and interventional tasks (e.g., energy selection for photoacoustic-guided surgeries) may also be improved with the proposed framework.

Ultrasound stimulation of the motor cortex during tonic muscle contraction

Transcranial ultrasound stimulation (tUS) shows potential as a noninvasive brain stimulation (NIBS) technique, offering increased spatial precision compared to other NIBS techniques. However, its reported effects on primary motor cortex (M1) are limited. We aimed to better understand tUS effects in human M1 by performing tUS of the hand area of M1 (M1hand) during tonic muscle contraction of the index finger. Stimulation during muscle contraction was chosen because of the transcranial magnetic stimulation-induced phenomenon known as cortical silent period (cSP), in which transcranial magnetic stimulation (TMS) of M1hand involuntarily suppresses voluntary motor activity. Since cSP is widely considered an inhibitory phenomenon, it presents an ideal parallel for tUS, which has often been proposed to preferentially influence inhibitory interneurons. Recording electromyography (EMG) of the first dorsal interosseous (FDI) muscle, we investigated effects on muscle activity both during and after tUS. We found no change in FDI EMG activity concurrent with tUS stimulation. Using single-pulse TMS, we found no difference in M1 excitability before versus after sparsely repetitive tUS exposure. Using acoustic simulations in models made from structural MRI of the participants that matched the experimental setups, we estimated in-brain pressures and generated an estimate of cumulative tUS exposure experienced by M1hand for each subject. We were unable to find any correlation between cumulative M1hand exposure and M1 excitability change. We also present data that suggest a TMS-induced MEP always preceded a near-threshold cSP.

SIMPA: an open-source toolkit for simulation and image processing for photonics and acoustics

SIGNIFICANCE

Optical and acoustic imaging techniques enable noninvasive visualisation of structural and functional properties of tissue. The quantification of measurements, however, remains challenging due to the inverse problems that must be solved. Emerging data-driven approaches are promising, but they rely heavily on the presence of high-quality simulations across a range of wavelengths due to the lack of ground truth knowledge of tissue acoustical and optical properties in realistic settings.

AIM

To facilitate this process, we present the open-source simulation and image processing for photonics and acoustics (SIMPA) Python toolkit. SIMPA is being developed according to modern software design standards.

APPROACH

SIMPA enables the use of computational forward models, data processing algorithms, and digital device twins to simulate realistic images within a single pipeline. SIMPA's module implementations can be seamlessly exchanged as SIMPA abstracts from the concrete implementation of each forward model and builds the simulation pipeline in a modular fashion. Furthermore, SIMPA provides comprehensive libraries of biological structures, such as vessels, as well as optical and acoustic properties and other functionalities for the generation of realistic tissue models.

RESULTS

To showcase the capabilities of SIMPA, we show examples in the context of photoacoustic imaging: the diversity of creatable tissue models, the customisability of a simulation pipeline, and the degree of realism of the simulations.

CONCLUSIONS

SIMPA is an open-source toolkit that can be used to simulate optical and acoustic imaging modalities. The code is available at: https://github.com/IMSY-DKFZ/simpa, and all of the examples and experiments in this paper can be reproduced using the code available at: https://github.com/IMSY-DKFZ/simpa_paper_experiments.

DeepUCT: Complex cascaded deep learning network for improved ultrasound tomography

Ultrasound computed tomography is an inexpensive and radiation-free medical imaging technique used to quantify the tissue acoustic properties for advanced clinical diagnosis. Image reconstruction in ultrasound tomography is often modeled as an optimization scheme solved by iterative methods like full-waveform inversion. These iterative methods are computationally expensive, while the optimization problem is ill-posed and nonlinear. To address this problem, we propose to use deep learning to overcome the computational burden and ill-posedness, and achieve near real-time image reconstruction in ultrasound tomography. We aim to directly learn the mapping from the recorded time-series sensor data to a spatial image of acoustical properties. To accomplish this, we develop a deep learning model using two cascaded convolutional neural networks with an encoder–decoder architecture. We achieve a good representation by first extracting the intermediate mapping-knowledge and later utilizing this knowledge to reconstruct the image. This approach is evaluated on synthetic phantoms where simulated ultrasound data are acquired from a ring of transducers surrounding the region of interest. The measurement data is acquired by forward modeling the wave equation using the k-wave toolbox. Our simulation results demonstrate that our proposed deep-learning method is robust to noise and significantly outperforms the state-of-the-art traditional iterative method both quantitatively and qualitatively. Furthermore, our model takes substantially less computational time than the conventional full-wave inversion method.

Evaluation of a pipeline for simulation, reconstruction, and classification in ultrasound-aided diffuse optical tomography of breast tumors

SIGNIFICANCE

Diffuse optical tomography is an ill-posed problem. Combination with ultrasound can improve the results of diffuse optical tomography applied to the diagnosis of breast cancer and allow for classification of lesions.

AIM

To provide a simulation pipeline for the assessment of reconstruction and classification methods for diffuse optical tomography with concurrent ultrasound information.

APPROACH

A set of breast digital phantoms with benign and malignant lesions was simulated building on the software VICTRE. Acoustic and optical properties were assigned to the phantoms for the generation of B-mode images and optical data. A reconstruction algorithm based on a two-region nonlinear fitting and incorporating the ultrasound information was tested. Machine learning classification methods were applied to the reconstructed values to discriminate lesions into benign and malignant after reconstruction.

RESULTS

The approach allowed us to generate realistic US and optical data and to test a two-region reconstruction method for a large number of realistic simulations. When information is extracted from ultrasound images, at least 75% of lesions are correctly classified. With ideal two-region separation, the accuracy is higher than 80%.

CONCLUSIONS

A pipeline for the generation of realistic ultrasound and diffuse optics data was implemented. Machine learning methods applied to a optical reconstruction with a nonlinear optical model and morphological information permit to discriminate malignant lesions from benign ones.

Quantitative Sound Speed Imaging of Cortical Bone and Soft Tissue: Results From Observational Data Sets

This work presents the first quantitative ultrasonic sound speed images of ex vivo limb cross-sections containing both soft tissue and bone using Full Waveform Inversion (FWI) with level set (LS) and travel time regularization. The estimated bulk sound speed of bone and soft tissue are within 10% and 1%, respectively, of ground truth estimates. The sound speed imagery shows muscle, connective tissue and bone features. Typically, ultrasound tomography (UST) using FWI is applied to imaging breast tissue properties (e.g. sound speed and density) that correlate with cancer. With further development, UST systems have the potential to deliver volumetric operator independent tissue property images of limbs with non-ionizing and portable hardware platforms. This work addresses the algorithmic challenges of imaging the sound speed of bone and soft tissue by combining FWI with LS regularization and travel time methods to recover soft tissue and bone sound speed with improved accuracy and reduced soft tissue artifacts when compared to conventional FWI. The value of leveraging LS and travel time methods is realized by evidence of improved bone geometry estimates as well as promising convergence properties and reduced risk of final model errors due to un-modeled shear wave propagation. Ex vivo bulk measurements of sound speed and MRI cross-sections validates the final inversion results.