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.     

Modeling of the acoustic response from contrast agent microbubbles near a rigid wall

In ultrasonic targeted imaging, specially designed encapsulated microbubbles are used, which are capable of selectively adhering to the target site in the body. A challenging problem is to distinguish the echoes from such adherent agents from echoes produced by freely circulating agents. In the present paper, an equation of radial oscillation for an encapsulated bubble near a plane rigid wall is derived. The equation is then used to simulate the echo from a layer of contrast agents localized on a wall. The echo spectrum of adherent microbubbles is compared to that of free, randomly distributed microbubbles inside a vessel, in order to examine differences between the acoustic responses of free and adherent agents. It is shown that the fundamental spectral component of adherent bubbles is perceptibly stronger than that of free bubbles. This increase is accounted for by a more coherent summation of echoes from adherent agents and the acoustic interaction between the agents and the wall. For cases tested, the increase of the fundamental component caused by the above two effects is on the order of 8-9 dB. Bubble aggregates, which are observed experimentally to form near a wall due to secondary Bjerknes forces, increase the intensity of the fundamental component only if they are formed by bubbles whose radii are well below the resonant radius. If the formation of aggregates contributes to the growth of the fundamental component, the increase can exceed 17 dB. Statistical analysis for the comparison between adhering and free bubbles, performed over random space bubble distributions, gives p-values much smaller than 0.05.     

High-speed observation of cavitation bubble clouds near a tissue boundary in high-intensity focused ultrasound fields

Cavitation bubble clouds generated near a tissue boundary by high-intensity focused ultrasound (HIFU) were studied using high-speed photography. In all, 171 image series were captured during the initial 100 ms of continuous HIFU exposure, which showed that cavitation bubble clouds at the tissue boundary organized into two structures - “cone-shape bubble cloud structure” recorded in 146 image series and “crown-shape bubble cloud structure” recorded in 18 image series. The remaining 7 image series showed the interchanging of these two structures. It was found that when cavitation bubbles first appeared at the tissue boundary, they developed to cone-shape bubble cloud. The cone-shape bubble cloud structure was characterized by a nearly fixed tip in front of the tissue boundary. When the cavitation bubbles initially appeared away from the tissue boundary they evolved into a crown-shape bubble cloud. Deformation of tissue boundary was shown in all the recorded image series.     

Scaling of the viscoelastic shell properties of phospholipid encapsulated microbubbles with ultrasound frequency

Phospholipid encapsulated microbubbles are widely employed as clinical diagnostic ultrasound contrast agents in the 1–5 MHz range, and are increasingly employed at higher ultrasound transmit frequencies. The stiffness and viscosity of the encapsulating “shells” have been shown to play a central role in determining both the linear and nonlinear response of microbubbles to ultrasound. At lower frequencies, recent studies have suggested that shell properties can be frequency dependent. At present, there is only limited knowledge of how the viscoelastic properties of phospholipid shells scale at higher frequencies. In this study, four batches of in-house phospholipid encapsulated microbubbles were fabricated with decreasing volume-weighted mean diameters of 3.20, 2.07, 1.82 and 1.61 μm. Attenuation experiments were conducted in order to assess the frequency-dependent response of each batch, resulting in resonant peaks in response at 4.2, 8.9, 12.6 and 19.5 MHz, respectively. With knowledge of the size measurements, the attenuation spectra were then fitted with a standard linearized bubble model in order to estimate the microbubble shell stiffness S_p and shell viscosity S_f, resulting in a slight increase in S_p (1.53–1.76 N/m) and a substantial decrease in S_f (0.29 × 10⁻⁶–0.08 × 10⁻⁶ kg/s) with increasing frequency. These results performed on a single phospholipid agent show that frequency dependent shell properties persist at high frequencies (up to 19.5 MHz).     

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.     

Acoustic scattering of a high-order Bessel beam by an elastic sphere

The exact analytical solution for the scattering of a generalized (or “hollow”) acoustic Bessel beam in water by an elastic sphere centered on the beam is presented. The far-field acoustic scattering field is expressed as a partial wave series involving the scattering angle relative to the beam axis and the half-conical angle of the wave vector components of the generalized Bessel beam. The sphere is assumed to have isotropic elastic material properties so that the nth partial wave amplitude for plane wave scattering is proportional to a known partial-wave coefficient. The transverse acoustic scattering field is investigated versus the dimensionless parameter ka(k is the wave vector, a radius of the sphere) as well as the polar angle θ for a specific dimensionless frequency and half-cone angle β. For higher-order generalized beams, the acoustic scattering vanishes in the backward (θ = π) and forward (θ = 0) directions along the beam axis. Moreover it is possible to suppress the excitation of certain resonances of an elastic sphere by appropriate selection of the generalized Bessel beam parameters.     

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

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

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.     

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

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

Anisotropic elastic properties of microtubules

