"Compression-only" behavior: a second-order nonlinear response of ultrasound contrast agent microbubbles

Oscillating phospholipid-coated ultrasound contrast agent microbubbles display a so-called "compression-only" behavior, where it is observed that the bubbles compress efficiently while their expansion is suppressed. Here, a theoretical understanding of the source of this nonlinear behavior is provided through a weakly nonlinear analysis of the shell buckling model proposed by Marmottant et al. [J. Acoust. Soc. Am. 118, 3499-3505 (2005)]. It is shown that the radial dynamics of the bubble can be considered as a superposition of a linear response at the fundamental driving frequency and a second-order nonlinear low-frequency response that describes the negative offset of the mean bubble radius. The analytical solution deduced from the weakly nonlinear analysis shows that the compression-only behavior results from a rapid change of the shell elasticity with bubble radius. In addition, the radial dynamics of single phospholipid-coated microbubbles was recorded as a function of both the amplitude and the frequency of the driving pressure pulse. The comparison between the experimental data and the theory shows that the magnitude of compression-only behavior is mainly determined by the initial phospholipids concentration on the bubble surface, which slightly varies from bubble to bubble.     

Modeling subharmonic response from contrast microbubbles as a function of ambient static pressure

Variation of subharmonic response from contrast microbubbles with ambient pressure is numerically investigated for non-invasive monitoring of organ-level blood pressure. Previously, several contrast microbubbles both in vitro and in vivo registered approximately linear (5-15 dB) subharmonic response reduction with 188 mm Hg change in ambient pressure. In contrast, simulated subharmonic response from a single microbubble is seen here to either increase or decrease with ambient pressure. This is shown using the code BUBBLESIM for encapsulated microbubbles, and then the underlying dynamics is investigated using a free bubble model. The ratio of the excitation frequency to the natural frequency of the bubble is the determining parameter--increasing ambient pressure increases natural frequency thereby changing this ratio. For frequency ratio below a lower critical value, increasing ambient pressure monotonically decreases subharmonic response. Above an upper critical value of the same ratio, increasing ambient pressure increases subharmonic response; in between, the subharmonic variation is non-monotonic. The precise values of frequency ratio for these three different trends depend on bubble radius and excitation amplitude. The modeled increase or decrease of subharmonic with ambient pressure, when one happens, is approximately linear only for certain range of excitation levels. Possible reasons for discrepancies between model and previous experiments are discussed.     

调频激励增强微气泡的次谐波理论和实验研究

超声造影剂的次谐波成像可以提高造影组织比，提供更好的图像质量. 提出一种利用调频脉冲激励以增强造影剂微气泡产生的次谐波新方法. 基于修正的Church方程，从理论上讨论了次谐波的产生与调频激励声压的关系及产生阈值，并且实验证实了优化调频信号的带宽及调频时间可以提高次谐波信号幅度及改善主瓣和旁瓣特性. 理论与实验表明，与传统脉冲信号激励相比，调频信号激励产生的次谐波幅度可提高约22dB. 关键词： 调频激励 超声造影剂 微气泡 次谐波     

Resonance frequencies of lipid-shelled microbubbles in the regime of nonlinear oscillations

Knowledge of resonant frequencies of contrast microbubbles is important for the optimization of ultrasound contrast imaging and therapeutic techniques. To date, however, there are estimates of resonance frequencies of contrast microbubbles only for the regime of linear oscillation. The present paper proposes an approach for evaluating resonance frequencies of contrast agent microbubbles in the regime of nonlinear oscillation. The approach is based on the calculation of the time-averaged oscillation power of the radial bubble oscillation. The proposed procedure was verified for free bubbles in the frequency range 1-4 MHz and then applied to lipid-shelled microbubbles insonified with a single 20-cycle acoustic pulse at two values of the acoustic pressure amplitude, 100 kPa and 200 kPa, and at four frequencies: 1.5, 2.0, 2.5, and 3.0 MHz. It is shown that, as the acoustic pressure amplitude is increased, the resonance frequency of a lipid-shelled microbubble tends to decrease in comparison with its linear resonance frequency. Analysis of existing shell models reveals that models that treat the lipid shell as a linear viscoelastic solid appear may be challenged to provide the observed tendency in the behavior of the resonance frequency at increasing acoustic pressure. The conclusion is drawn that the further development of shell models could be improved by the consideration of nonlinear rheological laws.     

