Thermophysical properties of 1D materials: transient characterization down to atomic level

Amin Karamati ,Shen Xu ,Huan Lin ,Mahya Rahbar ,and Xinwei Wang

1Department of Mechanical Engineering, Iowa State University, Ames 50011, USA;

2School of Mechanical and Automotive Engineering, Shanghai University of Engineering Science, Shanghai 201620, China;

3School of Environmental and Municipal Engineering, Qingdao University of Technology, Qingdao 266033, China

Abstract:The thermophysical properties of 1D micro/nanoscale materials could differ significantly from those of their bulk counterparts due to intensive energy carrier scattering by structures.This work provides an in-depth review of cuttingedge techniques employed for transient characterization of thermophysical properties at the micro/nanoscale scale.In terms of transient excitation,step Joule heating,step laser heating,pulsed laser heating,and frequency domain amplitudemodulated laser heating are covered.For thermal probing,electrical and Raman scattering-based physical principles are used.These techniques enable the measurement of thermal conductivity,thermal diffusivity,and specific heat from the sub-mm level down to the atomic level (single-atom thickness).This review emphasizes the advantages of these techniques over steady state techniques and their physics,challenges,and potential applications,highlighting their significance in unraveling the intricate thermal transport phenomena to the atomic level of 1D materials.

Keywords: 1D materials;thermophysical properties;transient characterization;Raman scattering;micro/nanoscale

1 Introduction

The past decades have witnessed significant growing interest in one-dimensional (1D) micro/nanomaterials,which possess unique structural and physical properties. These materials have high aspect ratios,with two dimensions significantly smaller than the third one.Examples of 1D micro/nanomaterials include carbon micro/nanofibers[1],carbon nanotubes(CNTs)[2,3],silver nanowires[4,5],and other metallic/semiconductive nanowires[6,7].With their specific function design,1D micro/nanomaterials have become highly attractive for a broad range of applications,including electronics[8,9],photonics[10,11],energy storage[12,13],and sensing[14].Specifically,in electronics,1D materials offer the potential to revolutionize the development of high-performance transistors,lightemitting diodes,and solar cells[15,16].In photonics,1D materials have been leveraged to create nanolasers,optical fibers,and photonic crystals[17,18].In energy storage,1D materials are being used to develop high-capacity batteries[19]and supercapacitors[20].Last not least,1D materials are also suitable for sensing applications,such as detecting gases[21,22],chemicals[23],and biomolecules[24].It is noteworthy that for these applications,the thermophysical properties of 1D micro/nanostructures,particularly thermal conductivity and thermal diffusivity,are critical to thermal design,lifetime,and stability in hostile thermal environments.However,given the novel nature of 1D micro/nanoscale materials,their thermal conductivity/diffusivity differs significantly from those of their bulk counterparts.

The unique physical properties of novel 1D micro/nanoscale materials stem from their thickness/radial dimension,which can be at or below certain characteristic lengths,such as the phonon mean free path,photon absorption depth,or electron mean free path,leading to a very strong size confinement effect[25].The intensive energy carrier scattering caused by nanostructures and small size in their reduced dimensions significantly reduces the energy carrier’s mean free path and thermal conductivity/diffusivity.It is crucial to understand the intrinsic thermal properties of these materials via appropriate thermophysical property measurements. However,the extremely small size of these materials presents great challenges in applying well-defined thermal excitation and characterizing the temperature differential across them[26,27].

Various methods have been developed to measure the thermophysical properties of 1D micro/nanomaterials,which can be broadly categorized into two types: steady state techniques (such as microbridge and conventional Raman thermometry) and transient techniques (including the 3ω technique,transient hot wire,and flash techniques).Steady state techniques are generally based on simple physical models but require a very careful assessment of boundary conditions and the establishment of a stable and measurable temperature gradient and heat flux.One commonly used steady state technique is the microbridge method[28,29],which involves the fabrication of a microdevice with two suspended membranes,between which a 1D sample is placed to form a bridge structure.One membrane serves as a heating element,while the other acts as a temperature sensor.A DC current is passed through the heating membrane to provide stable heat,which then flows through the sample/bridge to the sensing membrane,causing a small temperature rise in the latter.By measuring the temperature of the two membranes,the heat flow through the sample,and the sample’s dimensions,the thermal conductivity can be directly determined according to Fourier’s law.This technique offers high accuracy and precision and can do the measurement over a wide temperature range,but it also has limitations that must be carefully considered when interpreting results.Factors such as the thermal resistance of the membranes,the thermal contact resistance between the sample and the membranes,and the thermal radiation between the membranes can affect the accuracy of the technique,especially for samples with a low thermal conductivity[30,31].Advancements in eliminating the effect of thermal contact resistance between samples and microdevices can improve the accuracy of the measurement[32,33].In general,steady state techniques require precise determination of the heat transfer and temperature drop across micro/nanoscale materials,which is challenging to implement and subject to large experimental uncertainties and laborious experimental control and conduction.