The present study aims at providing insight into the acceleration mechanism of a bubble chain rising in shear-thinning viscoelastic fluids. The experimental investigation by the Particle Image Velocimetry (PIV), birefringence visualisation and rheological simulation shows that two aspects are central to bubble interactions in such media: the stress creation by the passage of bubbles, and their relaxation due to the fluids memory forming an evanescent corridor of reduced viscosity. Interactions between bubbles were taken into account mainly through a linear superposition of the stress evolution behind each bubble. An analytical approach together with the rheological consideration was developed to compute the rise velocity of a bubble chain in function of the injection period and bubble volume. The model predictions compare satisfactorily with the experimental investigation.     

A performance analysis of echographic ultrasonic techniques for non-invasive temperature estimation in hyperthermia range using phantoms with scatterers

Optimization of efficiency in hyperthermia requires a precise and non-invasive estimation of internal distribution of temperature. Although there are several research trends for ultrasonic temperature estimation, efficient equipments for its use in the clinical practice are not still available. The main objective of this work was to research about the limitations and potential improvements of previously reported signal processing options in order to identify research efforts to facilitate their future clinical use as a thermal estimator.In this document, we have a critical analysis of potential performance of previous ultrasonic research trends for temperature estimation inside materials, using different processing techniques proposed in frequency, time and phase domains. It was carried out in phantom with scatterers, assessing at their specific applicability, linearity and limitations in hyperthermia range. Three complementary evaluation indexes: technique robustness, Mat-lab processing time and temperature resolution, with specific application protocols, were defined and employed for a comparative quantification of the behavior of the techniques. The average increment per °C and mm was identified for each technique (3 KHz/°C in the frequency analysis, 0.02 rad/°C in the phase domain, while increments in the time domain of only 1.6 ns/°C were found). Their linearity with temperature rising was measured using linear and quadratic regressions and they were correlated with the obtained data.New improvements in time and frequency signal processing in order to reveal the potential thermal and spatial resolutions of these techniques are proposed and their subsequent improved estimation results are shown for simulated and measured A-scans registers. As an example of these processing novelties, an excellent potential resolution of 0.12 °C into hyperthermia range, with near-to-linear frequency dependence, could be achieved.Specifically defined “numerical” and physical multi-scatter phantoms are described, which mimic ultrasound velocity in tissues of about 1560 m/s @ 35 °C and have a quasi-uniform internal scattering structure designed to assure standard signal patterns adequate for processing comparisons in the same time and sound velocity conditions for all the techniques analyzed, and to obtain easily repeatable multi-pulse echo-patterns.A perfect lineal dependence (100% of correlation coefficient) between the unitary average increment measured by each technique and temperature rising was observed while working with simulated A-scan registers, where all the parameters are under an accurate control. Nevertheless a very small quadratic tendency appeared in the results obtained from experimental echo registers, which are more similar to a real tissues case. It would be an interesting future work to analyze the behavior of these techniques in real tissues in order to confirm or reject this light quadratic tendency.Finally, new methods were detailed and applied in order to precisely quantify the advantages of each estimation technique; their respective intrinsic limitations were also underlined.     

Subharmonic behavior of phospholipid-coated ultrasound contrast agent microbubbles

Coated microbubbles, unlike tissue are able to scatter sound subharmonically. Therefore, the subharmonic behavior of coated microbubbles can be used to enhance the contrast in ultrasound contrast imaging. Theoretically, a threshold amplitude of the driving pressure can be calculated above which subharmonic oscillations of microbubbles are initiated. Interestingly, earlier experimental studies on coated microbubbles demonstrated that the threshold for these bubbles is much lower than predicted by the traditional linear viscoelastic shell models. This paper presents an experimental study on the subharmonic behavior of differently sized individual phospholipid coated microbubbles. The radial subharmonic response of the microbubbles was recorded with the Brandaris ultra high-speed camera as a function of both the amplitude and the frequency of the driving pulse. Threshold pressures for subharmonic generation as low as 5 kPa were found near a driving frequency equal to twice the resonance frequency of the bubble. An explanation for this low threshold pressure is provided by the shell buckling model proposed by Marmottant et al. [J. Acoust. Soc. Am. 118, 3499-3505 (2005)]. It is shown that the change in the elasticity of the bubble shell as a function of bubble radius as proposed in this model, enhances the subharmonic behavior of the microbubbles.     

Effect of surfactant addition on removal of microbubbles using ultrasound

It is difficult to control the bubble in a liquid by the external operation, because the behavior of the bubble is controlled in buoyancy and flow of liquid. On the other hand, microbubbles, whose diameter is several decades μm, stably disperse in static liquid because of their small buoyancy and electrical repulsion. When an ultrasound, whose frequency was 2.4 MHz, was irradiated, the milky white microbubbles suspended solution became rapidly clear. In this study, the effects of surfactant addition on the removal of microbubbles from a liquid in an ultrasonic field were investigated. The efficiency of removal of microbubbles decreased with surfactant addition. Surfactant type influenced the size of agglomerated microbubbles, and the efficiency of removal of microbubbles changed. The surface of microbubble was modified by surfactant adsorption, and the steric inhibition influenced the removal of microbubbles.     

Nanometric laser trapping of microbubbles based on nanostructured substrates