多层膜磁性微泡的非线性声振动特性

超顺磁性氧化铁纳米粒子与造影剂微泡结合形成磁性微泡,用于产生多模态造影剂,以增强医学超声和磁共振成像.将装载有纳米磁性颗粒的微泡包膜层看作由磁流体膜与磷脂膜组合而成的双层膜结构,同时考虑磁性纳米颗粒体积分数a对膜密度及黏度的影响,从气泡动力学基本理论出发,构建多层膜结构磁性微泡非线性动力学方程.数值分析了驱动声压和频率等声场参数、颗粒体积分数、膜层厚度以及表面张力等膜壳参数对微泡声动力学行为的影响.结果表明,当磁性颗粒体积分数较小且a≤0.1时,磁性微泡声响应特性与普通包膜微泡相似,微泡的声频响应与其初始尺寸和驱动压有关;当驱动声场频率f为磁性微泡共振频率f0的2倍(f=2f0)时,微泡振动失稳临界声压最低;磁性颗粒的存在抑制了泡的膨胀和收缩但抑制效果非常有限;磁性微泡外膜层材料的表面张力参数K及膜层厚度d也会影响微泡的振动,当表面张力参数及膜厚取值分别为0.2—0.4 N/m及50—150 nm时,可观察到气泡存在不稳定振动响应区.     

超声造影剂微泡非线性动力学响应的机理及相关应用

随着生命科学及现代医学的发展, 一体化无创精准诊疗已经日益成为人们关注的焦点问题, 而关于超声造影剂微泡的非线性效应的相关机理、动力学建模及其在超声医学领域中的应用研究也得到了极大的推动. 本文对下列课题进行了总结和讨论, 包括: 1)基于Mie散射技术和流式细胞仪对造影剂微泡参数进行定征的一体化解决方案; 2)通过对微泡包膜的黏弹特性进行非线性修正, 构建新的包膜微泡动力学模型; 3)探索造影剂惯性空化阈值与其包膜参数之间的相关性; 以及4)研究超声联合造影剂微泡促进基因/药物转染效率并有效降低其生物毒性的相关机理.     

Comparison between maximum radial expansion of ultrasound contrast agents and experimental postexcitation signal results

Experimental postexcitation signal data of collapsing Definity microbubbles are compared with the Marmottant theoretical model for large amplitude oscillations of ultrasound contrast agents (UCAs). After taking into account the insonifying pulse characteristics and size distribution of the population of UCAs, a good comparison between simulated results and previously measured experimental data is obtained by determining a threshold maximum radial expansion (Rmax) to indicate the onset of postexcitation. This threshold Rmax is found to range from 3.4 to 8.0 times the initial bubble radius, R0, depending on insonification frequency. These values are well above the typical free bubble inertial cavitation threshold commonly chosen at 2R0. The close agreement between the experiment and models suggests that lipid-shelled UCAs behave as unshelled bubbles during most of a large amplitude cavitation cycle, as proposed in the Marmottant equation.     

Ultrasonic bubbles in medicine: influence of the shell

Ultrasound contrast agents consist of microscopically small bubbles encapsulated by an elastic shell. These microbubbles oscillate upon ultrasound insonification, and demonstrate highly nonlinear behavior, ameliorating their detectability. (Potential) medical applications involving the ultrasonic disruption of contrast agent microbubble shells include release-burst imaging, localized drug delivery, and noninvasive blood pressure measurement. To develop and enhance these techniques, predicting the cracking behavior of ultrasound-insonified encapsulated microbubbles has been of importance. In this paper, we explore microbubble behavior in an ultrasound field, with special attention to the influence of the bubble shell. A bubble in a sound field can be considered a forced damped harmonic oscillator. For encapsulated microbubbles, the presence of a shell has to be taken into account. In models, an extra damping parameter and a shell stiffness parameter have been included, assuming that Hooke's Law holds for the bubble shell. At high acoustic amplitudes, disruptive phenomena have been observed, such as microbubble fragmentation and ultrasonic cracking. We analyzed the occurrence of ultrasound contrast agent fragmentation, by simulating the oscillating behavior of encapsulated microbubbles with various sizes in a harmonic acoustic field. Fragmentation occurs exclusively during the collapse phase and occurs if the kinetic energy of the collapsing microbubble is greater than the instantaneous bubble surface energy, provided that surface instabilities have grown big enough to allow for break-up. From our simulations it follows that the Blake critical radius is not a good approximation for a fragmentation threshold. We demonstrated how the phase angle differences between a damped radially oscillating bubble and an incident sound field depend on shell parameters.     

