Photoacoustic monopole radiation in one, two, and three dimensions

Air-coupled ultrasound detection using capillary-based optical ring resonators

We experimentally demonstrate and theoretically analyze high Q-factor (~107) capillary-based optical ring resonators for non-contact detection of air-coupled ultrasound. Noise equivalent pressures in air as low as 215 mPa/√Hz and 41 mPa/√Hz at 50 kHz and 800 kHz in air, respectively, are achieved. Furthermore, non-contact detection of air-coupled photoacoustic pulses optically generated from a 200 nm thick Chromium film is demonstrated. The interaction of an acoustic pulse and the mechanical mode of the ring resonator is also studied. Significant improvement in detection bandwidth is demonstrated by encapsulating the ring resonator in a damping medium. Our work will enable compact and sensitive ultrasound detection in many applications, such as air-coupled non-destructive ultrasound testing, photoacoustic imaging, and remote sensing. It will also provide a model system for fundamental study of the mechanical modes in the ring resonator.

Photoacoustic resonance by spatial filtering of focused ultrasound transducers

The effect of a spherically focused ultrasound (US) transducer (TD) on photoacoustic (PA) measurements is analytically investigated using the concept of a virtual point detector. The derived analytical results indicate that the limited numerical aperture (NA) of the PA detector takes on the role of spatial filtering of the induced PA waves, which leads to the occurrence of a peak frequency in the PA spectrum. The mathematical description of this phenomenon is similar to the result of resonance peaks of light propagation in dielectrics. This PA resonance peak depends on the NA of the US TD and the absorption coefficient of the PA source. Experimentally measured PA spectra from an ink solution using a frequency-domain PA system verify the PA resonance effect. Finally, we discuss the possibility that previously reported PA resonances interpreted as PA source viscosity might be actually caused by the US TD's spatial filtering.

Detection, numerical simulation and approximate inversion of optoacoustic signals generated in multi-layered PVA hydrogel based tissue phantoms

Optoacoustic (OA) measurements can not only be used for imaging purposes but as a more general tool to "sense" physical characteristics of biological tissue, such as geometric features and intrinsic optical properties. In order to pave the way for a systematic model-guided analysis of complex objects we devised numerical simulations in accordance with the experimental measurements. We validate our computational approach with experimental results observed for layered polyvinyl alcohol hydrogel samples, using melanin as the absorbing agent. Experimentally, we characterize the acoustic signal observed by a piezoelectric detector in the acoustic far-field in backward mode and we discuss the implication of acoustic diffraction on our measurements. We further attempt an inversion of an OA signal in the far-field approximation.

Preclinical detection of liver fibrosis using dual-modality photoacoustic/ultrasound system

Liver fibrosis is a major cause for increasing mortality worldwide. Preclinical research using animal models is required for the discovery of new anti-fibrotic therapies, but currently relies on endpoint liver histology. In this study, we investigated a cost-effective and portable photoacoustic/ultrasound (PA/US) imaging system as a potential non-invasive alternative. Fibrosis was induced in mice using CCl₄ followed by liver imaging and histological analysis. Imaging showed significantly increased PA features with higher frequency signals in fibrotic livers versus healthy livers. This corresponds to more heterogeneous liver structure resulting from collagen deposition and angiogenesis. Importantly, PA response and its frequency were highly correlated with histological parameters. These results demonstrate the preclinical feasibility of the PA imaging approach and applicability of dual PA/US system.

Quantitative photoacoustic image reconstruction improves accuracy in deep tissue structures

Photoacoustic imaging (PAI) is emerging as a potentially powerful imaging tool with multiple applications. Image reconstruction for PAI has been relatively limited because of limited or no modeling of light delivery to deep tissues. This work demonstrates a numerical approach to quantitative photoacoustic image reconstruction that minimizes depth and spectrally derived artifacts. We present the first time-domain quantitative photoacoustic image reconstruction algorithm that models optical sources through acoustic data to create quantitative images of absorption coefficients. We demonstrate quantitative accuracy of less than 5% error in large 3 cm diameter 2D geometries with multiple targets and within 22% error in the largest size quantitative photoacoustic studies to date (6cm diameter). We extend the algorithm to spectral data, reconstructing 6 varying chromophores to within 17% of the true values. This quantitiative PA tomography method was able to improve considerably on filtered-back projection from the standpoint of image quality, absolute, and relative quantification in all our simulation geometries. We characterize the effects of time step size, initial guess, and source configuration on final accuracy. This work could help to generate accurate quantitative images from both endogenous absorbers and exogenous photoacoustic dyes in both preclinical and clinical work, thereby increasing the information content obtained especially from deep-tissue photoacoustic imaging studies.

Experimental observation of acoustic emissions generated by a pulsed proton beam from a hospital-based clinical cyclotron

PURPOSE

To measure the acoustic signal generated by a pulsed proton spill from a hospital-based clinical cyclotron.

METHODS

An electronic function generator modulated the IBA C230 isochronous cyclotron to create a pulsed proton beam. The acoustic emissions generated by the proton beam were measured in water using a hydrophone. The acoustic measurements were repeated with increasing proton current and increasing distance between detector and beam.

RESULTS

The cyclotron generated proton spills with rise times of 18 μs and a maximum measured instantaneous proton current of 790 nA. Acoustic emissions generated by the proton energy deposition were measured to be on the order of mPa. The origin of the acoustic wave was identified as the proton beam based on the correlation between acoustic emission arrival time and distance between the hydrophone and proton beam. The acoustic frequency spectrum peaked at 10 kHz, and the acoustic pressure amplitude increased monotonically with increasing proton current.