In contrast,transient techniques offer the advantages of fast implementation and far less strict experimental control.Instead of analyzing the temperature response at steady state,transient techniques analyze the response in the time and frequency domains.The main thermal property extracted from transient techniques is thermal diffusivity.Due to the higher signal-to-noise ratio in the transient temperature response,the accuracy of the measured results is guaranteed,even with far less strict experimental control.One commonly used transient technique is the 3ω technique[34],which belongs to the electrothermal measurement category.It feeds a sine wave current with a frequency of ω into the sample,and the generated Joule heat and resistance change in the sample respond in 2ω,while the monitored voltage is in 3ω.By using lock-in technology to pick up the 3ω voltage,the thermal conductivity and sometimes specific heat can be accurately determined based on the voltage-frequency relation[35].Although the thermophysical properties of coatings/films in the thickness direction are also related to energy transport at the micro/nanoscale,they are not the focus of this work.The 3ω,timedomain thermal reflectance (TDTR),laser flash,and photothermal techniques work well to measure the thermal conductivity/diffusivity in the thickness direction of coatings/films[36–45].The transient hot wire method is commonly used to measure the thermal conductivity of liquids and powders[46,47].Efforts have also been made to use it to measure coatings[48].

Instead of a general review,we intend to provide very indepth reviews of the progress in method development for 1D micro/nanomaterial characterization in Wang’s lab at Iowa State University.This review covers transient electrothermal(TET) and transient photoelectro-thermal (TPET) techniques in Section 2,including their experimental principles,physical models,experimental setups,and typical applications.While these two methods are very effective in measuring most 1D materials at the millimeter and micrometer scales,their limitations become apparent when samples are reduced to nanometer-scale thickness of μm length.Section 3 introduces the pulsed laser-assisted thermal relaxation (PLTR) technique,which overcomes the limitations of the previous two methods.Additionally,the mathematical and physical relationships among the three techniques are explained.Section 4 discusses the challenges in measuring 1D nanomaterials at the atomic level and introduces a differential concept for addressing these issues.Finally,Section 5 presents photon excitation and probing methods based on advanced Raman techniques that can measure materials of～nm thickness and a few μm length.

2 Thermal characterization under step thermal excitation

In this section,the TET and TPET techniques are discussed in detail.Both involve step heating and electrothermal sensing.The TET technique was invented by our lab in 2007[49].It has been shown to be very effective in measuring the in-plane thermal diffusivity of fiber-or film-like materials with high accuracy[50–53].A schematic of the TET technique is shown in Fig.1a. The to-be-measured sample is suspended over a trench and is secured using silver paste,which gives a good electrical and thermal connection between the sample and the electrodes.In the measurement,a transient voltage rise is induced on the sample by Joule heating using a step DC current.The transient temperature variation or,equivalently,the voltage rise (or drop) over the sample (Fig.1c) is then used to determine the thermal diffusivity.Due to the sample’s high length-to-thickness ratio,it is reasonable to assume that the heat conduction is one-dimensional along the sample.The 1D transient heat conduction along the sample (i.e.,x) is governed by the following equation[49]:

Fig.1. (a,b) Schematic of the TET technique and TPET technique.(c) Graphene fiber TET signal with time as an example.Reproduced with permission from Ref.[54].Copyright 2022.Elsevier.The inset is the suspended sample of 1.959 mm length.(d) Representation of the TET signal fitting: taking the natural of voltage subtracted by the steady state voltage.Reproduced with permission from Ref.[54].Copyright 2022.Elsevier.The difference among the fittings using three different initial data treatments.

whereρ,cp,andkare the sample’s density,specific heat,and thermal conductivity,respectively.Tis temperature,andin W?m-3refers to the heating induced by the electrical current.Moreover,the left end of the sample is taken as the starting point of thex-axis.The temperature distribution of the sample is then solved as[49]:

where α is the thermal diffusivity and the sample length isL.GX11is the Green function for pulsed heating in the sample at timeτand locationx′[49],

Eq.(2),as the transient temperature response over the sample,can ultimately be simplified to a normalized average temperature rise as[49]:

whereT?is the normalized temperature rise with respect to the steady state temperature rise caused by the same Joule heating for the sample.It should be noted that the temperature rise is averaged over the sample’s length.From the experimental data,the normalized temperature rise can be calculated as[54]:

whereV1is the voltage at the final steady state andV0is the voltage before the step heating begins.Defining the Fourier numberFo=αt/L2,Eq.(4) demonstrates that the normalized temperature rise versusFois always in the same pattern regardless of the dimension and the material of the sample[49].Several methods can be used to determine the α of the sample once the normalized temperature rise is determined by the experiment[49].Linear fitting at the initial stage of the experiment is one of the fitting methods.Knowing the length of the sample,it was shown that fitting the initialFo<0.001 experimental data will give a good α value[49].The second method is done by identifying a characteristic time (?tc),which is the time whenreaches 0.8665.After determining ?tcusing the～tcurve,α can be calculated as α=0.2026L2/?tc[49].The third method uses nonlinear global curve fitting,which uses different trial values of α to find the best fit based on the least square method.In addition to these methods,a linear fitting method was recently developed by our group[54].It was shown that after taking the natural log of the voltage variation over time for the TET experimental data,it perfectly resembles a line.Based on a detailed numerical calculation using Eq.(4),it is concluded that[54]:

