Super-shear evanescent waves for non-contact elastography of soft tissues

Dual-channel air-pulse optical coherence elastography for frequency-response analysis

Microliter air-pulse optical coherence elastography (OCE) has recently been proposed for the characterization of soft-tissue biomechanics using transient, sub-nanometer to micrometer-scale natural frequency oscillations. However, previous studies have not been able to provide real-time air-pulse monitoring during OCE natural frequency measurement, which could lead to inaccurate measurement results due to the unknown excitation spectrum. To address this issue, we introduce a dual-channel air-pulse OCE method, with one channel stimulating the sample and the other being simultaneously measured with a pressure sensor. This allows for more accurate natural frequency characterization using the frequency response function, as proven by a comprehensive comparison under different conditions with a diverse range of excitation spectra (from broad to narrow, clean to noisy) as well as a diverse set of sample response spectra. We also demonstrate the capability of the frequency-response analysis in distinguishing samples with different stiffness levels: the dominant natural frequencies increased with agar concentrations (181-359 Hz, concentrations: 1-2%, and maximum displacements: 0.12-0.47 µm) and intraocular pressures (IOPs) for the silicone cornea (333-412 Hz, IOP: 5-40 mmHg, and maximum displacements: 0.41-0.52 µm) under a 200 Pa stimulation pressure. These frequencies remained consistent across different air-pulse durations (3 ms to 35 ms). The dual-channel OCE approach that uses transient, low-pressure stimulation and high-resolution imaging holds the potential to advance our understanding of sample frequency responses, especially when investigating delicate tissues such as the human cornea in vivo.

Noninvasive identification of directionally-dependent elastic properties of soft tissues using full-field optical data

This paper introduces an innovative approach for elastic property characterization of soft tissues, having directional-dependent material behavior, via their vibration response measurement and interpretation. The full-field time-dependent surface displacements as a result of externally excited soft tissues are collected through digital image correlation (DIC). A developed analytical model, capturing the low-amplitude vibration behavior of anisotropic layered human skin with the incorporation of the influence of subcutaneous elasticity and inertia, is employed to accurately predict its resonant frequencies and pertaining displacement field images. An efficient solution approach for the model is implemented into an inverse algorithm to rapidly characterize the anisotropic elastic properties based on importing the vibration characteristics. To show the merit of the approach, a 3-D finite element (FE) simulation model was used to generate full-field data, detected and matched with a set of specific vibration modes via modal assurance criterion (MAC). The validity of the model implemented into the inverse characterization algorithm is demonstrated through a comparison of predicted vibration frequencies and mode-shapes simulated via the 3-D FE model for different cases with anisotropic elastic properties in different layers of the skin. It is shown that modes are influenced differently when anisotropic properties are introduced to the model. Thus, the established inverse characterization algorithm is capable of rapidly predicting the elastic material properties of anisotropic soft sheets with adequate accuracy.

Supershear Rayleigh wave imaging for quantitative assessment of biomechanical properties of brain using air-coupled optical coherence elastography

Recently, supershear Rayleigh waves (SRWs) have been proposed to characterize the biomechanical properties of soft tissues. The SRWs propagate along the surface of the medium, unlike surface Rayleigh waves, SRWs propagate faster than bulk shear waves. However, their behavior and application in biological tissues is still elusive. In brain tissue elastography, shear waves combined with magnetic resonance elastography or ultrasound elastography are generally used to quantify the shear modulus, but high spatial resolution elasticity assessment in 10 μm scale is still improving. Here, we develop an air-coupled ultrasonic transducer for noncontact excitation of SRWs and Rayleigh waves in brain tissue, use optical coherent elastography (OCE) to detect, and reconstruct the SRW propagation process; in combing with a derived theoretical model of SRWs on a free boundary surface, we quantify the shear modulus of brain tissue with high spatial resolution. We first complete validation experiments using a homogeneous isotropic agar phantom, and the experimental results clearly show the SRW is 1.9649 times faster than the bulk shear waves. Furthermore, the propagation velocity of SRWs in both the frontal and parietal lobe regions of the brain is all 1.87 times faster than the bulk shear wave velocity. Finally, we evaluated the anisotropy in different brain regions, and the medulla oblongata region had the highest anisotropy index. Our study shows that the OCE system using the SRW model is a new potential approach for high-resolution assessment of the biomechanical properties of brain tissue.

