Time-domain thermoreflectance

Time-domain thermoreflectance (TDTR) is a method by which the thermal properties of a material can be measured, most importantly thermal conductivity. This method can be applied most notably to thin film materials (up to hundreds of nanometers thick), which have properties that vary greatly when compared to the same materials in bulk. The idea behind this technique is that once a material is heated up, the change in the reflectance of the surface can be utilized to derive the thermal properties. The reflectivity is measured with respect to time, and the data received can be matched to a model with coefficients that correspond to thermal properties.

Experiment setup

The technique of this method is based on the monitoring of acoustic waves that are generated with a pulsed laser. Localized heating of a material will create a localized temperature increase, which induces thermal stress. This stress build in a localized region causes an acoustic strain pulse. At an interface, the pulse will be subjected to a transmittance/reflectance state, and the characteristics of the interface may be monitored with the reflected waves. A probe laser will detect the effects of the reflecting acoustic waves by sensing the piezo-optic effect.

The amount of strain is related to the optical laser pulse as follows. Take the localized temperature increase due to the laser,

${\displaystyle \Delta T(z)=(1-R){\frac {Q}{C(\zeta A)}}\exp(-z/\zeta ),}$

where R is the sample reflectivity, Q is the optical pulse energy, C is the specific heat (per unit volume), A is the optical spot area, ζ is the optical absorption length, and z is the distance into the sample (Ref A). This temperature increase results in a strain that can be estimated by multiplying it with the linear coefficient of thermal expansion of the film. Usually, a typical magnitude value of the acoustic pulse will be small, and for long propagation nonlinear effects could become important. But propagation of such short duration pulses will suffer acoustic attenuation if the temperature is not very low (Ref B). Thus, this method is most efficient with the utilization of surface acoustic waves, and studies on investigation of this method toward lateral structures are being conducted.

To sense the piezo-optic effect of the reflected waves, fast monitoring is required due to the travel time of the acoustic wave and heat flow. Acoustic waves travel a few nanometers in a picosecond, where heat flows about a hundred nanometers in a second.[1][2] Thus, lasers such as titanium sapphire (Ti:Al2O3) laser, with pulse width of ~200 fs, are used to monitor the characteristics of the interface. Other type of lasers include Yb:fiber, Yb:tungstate, Er:fiber, Nd:glass. Second-harmonic generation may be utilized to achieve frequency of double or higher.

The output of the laser is split into pump and probe beams by a half-wave plate followed by a polarizing beam splitter leading to a cross-polarized pump and probe. The pump beam is modulated on the order of a few megahertz by an acousto-optic or electro-optic modulator and focused onto the sample with a lens. The probe is directed into an optical delay line. The probe beam is then focused with a lens onto the same spot on the sample as the probe. Both pump and probe have a spot size on the order of 10–50 μm. The reflected probe light is input to a high bandwidth photodetector. The output is fed into a lock-in amplifier whose reference signal has the same frequency used to modulate the pump. The voltage output from the lock-in will be proportional to ΔR. Recording this signal as the optical delay line is changed provides a measurement of ΔR as a function of optical probe-pulse time delay.[3]

Modeling materials

The surface temperature of a single layer

The frequency domain solution for a semi-infinite solid which is heated by a point source with angular frequency ${\displaystyle \omega }$ can be expressed by the following equation.[4]

${\displaystyle g(r)={\frac {exp(-qr)}{(2\pi \Lambda r)}}}$ where ${\displaystyle q^{2}=(iw/d)}$ (1)

(Λ: thermal conductivity of the solid, D: thermal diffusivity of the solid, r: radial coordinate)

In a typical TDTR experiment the co-aligned laser beams have cylindrical symmetry, therefore the Hankel Transform can be used to simplify the computation of the convolution of equation (1) with the distributions of the laser intensities.

(The Hankel transform is an integral transform equivalent to a two-dimensional Fourier transform with a radially symmetric integral kernel)

