The effect of source backing materials and excitation pulse durations on laser-generated ultrasound waveforms.

In this article, it is shown experimentally that a planar laser-generated ultrasound source with a hard reflective backing will generate higher acoustic pressures than a comparable source with an acoustically matched backing when the stress confinement condition is not met. Furthermore, while the source with an acoustically matched backing will have a broader bandwidth when the laser pulse is short enough to ensure stress confinement, the bandwidths of both source types will converge as the laser pulse duration increases beyond stress confinement. The explanation of the results is supported by numerical simulations.

Modelling laser ultrasound waveforms: The effect of varying pulse duration and material properties.

Optical generation of ultrasound using nanosecond duration laser pulses has generated great interest both in industrial and biomedical applications. The availability of portable laser devices using semiconductor technology and optical fibres, as well as numerous source material types based on nanocomposites, has proliferated the applications of laser ultrasound. The nanocomposites can be deposited on the tip of optical fibres as well as planar hard and soft backing materials using various fabrication techniques, making devices suitable for a variety of applications. The ability to choose the acoustic material properties and the laser pulse duration gives considerable control over the ultrasound output. Here, an analytical time-domain solution is derived for the acoustic pressure waveform generated by a planar optical ultrasound source consisting of an optically absorbing layer on a backing. It is shown that by varying the optical attenuation coefficient, the thickness of the absorbing layer, the acoustic properties of the materials, and the laser pulse duration, a wide variety of pulse shapes and trains can be generated. It is shown that a source with a reflecting backing can generate pulses with higher amplitude than a source with an acoustically-matched backing in the same circumstances when stress-confinement has not been satisfied.

Measurement of the ultrasound attenuation and dispersion in 3D-printed photopolymer materials from 1 to 3.5 MHz.

Over the past decade, the range of applications in biomedical ultrasound exploiting 3D printing has rapidly expanded. For wavefront shaping specifically, 3D printing has enabled a diverse range of new, low-cost approaches for controlling acoustic fields. These methods rely on accurate knowledge of the bulk acoustic properties of the materials; however, to date, robust knowledge of these parameters is lacking for many materials that are commonly used. In this work, the acoustic properties of eight 3D-printed photopolymer materials were characterised over a frequency range from 1 to 3.5 MHz. The properties measured were the frequency-dependent phase velocity and attenuation, group velocity, signal velocity, and mass density. The materials were fabricated using two separate techniques [PolyJet and stereolithograph (SLA)], and included Agilus30, FLXA9960, FLXA9995, Formlabs Clear, RGDA8625, RGDA8630, VeroClear, and VeroWhite. The range of measured density values across all eight materials was 1120-1180 kg · m-3, while the sound speed values were between 2020 to 2630 m · s-1, and attenuation values typically in the range 3-9 dB · MHz-1· cm-1.

Laser generated ultrasound sources using carbon-polymer nanocomposites for high frequency metrology.

The characterization of ultrasound fields generated by diagnostic and therapeutic equipment is an essential requirement for performance validation and to demonstrate compliance against established safety limits. This requires hydrophones calibrated to a traceable standard. Currently, the upper calibration frequency range available to the user community is limited to 60 MHz. However, high frequencies are increasingly being used for both imaging and therapy necessitating calibration frequencies up to 100 MHz. The precise calibration of hydrophones requires a source of high amplitude, broadband, quasi-planar, and stable ultrasound fields. There are challenges to using conventional piezoelectric sources, and laser generated ultrasound sources offer a promising solution. In this study, various nanocomposites consisting of a bulk polymer matrix and multi-walled carbon nanotubes were fabricated and tested using pulsed laser of a few nanoseconds for their suitability as a source for high frequency calibration of hydrophones. The pressure amplitude and bandwidths were measured using a broadband hydrophone from 27 different nanocomposite sources. The effect of nonlinear propagation of high amplitude laser generated ultrasound on bandwidth and the effect of bandlimited sensitivity response on the deconvolved pressure waveform were numerically investigated. The stability of the nanocomposite sources under sustained laser pulse excitation was also examined.

Three dimensional photoacoustic tomography in Bayesian framework.

The image reconstruction problem (or inverse problem) in photoacoustic tomography is to resolve the initial pressure distribution from detected ultrasound waves generated within an object due to an illumination by a short light pulse. Recently, a Bayesian approach to photoacoustic image reconstruction with uncertainty quantification was proposed and studied with two dimensional numerical simulations. In this paper, the approach is extended to three spatial dimensions and, in addition to numerical simulations, experimental data are considered. The solution of the inverse problem is obtained by computing point estimates, i.e., maximum a posteriori estimate and posterior covariance. These are computed iteratively in a matrix-free form using a biconjugate gradient stabilized method utilizing the adjoint of the acoustic forward operator. The results show that the Bayesian approach can produce accurate estimates of the initial pressure distribution in realistic measurement geometries and that the reliability of these estimates can be assessed.