CONCLUSIONS

The authors report the first observation of acoustic emissions generated by a proton beam from a hospital-based clinical cyclotron. When modulated by an electronic function generator, the cyclotron is capable of creating proton spills with fast rise times (18 μs) and high instantaneous currents (790 nA). Measurements of the proton-generated acoustic emissions in a clinical setting may provide a method for in vivo proton range verification and patient monitoring.

Characterizing cellular morphology by photoacoustic spectrum analysis with an ultra-broadband optical ultrasonic detector

Photoacoustic spectrum analysis (PASA) has been demonstrated as a new method for quantitative tissue imaging and characterization. The ability of PASA in evaluating micro-size tissue features was limited by the bandwidth of detectors for photoacoustic (PA) signal acquisition. We improve upon such a limit, and report on developments of PASA facilitated by an optical ultrasonic detector based on micro-ring resonator. The detector's broad and flat frequency response significantly improves the performance of PASA and extents its characterization capability from the tissue level to cellular level. The performance of the system in characterizing cellular level (a few microns) stochastic objects was first shown via a study on size-controlled optically absorbing phantoms. As a further demonstration of PASA's potential clinical application, it was employed to characterize the morphological changes of red blood cells (RBCs) from a biconcave shape to a spherical shape as a result of aging. This work demonstrates that PASA equipped with the micro-ring ultrasonic detectors is an effective technique in characterizing cellular-level micro-features of biological samples.

Cost-effective imaging of optoacoustic pressure, ultrasonic scattering, and optical diffuse reflectance with improved resolution and speed

The idea of a method of cost-effective upgrades from an acoustic resolution photoacoustic microscope to a triple-modality imaging system is validated using phantoms. The newly developed experimental setup is based on a diode pumped solid state laser coupled to a fiber bundle with a spherically focused polyvinylidene fluoride detector integrated into the center of a ring shaped optical illuminator. Each laser pulse illuminating the sample performs two functions. While the photons absorbed by the sample provide a measurable optoacoustic (OA) signal, the photons absorbed by the detector provide the measurable diffuse reflectometry (DR) signal from the sample and the probing ultrasonic (US) pulse. At a 3 mm imaging depth, the axial resolution of the OA/US modalities is 38 μm/26 μm, while the lateral resolution of the DR/OA/US modalities is 3.5 mm/50 μm/35 μm. The maximum acquisition rate of the trimodal DR/OA/US A-scans is 2 kHz.

High resolution ultrasound and photoacoustic imaging of single cells

High resolution ultrasound and photoacoustic images of stained neutrophils, lymphocytes and monocytes from a blood smear were acquired using a combined acoustic/photoacoustic microscope. Photoacoustic images were created using a pulsed 532 nm laser that was coupled to a single mode fiber to produce output wavelengths from 532 nm to 620 nm via stimulated Raman scattering. The excitation wavelength was selected using optical filters and focused onto the sample using a 20× objective. A 1000 MHz transducer was co-aligned with the laser spot and used for ultrasound and photoacoustic images, enabling micrometer resolution with both modalities. The different cell types could be easily identified due to variations in contrast within the acoustic and photoacoustic images. This technique provides a new way of probing leukocyte structure with potential applications towards detecting cellular abnormalities and diseased cells at the single cell level.

Out-coupling of Longitudinal Photoacoustic Pulses by Mitigating the Phase Cancellation

Waves of any kinds, including sound waves and light waves, can interfere constructively or destructively when they are overlapped, allowing for myriad applications. However, unlike continuous waves of a single frequency, interference of photoacoustic pulses is often overlooked because of their broadband characteristics and short pulse durations. Here, we study cancellation of two symmetric photoacoustic pulses radiated in the opposite direction from the same photoacoustic sources near a free surface. The cancellation occurs when one of the two pulses is reflected with polarity reversal from the free surface and catches up with the other. The cancellation effect, responsible for reduced signal amplitudes, is systematically examined by implementing a thin transparent matching medium of the same acoustic impedance. By changing the thickness of the transparent layer, the overlap of the two symmetric pulses is controlled. For optimized matching layers, the cancellation effect can be significantly reduced, while the resulting output waveform remains unchanged. Similar to the planar absorber, different dimensional absorbers including cylinders and spheres also exhibit the cancellation between the outward and inward waves. This work could provide further understanding of photoacoustic generation and a simple strategy for increasing photoacoustic signal amplitudes.

Pulsed photoacoustic flow imaging with a handheld system

Flow imaging is an important technique in a range of disease areas, but estimating low flow speeds, especially near the walls of blood vessels, remains challenging. Pulsed photoacoustic flow imaging can be an alternative since there is little signal contamination from background tissue with photoacoustic imaging. We propose flow imaging using a clinical photoacoustic system that is both handheld and portable. The system integrates a linear array with 7.5 MHz central frequency in combination with a high-repetition-rate diode laser to allow high-speed photoacoustic imaging--ideal for this application. This work shows the flow imaging performance of the system in vitro using microparticles. Both two-dimensional (2-D) flow images and quantitative flow velocities from 12 to 75  mm/s were obtained. In a transparent bulk medium, flow estimation showed standard errors of ∼7% the estimated speed; in the presence of tissue-realistic optical scattering, the error increased to 40% due to limited signal-to-noise ratio. In the future, photoacoustic flow imaging can potentially be performed in vivo using fluorophore-filled vesicles or with an improved setup on whole blood.