Here g(r) is radially symmetric and by the definition of Hankel transform using Eq. (1),

${\displaystyle G(k)=2\pi \int _{0}^{\infty }g(r)J_{0}(2\pi kr)r\,dr={\frac {1}{\Lambda (4\pi ^{2}k^{2}+q^{2})^{1/2}}}}$ (2)

Since the pump and probe beams used here have Gaussian distribution, the ${\displaystyle 1/e^{2}}$ radius of the pump and probe beam are ${\displaystyle w_{0}}$ and ${\displaystyle w_{1}}$ respectively. The surface is heated by the pump laser beam with the intensity ${\displaystyle p(r)}$, i.e.

${\displaystyle p(r)={\frac {2A}{\pi w_{0}^{2}}}exp(-2r^{2}/w_{0}^{2})}$ (3)

where ${\displaystyle A}$ is the amplitude of the heat absorbed by the sample at frequency ${\displaystyle \omega }$. Then the Hankel transform of ${\displaystyle p(r)}$ is

${\displaystyle P(k)=Aexp(-\pi ^{2}k^{2}\omega _{0}^{2}/2)}$. (4)

Then the distributions of temperature oscillations at the surface ${\displaystyle \theta (r)}$ is the inverse Hankel transforms of the product ${\displaystyle G(k)}$ and ${\displaystyle P(k)}$, i.e.

${\displaystyle \theta (r)=2\pi \int _{0}^{\infty }P(k)G(k)J_{0}(2\pi kr)kdk}$ (5)

The surface temperatures are measured due to the change in the reflectivity ${\displaystyle R}$ with the temperature ${\displaystyle T}$, i.e. ${\displaystyle dR/dT}$, while this change is measured by the changes in the reflected intensity of a probe laser beam. The probe laser beam measures a weighted average of the temperature ${\displaystyle \theta (r)}$,i.e.

${\displaystyle \Delta T={\frac {4}{w_{1}^{2}}}\int _{0}^{\infty }\theta (r)exp\left(-{\frac {2r^{2}}{w_{1}^{2}}}\right)}$ (6a)

This last integral (6a) can be simplified to an integral over ${\displaystyle k}$:

${\displaystyle \Delta T=2\pi A\int _{0}^{\infty }G(k)exp(-\pi ^{2}k^{2}(\omega _{0}^{2}+\omega _{1}^{2})/2)kdk}$ (6b)

The surface temperature of a layered structure

In the similar way, frequency domain solution for the surface temperature of a layered structure can be acquired. Instead of Eq. (2), Eq. (7) will be used for a layered structure.

${\displaystyle G(k)=({\frac {B_{1}^{+}+B_{1}^{-}}{B_{1}^{-}-B_{1}^{+}}}){\frac {1}{\gamma _{1}}}}$ (7)

${\displaystyle \left({\begin{array}{c}B^{+}\\B^{-}\end{array}}\right)_{n}={\frac {1}{2\gamma _{n}}}\left({\begin{array}{cc}exp(-u_{n}L_{n})&0\\0&exp(u_{n}L_{n})\end{array}}\right)\left({\begin{array}{cc}\gamma _{n}+\gamma _{n+1}&\gamma _{n}-\gamma _{n+1}\\\gamma _{n}-\gamma _{n+1}&\gamma _{n}+\gamma _{n+1}\end{array}}\right)\left({\begin{array}{c}B^{+}\\B^{-}\end{array}}\right)_{n+1}}$

${\displaystyle u_{n}=(4\pi ^{2}k^{2}+q_{n}^{2})^{1/2},q_{n}^{2}={\frac {iw}{D_{n}}},\gamma _{n}=\Lambda _{n}u_{n}}$

(Λn: thermal conductivity of nth layer, Dn: thermal diffusivity of nth layer, Ln: thickness of nth layer) Using Eqs. (6) and (7), we can calculate the changes of temperature of a layered structure.