Resonance behaviors of encapsulated microbubbles oscillating nonlinearly with ultrasonic excitation

The resonance behaviors of a few lipid-coated microbubbles acoustically activated in viscoelastic media were comprehensively examined via radius response analysis. The size polydispersity and random spatial distribution of the interacting microbubbles, the rheological properties of the lipid shell and the viscoelasticity of the surrounding medium were considered simultaneously. The obtained radius response curves present a successive occurrence of linear resonances, nonlinear harmonic and sub-harmonic resonances with the acoustic pressure increasing. The microbubble resonance is radius-, pressure- and frequency-dependent. Specifically, the maximum bubble expansion ratio at the main resonance peak increases but the resonant radius decreases as the ultrasound pressure increases, while both of them decrease with the ultrasound frequency increasing. Moreover, compared to an isolated microbubble case, it is found that large microbubbles in close proximity prominently suppress the resonant oscillations while slightly increase the resonant radii for both harmonic and subharmonic resonances, even leading to the disappearance of the subharmonic resonance with the influences increasing to a certain degree. In addition, the results also suggest that both the encapsulating shell and surrounding medium can substantially dampen the harmonic and subharmonic resonances while increase the resonant radii, which seem to be affected by the medium viscoelasticity to a greater degree rather than the shell properties. This work offers valuable insights into the resonance behaviors of microbubbles oscillating in viscoelastic biological media, greatly contributing to further optimizing their biomedical applications.     

Nonlinear response of ultrasound contrast agent microbubbles:From fundamentals to applications

Modelling and biomedical applications of ultrasound contrast agent(UCA) microbubbles have attracted a great deal of attention. In this review, we summarize a series of researches done in our group, including(i) the development of an all-in-one solution of characterizing coated bubble parameters based on the light scattering technique and flow cytometry;(ii) a novel bubble dynamic model that takes into consideration both nonlinear shell elasticity and viscosity to eliminate the dependences of bubble shell parameters on bubble size;(iii) the evaluation of UCA inertial cavitation threshold and its relationship with shell parameters; and(iv) the investigations of transfection efficiency and the reduction of cytotoxicity in gene delivery facilitated by UCAs excited by ultrasound exposures.     

Excitation threshold for subharmonic generation from contrast microbubbles

Six models of contrast microbubbles are investigated to determine the excitation threshold for subharmonic generation. The models are applied to a commercial contrast agent; its characteristic parameters according to each model are determined using experimentally measured ultrasound attenuation. In contrast to the classical perturbative result, the minimum threshold for subharmonic generation is not always predicted at excitation with twice the resonance frequency; instead it occurs over a range of frequencies from resonance to twice the resonance frequency. The quantitative variation of the threshold with frequency depends on the model and the bubble radius. All models are transformed into a common interfacial rheological form, where the encapsulation is represented by two radius dependent surface properties-effective surface tension and surface dilatational viscosity. Variation of the effective surface tension with radius, specifically having an upper limit (resulting from strain softening or rupture of the encapsulation during expansion), plays a critical role. Without the upper limit, the predicted threshold is extremely large, especially near the resonance frequency. Having a lower limit on surface tension (e.g., zero surface tension in the buckled state) increases the threshold value at twice the resonance frequency, in some cases shifting the minimum threshold toward resonance.     

Optimal design and experimental investigation of surfactant encapsulated microbubbles

