Evaluation of acoustic-thermal simulations of in vivo magnetic resonance guided focused ultrasound ablative therapy.

PURPOSE: To evaluate numerical simulations of focused ultrasound (FUS) with a rabbit model, comparing simulated heating characteristics with magnetic resonance temperature imaging (MRTI) data collected during in vivo treatment. METHODS: A rabbit model was treated with FUS sonications in the biceps femoris with 3D MRTI collected. Acoustic and thermal properties of the rabbit muscle were determined experimentally. Numerical models of the rabbits were created, and tissue-type-specific properties were assigned. FUS simulations were performed using both the hybrid angular spectrum (HAS) method and k-Wave. Simulated power deposition patterns were converted to temperature maps using a Pennes' bioheat equation-based thermal solver. Agreement of pressure between the simulation techniques and temperature between the simulation and experimental heating was evaluated. Contributions of scattering and absorption attenuation were considered. RESULTS: Simulated peak pressures derived using the HAS method exceeded the simulated peak pressures from k-Wave by 1.6 ± 2.7%. The location and FWHM of the peak pressure calculated from HAS and k-Wave showed good agreement. When muscle acoustic absorption value in the simulations was adjusted to approximately 54% of the measured attenuation, the average root-mean-squared error between simulated and experimental spatial-average temperature profiles was 0.046 ± 0.019 °C/W. Mean distance between simulated and experimental COTMs was 3.25 ± 1.37 mm. Transverse FWHMs of simulated sonications were smaller than in in vivo sonications. Longitudinal FWHMs were similar. CONCLUSIONS: Presented results demonstrate agreement between HAS and k-Wave simulations and that FUS simulations can accurately predict focal position and heating for in vivo applications in soft tissue.

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.

Experimental validation of acoustic and thermal modeling in heterogeneous phantoms using the hybrid angular spectrum method.

PURPOSE: The aim was to quantitatively validate the hybrid angular spectrum (HAS) algorithm, a rapid wave propagation technique for heterogeneous media, with both pressure and temperature measurements. METHODS: Heterogeneous tissue-mimicking phantoms were used to evaluate the accuracy of the HAS acoustic modeling algorithm in predicting pressure and thermal patterns. Acoustic properties of the phantom components were measured by a through-transmission technique while thermal properties were measured with a commercial probe. Numerical models of each heterogeneous phantom were segmented from 3D MR images. Cylindrical phantoms 30-mm thick were placed in the pre-focal field of a focused ultrasound beam and 2D pressure measurements obtained with a scanning hydrophone. Peak pressure, full width at half maximum, and normalized root mean squared difference (RMSDn) between the measured and simulated patterns were compared. MR-guided sonications were performed on 150-mm phantoms to obtain MR temperature measurements. Using HAS-predicted power density patterns, temperature simulations were performed. Experimental and simulated temperature patterns were directly compared using peak and mean temperature plots, RMSDn metrics, and accuracy of heating localization. RESULTS: The average difference between simulated and hydrophone-measured peak pressures was 9.0% with an RMSDn of 11.4%. Comparison of the experimental MRI-derived and simulated temperature patterns showed RMSDn values of 10.2% and 11.1% and distance differences between the centers of thermal mass of 2.0 and 2.2 mm. CONCLUSIONS: These results show that the computationally rapid hybrid angular spectrum method can predict pressure and temperature patterns in heterogeneous models, including uncertainties in property values and other parameters, to within approximately 10%.

Development and characterization of a tissue mimicking psyllium husk gelatin phantom for ultrasound and magnetic resonance imaging.