Accurate simulation of transcranial ultrasound propagation for ultrasonic neuromodulation and stimulation.

Non-invasive, focal neurostimulation with ultrasound is a potentially powerful neuroscientific tool that requires effective transcranial focusing of ultrasound to develop. Time-reversal (TR) focusing using numerical simulations of transcranial ultrasound propagation can correct for the effect of the skull, but relies on accurate simulations. Here, focusing requirements for ultrasonic neurostimulation are established through a review of previously employed ultrasonic parameters, and consideration of deep brain targets. The specific limitations of finite-difference time domain (FDTD) and k-space corrected pseudospectral time domain (PSTD) schemes are tested numerically to establish the spatial points per wavelength and temporal points per period needed to achieve the desired accuracy while minimizing the computational burden. These criteria are confirmed through convergence testing of a fully simulated TR protocol using a virtual skull. The k-space PSTD scheme performed as well as, or better than, the widely used FDTD scheme across all individual error tests and in the convergence of large scale models, recommending it for use in simulated TR. Staircasing was shown to be the most serious source of error. Convergence testing indicated that higher sampling is required to achieve fine control of the pressure amplitude at the target than is needed for accurate spatial targeting.

Control of broadband optically generated ultrasound pulses using binary amplitude holograms.

In this work, the use of binary amplitude holography is investigated as a mechanism to focus broadband acoustic pulses generated by high peak-power pulsed lasers. Two algorithms are described for the calculation of the binary holograms; one using ray-tracing, and one using an optimization based on direct binary search. It is shown using numerical simulations that when a binary amplitude hologram is excited by a train of laser pulses at its design frequency, the acoustic field can be focused at a pre-determined distribution of points, including single and multiple focal points, and line and square foci. The numerical results are validated by acoustic field measurements from binary amplitude holograms, excited by a high peak-power laser.

A learned Born series for highly-scattering media.

A new method for solving the wave equation is presented, called the learned Born series (LBS), which is derived from a convergent Born series but its components are found through training. The LBS is shown to be significantly more accurate than the convergent Born series for the same number of iterations, in the presence of high contrast scatterers, while maintaining a comparable computational complexity. The LBS is able to generate a reasonable prediction of the global pressure field with a small number of iterations, and the errors decrease with the number of learned iterations.

open-UST: An Open-Source Ultrasound Tomography Transducer Array System.

Fast imaging methods are needed to promote clinical adoption of ultrasound tomography (UST), and more widely available UST hardware could support the experimental validation of new measurement configurations. In this work, an open-source 256-element transducer ring array was developed (morganjroberts.github. io/open-UST) and manufactured using rapid prototyping, for only £2k. Novel manufacturing techniques were used, resulting in a 1.17° mean beam axis skew angle, a [Formula: see text] mean element position error, and a [Formula: see text] deviation in matching layer thickness. The nominal acoustic performance was measured using hydrophone scans and watershot data, and the 61.2 dB signal-to-noise ratio (SNR), 55.4° opening angle, 10.2 mm beamwidth, and 54% transmit-receive bandwidth (-12 dB) were found to be similar to existing systems and compatible with state-of-the-art full-waveform-inversion image reconstruction methods. The interelement variation in acoustic performance was typically < 10% without using normalization, meaning that the elements can be modeled identically during image reconstruction, removing the need for individual source definitions based on hydrophone measurements. Finally, data from a phantom experiment were successfully reconstructed. These results demonstrate that the open-UST system is accessible for users and is suitable for UST imaging research.

The IPASC data format: A consensus data format for photoacoustic imaging.

Photoacoustic imaging (PAI) is an emerging modality that has shown promise for improving patient management in a range of applications. Unfortunately, the current lack of uniformity in PAI data formats compromises inter-user data exchange and comparison, which impedes: technological progress; effective research collaboration; and efforts to deliver multi-centre clinical trials. To overcome this challenge, the International Photoacoustic Standardisation Consortium (IPASC) has established a data format with a defined consensus metadata structure and developed an open-source software application programming interface (API) to enable conversion from proprietary file formats into the IPASC format. The format is based on Hierarchical Data Format 5 (HDF5) and designed to store photoacoustic raw time series data. Internal quality control mechanisms are included to ensure completeness and consistency of the converted data. By unifying the variety of proprietary data and metadata definitions into a consensus format, IPASC hopes to facilitate the exchange and comparison of PAI data.