The diagnostic capabilities of ultrasound imaging can be improved with contrast-specific nonlinear imaging modalities such as harmonic and subharmonic imaging. The nonlinear response of an encapsulated microbubble in an acoustic field is strongly influenced by the shell viscoelastic properties that are determined by the shell composition and thickness. In this paper, the subharmonic performance of a surfactant encapsulated microbubble was optimized by choosing the appropriate composition of shell material with the aid of theoretical model. To study the effects of viscoelastic properties of microbubble shell materials on the nonlinear scattered response of microbubbles, a theoretical model-modified Herring equation for the oscillation of encapsulated microbubbles in the ultrasound field was employed. Based on this model, a computer aided design system was developed to optimize and analyze the acoustic properties, particularly subharmonic responses, of microbubbles under different shell parameters. Furthermore, surfactant encapsulated microbubbles with different viscoelastic properties were prepared by changing the shell composition. Their shell viscoelastic behavior was measured indirectly as dilational modulus of monolayer film formed with surfactant molecular. Moreover, in vitro quantitative acoustic properties measurements of these microbubbles were carried out to evaluate their subharmonic performance. Both of the theoretical simulation and acoustic measurement showed that the surfactant encapsulated microbubbles with good subharmonic properties could be designed and prepared by adjusting the shell material composition with the guide of the computer aided design system.     

Ambient pressure dependence of the ultra-harmonic response from contrast microbubbles

Sub-harmonic response from ultrasound contrast agent microbubbles has been demonstrated to be an effective modality for noninvasive pressure measurement. In the present study, the dependence of ultra-harmonic response on the ambient overpressure was investigated by both experimental measurements and simulations. In the measurements, the microbubbles were exposed to Gaussian pulses with varied driving frequencies and pulse lengths, at an acoustic pressure of 0.3 MPa. The amplitudes of sub- and ultra-harmonic components were measured when the ambient overpressures varied from 0-25 kPa. At the driving frequency of 1.33 MHz, the ultra-harmonic energy decreased but the sub-harmonic energy increased with the increasing overpressure; while at the driving frequency of 4 MHz, both the sub- and ultra-harmonic components showed the same tendency that the corresponding energy decreased as the overpressure was increased. A 4-MHz Gaussian pulse with 64 cycles could provide an ultra-harmonic response with both good ambient pressure sensitivity and high linearity. Furthermore, the effects of shell parameters of a microbubble on the generation of ultra- and sub-harmonic responses were discussed based on simulations using Marmottant's model. This study suggests that the ultra-harmonic response from contrast microbubbles might be applicable for noninvasive pressure measurement.     

A model for the dynamics of gas bubbles in soft tissue

Understanding the behavior of cavitation bubbles driven by ultrasonic fields is an important problem in biomedical acoustics. Keller-Miksis equation, which can account for the large amplitude oscillations of bubbles, is rederived in this paper and combined with a viscoelastic model to account for the strain-stress relation. The viscoelastic model used in this study is the Voigt model. It is shown that only the viscous damping term in the original equation needs to be modified to account for the effect of elasticity. With experiment determined viscoelastic properties, the effects of elasticity on bubble oscillations are studied. Specifically, the inertial cavitation thresholds are determined using R(max)/R(0), and subharmonic signals from the emission of an oscillating bubble are estimated. The results show that the presence of the elasticity increases the threshold pressure for a bubble to oscillate inertially, and subharmonic signals may only be detectable in certain ranges of radius and pressure amplitude. These results should be easy to verify experimentally, and they may also be useful in cavitation detection and bubble-enhanced imaging.     

Numerical modeling of microbubble backscatter to optimize ultrasound particle image velocimetry imaging: initial studies

We have developed a promising non-invasive ultrasound-based method for performing particle image velocimetry (PIV) in vivo. This method, termed echo PIV, provides multi-component blood velocity data with good (2 ms) temporal resolution. The method takes advantage of the non-linear ultrasound backscatter characteristics of small gas-filled microbubbles (ultrasound contrast) that are seeded into the blood stream. In this study, we use a numerical model to explore potential areas to focus future work in echo PIV.