Based on Eq. (6),the slope of the linear fitting for ln|V-V1|～tis -π2α/L2,from which α can be easily obtained.This rigorous method not only gives the opportunity to take advantage of the simplicity of linear fitting but also results in uncertainty determination of the fitting,which is difficult to determine in nonlinear global data fitting.Since a similar relation can also be established for the TPET and PLTR techniques,their measurement fitting uncertainty can be determined rigorously in the same way. The final measurement uncertainty will also need to include the uncertainty in length measurement,which is mostly small.It should be noted that the data falling in therange of 0.1-0.8 are used for the fitting via this method,as it brings about the best accuracy[54].Fig.1d presents the linear fitting of the signal given in Fig.1c.The α of the sample (a micro-thick graphene fiber)shown in Fig.1c is determined to be (9.61×10-7±1.97×10-9)m2·s-1using the data range of 0.1<<0.8.When<0.1,some nonlinearity (even very small) arises in the data,making the fitting slightly less reliable,but the difference is still very small.As shown in Fig.1d,if we use the data of 0.05<<0.8 and<0.8 to do the fitting,theαvalue is determined to be (9.62×10-7±1.86×10-9)m2·s-1and (9.62×10-7±1.75×10-9)m2·s-1,very close to the fitting result using the data range of 0.1 <<0.8[54].

In the TET technique,heating and thermal sensing use the same step current.Sometimes this makes the experiment less controllable,especially for highly electrically conductive materials where a very high current is needed to achieve a sensible temperature/resistance/voltage change. Additionally,a high voltage has to be used for very poor electrically conductive samples,making the voltage too high to sense using an oscilloscope.The TPET technique is designed to overcome this problem,to provide a much faster rise time of heating,and to enable the measurement of samples of very highαand short length[55].In TPET,a step amplitude-modulated continuous wave (CW) laser is used to irradiate the whole sample(Fig.1b).Upon laser irradiation,the temperature of the suspended sample increases,and this transient temperature rise leads to a transient change in the sample’s electrical resistance as well as the voltage (Fig.1c).To detect this transient change in temperature,a small DC current is passed through the sample.It has been proven that the optimum amount of this DC current should be equal to 0. 79ITET[54]to cause the same level of transient temperature rise as TET but with the least accumulated heating.For uniform thermal excitation,the temperature (voltage) evolution is the same as that in the TET method.This voltage evolution is then used to obtain α using the processing methods mentioned earlier for the TET.

Instead of irradiating the whole sample in normal TPET,the laser beam can be resized to a narrow line-shaped beam oflwidth for localized heating and sensing.Integrating Eq.(2)along the length of the sample results in Eq.(4) for both the TET and the TPET with uniform heating.Now,if the laser beam with a width oflirradiates the sample at locationx,the normalized temperature rise would be[54]:

where

Using this localized heating,the α distribution along the sample’s length can be obtained.It was shown that if this localized laser beam irradiates the sample in a range ofx～(0.2L-0.8L),Eq.(6) can still be used to obtain α with negligible dependence on laser heating location if the sample has no spatial variation in its thermophysical properties[54].In fact,it has been proven that for localized laser heating,in a large portion of the data in the time domain,Eq.(6) is still applicable.One merit of this TPET scanning is that if the measured α has a certain distribution against the laser heating location,then such a distribution can be used to obtain the likely patterns for some of the properties (e.g.,temperature coefficient of resistivity) the sample may have.

The TET and TPET techniques are both very potent for measuring the thermal properties of a large variety of micro/nanoscale wires and films,whether they are conductive or nonconductive.However,they have advantages over each other from different viewpoints.The TPET,for instance,can be applied to extremely large or low-resistance samples.Additionally,the thermal excitation during TPET measurement has nothing to do with the electrical contact resistance.Another advantage of the TPET is that the step laser excitation could be much faster than the step electrical current used during TET measurement.The TET,on the other hand,has controlled and well-defined and known heating for a great variety of samples,which can simply be calculated asI2R.Ris the resistance of the sample,andIis the magnitude of the step current.Moreover,the effective thermal conductivity(keff) of the sample can be obtained using the sample’s steady state temperature rise askeff=I2RL/(12Ac?T)[56],whereAcis the sample’s cross-sectional area and ?Tis the average steady state temperature rise over the sample.This ?Tcan be evaluated as ?R/(dR/dT),where ?Rcan be calculated from the TET voltage change.Fitting calibrated (R-T) data points will obtain dR/dT.Oncekeffis calculated and αeffis obtained by fitting the TET signal,(ρcp)effof the sample can be determined.This is quite a unique capability since there are very few techniques that can measure the specific heat of extremely lightweight nanoscale materials.For instance,for a sample of 1 mm length,2 μm diameter,and 2000 kg?m-3density,its mass is only 6.2×10-9g.The specific heat of such a sample can be readily measured using the TET technique.Using the TET method,Zhu et al.[56]measured the thermal properties of a type of 3C-SiC nanocrystalline microfiber called LCVDFFF of 1495 μm length and 34.29 μm diameter coated with a 20 nm Ir layer.The αeffandkeffof the sample were measured to be 2.89 mm2·s-1and 5.59 W·m-1·K-1at room temperature,from which the (ρcp)effof the sample was determined as 1.94×106J·m-3·K-1.In another study by Deng et al.[57],three samples of carbon nanocoils (CNCs) were characterized via the TET technique.The range of the line diameter and total length of the CNCs are on the order of 300–500 nm and 100 μm,respectively. The corresponding αeffvalues of these samples were 1.58 mm2·s-1,2.7 mm2·s-1,and 4.11 mm2·s-1,and thekeffvalues were 1.67 W·m-1·K-1,2.77 W·m-1·K-1,and 3.20 W·m-1·K-1at 290 K.Having the αeffandkeffvalues of the samples and assuming ρ to be 2200 kg·m-3,the averagecpof the CNC was reported to be 910 J·kg-1·K-1at 290 K. DNA fiber (ρcp)effat 190 K was reported to be 7.3×105J·m-3·K-1calculated using αeff=0.605 mm2·s-1(from TET) andkeff=0.44 W·m-1·K-1[58]. The length and width of the fiber were 0.405 mm and 3.47 μm,respectively.The aforementioned studies reveal how powerful the TET technique is in thermophysical property characterization of low-dimensional fiber-or film-like materials.These fibers and films are of a wide variety.This variety includes polyacrylonitrile (PAN) wires (Fig.2a),carbon nanocoils (Fig.2b),3C crystalline SiC microwires (Fig.2c),carbon nanotubes(Fig.2d),and chemical vapor decomposition (CVD)-grown graphene supported on poly(methyl methacrylate)(PMMA)[59].Table 1 presents more details about these micro/nanoscale samples.

