## 1.

## Introduction

Accurate knowledge of the optical properties of turbid media, such as tissue, is important in many clinical applications that use light for diagnosis and treatment. The optical properties depend on the physiological state of the tissue and are influenced by many parameters. One such parameter is temperature, which can substantially alter measured optical properties. The effect can be quantified as a temperature coefficient, the rate of increase (or decrease) in optical absorption $\left(\alpha \right)$ and scattering $\left(\beta \right)$ coefficients with temperature.

Temperature-induced change in the optical properties can, for example, limit the application of photodynamic therapy.^{1} Other applications include the estimation of temperature from the change in optical, property measured. Temperature can be estimated from changes in the absorbance^{2} however, in turbid media, variation in scatter introduces an uncertainty in the path length, making it difficult to reliably estimate temperature from optical extinction coefficients. Uncoupling the effect of scattering on absorption may reduce this obscurity, offering the potential for accurate measurements of temperature in turbid solutions.

Kelly
^{3} measured the near-infrared (NIR) spectra of chicken, bovine, and porcine tissues from
$17\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}45\xb0\mathrm{C}$
with a spectrophotometer and reported that the measurements can predict tissue temperature with a standard error less than
$0.2\xb0\mathrm{C}$
. Laufer
^{4} measured the change in optical properties of human dermis and subdermis from
$25\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}40\xb0\mathrm{C}$
. They observed positive and negative temperature coefficients of scattering for dermis and subdermis, respectively, but no change in absorption with temperature.

A common phantom to simulate optical properties of tissue is Intralipid. Temperature-induced changes in the optical properties of Intralipid have been described previously. Kakuta
^{2} reported the temperature dependence of absorbance at
$1440\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
using a spectrophotometer and proposed a method to measure the temperature of turbid aqueous solutions, but they did not recover the absorption coefficient. McGlone
^{5} determined the temperature dependence of absorption and reduced scattering coefficients of Intralipid using continuous wave photon migration measurements from
$700\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}1000\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
. The authors relied on measurement of attenuation close to the source to separate scattering and absorption. The changes in absorption coefficients due to temperature were attributed to water,^{6, 7} the main absorbing species in Intralipid.^{8, 9} Kakuta found a discrepancy in absorption between heating and cooling cycle measurements and linked it to irreversible changes in the scattering properties of the micelles. McGlone observed a lack of repeatability, during heating and cooling, in temperature-induced changes in the reduced scattering properties.

The frequency-domain technique uncouples absorption from scattering^{10, 11} when determining the optical properties of turbid media, which potentially allows more accurate recovery of optical coefficients.^{12} Our work is a preliminary step toward monitoring these subtle changes. We have investigated the optical properties of a tissue-simulating phantom and present, for the first time, frequency-domain measurements of the temperature dependence of the absorption
$\left({\mu}_{a}\right)$
and reduced scattering
$\left({\mu}_{s}^{\prime}\right)$
coefficients of Intralipid in the wavelength range from
$710\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}850\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
. We have used a 1.8% (w/w) solution of Intralipid as the sample to simulate the reduced scattering properties of human brain and breast tissues.^{13} Measurements were made in the physiologically relevant temperature range of
$30\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}40\xb0\mathrm{C}$
while heating and cooling the liquid.

## 2.

## Method

Optical properties can be determined by measuring the amplitude attenuation and phase shift of intensity-modulated light passing through a turbid medium.^{14} The propagation of intensity-modulated light through turbid media has been extensively discussed in the literature.^{15, 16, 17, 18}

We employed the multidistance method, where the amplitude and phase of a photon density wave is measured at different separations from the source. An outline of our experimental setup is shown in Fig. 1
and has been described in detail previously.^{19} Using this configuration, we measured absorption
$\left(0.009\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}0.045\phantom{\rule{0.3em}{0ex}}{\mathrm{cm}}^{-1}\right)$
to an agreement with reported literature of better than 10% in the presence of reduced scattering coefficients between 10 and
$35\phantom{\rule{0.3em}{0ex}}{\mathrm{cm}}^{-1}$
. Briefly, the instrumentation employed a tunable Ti: Sapphire laser as the light source, which was intensity modulated by an acousto-optic modulator. A lock-in amplifier in conjunction with an avalanche photodiode (APD) module constituted a phase-sensitive detection system. Optical fibers, suspended from a translation stage in the middle of the Intralipid solution, were used to couple the light into and out of the medium.

