The photoacoustic or optoacoustic effect has seen two major renaissances since its discovery by Bell. The first transformation was facilitated by lasers and piezoelectric transducers to create an ultra-sensitive spectroscopy technique to characterize the absorptive properties of gaseous, liquid, and solid materials in the form of bulk, powders and surface films.[2–3] It has also been very helpful in characterizing photo-induced chemical reactions, and in chemical, thermal, and structural sensing, or in the characterization of two-photon absorption cross sections. Theoretical descriptions of various experimental situations were developed by a number of authors,[7–10] and summarized for instance in. More recently, a second renewed interest has developed motivated by biomedical imaging and image-guided therapy, although photoacoustic techniques have been used in the life sciences earlier. This renaissance is based on the development of tomographic techniques,[14–19] and molecular sensing strategies using highly absorbing nanostructures for biomedical imaging.[20–24] The strength of the photoacoustic technique for scattering or non-radiating samples lies in the fact that it is not light that is detected but sound, which has a much longer wavelength. Although nanoparticle based photoacoustics is not fundamentally different from the many forms of photoacoustics developed so far, the small volume of nanoparticles introduces new phenomena that reflect the dominance of surface effects in these systems. The goal of this paper is to show that in many situations the photoacoustic response is not from the nano-sized absorber but the very low-absorbing surrounding material and therefore photoacoustic imaging can provide information about the nanoparticle surface, and the immediate environment surrounding the particle, which is of high importance to molecular imaging. The information that can be gained includes the interfacial resistance, as well as the thermodynamic characteristic of the environment.
Nanoparticle-augmented photoacoustic imaging is a non-invasive diagnostic imaging modality capable of assessing tissue properties with good spatial resolution and a sufficient depth that is hard to achieve with common optical imaging techniques.[19, 17, 12] Photoacoustics uses a broadband ultrasound transducer to visualize biological tissue by measuring the photoacoustic transient – an ultrasound wave produced by localized thermal expansion of tissue absorbing the pulsed electromagnetic radiation. Using the duration of the radiation pulses shorter than the time needed for the generated heat and stress to propagate through the radiated area, the amplitude of the photoacoustic signal in optically homogeneous materials is determined by the optical absorption and thermal properties of the material, as well as the fluence of the impinging radiation. Therefore, the contrast in the photoacoustic image of media, which are homogeneous on a length scale much larger than the thermal diffusion length, is mainly attributed to the optical absorption of the materials.
The introduction of plasmonic nanoparticles as photoacoustic contrast agents and their enhanced optical absorption due to surface plasmon polaritons brought molecular sensitivity to photoacoustic imaging.[27, 20–21] Small size, ease of surface functionalization and the ability to tune the peak absorption wavelength away from tissue absorption allows nanoparticles to molecularly target imaged tissue.[20, 22, 27–29] However, with the presence of nanoparticles, the medium containing the nano-sized strong optical absorbers in a weakly absorbing tissue is no longer optically homogeneous – the medium is heterogeneous on the scale of the thermal diffusion length of the tissue. Unlike in systems that are homogeneous or where the heterogeneities are much larger than the thermal diffusion length and therefore the photoacoustic signal is produced directly by the illuminated absorber,[30–31] the origin of the photoacoustic signal in the case of plasmonic nanoparticles is expected to depend on the heat transfer from the nanoparticle to the environment, and has not been explored in detail so far. In this study, we estimate the relative contribution from a nanoparticle and its liquid environment based on continuums mechanics, and show experimentally with gold nanospheres that the photoacoustic signal has the thermal signature of the colloidal solvent. Additionally, silica coating the gold nanospheres is shown to modulate the signal by influencing the heat transfer from gold particle to liquid.
For an optically heterogeneous medium the photoacoustic signal results in general from both the particles and the surrounding medium. Their relative magnitudes depend on the size of the absorber and the size of the heated region in the surrounding, as well as the temporal profile of the heating in the absorber and the surrounding. In order to show that the photoacoustic signal for a plasmonic nanoparticle in a liquid environment with good thermal conductivity is generated mostly from the surrounding, we consider a metallic sphere of radius a heated by light absorption to have the temperature Ts(t) (Figure 1). We use a continuums approach and make the following assumptions for our estimation: (1) the problem has spherical symmetry, so that transversal wave components can be neglected; (2) the sphere is small enough and the absorption time scales are long so that the sphere can be considered to be heated homogeneously, an assumption made frequently and justified by the very high thermal conductivity of the gold and other noble metals we are interested in; (3) the temperature profile in the surrounding at time t and position r from the center of the sphere can be written as , with a characteristic decay length of r0. We then solve the equation of motion in each material in spherical coordinates (r, θ, ):
for the radial displacement ur and the diagonal stress tensor with components σii = −KβT + K (εrr + 2ε θθ) + 2 μ[εii − (εrr + 2ε θθ)] /3; here β is the volumetric expansion coefficient, K and μ are the compressive and shear modulus, respectively, and ε is the strain tensor. With the longitudinal speed of sound and ρ the density, equation (1) then reduces to:
Introducing the potential ϕ so that ur /t = ϕ and , gives the photoacoustic wave equation that has to be solved in both media:
Boundary conditions at the interface between the sphere and surrounding fluid are that the radial stresses and displacements are continuous. Additionally, the displacement at the sphere center needs to be finite (and zero), and no incoming wave (traveling from infinity) exists for the surrounding. With ψ = rϕ, the Laplace operator can be reduced to a simple second derivative in equation (3) and it can be solved in (temporal) Fourier space with F(r,ω) = F(ψ), Gs(ω) = F(Ts(t)), and Gm(ω) = F(Tm(t)), and F denotes the Fourier transform of the potential (ψ ), temperature in the particle (Ts), and temperature in the environment (Tm), respectively. The solution is given in the supplement. Since the pressure is measured in the medium far from the nanoparticle, the contributions from the solid sphere and the liquid environment can be separated by setting the liquid and the metal volume expansion coefficients, respectively, to zero. The ratio of the two pressure contributions is then approximately:
Here and are the pressure contribution at frequency ω from the sphere (s) and the medium (m) measured far away from the nanoparticle in the medium, and rn = r0/a. If no interfacial heat resistance exists between the nanosphere and the fluid, then spectral temperatures at the interface are equal and Gs(ω) = Gm(ω). If there is interfacial heat resistance at the interface, then Gs(ω) > Gm(ω), and the temperature undergoes a jump at the interface.
In order to emphasize the importance of the heat transfer from the solid sphere to the fluid environment a Gaussian temperature profile of width δ for both the sphere and the medium, such as T(t) = (2π)−1/2 δ−1 exp(−t2/2δ2) can be assumed. The ratio of the contributions to the pressure measured some distance away from the source is then
From equation (5), it is clear that the spatial decay of the temperature profile, the temporal broadening due to a finite heat transfer through the interface, and the characteristics of the interface influence the balance of the two contributions. It should be pointed out that equation (5) assumes a temporal temperature profile in the sphere; for the situation of constant fluence, however, the temperature in the sphere is higher with large interfacial resistance than with lower resistance. This effect is included in the Gaussian profile to the order of approximation used here. For a gold sphere in water and rn ≈ 1, little broadening of a nanosecond laser pulse δm ≈ δs, and the range of frequencies used in photoacoustic imaging with ultrasound transducers (f = 2πω = 1 – 50 MHz), so that and βs /βm ≈ 0.2 at 20° C in aqueous solvent, and equation (5) predicts the pressure ratio between the sphere and the water medium of about one percent. Therefore, the photoacoustic signal is generated mostly in the fluid surrounding the spherical particle.
Water is ideally suited to study the origin of the photoacoustic signal from a nanoparticle colloid because the thermal expansion of water at atmospheric pressure vanishes at 3.98 °C, and therefore any photoacoustic signal generated by the heating of water at this temperature vanishes, while sound generated by the sphere will still be transported. Experiments were performed with polyethylene glycol stabilized gold spheres, and silica coated gold spheres with silica thickness of 6, 18 and 38 nm (Figure 2a, and Figure 2b) using a flow capillary suspended in a thermostated bath, 5 ns pulses from a Nd:YAG pumped OPO laser, and imaging with a single element focused 7.5 MHz ultrasound transducer (Figure 2c and Figure 2d). The PEGylated gold nanospheres were produced from citrate stabilized gold nanospheres by exchanging citrate with the biocompatible mPEG-thiol, and silica coated gold nanospheres were prepared by coating silica directly on PEGylated gold nanospheres via a modified Stöber method.[33–34, 25] The peak optical density of each sample was adjusted to 1.5 to achieve equal absorption, and the laser was tuned to the peak wavelength of each sample. We verified with Mie theory that at the extinction peak the scattering contribution to the extinction cross-section due to the addition of a silica shell was well below 10% for silica shells between 0 nm and 40 nm thickness (see supplement for the experimental extinction spectra of gold nanosphere solutions). A laser fluence of 5 mJ/cm2 was chosen to prevent bubble formation in organic solvents, and to avoid damage to the nanoparticles. We have not found that the fluence used here causes any damage to the bare or coated nanospheres, and the PA signal remained stable over at least 128 laser pulses.
Figure 3a shows that the photoacoustic signal of PEGylated and silica-coated gold nanospheres in aqueous solution indeed falls below noise level around 4 °1C. The photoacoustic signal increases nearly linearly with temperature with a slope that depends on the silica thickness. The slope is lowest for the PEGylated gold nanospheres and is maximal for 18 nm silica. An approximately linear dependence of the photoacoustic signal with temperature owing to the raise of the Grüneisen parameter of water has been seen for gold nanospheres in water above 25 °C, and for homogeneous aqueous absorbers without nanoparticles, where the photoacoustic signal also disappeared at 4 °C. The photoacoustic signal therefore originates mostly from the surrounding water, and the signal from the absorber is below the noise level even if the size of the nanoparticle is increased 2-fold by the silica coating. The relative contributions of the nanoparticle, including shell, and the water can be estimated by fitting the measured photoacoustic signal to the equation Ptotout = PNP+bΓH2O(T), assuming that the nanoparticle contribution PNP, including the PEG-thiol capping and binding layer, is constant, and the Grüneisen parameter’s temperature dependence Γ(T) is dominated by that of the expansion coefficient, as demonstrated in Figure 3d. The proportionality factor b is a function of shell thickness, and the nanoparticle contribution to the signal PNP, which is below noise level for the fluence used here, also varies linearly with thickness (see Figure S1 in the supplement). The good fit to the thermal expansion of water, the magnitude of the nanoparticle contribution, and its linear dependence with thickness all confirm the trends of the theoretical estimate for a simple sphere with and without a shell (see supplement).
The increase in slope reflects the increase in amplitude of the photoacoustic signal at 20 °C, which we have recently reported for the silica-coated nanorods. The amplitude of the photoacoustic signal is determined by the Grüneisen parameter of the solvent, as well as the details of the profile of heat transferred from the nanoparticles to the liquid, which depends on the optical absorption of the nanoparticles, the laser fluence, and the interfacial heat conductivities between the two materials. We had concluded there that the signal increase due to silica coating was mainly due to changes in the interfacial heat resistance between the gold and water, and an enhanced efficiency of the heat transfer between the solid and water due to the coating, since all other parameters were approximately constant.
The effect of the heat transfer between nanoparticle and solvent on the photoacoustic signal was further investigated by replacing the aqueous solvent with two non-polar solvents with different thermal properties: polydimethylsiloxane (silicon oil) and toluene. The surface functionalization of the bare and silica-coated nanoparticles was changed to be able to suspend them in non-aqueous solvent. As shown in Figure 2b, the citrate of the citrate-stabilized nanospheres was replaced by octadecane-1-thiol, and the silica surface was functionalized with trimethoxy(octadecyl)silane by co-polymerization (see supplement for details).
Silicon oil (Figure 3b) shows three different features of the photoacoustic signal compared to water: a non-zero signal at 4 °C, a bi-phasic and mostly opposite slope of the temperature dependence, and a decrease of the signal amplitude with silicon coating thickness. The signal follows the bi-phasic profile of the Grüneisen parameter of pure silicon oil, measured using the intrinsic absorption of the oil without nanoparticles (Figure 3d), signaling that the oil is the major signal-generating medium despite the differences to water. However, the temperature dependence is more complex, as the signal is not proportional to the measured Grüneisen coefficient of oil, and there is a leveling off at high temperatures. One possible reason may be residues of water in the silica shell left from the solvent exchange. The fact that silica coating the nanospheres reduces the photoacoustic signal in the oil over most of the temperature range is in agreement with our recent findings at 20 °C.
Toluene is known to have a very high interfacial resistance with citrate coated gold nanospheres and could therefore be a good test of the role that thermal interfacial transfer has on the balance of the contributions from gold and fluid to the signal. Figure 3c shows the photoacoustic signal as a function of temperature in toluene. Between 15 and 45 °C the signal follows the trend of the temperature dependent Grüneisen parameter of toluene (Figure 3d) with a slight negative slope. A small signal amplification due to the addition of a silica shell can be seen, which differs from our finding for nanorods at 25°C. This difference may come in part from changes in the absorption cross section of gold spheres when a silica shell is added, which are not present for the longitudinal resonance of nanorods.
In conclusion, the photoacoustic signal of a solution of gold nanospheres in different solvents shows that the dominant photoacoustic signal is produced by its surrounding medium and hence carries the information of the thermal properties of the surrounding medium. This is despite the fact that the heat is generated by the absorption of light in the nanoparticle and therefore carries also the signature of the absorptivity of the gold nanoparticle and its interface to the solvent by determining the amount and time-dependence of the heat provided to the solvent. Photoacoustics can therefore be a new tool for investigating the thermal properties and transport through nanoparticle interfaces.