Table 1. Details of the micro/nanoscale sample characterized via TET.

Fig.2. (a) SEM images of single polyacrylonitrile (PAN) wires.Reproduced with permission from Ref.[60].Copyright 2007.Springer Nature.(b) SEM images of carbon nanocoils (CNCs).Reprinted with permission from Ref.[57].Copyright 2016.American Chemical Society.(c) SEM images of 3C crystalline silicon carbide (SiC) microwires.Reproduced with permission from Ref.[56].Copyright 2018. Elsevier.(d) SEM images of carbon nanotube(CNT) bundles.Reproduced with permission from Ref.[61].Copyright 2018.Elsevier.

In this section,the merits and drawbacks of the TET and TPET techniques over each other are discussed.In the next section,the PLTR technique as a potent technique to overcome some common limitations of the TET and TPET techniques will be discussed in detail.The αeffvalues summarized in Table 1 have an uncertainty of ±10% based on the global TET fitting.

3 Nanosecond pulsed thermal excitation and relaxation

While the physical principle underlying the TET/TPET techniques is applicable to a wide range of 1D micro/nanoscale materials,practical implementation requires electronic devices to supply a current and/or measure the sensing voltage.Unfortunately,the slow rise time (on the order of microseconds) of the current source and amplitude-modulated laser sources can severely limit the measurement capabilities of the TET/TPET techniques,as they rely on detecting the Joule/photon heating-induced voltage change.In particular,when the sample has a short thermal characteristic time (L2/α)falling within the same time span as the rise time,it can be difficult to differentiate the heating effect from the device’s rise time.Moreover,continuous heating of the sample in the TET/TPET techniques can cause extensive heat generation and conduction to the electrodes,leading to a continuous temperature rise and potential damage to the sample if the electrode support has poor heat dissipation capability or the sample has a high thermal conductivity that needs high electrical heating to induce a sensible resistance change.Currently,current and amplitude-modulated laser sources with faster rise times are becoming more available and can be used to overcome the above challenges.Here,the PLTR technique[62–66]overcomes these limitations by using an ultrafast heat source at the nanosecond or even picosecond timescale and focusing on detecting the relatively longer cooling process.The ultrashort duration of the heating period introduces minimal heat into the sample,significantly reducing the risk of sample damage.Additionally,this technique has no continuous heating effect and is applicable to a wide range of materials.

Fig.3a and b shows the experimental principle of the PLTR technique.Similar to those in the TPET technique,a 1D sample is suspended between two electrodes,and a small constant DC is passed through the sample to monitor changes in voltage/resistance that indicate the temperature rise.As mentioned above,in contrast to heating and voltage measurement simultaneously over a relatively long period of time,PLTR employs a nanosecond pulsed laser to introduce very short-time heating into the sample and records the subsequent cooling process.The small DC current provides a steady and constant temperature rise in the sample without affecting the temperature rise evolution θPLTRcaused by pulsed laser heatingMoreover,the cooling process in the length direction of the 1D sample is typically much slower than the nanosecond-scale heating duration.Consequently,the governing equation for 1D heat conduction is given by:

Fig.3. (a–c) The PLTR technique in the physical principle,schematics of the temperature response at the back side of the sample,and typical experimental data and fitting curve for in-plane heat conduction,respectively.Reproduced with permission from Ref.[65].Copyright 2008.World Scientific Publishing.(d) Typical experimental data and fitting for cross-plane heat conduction.Reproduced with permission from Ref.[63].Copyright 2017.Elsevier.

By normalizing θPLTRwith respect to the maximum temperature rise θPLTR,maximmediately after pulse heating,the resulting normalized temperature evolutionis given as:

Fig.3c shows a typical normalized temperature rise in the PLTR measurement.The thermal diffusivity can be determined by global data fitting based on Eq.(12) or characteristic point evaluation as α=0.2026L2/tc,wheretcis the time taken forto decrease to 0.1097[65,66].

Technically,heat conduction after laser pulse heating occurs simultaneously in both the thickness and length directions of the sample,as shown in Fig.3a.However,due to the high aspect ratio of the 1D sample,a thermal steady state is very quickly established in the thickness direction,which typically requires several orders of magnitude less time than cooling in the axial direction.In other words,the temperature at the back side of the sample will increase first according to the cross-plane heat conduction from the heated front surface and then decreases due to in-plane cooling to heat sinks(Fig.3b).As a result,it is possible to distinguish cooling in the thickness and length directions,as depicted in Fig.3a.Similar to the flash technique,the cross-plane thermal diffusivity can be determined as:

whereDis the thickness of the sample andt1/2is called the“half-rise time”,which is the time for the back side of the sample to reach half of the maximum temperature rise.Fig.3d shows the cross-plane thermal diffusivity measurement and the resulting thermal conductivity data for a 21.4 μm-thick graphene paper from 30 K to 295 K.For multilayered structures,a more complicated 1D physical model is needed to determine the thermal diffusivity.Usually,numerical heat conduction simulation is a good choice for such data processing[63].

The TET/TPET and PLTR techniques employ different functioned heat sources to introduce temperature changes in 1D materials,which are measured based on a constant resistance temperature coefficient (RTC) of metallic sensors.The thermal diffusivity of the sample is then determined using the corresponding physical model.Despite their differences,these methods are closely related from a deeper physical and mathematical standpoint,as shown in Fig.4.Specifically,they share the same governing equation (Eq.(9)) except for the heat source: a step function for TET/TPET and a delta function (i.e.,a pulse) for PLTR.

Fig.4.Mathematical relation behind the physical models of different types of heat sources.

For 1D suspended micro/nanoscale materials under heating and thermal sensing,if the step heat source is designated asH(t)=0 fort<τ andH(t)=fort≥τ,the temperature rise solution of θH(x,t) takes the form of the TET and TPET techniques (Eq.(4)).For the PLTR technique,the heat source(Θ) is a delta function (δ),which can be expressed as Θδ=?H/?t.Since the heat conduction equation is a linear equation,the resulting temperature rise solution θδ(x,t) can be expressed as θδ(x,t)=?θH(x,t)/?t.conversion makes it straightforward to construct the solution of heat conduction under an arbitrary heat source based on a simple temperature rise solution under step heating.Additionally,since the temperature rise under step heating can be tabulated for discrete moments (Δt,2Δt,3Δt,...) in the time domain,the temperature rise under an arbitrary heat source can be obtained by simple mathematical summation,leading to superior computational efficiency.

4 Differential concept for down to atomiclevel characterization

In the TET or TPET measurements,the samples need to be suspended between two electrodes.Due to their relatively weak structure,nanometer-thick materials cannot be directly suspended.To overcome this problem,the differential TET technique[58,73,74] has been developed to measure the in-plane thermal conductivity of nanostructures.In the experiment,a low-dimensional material with a low thermal conductivity is used as the substrate to support the to-be-measured ultrathin film,as shown in Fig.5a.The effective thermal diffusivity(αeff) of this composite structure can be characterized by using the TET technique.First,the substrate is deposited with one layer of film with a thicknessδc,and the effective thermal diffusivity of the thin structure-substrate system is measured asαeff,1.Then the same procedure is repeated (depositing another layer ofδcthickness and measuring αeff)ntimes.Fig.5b shows the results for theαeffof a glass fiber-Ir coating structure.From Fig.5b,it is clear that αeffchanges linearly againstR-1.Finally,it is possible to establish the relationship betweenαeff,nand the number of layers.The thermal conductivity of the thin structure can be determined by establishing heat conduction models for the thin structure-substrate system.If the thin structure cannot be easily coated multiple times,we will be only able to measureαeff,1and then subtract the influence of the substrate.More details can be found in Refs.[59,75].

Fig.5.(a) The cross-sectional schematic of a substrate with different nanofilm layers deposited.(b) Variation of the effective thermal diffusivity of an Ircoated glass fiber against the inverse of the electrical resistance for 6.4 nm Ir layers coated on the glass fiber.(c) Linear fitting of αeff ～L2 for 1.33-layered graphene on PMMA.(d) Graphene temperature evolution with time and the temperature difference between PMMA and graphene in the middle point of the length direction.(b) Reprinted with permission from Ref.[75].Copyright 2013.John Wiley and Sons,Inc.(c,d) Reproduced with permission from Ref.[59].Copyright 2016.The Royal Society of Chemistry.

When a step current is applied to the thin structuresubstrate system,heat is transferred from the thin structure to the substrate,but the temperature difference between them is negligible compared with the temperature rise of the sample and has a negligible effect on the TET measurement.Fig.5d shows the numerical simulation results of the graphene temperature variation with time and the temperature difference between PMMA and graphene in the middle of the axial direction.This is for a 1.54 mm long PMMA film of 790 nm thickness on which monolayer graphene is deposited.It can be seen that the overall temperature of the graphene increased by～ 5 °C.The temperature difference between PMMA and graphene is approximately 10-5°C.Therefore,the temperature of graphene and PMMA could be treated the same during Joule heating. In the case of a PMMA substrate(k=0.2 W·m-1·K-1) with a width of 1 mm,length of 2 mm,and thickness of 0.5 μm,we can estimate that its thermal resistance along the length direction is 2×107K·W-1,and the total thermal resistance between the graphene layer and substrate is 0.5 K·W-1for a thermal contact resistance of=10-6K·m2·W-1.Therefore,the thermal resistance along the length direction is much larger than that along the thickness direction.The temperature rise of the sample is mainly affected by the length direction heat conduction.In the thin structure-substrate system,when heat is transferred from the film to the substrate,the characteristic time to reach thermal equilibrium is ?t1=whereδcis the film’s thickness.The characteristic thermal diffusion time in the length direction is Δt2≈L2/α.In the case of a δc=10nm thick film on a PMMA substrate with a length of 1 mm and a thickness of 0.5 μm,=10-6K·m2·W-1,of 2×106J·m-3·K-1,and Δt1is 20 ns. The order ofαis approximately 10-6–10-3m2?s-1,so Δt2is estimated to be 1 ms to 1 s.It is obvious that Δt2is much longer than Δt1.Consequently,in the TET measurement,the time for the cross-sectional heat conduction to reach thermal equilibrium can be negligible.