Ultrasound backscatter from encapsulated microbubbles was modeled using a modified Rayleigh–Plesset equation (Church model, 1995), taking into account the protein/lipid shell layer as a thick, mass-conserving incompressible fluid surrounded by incompressible blood-like fluid. The equation of motion was solved numerically to characterize the fundamental and second harmonic components of the backscattered pressure. Results show a significant advantage in using the second harmonic component for echo PIV, especially for small bubble sizes less than 3 μm in diameter at 2.2 MHz frequency. The effect of the shell thickness ranging from 10 to 500 nm on the vibration amplitude of the bubble was examined and it is shown that the presence of the shell requires mechanical index (MI)>0.2 of incident pressure amplitude to improve bubble detectability. Analysis of the effect of pulse length shows a tradeoff between axial resolution (short pulse length) and bubble detectability (longer pulse length) will most likely be required. The effect of varying MI between 0.1 and 0.6 was also studied at a center frequency of 2.2 MHz and the results indicate that the resonance of the second harmonic is maximized for bubbles with diameter of approximately 2.75 μm. Bubble non-linearities at MI>0.2 induced a resonant frequency shift away from the integer multiple of the incident frequency in the second harmonic backscatter. For a given bubble size, there is a combination of optimal incident frequency and mechanical index range that maximizes the ratio of the second harmonic compared to the fundamental. This resonant frequency decreases with increasing bubble radius. Further, a narrow bandwidth pulse is shown to increase signal strength. Both of these effects may cause conflict with factors governing spatial resolution. Optimization of the incident frequency, microbubble size and mechanical index to enhance bubble detectability will depend on the particular clinical application. These theoretical predictions provide further understanding of the physics behind our echo PIV technique, and should be useful for guiding the design of echo PIV systems.     

Modeling of nonlinear viscous stress in encapsulating shells of lipid-coated contrast agent microbubbles

A general theoretical approach to the development of zero-thickness encapsulation models for contrast microbubbles is proposed. The approach describes a procedure that allows one to recast available rheological laws from the bulk form to a surface form which is used in a modified Rayleigh-Plesset equation governing the radial dynamics of a contrast microbubble. By the use of the proposed procedure, the testing of different rheological laws for encapsulation can be carried out. Challenges of existing shell models for lipid-encapsulated microbubbles, such as the dependence of shell parameters on the initial bubble radius and the “compression-only” behavior, are discussed. Analysis of the rheological behavior of lipid encapsulation is made by using experimental radius-time curves for lipid-coated microbubbles with radii in the range 1.2-2.5 μm. The curves were acquired for a research phospholipid-coated contrast agent insonified with a 20 cycle, 3.0 MHz, 100 kPa acoustic pulse. The fitting of the experimental data by a model which treats the shell as a viscoelastic solid gives the values of the shell surface viscosity increasing from 0.30 × 10⁻⁸ kg/s to 2.63 × 10⁻⁸ kg/s for the range of bubble radii, indicated above. The shell surface elastic modulus increases from 0.054 N/m to 0.37 N/m. It is proposed that this increase may be a result of the lipid coating possessing the properties of both a shear-thinning and a strain-softening material. We hypothesize that these complicated rheological properties do not allow the existing shell models to satisfactorily describe the dynamics of lipid encapsulation. In the existing shell models, the viscous and the elastic shell terms have the linear form which assumes that the viscous and the elastic stresses acting inside the lipid shell are proportional to the shell shear rate and the shell strain, respectively, with constant coefficients of proportionality. The analysis performed in the present paper suggests that a more general, nonlinear theory may be more appropriate. It is shown that the use of the nonlinear theory for shell viscosity allows one to model the “compression-only” behavior. As an example, the results of the simulation for a 2.03 μm radius bubble insonified with a 6 cycle, 1.8 MHz, 100 kPa acoustic pulse are given. These parameters correspond to the acoustic conditions under which the “compression-only” behavior was observed by de Jong et al. [Ultrasound Med. Biol. 33 (2007) 653-656]. It is also shown that the use of the Cross law for the modeling of the shear-thinning behavior of shell viscosity reduces the variance of experimentally estimated values of the shell viscosity and its dependence on the initial bubble radius.     

Harmonic responses and cavitation activity of encapsulated microbubbles coupled with magnetic nanoparticles