Purpose: To develop and characterize a tissue-mimicking phantom that enables the direct comparison of magnetic resonance (MR) and ultrasound (US) imaging techniques useful for monitoring high-intensity focused ultrasound (HIFU) treatments. With no additions, gelatin phantoms produce little if any scattering required for US imaging. This study characterizes the MR and US image characteristics as a function of psyllium husk concentration, which was added to increase US scattering.Methods: Gelatin phantoms were constructed with varying concentrations of psyllium husk. The effects of psyllium husk concentration on US B-mode and MR imaging were evaluated at nine different concentrations. T1, T2, and T2* MR maps were acquired. Acoustic properties (attenuation and speed of sound) were measured at frequencies of 0.6, 1.0, 1.8, and 3.0 MHz using a through-transmission technique. Phantom elastic properties were evaluated for both time and temperature dependence.Results: Ultrasound image echogenicity increased with increasing psyllium husk concentration while quality of gradient-recalled echo MR images decreased with increasing concentration. For all phantoms, the measured speed of sound ranged between 1567-1569 m/s and the attenuation ranged between 0.42-0.44 dB/(cm·MHz). Measured T1 ranged from 974-1051 ms. The T2 and T2* values ranged from 97-108 ms and 48-88 ms, respectively, with both showing a decreasing trend with increased psyllium husk concentration. Phantom stiffness, measured using US shear-wave speed measurements, increased with age and decreased with increasing temperature.Conclusions: The presented dual-use tissue-mimicking phantom is easy to manufacture and can be used to compare and evaluate US-guided and MR-guided HIFU imaging protocols.

Efficient shear wave elastography using transient acoustic radiation force excitations and MR displacement encoding.

PURPOSE: To present a novel MR shear wave elastography (MR-SWE) method that efficiently measures the speed of propagating wave packets generated using acoustic radiation force (ARF) impulses. METHODS: ARF impulses from a focused ultrasound (FUS) transducer were applied sequentially to a preselected set of positions and motion encoded MRI was used to acquire volumetric images of the propagating shear wavefront emanating from each point. The wavefront position at multiple propagation times was encoded in the MR phase image using a train of motion encoding gradient lobes. Generating a transient propagating wavefront at multiple spatial positions and sampling each at multiple time-points allowed for shear wave speed maps to be efficiently created. MR-SWE was evaluated in tissue mimicking phantoms and ex vivo bovine liver tissue before and after ablation. RESULTS: MR-SWE maps, covering an in-plane area of ~5 × 5 cm, were acquired in 12 s for a single slice and 144 s for a volumetric scan. MR-SWE detected inclusions of differing stiffness in a phantom experiment. In bovine liver, mean shear wave speed significantly increased from 1.65 ± 0.18 m/s in normal to 2.52 ± 0.18 m/s in ablated region (n = 581 pixels; P-value < 0.001). CONCLUSION: MR-SWE is an elastography technique that enables precise targeting and excitation of the desired tissue of interest. MR-SWE may be particularly well suited for treatment planning and endpoint assessment of MR-guided FUS procedures because the same device used for therapy can be used as an excitation source for tissue stiffness quantification.

Validation of hybrid angular spectrum acoustic and thermal modelling in phantoms.

In focused ultrasound (FUS) thermal ablation of diseased tissue, acoustic beam and thermal simulations enable treatment planning and optimization. In this study, a treatment-planning methodology that uses the hybrid angular spectrum (HAS) method and the Pennes' bioheat equation (PBHE) is experimentally validated in homogeneous tissue-mimicking phantoms. Simulated three-dimensional temperature profiles are compared to volumetric MR thermometry imaging (MRTI) of FUS sonications in the phantoms, whose acoustic and thermal properties are independently measured. Additionally, Monte Carlo (MC) uncertainty analysis is performed to quantify the effect of tissue property uncertainties on simulation results. The mean error between simulated and experimental spatiotemporal peak temperature rise was +0.33°C (+6.9%). Despite this error, the experimental temperature rise fell within the expected uncertainty of the simulation, as determined by the MC analysis. The average errors of the simulated transverse and longitudinal full width half maximum (FWHM) of the profiles were -1.9% and 7.5%, respectively. A linear regression and local sensitivity analysis revealed that simulated temperature amplitude is more sensitive to uncertainties in simulation inputs than in the profile width and shape. Acoustic power, acoustic attenuation and thermal conductivity had the greatest impact on peak temperature rise uncertainty; thermal conductivity and volumetric heat capacity had the greatest impact on FWHM uncertainty. This study validates that using the HAS and PBHE method can adequately predict temperature profiles from single sonications in homogeneous media. Further, it informs the need to accurately measure or predict patient-specific properties for improved treatment planning of ablative FUS surgeries.