αeffis an effective value that combines the effects of the substrate,the coating,and thermal radiation.Thus,it can be written as[49,58,75,76]:

where ε,σ,andδare the emissivity,Stefan-Boltzmann constant,and thickness,respectively.Subscripts “s” and “c” represent substrate and coating,respectively.LandT0are the sample length and initial temperature,respectively.M=16 for circular cross-section substrates andM=8 for substrates with rectangular cross-sections. Here,(ρcp)eff=[δc(ρc)c+δs(ρc)s]/(δc+δs).

To eliminate the influence of radiation,TET experiments need to be conducted using the same sample of different lengths.According to Eq.(14),we can see thatαeffand length square (L2) have a linear relationship.The intercept at theyaxis of theαeff～L2relationship is the real thermal diffusivity of the sample.To demonstrate this physics,theαeff～L2linear fitting for the graphene sample supported by PMMA (1.33 layers on average) of different lengths is shown in Fig.5c.With the knowledge of the fitting slope,PMMA substrate thickness,and temperature,the emissivity of the graphene samples can also be determined by dividing the slope by

For the glass fiber measured in Ref. [75],where the thermal conductivity of ultrathin Ir coating is determined,the impact of thermal radiation is negligible due to its short length.If the length of the glass fiber is long,the impact of thermal radiation cannot be neglected.It should be noted that if the substrate is very thin or has a low thermal conductivity,the thermal radiation effect will become much weaker when the Ir coating thickness increases.The emissivity of glass fiber is approximately 1 and much less than 0.1 for Ir.When layers of coating are added,the emissivity of the thin structure-substrate system will be reduced,resulting in a decrease in the thermal radiation effect.Therefore,as the film thickness increases,there may be a decrease inαeff,which could overshadow the effect of the coating.This problem can be solved by measuring samples of the same coating thickness but different lengths so as to subtract the influence of thermal radiation.

It is clear in Eq.(14) that if the substrate has a low thermal conductivity and thin thickness,the increase inαeffwill be more preeminent when the film thickness is increased.Therefore,in the differential TET technique,substrates of low thermal conductivity and thin thickness are preferred to maximize the measurement sensitivity.

5 Photon excitation and probing: Raman scattering

As mentioned earlier,in the TET,TPET,and PLTR techniques,the transient temperature rise of the sample leads to a transient change in the sample’s electrical resistance.Therefore,in all the abovementioned techniques,the voltage variation across the sample is used to measure the thermal properties of the sample[41,49,55,65,66,77–79].A linear relation is required between the changes in temperature and resistance.Since the needed temperature rise in measurement is small (a few to 10 degrees),such a linear relation usually holds true even for materials with a nonlinear temperature-resistance relation over a large temperature range (e.g.,300 K)[80].Such electrical sensing sometimes requires very difficult sample preparation,especially for nm-thick wires/tubes/films and very low electrical resistance samples.Additionally,for samples whose resistance has a very weak or no response to temperature change,these techniques become difficult/impossible to use.The Raman-based techniques reviewed in this section are designed to resolve the above problems (even not completely).

To date,Raman scattering has been used extensively for thermal probing and measuring the thermal properties of 2D materials such as graphene,in which measuring the exact temperature rise and the sample’s laser absorption during heating are required[81–84].Various studies have used very different optical absorption values,which is a major source of uncertainty in Raman-based techniques.In addition,for 1D nanoscale materials such as nanotubes,such measurements become almost impossible since precise determination of laser beam absorption in the material cannot be conducted with sound confidence.For example,CNT chirality,defect level,and whether they are single or in bundles strongly affect how much laser power is absorbed in the nanotube.Hsu et al.[85]measured the optical absorption of single-walled carbon nanotubes (SWCNTs).A large number of factors have to be considered and treated for evaluating the heat flux caused by laser beam absorption,and quite large uncertainties arise in the estimation.

To eliminate the need for an exact temperature rise and laser absorption evaluation in Raman-based experiments for 2D nanoscale materials,different Raman techniques have been developed in our lab,such as ET-Raman,frequencyresolved Raman (FR-Raman),and FET-Raman[52,70,71,86–89].

The ET-Raman technique uses two distinct energy transport states in both space and time domains to study the thermal response of the sample: the near zero-transport case using picosecond or nanosecond (ps or ns) lasers and the steady state case using a CW laser.The FR-Raman uses an amplitude-modulated laser with a square-wave form in a wide range of frequencies for both heating and Raman probing[70,89].The FET-Raman was developed based on the ETRaman and FR-Raman to measure the thermal properties of materials.It uses a CW laser and amplitude-modulated squarewave laser (FR) with a specific frequency,as shown in Fig.6a and b[52,71].Since it uses exactly the same configuration for both CW and FR situations,the laser absorption is exactly the same for both cases.In addition,the FR case makes the signal very stable and strong to collect.