Encapsulated microbubbles coupled with magnetic nanoparticles, one kind of hybrid agents that can integrate both ultrasound and magnetic resonance imaging/therapy functions, have attracted increasing interests in both research and clinic communities. However, there is a lack of comprehensive understanding of their dynamic behaviors generated in diagnostic and therapeutic applications. In the present work, a hybrid agent was synthesized by integrating superparamagnetic iron oxide nanoparticles (SPIOs) into albumin-shelled microbubbles (named as SPIO-albumin microbubbles). Then, both the stable and inertial cavitation thresholds of this hybrid agent were measured at varied SPIO concentrations and ultrasound parameters (e.g., frequency, pressure amplitude, and pulse length). The results show that, at a fixed acoustic driving frequency, both the stable and inertial cavitation thresholds of SPIO-albumin microbubble should decrease with the increasing SPIO concentration and acoustic driving pulse length. The inertial cavitation threshold of SPIO-albumin microbubbles also decreases with the raised driving frequency, while the minimum sub- and ultra-harmonic thresholds appear at twice and two thirds resonance frequency, respectively. It is also noticed that both the stable and inertial cavitation thresholds of SonoVue microbubbles are similar to those measured for hybrid microbubbles with a SPIO concentration of 114.7 μg/ml. The current work could provide better understanding on the impact of the integrated SPIOs on the dynamic responses (especially the cavitation activities) of hybrid microbubbles, and suggest the shell composition of hybrid agents should be appropriately designed to improve their clinical diagnostic and therapeutic performances of hybrid microbubble agents.     

A preliminary in vitro assessment of polymer-shelled microbubbles in contrast-enhanced ultrasound imaging

This paper focuses on the use of poly (vinyl alcohol)-shelled microbubbles as a contrast agent in ultrasound medical imaging. The objective was an in vitro assessment of the different working conditions and signal processing methods for the visual detection (especially in small vessels) of such microbubbles, while avoiding their destruction. Polymer-shelled microbubbles have recently been proposed as ultrasound contrast agents with some important advantages. The major drawback is a shell that is less elastic than that of the traditional lipidic microbubbles. Weaker echoes are expected, and their detection at low concentrations may be critical. In vitro experiments were performed with a commercial ultrasound scanner equipped with a dedicated acquisition board. A concentration of 100 bubbles/mm³, excitation pressure amplitudes from 120 kPa to 320 kPa, and a central frequency of 3 MHz or 4.5 MHz were used. Three multi-pulse techniques (i.e., pulse inversion, contrast pulse sequence based on three transmitted signals, and contrast pulse sequence in combination with the chirp pulse) were compared. The results confirmed that these microbubbles produce a weaker ultrasound response than lipidic bubbles with a reduced second-order nonlinear component. Nevertheless, these microbubbles can be detected by the contrast pulse sequence technique, especially when the chirp pulse is adopted. The best value of the contrast-to-tissue ratio was obtained at an excitation pressure amplitude of 230 kPa: although this pressure amplitude is higher than what is typically used for lipidic microbubbles, it does not cause the rupture of the polymeric contrast agent.     

Forced linear oscillations of microbubbles in blood capillaries

A theoretical investigation of the forced linear oscillations of a gas microbubble in a blood capillary, whose radius is comparable in size to the bubble radius is presented. The natural frequency of oscillation, the thermal and viscous damping coefficients, the amplitude resonance, the energy resonance, as well as the average energy absorbed by the system, bubble plus vessel, have been computed for different kinds of gas microbubbles, containing air, octafluropropane, and perflurobutane as a function of the bubble radius and applied frequency. It has been found that the bubble behavior is isothermal at low frequencies and for small bubbles and between isothermal and adiabatic for larger bubbles and higher frequencies, with the viscous damping dominating over the thermal damping. Furthermore, the width of the energy resonance is strongly dependent on the bubble size and the natural frequency of oscillation is affected by the presence of the vessel wall and position of the bubble in the vessel. Therefore, the presence of the blood vessel affects the way in which the bubble absorbs energy from the ultrasonic field. The motivation of this study lies in the possibility of using gas microbubbles as an aid to therapeutic focused ultrasound treatments.     

Microbubble spectroscopy of ultrasound contrast agents

A new optical characterization of the behavior of single ultrasound contrast bubbles is presented. The method consists of insonifying individual bubbles several times successively sweeping the applied frequency, and to record movies of the bubble response up to 25 million frames/s with an ultrahigh speed camera operated in a segmented mode. The method, termed microbubble spectroscopy, enables to reconstruct a resonance curve in a single run. The data is analyzed through a linearized model for coated bubbles. The results confirm the significant influence of the shell on the bubble dynamics: shell elasticity increases the resonance frequency by about 50%, and shell viscosity is responsible for about 70% of the total damping. The obtained value for shell elasticity is in quantative agreement with previously reported values. The shell viscosity increases significantly with the radius, revealing a new nonlinear behavior of the phospholipid coating.