The concentration of water in tissue reflects both physiological and pathological properties. Many techniques have been used to detect and measure water content in tissue. For instance, magnetic resonance imaging (MRI) was traditionally used to image water content in the brain,1, 2 bone,3, 4 and other tissues,5, 6 and diffusion-weight MRI (DW-MRI) has been used to map water mobility in brain and muscle tissues.7, 8 Recently, diffuse optical methods have been applied to measure the increased water concentration in tumors.9 The bound water fraction, which is generally associated with tissue pathophysiology,10 affects the molecular vibrational states of water. Diffuse optical spectroscopy (DOS) can measure the spectral changes due to bound water.11 Although DOS is accurate and sensitive to water state,11 its spatial location is poor.
As an emerging imaging technique, photoacoustic tomography (PAT) has been successfully implemented in biomedicine.12, 13, 14, 15 Combining high optical contrast and ultrasonic detection, PAT overcomes the depth limitation of other high-resolution optical imaging methods, and it is also free from speckle artifacts.16, 17, 18
PAT has been implemented to image various tissues19 from centimeter-large breast tumors to micrometer-large single red blood cells. Image contrast has been provided by many endogenous and exogenous optical absorbers, such as hemoglobin, melanoma, and various natural and artificial contrast agents. Since water absorbs electromagnetic (EM) waves in the radio frequency (rf) and microwave spectra, it has served as the absorber in rf-based photoacoustic (PA) detection, which is also referred to as thermoacoustic (TA) detection. In addition, the absorption spectrum of water has a peak absorption coefficient of around , with a FWHM of (Ref. 20). However, to our knowledge, water has never been detected by the PA imaging using a near-IR (NIR) light source. In this paper, we explore the potential of laser-based PA detection of water. From here on, the EM source is always a laser, unless otherwise noted.
The next section briefly introduces the theoretical basis. Then, we describe how the absorption spectra of water-ethanol mixtures at various concentrations were measured by using the PA technique. After that, both phantom and tissue experiments to image water are detailed. Finally, conclusions and a discussion are provided.
Theoretical Basisdenotes the acoustic pressure at position and time , is the speed of sound in the medium, denotes the thermal coefficient of volume expansion, denotes the isobaric specific heat capacity, and is the heating function defined as the thermal energy converted at and per unit volume and time. For optical absorption, the heating function generally equals , where is the percentage of energy that is converted into heat, is the optical absorption coefficient, and is the optical fluence rate.
In general, the initial pressure rise at immediately after illumination by a short laser pulse is given by23is defined as the Grueneisen coefficient (dimensionless), and is the optical fluence. In many cases, is approximately equal to 1. The initial pressure , serving as the source of the propagating PA pressure wave, depends on not only the illumination fluence and absorption coefficient, but also on the Grueneisen coefficient. In addition, the absorption coefficient is generally spectrally dependent.
Experimental and Theoretical Methods
Figure 1 shows the experimental setup of our absorption spectrum measurement system. The laser source was a tunable laser [Vibrant (HE) 355 I, OPOTEK], which provided laser pulses at a repetition rate of with a wavelength-tuning range of . Water-ethanol mixtures with the same volume in a thin-walled plastic container were placed in front of a unfocused ultrasonic transducer (V323, Panametrics), which was immersed in a tank of mineral oil.24 The distance from the bottom surface of the mixture to the front surface of the transducer was . Since the detected signal amplitude is proportional to the initial pressure, here we used the peak-to-peak PA signal value to represent . The signal from the transducer was amplified by by a pulse amplifier (5072 PR, Panametrics) and then recorded by a data acquisition card (CompuScope 14100, GaGe, Lockport, Illinois) at a sampling rate of . Both the mixture and the container generated PA signals, but the discontinuity at their interface helped to clearly distinguish between the two signals according to their different arrival times.
Six mixtures of water-ethanol with different volume concentrations of water were made: 100 (pure deionized water), 80, 60, 40, 20, and 0% (ethanol only). Since the ethanol used was 95% pure (Pharmco-Aaper, 190 proof) with the remaining 5% being water, the actual volume concentrations of water were 100, 81, 62, 43, 24, and 5%, respectively. The total volume of each mixture was set to . PA measurement can quantify the relative initial pressure as a function of the optical wavelength at each of these concentrations. We calibrated the power of the laser in the spectral range from at intervals. The laser pulse energy was recorded by a powermeter (Melles Griot broadband power/energy meter 13PEM001), and 20 laser pulses were averaged at each wavelength.
We calculated dependence of the Grueneisen coefficient on the mixture composition. The speed of sound in the water-ethanol mixtures with different water concentrations were obtained from Ref. 25. The thermal coefficient of volume expansion and specific heat capacity of deionized water and pure ethanol were taken from Ref. 26. It was assumed that the thermal coefficient of volume expansion is additive with volume. That is, , where and denote the volumes of water and ethanol, respectively; is the thermal coefficient of the mixture; and and are the thermal coefficients of water and ethanol, respectively. We further assumed that is also additive with mass. So , where and denote the mass concentrations of water and ethanol, respectively. The mass could be calculated with the volume and density shown in Ref. 27. The specific heat capacities of the mixture, water, and ethanol are , , and , respectively. Then, the Grueneisen coefficient of the mixture was calculated using .
Figure 1 shows the experimental setup of our PAT system, which was used to image water at various concentrations on the basis of the dependence of the PA signal on the water concentration. The tunable laser operated at the wavelength. A pulse energy was used to excite PA signals. An ultrasound transducer (842-300, GE) with a center frequency of was immersed in the mineral oil in the tank. A step motor controlled by a computer drove the transducer to scan around the phantom through 120 positions evenly distributed along the horizontal scanning trajectory, while the PA signals were detected by the transducer. The same amplifier and data acquisition card shown in Fig. 1 were used. The image reconstruction was implemented in the computer after data acquisition using the image reconstruction algorithm described in Ref. 28.
Two samples were imaged using this PAT system. For sample 1, a gel cylinder made of an agar-water mixture29 (4% agar) with three cylindrical holes was used. The diameter of each hole was . One hole was filled with deionized water, another with a 40% ethanol and 60% deionized water mixture, and the third hole was left in air [Fig. 5 in Sec. 4]. For sample 2, a piece of porcine fat was used. A hole was drilled in the fat and filled with 2% agar gel [Fig. 6 in Sec. 4].
Spectral Measurement of Water-Ethanol Mixtures
Figure 2 shows the estimated Grueneisen coefficients at various water concentrations of the mixtures. A line, fitted with the coefficient of simple determination , was used to provide the Grueneisen coefficient values in estimating the optical absorption coefficient of each mixture.
By assuming that the absorption coefficient is additive on the water concentration and using the linear fit of the Grueneisen coefficient, we can get a quadratic dependence of the PA signal on the water concentration, as seen in Fig. 2. The figure plots the normalized PA signal to that of the pure deionized water. This quadratic relation is spectrally independent. The nonmonotonic curve in Fig. 2 implies that the water fraction of a water-ethanol mixture might not be determined uniquely according to the PA measurement alone. In general, one PA signal corresponds to two possible water concentrations. We measured relative PA signals at the wavelength for five mixtures, and estimated the water concentration in Table 1 . For the mixture with 50% water concentration, the ratio of its PA signal to that of the deionized water is greater than the maximum in the curve shown in Fig. 2. As a result, the solutions for the water concentration from multiple PA measurements were complex with a constant real part but a varying imaginary part. Only the real part is used as the photoacoustically measured water concentration, which matches the concentration for the peak of the curve in Fig. 2.
Comparison of water concentrations with corresponding PA-derived fractional values (five means; for each mean, 100 PA signals were used).
|Preset water concentration||10%||30%||50%||70%||90%|
|Photoacoustically measured waterconcentration (solution 1)||50.6%|
|Photoacoustically measured waterconcentration (solution 2)||50.6%|
The relative absorption spectra were computed with Eq. 2 and then converted to absolute values [Fig. 3 ] by calibrating with the absorption coefficient of deionized water at measured by a spectrophotometer (Cary 5E, Varian, Walnut Greek, California). For Eq. 2, we used the Grueneisen coefficient shown in Fig. 2 and the calibrated optical fluence; the latter is assumed to be uniformly distributed owing to the ground glass shown in Fig. 1. For comparison, the absorption spectra were also measured with the spectrophotometer [Fig. 3]. The relative errors between the PA and spectrophotometric data are plotted in Fig. 3, which shows good accordance except at , where the absorption coefficient is small. In addition, the PA estimations for mixtures with less water were less accurate, possibly because of the inaccuracy in estimating the Grueneisen coefficient.
Because ethanol has negligible absorption in comparison to water at the spectral region of interest and the absorption coefficient is the product of the absorption cross section and the number density of the absorbers,30 the absorption coefficient of each mixture is proportional to the water fraction. Therefore, we can use the absorption coefficients of deionized water and the ethanol stock solution (5% water) derived from the PA measurement to estimate the absorption coefficient of the known ratio of water and ethanol. The estimates and the PA experimental results are plotted in Fig. 4 , and the relative errors are shown in Fig. 4. We can see that the estimations were accurate, especially at the wavelength of and at high water concentrations, where the absorption coefficient was relatively large. Yet the estimation for the 24% water concentration was less accurate (12.5% error at the wavelength) than that for other mixtures (less than 5% error at the wavelength).
Imaging Water Content in Phantoms and Tissue
Figure 5 shows a PAT image of objects with different water concentrations at . The negative value is partially due to the limited bandwidth of the ultrasonic detection system. The reconstructed signal was integrated within the square marked in Fig. 5. The 40% water to 60% ethanol mixture yielded a higher PA signal than deionized water; about 1.65 times as great, which agrees with the ratio of 1.84 predicted from Fig. 2.
The tissue experiment [Fig. 6 ] indicated that water generates nine times stronger PA signals than fat at the wavelength [Fig. 6]. This difference is possibly due to differences in both the absorption coefficient and the Grueneisen coefficient between water and fat. Thus, PAT has a great potential to detect localized water in fat.
We applied an NIR laser-based PA technique for the first time to image water. Absorption spectral measurements demonstrated that we could distinguish water signals of various concentrations. A tissue experiment showed the potential of water imaging in vivo. Different from previous applications in medical imaging, water imaging reflects not only the absorption coefficients but also the Grueneisen coefficients. In other words, different water concentrations lead to differences in both the optical absorption coefficient and the Grueneisen coefficient. PAT can image the differences in water content, especially in tissue with low water concentration such as fat, although the impact of absorption of skin and blood should be considered in animal experiments in the future.
This work was sponsored in part by National Institutes of Health Grants R01 EB000712, EB000712A2S1, R01 EB00071207S2, R01 EB008085, R01 CA113453901, U54 CA136398, and 5P60 DK02057933. L.W. has a financial interest in Microphotoacoustics, Inc., and Endra, Inc., which, however, did not support this work.