Experimental assessment of phase aberration correction for breast MRgFUS therapy.

PURPOSE: This study validates that phase aberrations in breast magnetic resonance-guided focussed ultrasound (MRgFUS) therapies can be corrected in a clinically relevant time frame to generate more intense, smaller and more spatially accurate foci. MATERIALS AND METHODS: Hybrid angular spectrum (HAS) ultrasound calculations in an magnetic resonance imaging (MRI)-based tissue model, were used to compute phase aberration corrections for improved experimental MRgFUS heating in four heterogeneous breast-mimicking phantoms (n = 18 total locations). Magnetic resonance(MR) temperature imaging was used to evaluate the maximum temperature rise, focus volume and focus accuracy for uncorrected and phase aberration-corrected sonications. Thermal simulations assessed the effectiveness of the phase aberration correction implementation. RESULTS: In 13 of 18 locations, the maximum temperature rise increased by an average of 30%, focus volume was reduced by 40% and focus accuracy improved from 4.6 to 3.6 mm. Mixed results were observed in five of the 18 locations, with focus accuracy improving from 6.1 to 2.5 mm and the maximum temperature rise decreasing by 8% and focus volume increasing by 10%. Overall, the study demonstrated significant improvements (p < 0.005) in maximum temperature rise, focus volume and focus accuracy. Simulations predicted greater improvements than observed experimentally, suggesting potential for improvement in implementing the technique. The complete phase aberration correction procedure, including model generation, segmentation and phase aberration computations, required less than 45 min per sonication location. CONCLUSION: The significant improvements demonstrated in this study i.e., focus intensity, size and accuracy from phase aberration correction have the potential to improve the efficacy, time-efficiency and safety of breast MRgFUS therapies.

Development and validation of a MRgHIFU non-invasive tissue acoustic property estimation technique.

MR-guided high-intensity focussed ultrasound (MRgHIFU) non-invasive ablative surgeries have advanced into clinical trials for treating many pathologies and cancers. A remaining challenge of these surgeries is accurately planning and monitoring tissue heating in the face of patient-specific and dynamic acoustic properties of tissues. Currently, non-invasive measurements of acoustic properties have not been implemented in MRgHIFU treatment planning and monitoring procedures. This methods-driven study presents a technique using MR temperature imaging (MRTI) during low-temperature HIFU sonications to non-invasively estimate sample-specific acoustic absorption and speed of sound values in tissue-mimicking phantoms. Using measured thermal properties, specific absorption rate (SAR) patterns are calculated from the MRTI data and compared to simulated SAR patterns iteratively generated via the Hybrid Angular Spectrum (HAS) method. Once the error between the simulated and measured patterns is minimised, the estimated acoustic property values are compared to the true phantom values obtained via an independent technique. The estimated values are then used to simulate temperature profiles in the phantoms, and compared to experimental temperature profiles. This study demonstrates that trends in acoustic absorption and speed of sound can be non-invasively estimated with average errors of 21% and 1%, respectively. Additionally, temperature predictions using the estimated properties on average match within 1.2 °C of the experimental peak temperature rises in the phantoms. The positive results achieved in tissue-mimicking phantoms presented in this study indicate that this technique may be extended to in vivo applications, improving HIFU sonication temperature rise predictions and treatment assessment.

Kinetics of Ultrasonic Drug Delivery from Targeted Micelles.

To minimize the adverse side effects of conventional chemotherapy, a targeted micellar drug carrier was investigated that retains hydrophobic drugs in its core and then releases the drug via ultrasonic activation. This paper compares the percent drug release from folated versus non-folated micelles by insonation at 70 kHz and different acoustic power densities. The encapsulated drug is Doxoru- bicin (Dox). A physical model of zero-order release with first-order re-encapsulation was used to fit the experimental kinetic data. Additionally, the acoustic activation power density and Gibbs free energy were introduced and calculated for folated and non-targeted micelles. The data suggests an important role of inertial cavitation in drug release and the presence of a power density threshold for inertial cavitation.