Fig.6.Physical principle of the FET-Raman technique.(a) Temperature rise of the suspended sample under CW laser heating.(b) Amplitude-modulated laser heating.

Herein,we present the FET-Raman technique used to measure the thermal diffusivity of suspended well-aligned 1D SWCNTs over a silicon substrate trench using a 50× objective lens at 257 K.The sample consists of SWCNT bundles with a length of 11 μm and diameters in the range of 1.61 nm to 6.70 nm,which are composed of a few individual SWCNTs.The diameter of the laser spot under a 50× objective lens is 1.4 μm,covering several SWCNT bundles.In both CW and FR laser states,the laser irradiates the sample,causes the sample’s temperature to rise,and measures its Raman shift simultaneously under different laser powers (P).

Afterward,ψ=?ω/?P,termed the Raman shift power coefficient (RSC),is determined for each heating state.Here,ωis the Raman wavenumber.ψCWdepends on several factors,including the laser absorption coefficient (αl),temperature coefficient of Raman shift (?ω/?T),and sample thermal conductivity (k) as ψCW∝αl(?ω/?T)·f1/k.On the other hand,ψFRdepends on αl,?ω/?T,k,and volumetric heat capacity (ρcp) as ψFR∝αl(?ω/?T)·f2(α)/(ρcp).RSC normalization as Θ=ψFR/ψCWeliminates the effect of Raman temperature coefficients and laser absorption.Therefore,the value of Θ is only dependent on the sample’s thermal diffusivity (α),which will be used to determine α.

In the CW laser state,the energy transport governing equation is[90,91]

in whichTCWis the temperature rise under CW laser heating andis the Gaussian beam heating,which is defined below[88]:

In the FR laser state,the energy transport governing equation is[92]

in whichTFRis the temperature rise in the FR laser state,andis the same as in Eq.(16).However,at frequencyf,the square-modulating duty cycle is 50%.

Solving Eqs.(15) and (17) will yield the theoretical ratio of the sample’s temperature rise of the FR case to the CW case.The relationship between this theoretical ratio and the sample’s thermal diffusivity can be established.By comparing theoretical and experimental ratios,the SWCNT sample’s thermal diffusivity can be determined.Note that while simulating the 1D model,we rigorously consider the Raman intensity weighted average temperature rise over spacein the CW case and the Raman intensity weighted average temperature over space and timein the FR case.

Since the sample might be damaged and we might not obtain a linear ω ～Prelationship at high laser powers,we keep the laser power as low as possible for both laser heating states as 1.79–17.9 mW using a 50× objective lens.To achieve the highest sensitivity of the thermal diffusivity measurement,it is important to keep Θ as close to 0.75 as possible[71].Fig.7a shows the 2D contour plot of the Raman shift of the G peak at～1582 cm-1at 257 K under different laser powers.For both CW and FR laser states,the Raman shift of the G peak is linearly related to the laser power at 257 K with ψCW=(-0.161±0.007)cm-1·mW-1and ψFR=(-0.118±0.006)cm-1·mW-1,which gives a Θ of 0.73±0.02.It is important to note that the laser power is measured during the laser-on time,not as an average over the laser-on and off times.Fig.7b shows the modeling result of theoretical values of Θ at 1 MHz versus thermal diffusivityα.Using this curve and the experimental Θ of 0.73,the sample’s thermal diffusivity is determined to be (6.54±0.03)×10-5m2·s-1at 257 K.

Fig.7.(a) 2D contour plot of the Raman shift for the G peak against laser power at 257 K for SWCNT bundles.Left side: for the CW case,the slope of linear fitting,as indicated with a solid black line,is ψCW=(-0.161±0.007)cm-1·mW-1.Right side: for the FR case,the slope of linear fitting,as indicated with a solid black line,is ψFR=(-0.118±0.006)cm-1·mW-1.(b) Theoretical values of the temperature rise ratio at 1 MHz from numerical modeling against the thermal diffusivity.Using the normalized RSC (Θ) from the experiment,the thermal diffusivity of SWCNTs can be interpolated,as indicated by the blue dashed line.

The ET-Raman can also be used to measure 1D micro/nanoscale materials.Because a pulsed laser provides very intense photon heating in a very short time and the pulsepulse interference is low since the temperature rise caused by one heating pulse can be quite relaxed before the next pulse comes in,the physical sensitivity of the ET-Raman is higher than that of the FET-Raman.However,it might damage the sample due to intense pulse heating.Furthermore,involving more intense laser heating during the very short time in the ET-Raman could exceed the absorption saturation level,meaning that the sample will not absorb any more laser,which is not preferred.The advantage of ET-Raman is that we can go to an extremely small time domain (ps) if a picosecond laser is used for heating[86–88].Due to the same laser-on and off times in FET-Raman,the sample may not completely cool down after one laser heating before the next heating interval,resulting in decreased sensitivity.

6 Concluding remarks and perspectives