Three-Dimensional Optoacoustic and Laser-Induced Ultrasound Tomography System for Preclinical Research in Mice: Design and Phantom Validation

In this work, we introduce a novel three-dimensional imaging system for in vivo high-resolution anatomical and functional whole-body visualization of small animal models developed for preclinical and other type of biomedical research. The system (LOUIS-3DM) combines a multiwavelength optoacoustic tomography (OAT) and laser-induced ultrasound tomography (LUT) to obtain coregistered maps of tissue optical absorption and speed of sound, displayed within the skin outline of the studied animal. The most promising applications of the LOUIS-3DM include 3D angiography, cancer research, and longitudinal studies of biological distributions of optoacoustic contrast agents.

Synthesis and Testing of Modular Dual-Modality Nanoparticles for Magnetic Resonance and Multispectral Photoacoustic Imaging

Magnetic resonance (MR) and photoacoustic (PA) imaging are currently being investigated as complementing strategies for applications requiring sensitive detection of cells in vivo. While combined MR/PAI detection of cells requires biocompatible cell labeling probes, water-based synthesis of dual-modality MR/PAI probes presents significant technical challenges. Here we describe facile synthesis and characterization of hybrid modular dextran-stabilized gold/iron oxide (Au-IO) multimetallic nanoparticles (NP) enabling multimodal imaging of cells. The stable association between the IO and gold NP was achieved by priming the surface of dextran-coated IO with silver NP resulting from silver(I) reduction by aldehyde groups, which are naturally present within the dextran coating of IO at the level of 19-23 groups/particle. The Au-IO NP formed in the presence of silver-primed Au-IO were stabilized by using partially thiolated MPEG5-gPLL graft copolymer carrying residual amino groups. This stabilizer served as a carrier of near-infrared fluorophores (e.g., IRDye 800RS) for multispectral PA imaging. Dual modality imaging experiments performed in capillary phantoms of purified Au-IO-800RS NPs showed that these NPs were detectible using 3T MRI at a concentration of 25 μM iron. PA imaging achieved approximately 2.5-times higher detection sensitivity due to strong PA signal emissions at 530 and 770 nm, corresponding to gold plasmons and IRDye integrated into the coating of the hybrid NPs, respectively, with no "bleaching" of PA signal. MDA-MB-231 cells prelabeled with Au-IO-800RS retained plasma membrane integrity and were detectable by using both MR and dual-wavelength PA at 49 ± 3 cells/imaging voxel. We believe that modular assembly of multimetallic NPs shows promise for imaging analysis of engineered cells and tissues with high resolution and sensitivity.

The Effect of Acoustic Impedance on Subsurface Absorber Geometry Reconstruction using 1D Frequency-Domain Photoacoustics

This paper considers the effect of an impedance mismatch between the absorber and its surroundings on the aborber reconstructions from the photoacoustic signal profile, in particular when a non-delta input pulse is used. A transfer function approach is taken, demonstrating in the case of impedance mismatch how the total response can be modeled using the sum of the mismatch-free response and its time-delayed, time-reversed replicas, which may or may not overlap. It is shown how this approach can be exploited to accommodate the effects of non-delta pulses and/or pulse-equivalent waveforms such as linear-frequency-modulated (LFM) chirps, and impedance mismatches in any inversion algorithms, even in the presence of large reflection coefficients. As a consequence, for simple-absorber reconstruction algorithms that assume impulses or 'short enough' pulses, the compressive portion of the measured response may be used in reconstruction formulas that do not model the impedance mismatch, regardless of the size of the mismatch. For longer-duration input waveforms, it is demonstrated how existing reconstruction methods can be successfully adapted to include the effect of the impedance mismatch. Simulations are used to illustrate these ideas. The gained physical insight into how components of the generated pressure wave carry absorber information is then exploited for signal inversion and absorber reconstruction in the frequency domain when multi-frequency modulation chirps are used for photoacoustic radar pressure measurements. The foundational theoretical developments ultimately address impendance mismatch issues germane to the major photoacoustic frequency-domain imaging modality to-date, which is the photoacoustic radar.

Nano-structural characteristics of carbon nanotube-polymer composite films for high-amplitude optoacoustic generation

We demonstrate nano-structural characteristics of carbon nanotube (CNT)-polydimethylsiloxane (PDMS) composite films that can be used as highly efficient and robust ultrasound transmitters for diagnostic and therapeutic applications. An inherent architecture of the nano-composite provides unique thermal, optical, and mechanical properties that are accommodated not just for efficient energy conversion but also for extraordinary robustness against pulsed laser ablation. First, we explain a thermoacoustic transfer mechanism within the nano-composite. CNT morphologies are examined to determine a suitable arrangement for heat transfer to the surrounding PDMS. Next, we introduce an approach to enhance optical extinction of the composite films, which uses shadowed deposition of a thin Au layer through an as-grown CNT network. Finally, the transmitter robustness is quantified in terms of laser-induced damage threshold. This reveals that the CNT-PDMS films can withstand an order-of-magnitude higher optical fluence (and extinction) than a Cr film used as a reference. Such robustness is crucial to increase the maximum-available optical energy for optoacoustic excitation and pressure generation. All of these structure-originated characteristics manifest the CNT-PDMS composite films as excellent optoacoustic transmitters for high-amplitude and high-frequency ultrasound generation.

Spatiospectral denoising framework for multispectral optoacoustic imaging based on sparse signal representation

PURPOSE