Ultra-wideband optical coherence elastography from acoustic to ultrasonic frequencies

Visualizing viscoelastic waves in materials and tissues through noninvasive imaging is valuable for analyzing their mechanical properties and detecting internal anomalies. However, traditional elastography techniques have been limited by a maximum wave frequency below 1-10 kHz, which hampers temporal and spatial resolution. Here, we introduce an optical coherence elastography technique that overcomes the limitation by extending the frequency range to MHz. Our system can measure the stiffness of hard materials including bones and extract viscoelastic shear moduli for polymers and hydrogels in conventionally inaccessible ranges between 100 Hz and 1 MHz. The dispersion of Rayleigh surface waves across the ultrawide band allowed us to profile depth-dependent shear modulus in cartilages ex vivo and human skin in vivo with sub-mm anatomical resolution. This technique holds immense potential as a noninvasive measurement tool for material sciences, tissue engineering, and medical diagnostics.

Recent advances in optical elastography and emerging opportunities in the basic sciences and translational medicine [Invited]

Optical elastography offers a rich body of imaging capabilities that can serve as a bridge between organ-level medical elastography and single-molecule biophysics. We review the methodologies and recent developments in optical coherence elastography, Brillouin microscopy, optical microrheology, and photoacoustic elastography. With an outlook toward maximizing the basic science and translational clinical impact of optical elastography technologies, we discuss potential ways that these techniques can integrate not only with each other, but also with supporting technologies and capabilities in other biomedical fields. By embracing cross-modality and cross-disciplinary interactions with these parallel fields, optical elastography can greatly increase its potential to drive new discoveries in the biomedical sciences as well as the development of novel biomechanics-based clinical diagnostics and therapeutics.

Supershear surface waves reveal prestress and anisotropy of soft materials

Surface waves play important roles in many fundamental and applied areas from seismic detection to material characterizations. Supershear surface waves with propagation speeds greater than bulk shear waves have recently been reported, but their properties are not well understood. Here we describe theoretical and experimental results on supershear surface waves in rubbery materials. We find that supershear surface waves can be supported in viscoelastic materials with no restriction on the shear quality factor. Interestingly, the effect of prestress on the speed of the supershear surface wave is opposite to that of the Rayleigh surface wave. Furthermore, anisotropy of material affects the supershear wave much more strongly than the Rayleigh surface wave. We offer heuristic interpretation as well as theoretical verification of our experimental observations. Our work points to the potential applications of supershear waves for characterizing the bulk mechanical properties of soft solid from the free surface.

Spatial resolution in optical coherence elastography of bounded media

Dynamic optical coherence elastography (OCE) tracks mechanical wave propagation in the subsurface region of tissue to image its shear modulus. For bulk shear waves, the lateral resolution of the reconstructed modulus map (i.e., elastographic resolution) can approach that of optical coherence tomography (OCT), typically a few tens of microns. Here we perform comprehensive numerical simulations and acoustic micro-tapping OCE experiments to show that for the typical situation of guided wave propagation in bounded media, such as cornea, the elastographic resolution cannot reach the OCT resolution and is mainly defined by the thickness of the bounded tissue layer. We considered the excitation of both broadband and quasi-harmonic guided waves in a bounded, isotropic medium. Leveraging the properties of broadband pulses, a robust method for modulus reconstruction with minimum artifacts at interfaces is demonstrated. In contrast, tissue bounding creates large instabilities in the phase of harmonic waves, leading to serious artifacts in modulus reconstructions.

