Phase-dependent dual-frequency contrast imaging at sub-harmonic frequency

GPU accelerated study of a dual-frequency driven single bubble in a 6-dimensional parameter space: The active cavitation threshold

The active cavitation threshold of a dual-frequency driven single spherical gas bubble is studied numerically. This threshold is defined as the minimum intensity required to generate a given relative expansion (R_max-R_E)/R_E, where R_E is the equilibrium size of the bubble and R_max is the maximum bubble radius during its oscillation. The model employed is the Keller-Miksis equation that is a second order ordinary differential equation. The parameter space investigated is composed by the pressure amplitudes, excitation frequencies, phase shift between the two harmonic components and by the equilibrium bubble radius (bubble size). Due to the large 6-dimensional parameter space, the number of the parameter combinations investigated is approximately two billion. Therefore, the high performance of graphics processing units is exploited; our in-house code is written in C++ and CUDA C software environments. The results show that for (R_max-R_E)/R_E=2, the best choice of the frequency pairs depends on the bubble size. For small bubbles, below 3μm, the best option is to use just a single frequency of a low value in the giant response region. For medium sized bubbles, between 3μm and 6μm, the optimal choice is the mixture of low frequency (giant response) and main resonance frequency. For large bubbles, above 6μm, the main resonance dominates the active cavitation threshold. Increasing the prescribed relative expansion value to (R_max-R_E)/R_E=3, the optimal choice is always single frequency driving with the lowest value (20kHz here). Thus, in this case, the giant response always dominates the active cavitation threshold. The phase shift between the harmonic components of the dual-frequency driving (different frequency values) has no effect on the threshold.

Harmonic generation with a dual frequency pulse

Nonlinear imaging was implemented in commercial ultrasound systems over the last 15 years offering major advantages in many clinical applications. In this work, pulsing schemes coupled with a dual frequency pulse are presented. The pulsing schemes considered were pulse inversion, power modulation, and power modulated pulse inversion. The pulse contains a fundamental frequency f and a specified amount of its second harmonic 2f. The advantages and limitations of this method were evaluated with both acoustic measurements of harmonic generation and theoretical simulations based on the KZK equation. The use of two frequencies in a pulse results in the generation of the sum and difference frequency components in addition to the other harmonic components. While with single frequency pulses, only power modulation and power modulated pulse inversion contained odd harmonic components, with the dual frequency pulse, pulse inversion now also contains odd harmonic components.

The effect of amplitude modulation on subharmonic imaging with chirp excitation

Subharmonic generation from ultrasound contrast agents depends on the spectral and temporal properties of the excitation signal. The subharmonic response can be improved by using wideband and long-duration signals. However, for sinusoidal tone-burst excitation, the effective bandwidth of the signal is inversely proportional to the signal duration. Linear frequency-modulated (LFM) and nonlinear frequency-modulated (NLFM) chirp excitations allow independent control over the signal bandwidth and duration; therefore, in this study LFM and NLFM signals were used for the insonation of microbubble populations. The amplitude modulation of the excitation waveform was achieved by applying different window functions. A customized window was designed for the NLFM chirp excitation by focusing on reducing the spectral leakage at the subharmonic frequency and increasing the subharmonic generation from microbubbles. Subharmonic scattering from a microbubble population was measured for various excitation signals and window functions. At a peak negative pressure of 600 kPa, the generated subharmonic energy by ultrasound contrast agents was 15.4 dB more for NLFM chirp excitation with 40% fractional bandwidth when compared with tone-burst excitation. For this reason, the NLFM chirp with a customized window was used as an excitation signal to perform subharmonic imaging in an ultrasound flow phantom. Results showed that the NLFM waveform with a customized window improved the subharmonic contrast by 4.35 ± 0.42 dB on average over a Hann-windowed LFM excitation.

Second harmonic and subharmonic for non-linear wideband contrast imaging using a capacitive micromachined ultrasonic transducer array

When insonified with suitable ultrasound excitation, contrast microbubbles generate various non-linear scattered components, such as the second harmonic (2H) and the subharmonic (SH). In this study, we exploit the wide frequency bandwidth of capacitive micromachined ultrasonic transducers (CMUTs) to enhance the response from ultrasound contrast agents by selective imaging of both the 2H and SH components simultaneously. To this end, contrast images using the pulse inversion method were recorded with a 64-element CMUT linear array connected to an open scanner. In comparison to imaging at 2H alone, the wideband imaging including both the 2H and SH contributions provided up to 130% and 180% increases in the signal-to-noise and contrast-to-tissue ratios, respectively. The wide-frequency band of CMUTs offers new opportunities for improved ultrasound contrast agent imaging.

In vivo characterization of ultrasound contrast agents: microbubble spectroscopy in a chicken embryo