One of the major challenges in dynamic multispectral optoacoustic imaging is its relatively low signal-to-noise ratio which often requires repetitive signal acquisition and averaging, thus limiting imaging rate. The development of denoising methods which prevent the need for signal averaging in time presents an important goal for advancing the dynamic capabilities of the technology.

METHODS

In this paper, a denoising method is developed for multispectral optoacoustic imaging which exploits the implicit sparsity of multispectral optoacoustic signals both in space and in spectrum. Noise suppression is achieved by applying thresholding on a combined wavelet-Karhunen-Loève representation in which multispectral optoacoustic signals appear particularly sparse. The method is based on inherent characteristics of multispectral optoacoustic signals of tissues, offering promise for general application in different incarnations of multispectral optoacoustic systems.

RESULTS

The performance of the proposed method is demonstrated on mouse images acquired in vivo for two common additive noise sources: time-varying parasitic signals and white noise. In both cases, the proposed method shows considerable improvement in image quality in comparison to previously published denoising strategies that do not consider multispectral information.

CONCLUSIONS

The suggested denoising methodology can achieve noise suppression with minimal signal loss and considerably outperforms previously proposed denoising strategies, holding promise for advancing the dynamic capabilities of multispectral optoacoustic imaging while retaining image quality.

Projection Onto Convex Sets (POCS) method for photoacoustic tomography with a non negative constraint

Photoacoustic imaging is a biomedical imaging modality capable of early cancer detection. In this paper, we proposed a novel iterative Projections Onto Convex Sets (POCS) method for improving photoacoustic reconstruction. This method aims to obtain a non negative pressure distribution satisfying the measured signals. This POCS method is performed in the Fourier Bessel space avoiding matrix inversions in the projections, speeding up projections and is capable of handling the large data sets present in photoacoustic imaging. The numerical experiments performed showed that improved reconstruction was obtained with a few iterations together with the recovery of some lost information.

Oscillatory Dynamics and In Vivo Photoacoustic Imaging Performance of Plasmonic Nanoparticle-Coated Microbubbles

Microbubbles bearing plasmonic nanoparticles on their surface provide contrast enhancement for both photoacoustic and ultrasound imaging. In this work, the responses of microbubbles with surface-bound gold nanorods-termed AuMBs-to nanosecond pulsed laser excitation are studied using high-speed microscopy, photoacoustic imaging, and numerical modeling. In response to laser fluences below 5 mJ cm(-2) , AuMBs produce weak photoacoustic emissions and exhibit negligible microbubble wall motion. However, in reponse to fluences above 5 mJ cm(-2) , AuMBs undergo dramatically increased thermal expansion and emit nonlinear photoacoustic waves of over 10-fold greater amplitude than would be expected from freely dispersed gold nanorods. Numerical modeling suggests that AuMB photoacoustic responses to low laser fluences result from conductive heat transfer from the surface-bound nanorods to the microbubble gas core, whereas at higher fluences, explosive boiling may occur at the nanorod surface, producing vapor nanobubbles that contribute to rapid AuMB expansion. The results of this study indicate that AuMBs are capable of producing acoustic emissions of significantly higher amplitude than those produced by conventional sources of photoacoustic contrast. In vivo imaging performance of AuMBs in a murine kidney model suggests that AuMBs may be an effective alternative to existing contrast agents for noninvasive photoacoustic and ultrasound imaging applications.

Photoacoustic transients generated by laser irradiation of thin films

Irradiation of an optically thin layer immersed in a transparent fluid with pulsed laser radiation can generate photoacoustic waves through two mechanisms. The first of these is the conventional optical heating of the layer followed by thermal expansion, in which the mechanical motion of the expansion launches a pair of oppositely directed sound waves. A second, recently reported mechanism, is operative when heat is conducted to the transparent medium raising its temperature, while at the same time reducing the temperature in the absorbing body. The latter mechanism has been shown to result in compressive transients at the leading edges of the photoacoustic waveforms. Here the photoacoustic effect produced by irradiating thin metal films which undergo negligible thermal expansion under optical irradiation, but which generate sound solely by the heat transfer mechanism is investigated. Solution to the wave equation for the photoacoustic effect from the heat transfer mechanism is given and compared with the results of experiments using nanosecond laser pulses to irradiate thin metal films.

Quantitative imaging of microvasculature in deep tissue with a spectrum-based photo-acoustic microscopy

We analyze photo-acoustic signals from capillaries and theoretically demonstrate the quantitative relationship between vascular diameter and spectral slope in a low-frequency band. Phantom experiments validate the theoretical analysis. Based on this finding, spectral slope is proposed as the imaging parameter of a photo-acoustic microscopy. This system effectively quantifies the microvasculature with diameters of 60 and 150 μm, which are smaller than the wavelength 342 μm at the central frequency 4.39 MHz of ultrasound transducer. The low frequency also guarantees the imaging depth in the order of centimeters. The proposed scheme could be potential for noninvasive diagnosis of diseases related to abnormal vasoconstriction or angiectasis.

Photoacoustic imaging of nanoparticle- containing cells using single-element focused transducer: a simulation study

A new theoretical approach for photoacoustic (PA) image simulation of an ensemble of cells with endocytosed gold nanoparticles is presented. Each cell was approximated as a fluid sphere and suspended in a nonabsorbing fluid medium. It was assumed that the cellular optical absorption coefficient changed greatly because of endocytosis of nanoparticles; however, thermophysical parameters remained unchanged because nanoparticles occupied negligible intracellular volume. A frequency-domain method was used to obtain a PA signal from a single cell and resultant signal detected by a focused single-element transducer was evaluated by convolving signals from many cells with the spatial impulse response function of the receiver. The proposed model was explored to simulate PA images of numerical phantoms. It was observed that features of the phantoms are retained precisely in those simulated images. Also, speckles in PA images are significantly suppressed because of strong boundary buildup when cells are bounded to a region. Nevertheless, speckle visibility increases when cells are not bounded to a region. This approach may be developed as a realistic simulation tool for PA imaging of tissue medium utilizing its cellular feature.