The temperature of the medium was maintained using a closed-loop control system incorporating a heated plate (VELP Scintifica, Italy) equipped with a magnetic stirrer and a semiconductor temperature sensor (LM35, National Semiconductor, Santa Clara, California) immersed in the solution. The temperature sensor was located so that it would not affect the photon density wave. The temperature sensor and the solution heater were connected via a data acquisition device (U 12, LabJack Corporation, Lakewood, Colorado) to a computer and controlled by a Labview (National Instruments, Austin, Texas) program. A loose polythene cover minimized water evaporation during the experiment. A 1.8% (w/w) solution of Intralipid was obtained by diluting $500\phantom{\rule{0.3em}{0ex}}\mathrm{mL}$ of Intralipid-20% (Pharmaco, New Zealand) in $5\phantom{\rule{0.3em}{0ex}}\mathrm{L}$ of distilled water.

Amplitude and phase measurements were initially collected, at $10\text{-}\mathrm{nm}$ intervals, between 710 and $850\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ at three different set-point temperatures (30, 35, and $40\pm 0.5\xb0\mathrm{C}$ —the “low-temperature resolution” data) as the distance between the source and detector fibers was increased from $10\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}30\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ ( $0.8\text{-}\mathrm{mm}$ steps). At each stage-position, five amplitude and phase measurements were recorded to estimate the errors in our estimated absorption and reduced scattering coefficients. The temperature control system was allowed to stabilize for at least $30\phantom{\rule{0.3em}{0ex}}\mathrm{min}$ before collecting each set of measurements. Additional amplitude and phase measurements—the “high-resolution” data—were made at 740, 800, and $840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ at room temperature $(22\xb0\mathrm{C})$ , then as the solution was first slowly heated (over $2\phantom{\rule{0.3em}{0ex}}\mathrm{h}$ ) from $30\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}40\xb0\mathrm{C}$ , then again as it cooled (over $4\phantom{\rule{0.3em}{0ex}}\mathrm{h}$ ), and finally again at room temperature. These measurements were collected every $1.1\pm 0.1\phantom{\rule{0.2em}{0ex}}\xb0\mathrm{C}$ between $30\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}40\phantom{\rule{0.2em}{0ex}}\xb0\mathrm{C}$ over the heating and cooling regimes. Replicate heating and cooling measurements were collected at $740\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ and $840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ with similar solutions prepared, on different days; a single set of measurements was recorded at $800\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ . Last, optical measurements were recorded when the solution reached $22\xb0\mathrm{C}$ after a cooling cycle, then again at $22\xb0\mathrm{C}$ the next day.

As Intralipid is known to decay over time,^{20} a fresh solution was prepared for each set of measurements.

## 3.

## Results

Figure 2
shows the absorption coefficient of Intralipid, calculated from the amplitude and phase measurements, over the wavelength range
$710\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}850\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
measured at the three different temperatures. The absorption coefficient from previous measurements made by us at
$22\pm 0.5\xb0\mathrm{C}$
is included.^{19} The error bars at each measurement illustrate the measurement uncertainty (
$\pm 1$
standard deviation). The uncertainty in the optical coefficients was estimated from the replicate error in phase and amplitude measurements through error propagation.^{21} Arrows on the plot show the trend in absorption coefficient with temperature. The absorption coefficient increases around 740 and
$840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
and decreases around
$800\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
.

The reduced scattering coefficient, again calculated from the amplitude and phase measurements, is illustrated in Fig. 3
for the three temperatures measured in this study and at
$22\xb0\mathrm{C}$
(Ref. 19). As predicted by Mie theory,^{22} the reduced scattering coefficient decreases, nearly linearly, as wavelength increases. Moreover, these measurements show that the reduced scattering coefficient decreases as temperature increases.

The temperature trend is more clearly illustrated for the absorption and reduced scattering coefficients using the data measured continuously as the sample was heated and cooled (Fig. 4
and Fig. 5
, respectively). Each plot illustrates an approximately linear relationship between the optical coefficient and temperature at 740 and
$840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
for two replicate measurements during heating and cooling of the sample. Results at
$22\pm 0.4\xb0\mathrm{C}$
are included for comparison from earlier work.^{19} Each data point between 30 and
$40\xb0\mathrm{C}$
represents the average of five measurements with one standard deviation error bars estimated by propagating the noise in the amplitude and phase measurements through the calculations for the optical coefficients.

As Fig. 4 shows, the absorption coefficient increases with temperature at both 740 and
$840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
. This is consistent with the change in absorption coefficient of pure water,^{7} the main absorbing species in Intralipid.^{8} Lipid, the other absorbing species in Intralipid, has an absorption coefficient of less than
$0.0005\phantom{\rule{0.3em}{0ex}}{\mathrm{cm}}^{-1}$
but is present only at 1.8% by volume, so its contribution will be negligible.^{9} The absorption coefficients over the heating and cooling period match well at
$740\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
. However, at
$840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
, absorption was generally higher during the cooling period. This may be due to a variation in dilution across the different batches of Intralipid.