In vivo stiffness measurement of epidermis, dermis, and hypodermis using broadband Rayleigh-wave optical coherence elastography

Traveling-wave optical coherence elastography (OCE) is a promising technique to measure the stiffness of biological tissues. While OCE has been applied to relatively homogeneous samples, tissues with significantly varying elasticity through depth pose a challenge, requiring depth-resolved measurement with sufficient resolution and accuracy. Here, we develop a broadband Rayleigh-wave OCE technique capable of measuring the elastic moduli of the 3 major skin layers (epidermis, dermis, and hypodermis) reliably by analyzing the dispersion of leaky Rayleigh surface waves over a wide frequency range of 0.1-10 kHz. We show that a previously unexplored, high frequency range of 4-10 kHz is critical to resolve the thin epidermis, while a low frequency range of 0.2-1 kHz is adequate to probe the dermis and deeper hypodermis. We develop a dual bilayer-based inverse model to determine the elastic moduli in all 3 layers and verify its high accuracy with finite element analysis and skin-mimicking phantoms. Finally, the technique is applied to measure the forearm skin of healthy volunteers. The Young's modulus of the epidermis (including the stratum corneum) is measured to be ∼ 4 MPa at 4-10 kHz, whereas Young's moduli of the dermis and hypodermis are about 40 and 15 kPa, respectively, at 0.2-1 kHz. Besides dermatologic applications, this method may be useful for the mechanical analysis of various other layered tissues with sub-mm depth resolution. STATEMENT OF SIGNIFICANCE: To our knowledge, this is the first study that resolves the stiffness of the thin epidermis from the dermis and hypodermis, made possible by using high-frequency (4 - 10 kHz) elastic waves and optical coherence elastography. Beyond the skin, this technique may be useful for mechanical characterizations of various layered biomaterials and tissues.

Simplified Green's function for surface waves in quasi-incompressible elastic plates with application to elastography

Surface wave elastography is a growing method to estimate the elasticity in soft solids. It is particularly useful in the case of agrifoods like meat, cheese, or fruits because it does not require major infrastructure or large equipment and could be developed in portable devices. However, estimating the shear elastic properties from surface wave measurements is not straightforward. The shear wavelength in those materials is cm sized for the excitation frequencies usually employed in elastography (∼10² Hz), and the size of samples is comparable to it. Thus, the surface wave speed is frequency dependent with no direct relation to the shear wave speed. In this work we propose a simplified Green's function for soft solid elastic plates which allows to retrieve the shear elasticity from near field measurements. The model is compared with experimental results obtained in agar-gelatin phantoms and food samples (cheese and bovine liver). The results show a good overall agreement although improvements can be achieved by incorporating diffraction and viscosity to the model.

Spatial Assessment of Heterogeneous Tissue Natural Frequency Using Micro-Force Optical Coherence Elastography

Analysis of corneal tissue natural frequency was recently proposed as a biomarker for corneal biomechanics and has been performed using high-resolution optical coherence tomography (OCT)-based elastography (OCE). However, it remains unknown whether natural frequency analysis can resolve local variations in tissue structure. We measured heterogeneous samples to evaluate the correspondence between natural frequency distributions and regional structural variations. Sub-micrometer sample oscillations were induced point-wise by microliter air pulses (60–85 Pa, 3 ms) and detected correspondingly at each point using a 1,300 nm spectral domain common path OCT system with 0.44 nm phase detection sensitivity. The resulting oscillation frequency features were analyzed via fast Fourier transform and natural frequency was characterized using a single degree of freedom (SDOF) model. Oscillation features at each measurement point showed a complex frequency response with multiple frequency components that corresponded with global structural features; while the variation of frequency magnitude at each location reflected the local sample features. Silicone blocks (255.1 ± 11.0 Hz and 249.0 ± 4.6 Hz) embedded in an agar base (355.6 ± 0.8 Hz and 361.3 ± 5.5 Hz) were clearly distinguishable by natural frequency. In a beef shank sample, central fat and connective tissues had lower natural frequencies (91.7 ± 58.2 Hz) than muscle tissue (left side: 252.6 ± 52.3 Hz; right side: 161.5 ± 35.8 Hz). As a first step, we have shown the possibility of natural frequency OCE methods to characterize global and local features of heterogeneous samples. This method can provide additional information on corneal properties, complementary to current clinical biomechanical assessments, and could become a useful tool for clinical detection of ocular disease and evaluation of medical or surgical treatment outcomes.