Simultaneous photoacoustic and optically mediated ultrasound microscopy: an in vivo study

We propose the use of thermoelastic (TE) excitation of an ultrasonic (US) detector by backscattered laser radiation as a means of upgrading a single-modality photoacoustic (PA) microscope to dual-modality PA/US imaging at minimal cost. The upgraded scanning head of our dual-modality microscope consists of a fiber bundle with 14 output arms and a 32MHz polyvinylidene difluoride (PVDF) detector with a 34 MHz bandwidth (-6 dB level), 12.7 mm focal length, and a 0.25 numerical aperture. A single optical pulse delivered through the fiber bundle to the biotissue being investigated, in combination with a metalized surface on the PVDF detector allows us to obtain both PA and US A-scans. To demonstrate the in vivo capabilities of the proposed method we present the results of bimodal imaging of the brain of a newborn rat, a mouse tail and a mouse tumor.

Ultrasonic-heating-encoded photoacoustic tomography with virtually augmented detection view

Photoacoustic (PA) imaging of arbitrarily-shaped or oriented objects may miss important features because PA waves propagate normal to structure boundaries and may miss the acoustic detectors when the detection view has a limited angular range. To overcome this long-standing problem, we present an ultrasonic thermal encoding approach that is universally applicable. We exploit the temperature dependence of the Grueneisen parameter and encode a confined [[What does confined mean here?]] voxel using heat generated by a focused ultrasonic transducer. The PA amplitude from the encoded voxel is increased while those from the neighboring voxels are unchanged. Consequently, the amplitude-increased PA waves propagate in all directions due to the round cross-section of the encoded region and thus can be received at any viewing angle on the cross-sectional plane [[Please check throughout the manuscript for similar places.]]. We built a mathematical model for the thermally encoded PA tomography, performed a numerical simulation, and experimentally validated the ultrasonic thermal encoding efficiency. As a proof of concept, we demonstrate full-view in vivo vascular imaging and compare it to the original linear-array PA tomography system, showing dramatically enhanced imaging of arbitrarily oriented blood vessels. Since ultrasonic heating can be focused deeply, this method can be applied to deep tissue imaging and is promising for full-view imaging of other features of biomedical interest, such as tumor margins.

Photoacoustic imaging of breast microcalcifications: a validation study with 3-dimensional ex vivo data and spectrophotometric measurement

This paper investigates whether photoacoustic imaging (PAI) can provide the visualization of microcalcifications in breast tissue. For this, the geometrical correlation between the 3-D PA images of breast microcalcifications within ex vivo specimens and the corresponding mammograms was ascertained. Also, the optical absorbance of the calcification compositions (i.e., hydroxyapatite and calcium oxalate) was measured and compared with the PA responses of the microcalcifications. The experimental results demonstrated that the PA images discriminated between the microcalcifications and the surrounding tissue, and their locations in PA images reasonably meshed with those of the microcalcifications appeared in the mammograms. Also, the change in PA signal amplitude along the laser wavelength agreed with the absorbance of hydroxyapatite associated with the relatively high potential of malignant cancers, but not calcium oxalate with only benign cases.

Proton beam characterization by proton-induced acoustic emission: simulation studies

Due to their Bragg peak, proton beams are capable of delivering a targeted dose of radiation to a narrow volume, but range uncertainties currently limit their accuracy. One promising beam characterization technique, protoacoustic range verification, measures the acoustic emission generated by the proton beam. We simulated the pressure waves generated by proton radiation passing through water. We observed that the proton-induced acoustic signal consists of two peaks, labeled α and γ, with two originating sources. The α acoustic peak is generated by the pre-Bragg peak heated region whereas the source of the γ acoustic peak is the proton Bragg peak. The arrival time of the α and γ peaks at a transducer reveals the distance from the beam propagation axis and Bragg peak center, respectively. The maximum pressure is not observed directly above the Bragg peak due to interference of the acoustic signals. Range verification based on the arrival times is shown to be more effective than determining the Bragg peak position based on pressure amplitudes. The temporal width of the α and γ peaks are linearly proportional to the beam diameter and Bragg peak width, respectively. The temporal separation between compression and rarefaction peaks is proportional to the spill time width. The pressure wave expected from a spread out Bragg peak dose is characterized. The simulations also show that acoustic monitoring can verify the proton beam dose distribution and range by characterizing the Bragg peak position to within ~1 mm.

Computational modeling of photoacoustic signals from mixtures of melanoma and red blood cells

A theoretical approach to model photoacoustic (PA) signals from mixtures of melanoma cells (MCs) and red blood cells (RBCs) is discussed. The PA signal from a cell approximated as a fluid sphere was evaluated using a frequency domain method. The tiny signals from individual cells were summed up obtaining the resultant PA signal. The local signal to noise ratio for a MC was about 5.32 and 5.40 for 639 and 822 nm illuminations, respectively. The PA amplitude exhibited a monotonic rise with increasing number of MCs for each incident radiation. The power spectral lines also demonstrated similar variations over a large frequency range (5-200 MHz). For instance, spectral intensity was observed to be 5.5 and 4.0 dB greater at 7.5 MHz for a diseased sample containing 1 MC and 22,952 RBCs than a normal sample composed of 22,958 RBCs at those irradiations, respectively. The envelope histograms generated from PA signals for mixtures of small numbers of MCs and large numbers of RBCs seemed to obey pre-Rayleigh statistics. The generalized gamma distribution found to facilitate better fits to the histograms than the Rayleigh and Nakagami distributions. The model provides a means to study PAs from mixtures of different populations of absorbers.