The transient techniques developed in Wang’s lab offer powerful means of measuring the thermophysical properties of various micro/nanoscale wires and films,with each technique having its own advantages and disadvantages.The TET technique provides controlled and well-defined heating,making it suitable for a wide range of samples.Samples with extremely high or low electrical resistance or a very fast thermal characteristic time will be very difficult to measure using the TET technique.On the other hand,the TPET technique overcomes the limitations of high electrical current/voltage that might be needed in TET and enables measurements of extremely large and low-resistance samples. However,both techniques rely on detecting the step heating-induced voltage change,which can be limited by the slow rise time of the current/laser sources.To address this limitation,the PLTR technique utilizes ultrafast heat sources (i.e.,laser) at the nanosecond or picosecond timescale and probes the cooling process.The PLTR technique offers the capability of measuring very fast heat transfer processes and reduces the risk of continuous heat generation and potential damage to the sample.

TET or TPET measurements of nanometer-thick materials can be problematic because they cannot be directly suspended due to their weak structure.The differential TET technique overcomes this problem by utilizing low-dimensional materials as supportive substrates for ultrathin films,enabling the measurement of in-plane thermal conductivity in nanostructures.The TET,TPET,and PLTR techniques rely on a good electrical resistance response to temperature changes.For materials whose resistance has a weak or no response to temperature change,Raman-based techniques can be used and have been widely applied to measure 1D and 2D nanoscale materials.Wang’s lab has developed Raman-based techniques,such as ET-Raman,FR-Raman,and FET-Raman,eliminating the need for precise temperature rise and laser absorption evaluation.On the other hand,these Raman techniques will have low signals when the temperature is high.For samples with easy photofluorescence excitation,Ramanbased techniques are difficult to use since the fluorescence signal can easily overshadow the Raman signal.For wide bandgap materials with a very strong Raman signal but negligible laser absorption,Raman-based techniques cannot be used directly since the Raman laser will not cause simultaneous heating.For such a scenario,a very thin material (e.g.,few-layered graphene) can be coated on the sample for both excellent Raman excitation and laser absorption/heating purposes.

Thermophysical property measurements at the micro/nanoscale still face significant challenges caused by extremely high or low thermal conductivity,material geometry,and material size.Additionally,studying the variation of thermophysical properties against temperature from very low (a few K) to very high levels (1000 K or higher) requires additional and rigorous technical and scientific considerations and design.Such a temperature effect is of great engineering and science interest in device and material design and optimization.Generally speaking,transient techniques offer great advantages over steady state techniques for several reasons.Steady state techniques usually require precise control and determination of the heat flux and temperature difference across the sample.Such experimental control is often very challenging and subject to large uncertainties.Transient techniques,however,have no such limitation.Most times,it only needs size information about the sample and its relative temperature response,which can be readily obtained.We foresee that transient Raman techniques,including the ET-Raman and FET-Raman,will see broad adoption and provide unprecedented measurement capability and physics knowledge about 1D micro/nanoscale materials.

Acknowledgements

This work was supported by the National Key Research and Development Program (2019YFE0119900),US National Science Foundation (CBET1930866,CMMI2032464),National Natural Science Foundation of China (52106220),and Natural Science Foundation of Shandong Province(ZR2020ME183).

Conflict of interest

The authors declare that they have no conflict of interest.

Biographies

Amin Karamatiis pursuing his Ph.D.degree in Mechanical Engineering at Iowa State University since 2021.He received his M.S.and B.S.in Mechanical Engineering from K.N.Toosi University of Technology,Iran,in 2016 and 2013,respectively.Thermal transport in 1D scale has been his research focus.Currently,he is conducting research on the thermal characterization of micro/nanowires using Raman Spectroscopy,TET,and TPET techniques.

Shen Xuis a Professor of Shanghai University of Engineering Science.She received her B.S.from East China University of Science and Technology in 2008,M.S.from Fudan University in 2011,and her Ph.D.degree from Iowa State University in 2015.Her research interests cover energy transfer at the micro/nanoscale,heat transfer behavior of new materials and complex nanostructures,and multiphysics in the coupling of optical and temperature fields.She has developed the first time-domain differential Raman (TD-Raman) and frequency-resolved Raman (FRRaman) techniques for high accuracy measurement of thermal conductivity and diffusivity of atomic-scale materials.

Huan Linis a Professor in Engineering Thermophysics at Qingdao University of Technology.She obtained her Ph.D.degree from the College of Engineering,Ocean University of China in 2014,M.S.(2008)from Beijing Institute of Technology and B.S.(2006) from Shandong University of Technology.Her current research focuses on thermal characterization of micro and nano materials.She has developed the differential TET technique for measuring thermal conductivity of fibers and films down to～nm scale,and conducted the first characterization of in-plane structure thermal domain size of 2D materials.

Mahya Rahbaris pursuing her Ph.D.degree in Mechanical Engineering at Iowa State University in 2022.She received her M.S.and B.S.in Chemical Engineering from the University of Tehran,Iran,in 2018 and 2021.Experimental investigation of thermal transport in micro/nanoscale materials has been her research focus.Currently,she is conducting research on optical-acoustic phonon temperature difference of 2D materials using Raman Spectroscopy.

Xinwei Wangis the Anson Marston Distinguished Professor and Wilkinson Professor in Interdisciplinary Engineering at Iowa State University.He obtained his Ph.D.degree from the School of Mechanical Engineering,Purdue University in 2001,M.S.(1996) and B.S.(1994) from the University of Science and Technology of China.Over the past 20 years,he has led his laboratory to develop new techniques for characterizing energy transport at the micro/nanoscale.His current research focuses on conjugated energy transport in 2D materials and interfaces.