A temperature coefficient,
$\alpha $
, has been estimated using a linear, least-squares fit to the measurements made in the heating and cooling periods. The discrepancy between heating and cooling at
$840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
had no significant effect on the slope, so data from the two periods has been combined in all cases. These measurements are in good agreement with results collected at
$22\xb0\mathrm{C}$
both during this work and in previous work.^{19}

Generally, the reduced scattering coefficient decreased as temperature increased at 740 and
$840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
(Fig. 5). However, at
$740\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
, a hysteresis was observed whereby the reduced scattering coefficient decreased during heating yet remained reasonably constant during cooling, even down to
$22\xb0\mathrm{C}$
. This behavior was observed in both replicates but was not a result of the self cooling, as the measurement time for a single set (a detector fiber scan from
$10\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}30\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$
) is about
$7\phantom{\rule{0.3em}{0ex}}\mathrm{min}$
and the change in temperature in this interval is only of the order of
$0.2\xb0\mathrm{C}$
. Earlier data at
$22\xb0\mathrm{C}$
and a subsequent measurement at
$22\xb0\mathrm{C}$
(green square) on the following day (
$15\phantom{\rule{0.3em}{0ex}}\mathrm{h}$
later), indicate that the scattering coefficient recovers to its original value eventually. This behavior was observed in both replicates at
$740\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
, but at
$840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
, although the slope is slightly lower during the cooling period, no significant difference was observed between the heating and cooling periods. Kakuta
^{2} and McGlone
^{5} also observed different behavior in their measurements during heating and cooling. Kakuta observed discrepancy in their absorbance measurement and postulated that it may be due to the permanent structural changes in the micelles during heating. McGlone also observed different scattering behavior during heating and cooling. Because of this hysteresis, the temperature coefficient,
$\beta $
, for the scattering coefficient was estimated using only data from the heating period, with a linear least-squares fit.

## 4.

## Discussion

A positive temperature coefficient of absorption
$\left(\alpha \right)$
of
$(2.0\pm 0.02)\times {10}^{-4}\phantom{\rule{0.3em}{0ex}}{\mathrm{cm}}^{-1}\phantom{\rule{0.2em}{0ex}}{\mathrm{K}}^{-1}$
at
$740\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
and
$(2.3\pm 0.3)\times {10}^{-4}\phantom{\rule{0.3em}{0ex}}{\mathrm{cm}}^{-1}\phantom{\rule{0.2em}{0ex}}{\mathrm{K}}^{-1}$
at
$840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
was calculated from the data plotted in Fig. 4. At
$800\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
(data not shown), we calculated a temperature coefficient of
$(-7\pm 1)\times {10}^{-5}\phantom{\rule{0.3em}{0ex}}{\mathrm{cm}}^{-1}\phantom{\rule{0.2em}{0ex}}{\mathrm{K}}^{-1}$
. The near-infrared absorption features arise from overtones and combination bands of the vibrations of the water molecule.^{23} At
$837\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
, the main contribution is a combination of symmetric and asymmetric stretch and bending vibrations. The
$739\text{-}\mathrm{nm}$
absorption peak is caused by symmetric and asymmetric stretch vibrations. The sensitivity of water absorption to temperature arises from changes to the strength of the hydrogen bonds and microscopic changes in the structure of the water.^{7, 24}

There is some disagreement in the literature over the precise value of water’s near-infrared temperature coefficient. Our results at 740 and
$840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
are in good agreement with those reported by Hollis^{25, 26} for water but are 40% higher than the temperature coefficient of water reported by Langford (Fig. 6
).^{7} The temperature coefficient of absorption for Intralipid measured by McGlone
^{5} are about 50% lower than our results. At
$800\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$
, our result is more consistent with the data reported by Langford
^{7} and McGlone, ^{5} underestimating Hollis’ observation by 75%.

The temperature coefficient of Intralipid absorption was estimated in $10\text{-}\mathrm{nm}$ steps between 710 and $850\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ using the data plotted in Fig. 3. At most wavelengths, the temperature coefficient was estimated from absorption measured at 22, 30, 35, and $40\xb0\mathrm{C}$ (squares). At 740, 800, and $840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ , the results, obtained from high-temperature resolution measurements, discussed earlier are plotted (triangles). The error bars indicate 68% confidence interval of the regression values. Large errors—at $820\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ , for example—occur when there is a significant variation in laser power during a measurement. Overall, our results agree in magnitude and spectral shape with values reported in the literature.