The accuracy and precision of two non-invasive, magnetic resonance-guided focused ultrasound-based thermal diffusivity estimation methods.

PURPOSE: The use of correct tissue thermal diffusivity values is necessary for making accurate thermal modelling predictions during magnetic resonance-guided focused ultrasound (MRgFUS) treatment planning. This study evaluates the accuracy and precision of two non-invasive thermal diffusivity estimation methods, a Gaussian temperature method and a Gaussian specific absorption rate (SAR) method. MATERIALS AND METHODS: Both methods utilise MRgFUS temperature data obtained during cooling following a short (<25 s) heating pulse. The Gaussian SAR method can also use temperatures obtained during heating. Experiments were performed at low heating levels (ΔT∼10 °C) in ex vivo pork muscle and in vivo rabbit back muscle. The non-invasive MRgFUS thermal diffusivity estimates were compared with measurements from two standard invasive methods. RESULTS: Both non-invasive methods accurately estimated thermal diffusivity when using MR temperature cooling data (overall ex vivo error <6%, in vivo <12%). Including heating data in the Gaussian SAR method further reduced errors (ex vivo error <2%, in vivo <3%). The significantly lower standard deviation values (p < 0.03) of the Gaussian SAR method indicated that it had better precision than the Gaussian temperature method. CONCLUSIONS: With repeated sonications, either MR-based method could provide accurate thermal diffusivity values for MRgFUS therapies. Fitting to more data simultaneously likely made the Gaussian SAR method less susceptible to noise, and using heating data helped it converge more consistently to the FUS fitting parameters and thermal diffusivity. These effects led to the improved precision of the Gaussian SAR method.

Adaptive model-predictive controller for magnetic resonance guided focused ultrasound therapy.

PURPOSE: Minimising treatment time and protecting healthy tissues are conflicting goals that play major roles in making magnetic resonance image-guided focused ultrasound (MRgFUS) therapies clinically practical. We have developed and tested in vivo an adaptive model-predictive controller (AMPC) that reduces treatment time, ensures safety and efficacy, and provides flexibility in treatment set-up. MATERIALS AND METHODS: The controller realises time savings by modelling the heated treatment cell's future temperatures and thermal dose accumulation in order to anticipate the optimal time to switch to the next cell. Selected tissues are safeguarded by a configurable temperature constraint. Simulations quantified the time savings realised by each controller feature as well as the trade-offs between competing safety and treatment time parameters. In vivo experiments in rabbit thighs established the controller's effectiveness and reliability. RESULTS: In all in vivo experiments the target thermal dose of at least 240 CEM43 was delivered everywhere in the treatment volume. The controller's temperature safety limit reliably activated and constrained all protected tissues to <9 CEM43. Simulations demonstrated the path independence of the controller, and that a path which successively proceeds to the hottest untreated neighbouring cell leads to significant time savings, e.g. when compared to a concentric spiral path. Use of the AMPC produced a compounding time-saving effect; reducing the treatment cells' heating times concurrently reduced heating of normal tissues, which eliminated cooling periods. CONCLUSIONS: Adaptive model-predictive control can automatically deliver safe, effective MRgFUS treatments while significantly reducing treatment times.

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

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

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