100 MHz bandwidth planar laser-generated ultrasound source for hydrophone calibration.

High-frequency calibration of hydrophones is becoming increasingly important, both for clinical and scientific applications of ultrasound, and user safety. At present, the calibrations available routinely to the user community extend to 60 MHz. However, hydrophones that can measure beyond this are available, and ultrasonic fields often contain energy at higher frequencies, e.g., generated through nonlinear propagation of high-amplitude ultrasound used for therapeutic applications, and the increasing use of higher frequencies in imaging. Therefore, there is a need for calibrations up to at least 100 MHz, to allow ultrasonic fields to be accurately characterized, and the risk of harmful bioeffects to be properly assessed. Currently, sets of focused piezoelectric transducers are used to meet the pressure amplitude and bandwidth requirements of Primary Standard calibration facilities. However, when the frequency is high enough such that the size of the ultrasound focus becomes less than the hydrophone element's diameter, the uncertainty due to spatial averaging becomes significant, and can be as high as 20% at 100 MHz. As an alternate to piezoelectric transducers, a laser-generated ultrasound calibration source was designed, fabricated, and characterized. The source consists of an optically absorbing carbon-polymer nanocomposite excited by a large-diameter 1064 nm laser pulse of 2.6 ns duration. Peak pressure amplitudes of several Mega-Pascal were readily achievable, and the signal contained measurable frequency components up to 100 MHz. The variation in the pressure amplitudes was less than 2% from its mean over a three-hour test period. The ultrasound beam was sufficiently broad that the uncertainties due to spatial averaging were negligible.

The Effect of Curing Temperature and Time on the Acoustic and Optical Properties of PVCP.

Polyvinyl chloride plastisol (PVCP) has been increasingly used as a phantom material for photoacoustic and ultrasound imaging. As one of the most useful polymeric materials for industrial applications, its mechanical properties and behavior are well-known. Although the acoustic and optical properties of several formulations have previously been investigated, it is still unknown how these are affected by varying the fabrication method. Here, an improved and straightforward fabrication method is presented, and the effect of curing temperature and curing time on the PVCP acoustic and optical properties, as well as their stability over time, is investigated. The speed of sound and attenuation were determined over a frequency range from 2 to 15 MHz, while the optical attenuation spectra of samples were measured over a wavelength range from 500 to 2200 nm. The results indicate that the optimum properties are achieved at curing temperatures between 160 °C and 180 °C, while the required curing time decreases with increasing temperature. The properties of the fabricated phantoms were highly repeatable, meaning that the phantoms are not sensitive to the manufacturing conditions provided that the curing temperature and time are within the range of complete gelation-fusion (samples are optically clear) and below the limit of thermal degradation (indicated by the yellowish appearance of the sample). The samples' long-term stability was assessed over 16 weeks, and no significant change was observed in the measured acoustic and optical properties.

Analysis of the Directivity of Glass-Etalon Fabry-Pérot Ultrasound Sensors.

Planar glass-etalon Fabry-Pérot (FP) optical ultrasound sensors offer an alternative to piezoelectric sensors for the measurements of high-intensity focused ultrasound (HIFU) fields and other metrological applications. In this work, a model of the frequency-dependent directional response of the FP sensor was developed using the global matrix method, treating the sensor as a multilayered elastic structure. The model was validated against the experimentally measured directional response of an air-backed cover-slip FP sensor with well-known material properties. In addition, the model was compared with the measurements of an all-hard-dielectric sensor suitable for HIFU measurements. The model was then used to calculate modal dispersion curves for both glass-etalon sensors, allowing the features of the directional response to be linked to specific wave phenomena. The features in the directivity of the air-backed cover-slip sensor are due to guided Lamb waves. Symmetric Lamb modes give rise to regions of high sensitivity, whereas anti-symmetric modes cause regions of low sensitivity. For the all-hard-dielectric sensor, two features correspond to the water-substrate and water-spacer compressional and shear critical angles. A region of high sensitivity close to the shear critical angle is associated with a leaky-Rayleigh wave, which has a frequency-dependent phase speed. At higher frequencies, this feature is counteracted by a region of low sensitivity, which occurs when there is no difference in the vertical displacement of the mirrors forming the FP cavity. The model may be used to improve and optimize the design of FP sensors or could be used to assist with the accurate deconvolution of the directional response from array measurements in metrological and imaging applications.