Analytic theory of photoacoustic wave generation from a spheroidal droplet

In this paper, we develop an analytic theory for describing the photoacoustic wave generation from a spheroidal droplet and derive the first complete analytic solution. Our derivation is based on solving the photoacoustic Helmholtz equation in spheroidal coordinates with the separation-of-variables method. As the verification, besides carrying out the asymptotic analyses which recover the standard solutions for a sphere, an infinite cylinder and an infinite layer, we also confirm that the partial transmission and reflection model previously demonstrated for these three geometries still stands. We expect that this analytic solution will find broad practical uses in interpreting experiment results, considering that its building blocks, the spheroidal wave functions (SWFs), can be numerically calculated by the existing computer programs.

Photoacoustic radiation force on a microbubble

We investigate the radiation force on a microbubble due to the photoacoustic wave which is generated by using a pulsed laser. In particular, we focus on the dependence of pulsed laser parameters on the radiation force. In order to do so, we first obtain a new and comprehensive analytical solution to the photoacoustic wave equation based on the Fourier transform for various absorption profiles. Then, we write an expression of the radiation force containing explicit laser parameters, pulse duration, and beamwidth of the laser. Furthermore, we calculate the primary radiation force acting on a microbubble. We show that laser parameters and the position of the microbubble relative to a photoacoustic source have a considerable effect on the primary radiation force. By means of recent developments in laser technologies that render tunability of pulse duration and repetition frequency possible, an adjustable radiation force can be applied to microbubbles. High spatial control of applied force is ensured on account of smaller focal spots achievable by focused optics. In this context, conventional piezoelectric acoustic source applications could be surpassed. In addition, it is possible to increase the radiation force by making source wavelength with the absorption peak of absorber concurrent. The application of photoacoustic radiation force can open a cache of opportunities such as manipulation of microbubbles used as contrast agents and as carrier vehicles for drugs and genes with a desired force along with in vivo applications.

Photoacoustic microscopy and computed tomography: from bench to bedside

Photoacoustic imaging (PAI) of biological tissue has seen immense growth in the past decade, providing unprecedented spatial resolution and functional information at depths in the optical diffusive regime. PAI uniquely combines the advantages of optical excitation and those of acoustic detection. The hybrid imaging modality features high sensitivity to optical absorption and wide scalability of spatial resolution with the desired imaging depth. Here we first summarize the fundamental principles underpinning the technology, then highlight its practical implementation, and finally discuss recent advances toward clinical translation.

Multiscale multispectral optoacoustic tomography by a stationary wavelet transform prior to unmixing

Multispectral optoacoustic tomography (MSOT) utilizes broadband ultrasound detection for imaging biologically-relevant optical absorption features at a range of scales. Due to the multiscale and multispectral features of the technology, MSOT comes with distinct requirements in implementation and data analysis. In this work, we investigate the interplay between scale, which depends on ultrasonic detection frequency, and optical multispectral spectral analysis, two dimensions that are unique to MSOT and represent a previously unexplored challenge. We show that ultrasound frequency-dependent artifacts suppress multispectral features and complicate spectral analysis. In response, we employ a wavelet decomposition to perform spectral unmixing on a per-scale basis (or per ultrasound frequency band) and showcase imaging of fine-scale features otherwise hidden by low frequency components. We explain the proposed algorithm by means of simple simulations and demonstrate improved performance in imaging data of blood vessels in human subjects.

Characterization of the spatio-temporal response of optical fiber sensors to incident spherical waves

In this study a theoretical framework for calculating the acoustic response of optical fiber-based ultrasound sensors is presented. The acoustic response is evaluated for optical fibers with several layers of coating assuming a harmonic point source with arbitrary position and frequency. First, the fiber is acoustically modeled by a layered cylinder on which spherical waves are impinged. The scattering of the acoustic waves is calculated analytically and used to find the normal components of the strains on the fiber axis. Then, a strain-optic model is used to calculate the phase shift experienced by the guided mode in the fiber owing to the induced strains. The framework is showcased for a silica fiber with two layers of coating for frequencies in the megahertz regime, commonly used in medical imaging applications. The theoretical results are compared to experimental data obtained with a sensing element based on a pi-phase-shifted fiber Bragg grating and with photoacoustically generated ultrasonic signals.

24-MHz scanner for optoacoustic imaging of skin and burn

Optoacoustic (photoacoustic) imaging uniquely visualizes optical contrast in high resolution and comes with very attractive characteristics for clinical imaging applications. In this paper, we showcase the performance of a scanner based on a 24 MHz center-frequency 128 element array, developed for applications in dermatology. We perform system characterization to examine the imaging performance achieved. We then showcase its imaging ability on healthy tissue and cancer. Finally, we image burns and human lesions in vivo and gain insights on the benefits and challenges of this approach as it is considered for diagnostic and treatment follow-up applications in dermatology and beyond.

Model-based optoacoustic image reconstruction of large three-dimensional tomographic datasets acquired with an array of directional detectors