Laser trapping near the surface of a nanostructured substrate is demonstrated. Stable microbubbles with radii of 1-20 μm have been created and manipulated with sub-micron precision by a focused laser beam in an immersion oil covering arrays of pairs of gold nanopillars deposited on a glass substrate. The threshold for bubble creation and trapping characteristics depended on near-field coupling of nanopillars. The nanometric laser tweezers showed giant trapping efficiency of Q ∼ 50 for the trapped microbubbles.     

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.     

"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.     

Influence of nesting shell size on brightness longevity and resistance to ultrasound-induced dissolution during enhanced B-mode contrast imaging

This study aims to bridge the gap between transport mechanisms of an improved ultrasound contrast agent (UCA) and its resulting behavior in a clinical imaging study. Phospholipid-shelled microbubbles nested within the aqueous core of a polymer microcapsule are examined for their use and feasibility as an improved UCA. The nested formulation provides contrast comparable to traditional formulations, specifically an SF₆ microbubble coated by a DSPC PEG-3000 monolayer, with the advantage that contrast persists at least nine times longer in a mock clinical, in vitro setting. The effectiveness of the sample was measured using a contrast ratio in units of decibels (dB) which compares the brightness of the nested microbubbles to a reference value of a phantom tissue mimic. During a 40 min imaging study, six nesting formulations with average outer capsule diameters of 1.95, 2.53, 5.55, 9.95, 14.95, and 20.51 μm reached final contrast ratio values of 0.25, 2.35, 3.68, 4.51, 5.93, and 8.00 dB, respectively. The starting contrast ratio in each case was approximately 8 dB and accounts for the brightness attributed to the nesting shell. As compared with empty microcapsules (no microbubbles nested within), enhancement of the initial contrast ratio increased systematically with decreasing microcapsule size. The time required to reach a steady state in the temporal contrast ratio profile also varied with microcapsule diameter and was found to be 420 s for each of the four smallest shell diameters and 210 s and 150 s, respectively, for the largest two shell diameters. All nested formulations were longer-lived and gave higher final contrast ratios than a control sample comprising un-nested, but otherwise equivalent, microbubbles. Specifically, the contrast ratio of the un-nested microbubbles decreased to a negative value after 4 min of continuous ultrasound exposure with complete disappearance of the microbubbles after 15 min whereas all nested formulations maintained positive contrast ratio values for the duration of the 40 min trial. The results are consistent with two distinct stages of gas transport: in the first stage, passive diffusion occurs under ambient conditions across the microbubble monolayer within the first few minutes after formulation until the aqueous interior of the microcapsule is saturated with gas; in the second stage ultrasound drives additional gas dissolution even further due to pressure modulation. It is important to understand the chemistry and transport mechanisms of this contrast agent under the influence of ultrasound to attain better perspicacity for enhanced applications in imaging. Results from this study will facilitate future preclinical studies and clinical applications of nested microbubbles for therapeutic and diagnostic imaging.     

调频超声脉冲驱动微气泡运动偏移的研究

提出调频超声辐射力技术驱动微泡群,以加强微泡的吸附效率.基于改进的RP方程及粒子轨迹方程研究了微泡群整体的运动位移与调频信号的中心频率、调频范围、信号声压,以及微泡半径分布关系.研究结果表明调频信号在驱动半径具有宽泛分布的气泡群,以及半径分布远离谐振半径的气泡群时,作用效果好于传统正弦波信号.例如中心频率1 MHz、调频范围0.75 MHz的调频脉冲作用高斯分布（平均半径3.5 μm、均方差为1）的微泡群200 μs,可比同等声压的正弦波多约12%的微气泡产生位移30 μm. 关键词： 超声辐射力 调频波 高斯分布     

Physical forces exerted by microbubbles on a surface in a traveling wave field

The effect of a wave with a varying traveling component on the bubble activity as well as the physical force generated by microbubbles on a surface has been studied. The acoustic emission from a collection of bubbles is measured in a 928 kHz sound field. Particle removal tests on a surface, which actually measures the applied physical force by the bubbles on that surface, indicate a very strong dependence on the angle of incidence. In other words, when the traveling wave component is maximized, the average physical force applied by microbubbles reaches a maximum. Almost complete particle removal for 78 nm silica particles was obtained for a traveling wave, while particle removal efficiency was reduced to only a few percent when a standing wave was applied. This increase in particle removal for a traveling wave is probably caused by a decrease in bubble trapping at nodes and antinodes in a standing wave field.     

微泡对高强度聚焦超声焦域影响的研究

微泡对高强度聚焦超声(HIFU)治疗具有增效作用,而HIFU治疗中不同声学条件下微泡对HIFU治疗焦域的影响尚不清楚。本文基于声传播方程、Yang-Church气泡运动方程、生物热传导方程、时域有限差分法(FDTD)、龙格-库塔(RK)法数值仿真研究输入功率、激励频率和气泡初始半径对HIFU在含气泡体模中形成焦域的影响,并利用含Sono Vue造影剂的仿组织体模研究进行实验验证。结果表明,增大输入功率、气泡初始半径和升高激励频率均可增大焦域,随着输入功率的增大,焦域形状可能发生变化,而随着激励频率升高和气泡初始半径的增大,焦域会向远离换能器的方向移动。