BACKGROUND: Ultrasound beam aberration correction is vital when focusing ultrasound through the skull bone in transcranial magnetic resonance-guided focused ultrasound (tcMRgFUS) applications. Current methods make transducer element phase adjustments to compensate for the variation in skull properties (shape, thickness, and acoustic properties), but do not account for variations in the internal brain anatomy. PURPOSE: Our objective is to investigate the effect of cerebrospinal fluid (CSF) and brain anatomy on beam focusing in tcMRgFUS treatments. METHODS: Simulations were conducted with imaging data from 20 patients previously treated with focused ultrasound for disabling tremor. The Hybrid Angular Spectrum (HAS) method was used to test the effect of including cerebral spinal fluid (CSF) and brain anatomy in determining the element phases used for aberration correction and beam focusing. Computer tomography (CT) and magnetic resonance imaging (MRI) images from patient treatments were used to construct a segmented model of each patient's head. The segmented model for treatment simulation consisted of water, skin, fat, brain, CSF, diploë, and cortical bone. Transducer element phases used for treatment simulation were determined using time reversal from the desired focus, generating a set of phases assuming a homogeneous brain in the intracranial volume, and a second set of phases assigning CSF acoustic properties to regions of CSF. In addition, for three patients, the relative effect of separately including CSF speed of sound values compared to CSF attenuation values was found. RESULTS: We found that including CSF acoustic properties (speed of sound and attenuation) during phase planning compared to phase correction without considering CSF increased the absorbed ultrasound power density ratios at the focus over a range of 1.06 to 1.29 (mean of 17% ± 6%) for 20 patients. Separately considering the CSF speed of sound and CSF attenuation showed that the increase was due almost entirely to including the CSF speed of sound; considering only the CSF attenuation had a negligible effect. CONCLUSIONS: Based on HAS simulations, treatment planning phase determination using morphologically realistic CSF and brain anatomy yielded an increase of up to 29% in the ultrasound focal absorbed power density. Future work will be required to validate the CSF simulations.

Thermal sensitivity of endothelial cells on synthetic vascular graft material.

PURPOSE: The goal is to identify thermal exposures capable of reducing or eliminating cell survival on expanded polytetrafluoroethylene (ePTFE), in an effort to develop a mild hyperthermia treatment of neointimal hyperplasia in ePTFE vascular grafts. MATERIALS AND METHODS: Viable and dead bovine aortic endothelial cells were quantified following different thermal exposure conditions: cells on collagen-coated ePTFE sheets or tissue culture polystyrene dishes were heated at 42° and 45°C to determine their thermal sensitivity on different surfaces, and cells cultured on collagen-coated ePTFE sheets were heated at 43-50°C for various durations, followed by incubation at 37°C for 0 and 20 h, respectively. Significant cell death was set to be 50%. Two types of cell death, apoptosis and necrosis, were distinguished by cell morphology and membrane integrity assessments. RESULTS: The attachment and survival of cells on ePTFE sheets were more sensitive to inhibition by mild heating than those on tissue culture dishes. Exposure to 45°C for 90 min and 50°C for 30 min caused significant necrotic cell death on ePTFE (65% and 75%, respectively). A 37°C/20-h incubation following 30-min exposures at 47° and 50°C increased total cell death (necrosis + apoptosis) from 20% to 50% and 75% to 100%, respectively. CONCLUSION: Cells grown on ePTFE were more susceptible to mild hyperthermia-induced death, compared to those on tissue culture dishes. Significant cell death on ePTFE mainly via apoptosis can be achieved by optimising temperature and duration of exposure.

The effect of electronically steering a phased array ultrasound transducer on near-field tissue heating.

PURPOSE: This study presents the results obtained from both simulation and experimental techniques that show the effect of mechanically or electronically steering a phased array transducer on proximal tissue heating. METHODS: The thermal response of a nine-position raster and a 16-mm diameter circle scanning trajectory executed through both electronic and mechanical scanning was evaluated in computer simulations and experimentally in a homogeneous tissue-mimicking phantom. Simulations were performed using power deposition maps obtained from the hybrid angular spectrum (HAS) method and applying a finite-difference approximation of the Pennes' bioheat transfer equation for the experimentally used transducer and also for a fully sampled transducer to demonstrate the effect of acoustic window, ultrasound beam overlap and grating lobe clutter on near-field heating. RESULTS: Both simulation and experimental results show that electronically steering the ultrasound beam for the two trajectories using the 256-element phased array significantly increases the thermal dose deposited in the near-field tissues when compared with the same treatment executed through mechanical steering only. In addition, the individual contributions of both beam overlap and grating lobe clutter to the near-field thermal effects were determined through comparing the simulated ultrasound beam patterns and resulting temperature fields from mechanically and electronically steered trajectories using the 256-randomized element phased array transducer to an electronically steered trajectory using a fully sampled transducer with 40 401 phase-adjusted sample points. CONCLUSIONS: Three distinctly different three distinctly different transducers were simulated to analyze the tradeoffs of selected transducer design parameters on near-field heating. Careful consideration of design tradeoffs and accurate patient treatment planning combined with thorough monitoring of the near-field tissue temperature will help to ensure patient safety during an MRgHIFU treatment.