The prediction intervals in Fig. 4 give an indication of the temperature sensitivity of our setup. The minimum temperature change observable within the error limits is $4\xb0\mathrm{C}$ and $6\xb0\mathrm{C}$ at 740 and $840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ , respectively. This is probably insufficient for applications requiring precision measurement. Moreover, the presence of other absorbing species, such as hemoglobin, in biological samples will likely demand even greater accuracy. However, we believe that the largest source of error lies in the precision of the phase measurements and the short-term stability of the source. We believe that more precise temperature measurements can be achieved by including active feedback to control the amplitude stability of our laser source.

The temperature coefficient of the reduced scattering coefficient $\left(\beta \right)$ was estimated, as for absorption, using linear regression. However, as we found a significant difference in scattering behavior during heating and cooling, only heating data was used. The results are plotted in Fig. 7 . At most wavelengths, the temperature coefficient was estimated from reduced scattering measured at 22, 30, 35, and $40\xb0\mathrm{C}$ (squares); at 740, 800, and $840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ , approximately 10 measurements between 30 and $40\xb0\mathrm{C}$ , as well as $22\xb0\mathrm{C}$ were used (triangles). A linear regression line through the reduced scattering temperature coefficient shows a general decrease in temperature sensitivity as wavelength increases.

The scattering coefficient of a turbid media is affected by temperature in two ways. First, as temperature increases, the volume of the medium expands, diluting the scattering effect.^{5} However, more significantly, the refractive index of water also decreases with temperature,^{27} leading to a greater mismatch between the scattering particles and the absorbing medium. We have modeled the temperature coefficient of the reduced scattering coefficient of Intralipid based on Mie scattering and changes of the refractive index with temperature. We used published particle-size distribution data for Intralipid^{9} and assumed a refractive index temperature sensitivity of
$3.54\times {10}^{-4}\phantom{\rule{0.3em}{0ex}}{\mathrm{K}}^{-1}$
, as that value is typical of the major fatty acids components of soya bean oil and similar derivatives.^{28} The result is included in Fig. 7 as a green dashed line for comparison. McGlone measured temperature coefficients for scattering in a 1.2% Intralipid solution by continuous wave measurements. Their data, scaled to 1.8%, are illustrated as well. Within the large uncertainty, our measurements agree with theory and other measurements.

## 5.

## Conclusion

We have used frequency-domain, photon-migration spectroscopy to separately measure absorption and reduced scattering temperature coefficients in a turbid medium for the first time. The optical temperature coefficients of the liquid, tissue-simulating phantom 1.8% Intralipid were estimated from $710\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}850\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ in the physiologically relevant $30\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}40\xb0\mathrm{C}$ temperature range. At 740, 800, and $840\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ (key wavelengths in the water absorption temperature coefficient), we measured absorption temperature coefficients of $(2.0\pm 0.02)\times {10}^{-4}$ , $(-7\pm 1)\times {10}^{-5}$ , and $(2.32\pm 0.03)\times {10}^{-4}\phantom{\rule{0.3em}{0ex}}{\mathrm{cm}}^{-1}\phantom{\rule{0.2em}{0ex}}{\mathrm{K}}^{-1}$ and reduced scattering temperature coefficients of $-0.18\pm 0.02$ , $-0.09\pm 0.01$ , and $-0.09\pm 0.02\phantom{\rule{0.3em}{0ex}}{\mathrm{cm}}^{-1}\phantom{\rule{0.2em}{0ex}}{\mathrm{K}}^{-1}$ , respectively. The temperature coefficients observed for absorption closely follow the published results for water, as expected. The temperature coefficients observed for reduced scattering exhibited a general decrease as wavelength increased, consistent with predictions made by Mie theory and literature.

We observed hysteresis behavior in the reduced scattering measurements made at $740\phantom{\rule{0.3em}{0ex}}\mathrm{nm}$ that did not appear to occur at higher wavelengths. The reduced scattering coefficient decreased with temperature but did not increase again as the solution was cooled until sometime later. This hysteresis did not affect the absorption measurements. This result suggests that similar observations of this behavior reported in extinction measurements are caused by changes in the scattering micelles.

These results suggest that temperature measurement, with a precision of $\pm 4\xb0\mathrm{C}$ , in turbid media is feasible based on the absorption coefficient’s temperature coefficient. However, an improvement in precision will be required for many applications.

## Acknowledgments

We acknowledge financial support from the University of Waikato and the Foundation for Research Science and Technology.

## References

*ex vivo*human dermis and subdermis,” Phys. Med. Biol.0031-9155 43(9), 2479–2489 (1998). 10.1088/0031-9155/43/9/004 Google Scholar