The dynamics of coated microbubbles was studied in an in vivo model. Biotinylated lipid-coated microbubbles were prepared in-house and were injected into a chick embryo chorioallantoic membrane (CAM) model on the fifth day of incubation. The microbubbles, ranging between 1.0 and 3.5 μm in diameter, were insonified in the frequency range of 4-7 MHz. Two amplitudes of acoustic pressure were applied: 300 kPa and 400 kPa. The fundamental and subharmonic responses were recorded optically with an ultra-fast camera (Brandaris 128) at 20 million frames per second. A subharmonic response was observed for 44% of the studied bubbles. From the data the frequency of the maximum fundamental and subharmonic response was derived for each individual bubble and resulted in the resonance curves of the microbubbles. All the bubbles showed shell (strain) hardening behavior for a higher acoustic pressure. We conclude that the subharmonic oscillations observed in this study belonged to the transmit at resonance (TR) regime.

Processing of subharmonic signals from ultrasound contrast agents to determine ambient pressures

Subharmonic-aided pressure estimation (SHAPE) is a technique that utilizes the subharmonic emissions, occurring at half the insonation frequency, from ultrasound contrast agents to estimate ambient pressures. The purpose of this work was to compare the performance of different processing techniques for the raw radiofrequency (rf) data acquired for SHAPE. A closed loop flow system was implemented circulating reconstituted Sonazoid (GE Healthcare, Oslo, Norway; 0.2 ml for 750 ml diluent) and the beam-formed unprocessed rf data were obtained from a 4 mm diameter lumen of a Doppler flow phantom (ATS Laboratories, Inc., Bridgeport, CT) using a SonixRP scanner (Ultrasonix, Richmond, BC, Canada). The transmit frequency and incident acoustic pressures were set to 2.5 MHz and 0.22 MPa, respectively, in order to elicit Sonazoid subharmonic emissions that are ambient-pressure sensitive. The time-varying ambient pressures within the flow phantom were recorded by a Millar pressure catheter. Four techniques for extracting the subharmonic amplitude from the rf data were tested along with two noise filtering techniques to process this data. Five filter orders were tested for the noise removing filters. The performance was evaluated based on the least root-mean-square errors reported after linear least-square regression analyses of the subharmonic data and the pressure catheter data and compared using a repeated ANOVA. When the subharmonic amplitudes were extracted as the mean value within a 0.2 MHz bandwidth about 1.25 MHz and when the resulting temporally-varying subharmonic signal was median filtered with an order of 500, the filtered subharmonic signal significantly predicted the ambient pressures (r2 = 0.90; p < 0.001) with the least error. The resulting root mean square and mean absolute errors were 8.16 +/- 0.26 mmHg and 6.70 +/- 0.17 mmHg, respectively. Thus, median processing the subharmonic data extracted as the mean value within a 0.2 MHz bandwidth about the theoretical subharmonic frequency turned out to be the best technique to process acoustic data for SHAPE. The implementation of this technique on ultrasound scanners may permit real-time SHAPE applications.

Phantom investigation of phase-inversion-based dual-frequency excitation imaging for improved contrast display

OBJECTIVE

The goal of this work is to examine the effects of pulse-inversion (PI) technique in combination with dual-frequency (DF) excitation method to separate the high-order nonlinear responses from microbubble contrast agents for improvement of image contrast. DF excitation method has been previously developed to induce the low-frequency ultrasound nonlinear responses from bubbles by using the composition of two high-frequency sinusoids (f(1) and f(2)).

MOTIVATION

Although the simple filtering was conventionally utilized to provide signal separation, the PI approach is better in the sense that it minimizes the mutual interferences among these high-order nonlinear responses in the presence of spectral overlap. The novelty of the work is that, in addition to the common PI summation, the PI subtraction was also applied in DF excitation method.

METHODS

DF excitation pulses having an envelope frequency of 3MHz (i.e., f(1)=8.5MHz and f(2)=11.5MHz) with pulse lengths of 3-10μs and the pressure amplitudes from 0.5 to 1.5MPa were used to interrogate the nonlinear responses of SonoVue™ microbubbles in the phantom experiments. The high-order nonlinear responses in the DF excitation were extracted for contrast imaging using PI summation for even-order nonlinear components or PI subtraction for odd-order nonlinear ones.

RESULTS

Our results indicated that, as compared to the conventional filtering technique, the PI processing effectively increases the contrast-to-tissue ratio (CTR) of the third-order nonlinear response at 5.5MHz and the fourth-order nonlinear response at 6MHz by 2-5dB. For these high-order nonlinear components, the CTR increase varies with the transmission pressures from 0.5 to 1.5MPa due to the microbubbles' displacement induced by the radiation force of DF excitation.

CONCLUSIONS

For DF excitation technique, the PI processing can help to extract either the odd-order or the even-order nonlinear components for higher CTR estimates.