Modeling of data acquired in TDTR

The acquired data from TDTR experiments are required to be compared with the model.

${\displaystyle Re[\Delta RM(t)]={\frac {dR}{dT}}\sum _{m=-M}^{M}(\Delta T(m/\tau +f)+\Delta T(m/\tau -f))exp(i2\pi mt/\tau )}$ (8)

${\displaystyle Im[\Delta RM(t)]=-i{\frac {dR}{dT}}\sum _{m=-M}^{M}(\Delta T(m/\tau +f)-\Delta T(m/\tau -f))exp(i2\pi mt/\tau )}$ (9)

${\displaystyle {\frac {V_{f}(t)}{V_{0}}}={\frac {Q}{\sqrt {2}}}{\frac {\Delta R(t)}{R}}}$ (10)

(Q: quality factor of the resonant circuit) This calculated Vf/V0 would be compared with the measured one.

Application

Through this process of TDTR, the thermal properties of many materials can be obtained. Common test setups include having multiple metal blocks connected together in a diffusion multiple, where once subjected to high temperatures various compounds can be created as a result of the diffusion of two adjacent metal blocks. An example would be a Ni-Cr-Pd-Pt-Rh-Ru diffusion multiple which would have diffusion zones of Ni-Cr, Ni-Pd, Ni-Pt and so on. In this way, many different materials can be tested at the same time.[5] Lowest thermal conductivity for a thin film of solid, fully dense material (i.e. not porous) was also recently reported with measurements using this method.[6]

Once this test sample is obtained, TDTR measurements can take place, with laser pulses of very short duration for both the pump and the probe lasers (<1 ps). The thermoreflected signal is then measured by a photodiode which is connected to a RF lock-in amplifier. The signals that come out of the amplifier consist of an in phase and out of phase component, and the ratio of these allow thermal conductivity data to be measured for a specific delay time.

The data received from this process can then be compared to a thermal model, and the thermal conductivity and thermal conductance can then be derived. It is found that these two parameters can be derived independently based on the delay times, with short delay times (0.1 - .5 ns) resulting in the thermal conductivity and longer delay times (> 2ns) resulting in the thermal conductance.

There is much room for error involved due to phase errors in the RF amplifier in addition to noise from the lasers. Typically, however, accuracy can be found to be within 8%.

References

1. G. Andrew Antonelli, Bernard Perrin, Brian C. Daly, and David G. Cahill, "Characterization of mechanical and thermal properties using ultrafast optical metrology", MRS Bulletin, August 2006.
2. Scott Huxtable, David G. Cahill, Vincent Fauconnier, Jeffrey O. White, and Ji-Cheng Zhao, "Thermal conductivity imaging at micrometre-scale resolution for combinatorial studies of materials", Nature Materials 3 298-301 (2004), doi:10.1038/nmat1114.
3. David G. Cahill, Wayne K. Ford, Kenneth E. Goodson, Gerald D. Mahan, Arun Majudar, Humphrey J. Maris, Roberto Merlin, and Simon R. Phillpot. "Nanoscale thermal transport", J. Appl. Phys. 93, 793 (2003), doi:10.1063/1.1524305.
4. Cahill, DG "Analysis of heat flow in layered structures for time-domain thermoreflectance" Rev Sci Instrum 2007;75:5119, doi:10.1063/1.1819431
5. X. Zheng, D. G. Cahill, P. Krasnochtchekov, R. S. Averback, and J.-C. Zhao, "High-throughput thermal conductivity measurements of nickel solid solutions and the applicability of the Wiedemann–Franz law", Acta Materialia 55, 5177-5185 (2007)
6. Catalin Chiritescu, David G. Cahill, Ngoc Nguyen, David Johnson, Arun Bodapati, Pawel Keblinski, and Paul Zschack, "Ultralow Thermal Conductivity in Disordered, Layered WSe2 Crystals" Science 315, 351-353 (2007) doi:10.1126/science.1136494