Image quality in 3-D optoacoustic (photoacoustic) tomography is greatly influenced by both the measurement system, in particular the number and spatial arrangement of ultrasound sensors, and the ability to account for the spatio-temporal response of the sensor element(s) in the reconstruction algorithm. Herein we present a reconstruction procedure based on the inversion of a time-domain forward model incorporating the spatial impulse response due to the shape of the transducer, which is subsequently applied in a tomographic system based on a translation-rotation scan of a linear detector array. The proposed method was also adapted to cope with the data-intensive requirements of high-resolution volumetric optoacoustic imaging. The processing of 2 · 10 (4) individual signals resulted in well-resolved images of both ~ 200 μm absorbers in phantoms and complex vascular structures in biological tissue. The results reported herein demonstrate that the introduced model-based methodology exhibits a better contrast and resolution than standard back-projection and model-based algorithms that assume point detectors. Moreover, the capability of handling large datasets anticipates that model-based methods incorporating the sensor properties can become standard practice in volumetric opto acoustic image formation.

Photoacoustic pulse wave forming along the rotation axis of an ellipsoid droplet: a geometric calculation study

A geometric calculation method is developed to study the pulsed photoacoustic wave forming of an arbitrarily shaped droplet. It is found that for an ellipsoid droplet, either a prolate ellipsoid or an oblate ellipsoid, strict analytical formulas for describing the wave profile developed along the rotation axis can be derived. The results show intriguing differences compared to those of a sphere droplet in terms of the multiple geometric parameters being in effect, the pulse wave profile variant, and the existing of unlimited points of infinite tensile pressure.

Raster-scan optoacoustic mesoscopy in the 25-125 MHz range

We developed a raster-scan acoustic resolution broadband optoacoustic mesoscopy system and investigated the imaging performance using ultrasonic frequencies up to 125 MHz. The developed system achieves 7 μm axial resolution and transverse resolution of 30 μm reaching depths of at least 5 mm. This unprecedented performance is achieved by operating at out-of-focus ultrasonic detection and tomographic reconstruction. We demonstrate the limits reached due to the width of the laser pulse employed and showcase the technique on drosophila fly and drosophila pupae ex vivo.

A theoretical investigation of photoacoustic contrast agents

Photoacoustic imaging offers significant potential as a biomedical imaging modality. For some applications, however, there is a need for contrast enhancement. In this paper, a theoretical comparison is presented of the efficacy of three different designs for photoacoustic contrast agents (PACAs), specifically, a droplet of dye, a bubble filled with gas coated by a dye loaded shell, and a droplet of volatile dye. For each case, the governing equations describing the dynamics of a single PACA in a homogenous incompressible fluid are derived. The coupled sets of equations describing the bubble oscillation and resulting radiated pressure, the photo-acoustic energy equation, and the equation of state are then solved numerically. The numerical results predict a stronger radiated acoustic signal for the same optical source energy density in the case of the volatile dye droplet by a factor of up to two orders of magnitude compared with the other two types of agent.

Optoacoustic determination of spatio-temporal responses of ultrasound sensors

The characterization of the spatial and frequency response of acoustic detectors is important for enabling accurate optoacoustic imaging. In this work, we developed a hybrid method for the characterization of the spatially dependent response of ultrasound detectors. The method is based on the experimental determination of the receive-mode electrical impulse response (EIR) of the sensor, which is subsequently convolved with the corresponding spatial impulse response (SIR), computed numerically. The hybrid method is shown to have superior performance over purely experimental techniques in terms of accurate determination of the spatial and temporal responses of ultrasonic detectors, in high as well as low sensitivity regions of the sensor.

Co-registered pulse-echo/photoacoustic transvaginal probe for real time imaging of ovarian tissue

We present the design and construction of a prototype imaging probe capable of co-registered pulse-echo ultrasound and photoacoustic (optoacoustic) imaging in real time. The probe consists of 36 fibers of 200 micron core diameter each that are distributed around a commercial transvaginal ultrasound transducer, and housed in a protective shield. Its performance was demonstrated by two sets of experiments. The first set involved imaging of blood flowing through a tube mimicking a blood vessel, the second set involved imaging of human ovaries ex vivo. The results suggest that the system along with the probe has great potential for imaging and characterizing of ovarian tissue in vivo.

Carbon-nanotube optoacoustic lens for focused ultrasound generation and high-precision targeted therapy

We demonstrate a new optical approach to generate high-frequency (>15 MHz) and high-amplitude focused ultrasound, which can be used for non-invasive ultrasound therapy. A nano-composite film of carbon nanotubes (CNTs) and elastomeric polymer is formed on concave lenses, and used as an efficient optoacoustic source due to the high optical absorption of the CNTs and rapid heat transfer to the polymer upon excitation by pulsed laser irradiation. The CNT-coated lenses can generate unprecedented optoacoustic pressures of >50 MPa in peak positive on a tight focal spot of 75 μm in lateral and 400 μm in axial widths. This pressure amplitude is remarkably high in this frequency regime, producing pronounced shock effects and non-thermal pulsed cavitation at the focal zone. We demonstrate that the optoacoustic lens can be used for micro-scale ultrasonic fragmentation of solid materials and a single-cell surgery in terms of removing the cells from substrates and neighboring cells.

Ultrashort microwave-induced thermoacoustic imaging: a breakthrough in excitation efficiency and spatial resolution

With theoretical prediction and experimental validation, we propose a novel approach to significantly enhance the conversion efficiency of thermoacoustic (TA) imaging by using an ultrashort microwave pulse. The implementation of the ultrashort microwave pulse leads to orders of magnitude enhancement in excitation efficiency and spatial resolution, compared to that from existing TA imaging techniques. This allows high-resolution (~ 100 micron resolution) TA imaging to be acquired noninvasively. The present work represents a major breakthrough in the conversion efficiency of the TA effect and the resolution of TA imaging, which can potentially be used for clinical imaging.