Modelling ultrasound-induced mild hyperthermia of hyperplasia in vascular grafts.

BACKGROUND: Expanded polytetrafluoroethylene (ePTFE) vascular grafts frequently develop occlusive neointimal hyperplasia as a result of myofibroblast over-growth, leading to graft failure. ePTFE exhibits higher ultrasound attenuation than native soft tissues. We modelled the selective absorption of ultrasound by ePTFE, and explored the feasibility of preventing hyperplasia in ePTFE grafts by ultrasound heating. Specifically, we simulated the temperature profiles of implanted grafts and nearby soft tissues and blood under ultrasound exposure. The goal was to determine whether ultrasound exposure of an ePTFE graft can generate temperatures sufficient to prevent cell growth on the graft without damaging nearby soft tissues and blood. METHODS: Ultrasound beams from two transducers (1.5 and 3.2 MHz) were simulated in two graft/tissue models, with and without an intra-graft cellular layer mimicking hyperplasia, using the finite-difference time-domain (FDTD) method. The resulting power deposition patterns were used as a heat source for the Pennes bioheat equation in a COMSOL(®) Multiphysics heat transfer model. 50°C is known to cause cell death and therefore the transducer powers were adjusted to produce a 13°C temperature rise from 37°C in the ePTFE. RESULTS: Simulations showed that both the frequency of the transducers and the presence of hyperplasia significantly affect the power deposition patterns and subsequent temperature profiles on the grafts and nearby tissues. While neither transducer significantly raised the temperature of the blood, the 1.5-MHz transducer was less focused and heated larger volumes of the graft and nearby soft tissues than the 3.2-MHz transducer. The presence of hyperplasia had little effect on the blood's temperature, but further increased the temperature of the graft and nearby soft tissues in response to either transducer. Skin cooling and blood flow play a significant role in preventing overheating of the native tissues. CONCLUSIONS: Modelling shows that ultrasound can selectively heat ePTFE grafts and produce temperatures that cause cell death on the graft. The temperature increase in blood is negligible and that in the adjacent soft tissues may be minimized by skin cooling and using appropriate transducers. Therefore, ultrasound heating may have the potential to reduce neointimal hyperplasia and failure of ePTFE vascular grafts.

Extension of the angular spectrum method to calculate pressure from a spherically curved acoustic source.

The angular spectrum method is an accurate and computationally efficient method for modeling acoustic wave propagation. The use of the typical 2D fast Fourier transform algorithm makes this a fast technique but it requires that the source pressure (or velocity) be specified on a plane. Here the angular spectrum method is extended to calculate pressure from a spherical transducer-as used extensively in applications such as magnetic resonance-guided focused ultrasound surgery-to a plane. The approach, called the Ring-Bessel technique, decomposes the curved source into circular rings of increasing radii, each ring a different distance from the intermediate plane, and calculates the angular spectrum of each ring using a Fourier series. Each angular spectrum is then propagated to the intermediate plane where all the propagated angular spectra are summed to obtain the pressure on the plane; subsequent plane-to-plane propagation can be achieved using the traditional angular spectrum method. Since the Ring-Bessel calculations are carried out in the frequency domain, it reduces calculation times by a factor of approximately 24 compared to the Rayleigh-Sommerfeld method and about 82 compared to the Field II technique, while maintaining accuracies of better than 96% as judged by those methods for cases of both solid and phased-array transducers.