Direct Estimation of Optical Parameters From Photoacoustic Time Series in Quantitative Photoacoustic Tomography.

Estimation of optical absorption and scattering of a target is an inverse problem associated with quantitative photoacoustic tomography. Conventionally, the problem is expressed as two folded. First, images of initial pressure distribution created by absorption of a light pulse are formed based on acoustic boundary measurements. Then, the optical properties are determined based on these photoacoustic images. The optical stage of the inverse problem can thus suffer from, for example, artefacts caused by the acoustic stage. These could be caused by imperfections in the acoustic measurement setting, of which an example is a limited view acoustic measurement geometry. In this work, the forward model of quantitative photoacoustic tomography is treated as a coupled acoustic and optical model and the inverse problem is solved by using a Bayesian approach. Spatial distribution of the optical properties of the imaged target are estimated directly from the photoacoustic time series in varying acoustic detection and optical illumination configurations. It is numerically demonstrated, that estimation of optical properties of the imaged target is feasible in limited view acoustic detection setting.

Reconstruction-classification method for quantitative photoacoustic tomography.

We propose a combined reconstruction-classification method for simultaneously recovering absorption and scattering in turbid media from images of absorbed optical energy. This method exploits knowledge that optical parameters are determined by a limited number of classes to iteratively improve their estimate. Numerical experiments show that the proposed approach allows for accurate recovery of absorption and scattering in two and three dimensions, and delivers superior image quality with respect to traditional reconstruction-only approaches.

Quantitative photoacoustic tomography using illuminations from a single direction.

Quantitative photoacoustic tomography is an emerging imaging technique aimed at estimating optical parameters inside tissues from photoacoustic images, which are formed by combining optical information and ultrasonic propagation. This optical parameter estimation problem is ill-posed and needs to be approached within the framework of inverse problems. It has been shown that, in general, estimating the spatial distribution of more than one optical parameter is a nonunique problem unless more than one illumination pattern is used. Generally, this is overcome by illuminating the target from various directions. However, in some cases, for example when thick samples are investigated, illuminating the target from different directions may not be possible. In this work, the use of spatially modulated illumination patterns at one side of the target is investigated with simulations. The results show that the spatially modulated illumination patterns from a single direction could be used to provide multiple illuminations for quantitative photoacoustic tomography. Furthermore, the results show that the approach can be used to distinguish absorption and scattering inclusions located near the surface of the target. However, when compared to a full multidirection illumination setup, the approach cannot be used to image as deep inside tissues.

Bayesian Image Reconstruction in Quantitative Photoacoustic Tomography.

Quantitative photoacoustic tomography is an emerging imaging technique aimed at estimating chromophore concentrations inside tissues from photoacoustic images, which are formed by combining optical information and ultrasonic propagation. This is a hybrid imaging problem in which the solution of one inverse problem acts as the data for another ill-posed inverse problem. In the optical reconstruction of quantitative photoacoustic tomography, the data is obtained as a solution of an acoustic inverse initial value problem. Thus, both the data and the noise are affected by the method applied to solve the acoustic inverse problem. In this paper, the noise of optical data is modelled as Gaussian distributed with mean and covariance approximated by solving several acoustic inverse initial value problems using acoustic noise samples as data. Furthermore, Bayesian approximation error modelling is applied to compensate for the modelling errors in the optical data caused by the acoustic solver. The results show that modelling of the noise statistics and the approximation errors can improve the optical reconstructions.

Measurement of the ultrasound attenuation and dispersion in whole human blood and its components from 0-70 MHz.

The ultrasound attenuation coefficient and dispersion from 0-70 MHz in whole human blood and its components (red blood cells and plasma) at 37°C is reported. The measurements are made using a fixed path substitution technique that exploits optical mechanisms for the generation and detection of ultrasound. This allows the measurements to cover a broad frequency range with a single source and receiver. The measured attenuation coefficient and dispersion in solutions of red blood cells and physiological saline for total haemoglobin concentrations of 10, 15 and 20 g/dL are presented. The attenuation coefficient and dispersion in whole human blood taken from four healthy volunteers by venipuncture is also reported. The power law dependence of the attenuation coefficient is shown to vary across the measured frequency range. This is due to the varying frequency dependence of the different mechanisms responsible for the attenuation. The attenuation coefficient measured at high frequencies is found to be significantly higher than that predicted by historical power law parameters. A review of the attenuation mechanisms in blood along with previously reported experimental measurements is given. Values for the sound speed and density in the tested samples are also presented.