Frequency-radial duality based photoacoustic image reconstruction

Photoacoustic image reconstruction algorithms are usually slow due to the large sizes of data that are processed. This paper proposes a method for exact photoacoustic reconstruction for the spherical geometry in the limiting case of a continuous aperture and infinite measurement bandwidth that is faster than existing methods namely (1) backprojection method and (2) the Norton-Linzer method [S. J. Norton and M. Linzer, "Ultrasonic reflectivity imaging in three dimensions: Exact inverse scattering solution for plane, cylindrical and spherical apertures," Biomedical Engineering, IEEE Trans. BME 28, 202-220 (1981)]. The initial pressure distribution is expanded using a spherical Fourier Bessel series. The proposed method estimates the Fourier Bessel coefficients and subsequently recovers the pressure distribution. A concept of frequency-radial duality is introduced that separates the information from the different radial basis functions by using frequencies corresponding to the Bessel zeros. This approach provides a means to analyze the information obtained given a measurement bandwidth. Using order analysis and numerical experiments, the proposed method is shown to be faster than both the backprojection and the Norton-Linzer methods. Further, the reconstructed images using the proposed methodology were of similar quality to the Norton-Linzer method and were better than the approximate backprojection method.

Mathieu function solutions for photoacoustic waves in sinusoidal one-dimensional structures

The photoacoustic effect for a one-dimensional structure, the sound speed of which varies sinusoidally in space, is shown to be governed by an inhomogeneous Mathieu equation with the forcing term dependent on the spatial and temporal properties of the exciting optical radiation. New orthogonality relations, traveling wave Mathieu functions, and solutions to the inhomogeneous Mathieu equation are found, which are used to determine the character of photoacoustic waves in infinite and finite length phononic structures. Floquet solutions to the Mathieu equation give the positions of the band gaps, the damping of the acoustic waves within the band gaps, and the dispersion relation for photoacoustic waves. The solutions to the Mathieu equation give the photoacoustic response of the structure, show the space equivalent of subharmonic generation and acoustic confinement when waves are excited within band gaps.

Determination of optical absorption coefficient with focusing photoacoustic imaging

Absorption coefficient of biological tissue is an important factor for photothermal therapy and photoacoustic imaging. However, its determination remains a challenge. In this paper, we propose a method using focusing photoacoustic imaging technique to quantify the target optical absorption coefficient. It utilizes the ratio of the amplitude of the peak signal from the top boundary of the target to that from the bottom boundary based on wavelet transform. This method is self-calibrating. Factors, such as absolute optical fluence, ultrasound parameters, and Grüneisen parameter, can be canceled by dividing the amplitudes of the two peaks. To demonstrate this method, we quantified the optical absorption coefficient of a target with various concentrations of an absorbing dye. This method is particularly useful to provide accurate absorption coefficient for predicting the outcomes of photothermal interaction for cancer treatment with absorption enhancement.

Dependence of photoacoustic speckles on boundary roughness

Speckles have been considered ubiquitous in all scattering-based coherent imaging technologies. However, as an optical-absorption-based coherent imaging technology, photoacoustic (PA) tomography (PAT) suppresses speckles by building up prominent boundary signals. We theoretically study the dependence of PAT speckles on the boundary roughness, which is quantified by the root-mean-squared value and the correlation length of the boundary height. Both the speckle visibility and the correlation coefficient between the reconstructed and actual boundaries are quantified. If the root-mean-squared height fluctuation is much greater than, and the height correlation length is much smaller than the imaging resolution, the reconstructed boundaries become fully developed speckles. In other words, speckle formation requires large uncorrelated height fluctuations within the resolution cell. The first- and second-order statistics of PAT speckles are also studied experimentally. While the amplitude of the speckles follows a Gaussian distribution, the autocorrelation of the speckle patterns tracks that of the system point spread function.

Model-based reconstruction integrated with fluence compensation for photoacoustic tomography

Photoacoustic (PA) tomography (PAT) is a rapidly developing imaging modality that can provide high contrast and spatial-resolution images of light-absorption distribution in tissue. However, reconstruction of the absorption distribution is affected by nonuniform light fluence. This paper introduces a reconstruction method for reducing amplification of noise and artifacts in low-fluence regions. In this method, fluence compensation is integrated into model-based reconstruction, and the absorption distribution is iteratively updated. At each iteration, we calculate the residual between detected PA signals and the signals computed by a forward model using the initial pressure, which is the product of estimated voxel value and light fluence. By minimizing the residual, the reconstructed values converge to the true absorption distribution. In addition, we developed a matrix compression method for reducing memory requirements and accelerating reconstruction speed. The results of simulation and phantom experiments indicate that the proposed method provides a better contrast-to-noise ratio (CNR) in low-fluence regions. We expect that the capability of increasing imaging depth will broaden the clinical applications of PAT.

Environment-dependent generation of photoacoustic waves from plasmonic nanoparticles

Nanoparticle-augmented photoacoustics is an emerging technique for molecular imaging. This study investigates the fundamental process of the photoacoustic signal generation by plasmonic nanoparticles suspended in a weakly absorbing fluid. The photoacoustic signal of gold nanospheres with varying silica shell thicknesses is shown to be dominated by the heat transfer between the nanoparticles and the surrounding environment.