Magnetic resonance shear wave elastography using transient acoustic radiation force excitations and sinusoidal displacement encoding.

A magnetic resonance (MR) shear wave elastography technique that uses transient acoustic radiation force impulses from a focused ultrasound (FUS) transducer and a sinusoidal-shaped MR displacement encoding strategy is presented. Using this encoding strategy, an analytic expression for calculating the shear wave speed in a heterogeneous medium was derived. Green's function-based simulations were used to evaluate the feasibility of calculating shear wave speed maps using the analytic expression. Accuracy of simulation technique was confirmed experimentally in a homogeneous gelatin phantom. The elastography measurement was compared to harmonic MR elastography in a homogeneous phantom experiment and the measured shear wave speed values differed by less than 14%. This new transient elastography approach was able to map the position and shape of inclusions sized from 8.5 to 14 mm in an inclusion phantom experiment. These preliminary results demonstrate the feasibility of using a straightforward analytic expression to generate shear wave speed maps from MR images where sinusoidal-shaped motion encoding gradients are used to encode the displacement-time history of a transiently propagating wave-packet. This new measurement technique may be particularly well suited for performing elastography before, during, and after MR-guided FUS therapies since the same device used for therapy is also used as an excitation source for elastography.

In-plane parallel scanning: a microarray technology for point-of-care testing.

A new microarray technology is described for rapid, inexpensive, multiplex diagnostics assays. Referred to as "in-plane parallel scanning" (IPPS), this technology replaces expensive laser scanning with a grid of 100-µm-wide waveguides embedded in the chip's substrate, enabling real-time quantification of molecular complex formation on the chip's surface. Compared to conventional microarray technology, IPPS has advantages of shorter assay time and lower instrument cost and complexity so that the platform can potentially be used in point-of-care (POC) settings. Two different chip formats are described: a low-density microarray with 10 sensing wells (IPPS-10) and a medium-density one with 100 sensing wells (IPPS-100). Performance was evaluated in two different proof-of-principle immunoassays: interleukin-1ß (IL-1ß) and Clostridium difficile toxin A. The two assays gave similar limits of detection of 0.67 and 0.94 pM, respectively. A saturation kinetics model described the sensor response with apparent dissociation constants of 511 pM for IL-1ß and 6.47 nM for C. difficile toxin A toxoid. The multiplexing capabilities of the IPPS technology were also demonstrated in a multiplex assay for both analytes on the same IPPS-10 chip. Based on these results, the IPPS technology holds promise for translating diagnostic microarrays into near-patient environments.

Minimisation of HIFU pulse heating and interpulse cooling times.

This study presents results from a new optimisation technique that reduces HIFU treatment times by minimising individual heating and interpulse cooling times while adhering to normal tissue constraint limits at each sonication position. The potential clinical usefulness of this technique is demonstrated through its implementation in three dimentsional (3D) simulations of HIFU treatments for a range of tumour geometries, normal tissue constraint values, tissue perfusion levels and focal zone scanning path trajectories, all studied as a function of the applied power magnitude. When compared to typical open loop values the optimised treatment times were lower for all conditions studied, including when treatment-limiting normal tissue thermal build-up was present. While use of this technique guarantees minimum pulse heating and interpulse cooling times for each pulse, the total treatment time gains realised depend on the individual clinical treatment configuration. In combination with a judiciously selected scan path, use of the pulse time optimisation procedure reduced treatment times in a small, superficial tumour by 85%. In addition, in all cases studied the use of an increased applied power always decreased the treatment time, including cases when significant normal tissue thermal build-up was present. Importantly, the power maximisation and pulse time minimisation procedures can be applied independently of the optimisation of the focal zone's scan path, size and shape. Given the basic nature, universal applicability and ready clinical adaptability for use in real time model predictive control, the pulse time minimisation and power maximisation approaches have significant clinical promise for reducing HIFU treatment times.