All-optical noncontact phase-domain photoacoustic elastography

Mechanical properties such as elasticity are important indicators of tissue functions that can be used for clinical diagnosis and disease monitoring. However, most current elastography techniques are limited in their ability to distinguish localized microstructural mechanical variations due to employing elastic wave velocity measurement. In addition, their contact-based measurement manner is not favored and may even be prohibited in many applications. In this Letter, we propose all-optical noncontact phase-domain photoacoustic elastography (NPD-PAE), leveraging the temporal response characteristics of laser-induced thermoelastic displacement using optical interferometric detection to calculate the elastic modulus. The all-optical pump-probe method allows the capture of the initial displacement profiles generated at the origin, thus enabling the extraction of in situ elasticity. The feasibility of the method was verified using a tissue-mimicking phantom. The capability to map the mechanical contrast was demonstrated on an ex vivo biological tissue. NPD-PAE opens a new avenue for development of a noncontact elastography technique, holding great potential in the biomedical field and materials science.

Surface wave elastography is a reliable method to correlate muscle elasticity, torque, and electromyography activity level

The shear elastic modulus is one of the most important parameters to characterize the mechanical behavior of soft tissues. In biomechanics, ultrasound elastography is the gold standard for measuring and mapping it locally in skeletal muscle in vivo. However, their applications are limited to the laboratory or clinic. Thus, low-frequency elastography methods have recently emerged as a novel alternative to ultrasound elastography. Avoiding the use of high frequencies, these methods allow obtaining a mean value of bulk shear elasticity. However, they are frequently susceptible to diffraction, guided waves, and near field effects, which introduces biases in the estimates. The goal of this work is to test the performance of the non-ultrasound surface wave elastography (NU-SWE), which is portable and is based on new algorithms designed to correct the incidence of such effects. Thus, we show its first application to muscle biomechanics. We performed two experiments to assess the relationships of muscle shear elasticity versus joint torque (experiment 1) and the electromyographic activity level (experiment 2). Our results were comparable regarding previous works using the reference ultrasonic methods. Thus, the NU-SWE showed its potentiality to get wide the biomechanical applications of elastography in many areas of health and sports sciences.

Characterization of natural frequencies from nanoscale tissue oscillations using dynamic optical coherence elastography

We demonstrate the use of OCT-based elastography for soft-tissue characterization using natural frequency oscillations. Sub-micrometer to sub-nanometer oscillations were induced in tissue phantoms and human cornea in vivo by perpendicular air-pulse stimulation and observed by common-path OCT imaging (sensitivity: 0.24 nm). Natural frequency and damping ratio were acquired in temporal and frequency domains using a single degree of freedom method. The dominant natural frequency was constant for different stimulation pressures (4-32 Pa) and measured distances (0.3-5.3 mm), and decreased as the sample thickness increased. The dominant natural frequencies of 0.75-2% agar phantoms were 127-774 Hz (mean coefficient of variation [CV]: 0.9%), and correlated with the square root of Young's moduli (16.5-117.8 kPa, mean CV: 5.8%). These preliminary studies show repeatable in vivo corneal natural frequency measurements (259 Hz, CV: 1.9%). This novel OCE approach can distinguish tissues and materials with different mechanical properties using the small-amplitude tissue oscillation features, and is suitable for characterizing delicate tissues in vivo such as the eye.