Next Article in Journal
Exploring the Antibacterial Potential of Green-Synthesized MgO and ZnO Nanoparticles from Two Plant Root Extracts
Next Article in Special Issue
Electrical and Thermal Conductivities of Single CuxO Nanowires
Previous Article in Journal
Synthesis of Zinc Oxide Nanorods from Zinc Borate Precursor and Characterization of Supercapacitor Properties
Previous Article in Special Issue
Scan-Free GEXRF in the Soft X-ray Range for the Investigation of Structured Nanosamples
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Measurement of Thermal Conductivity by SThM Technique: Measurements, Calibration Protocols and Uncertainty Evaluation

1
Laboratoire National de Métrologie et d’Essais (LNE), 29, Avenue Roger Hennequin, 78197 Trappes, France
2
Detection, Sensors and Measurements Laboratory, Ifremer, 1625 Route de Sainte-Anne, 29280 Plouzané, France
3
CETHIL UMR5008, CNRS, INSA-Lyon, Université Claude Bernard Lyon 1, 69621 Villeurbanne, France
*
Author to whom correspondence should be addressed.
Nanomaterials 2023, 13(17), 2424; https://0-doi-org.brum.beds.ac.uk/10.3390/nano13172424
Submission received: 29 July 2023 / Revised: 18 August 2023 / Accepted: 21 August 2023 / Published: 25 August 2023
(This article belongs to the Special Issue Metrology for Energy Nanomaterials)

Abstract

:
Thermal management is a key issue for the downsizing of electronic components in order to optimise their performance. These devices incorporate more and more nanostructured materials, such as thin films or nanowires, requiring measurement techniques suitable to characterise thermal properties at the nanoscale, such as Scanning Thermal Microscopy (SThM). In active mode, a hot thermoresistive probe scans the sample surface, and its electrical resistance R changes as a function of heat transfers between the probe and sample. This paper presents the measurement and calibration protocols developed to perform quantitative and traceable measurements of thermal conductivity k using the SThM technique, provided that the heat transfer conditions between calibration and measurement are identical, i.e., diffusive thermal regime for this study. Calibration samples with a known k measured at the macroscale are used to establish the calibration curve linking the variation of R to k. A complete assessment of uncertainty (influencing factors and computational techniques) is detailed for both the calibration parameters and the estimated k value. Outcome analysis shows that quantitative measurements of thermal conductivity with SThM (with an uncertainty value of 10 % ) are limited to materials with low thermal conductivity ( k < 10 W   m 1   K 1 ).

1. Introduction

The issues of thermal management in electronic devices (cell phones, laptops, batteries…) became more and more important with the progressive miniaturisation of components observed during the last decades. Understanding heat transfer processes in the corresponding nanostructured materials and devices requires accurate knowledge of the thermal properties of the materials used at the nanoscale. To meet these needs, several optical and near-field techniques have been developed to obtain local thermal information on a small scale [1]. Among these measurement methods, the Scanning Thermal Microscopy (SThM) [2,3] is the most used technique as it enables to reach a lateral resolution less than 100 nm with appropriate probes [4,5], whereas optical approaches such as photothermal radiometry [6], thermo-reflectance [7], and photo-reflectance [8] are limited by light diffraction and have thus higher lateral spatial resolution. Based on conventional atomic force microscopy (AFM) equipped with a miniaturised thermal sensor, SThM devices have been developed actively since the 1990s in order to operate either in passive mode for surface temperature measurements [9] or in active mode for phase transition detection [10], thermal contact resistance [11], and thermal conductivity contrast imaging [12,13,14].
SThM is, in addition, a promising technique for performing quantitative measurements of thermal conductivity at the nanoscale [15]. During the last decade, abundance experimental measurements have been performed with the SThM technique to quantify the thermal properties of nanoscale structures. For example, the thermal conductivity of nanowires embedded in a matrix was estimated [16,17]. SThM has also been used to measure the thermal conductivity of graphene oxide nanosheets [18] or ZnO thin films [14,19]. Unfortunately, most of the results in the literature are given without associated uncertainty or with uncertainties corresponding only to the range or the standard deviation [20] of multiple measurements, and are not traceable to the International System of Units (SI).
The objective of this paper is to provide complete guidance for the calibration of SThM, for traceable estimation of the thermal conductivity k at nanoscale and the evaluation of its associated uncertainty using the SThM apparatus, which is a new contribution in this field. This assessment is not trivial because the measurement of the thermal conductivity k of a sample by SThM technique is an indirect measurement method [21] which can be represented with the following implicit measurement model [21]
h ( k , Y , Q 1 , , Q n ) = 0 ,
where h is a function that relates the measurand k (the quantity intended to be measured [22]) to an intermediate measurand Y linked to the variation of the thermo-electrical response of the probe (such as the electrical resistance, the temperature, and the electromotive force (EMF), depending on the type of probe used), and to other quantities Q i , i = 1 , , N involved in the measurement. The elements h, Y and Q i of (1) are described thereafter.
In SThM active mode, the probe acts both as a heater and as a sensor. In steady state, the hot probe exchanges a constant heat flow corresponding to a G T o t thermal conductance (reciprocal of the thermal resistance) with its surroundings. In air, when the probe is far from the sample, the heat flux dissipates only through the cantilever and the air. When the probe is in contact with a sample, an additional channel is opened for the heat flux to the sample. This thermal flux is characterised by a thermal conductance G s a m p l e and is function of the thermal conductivity k of the sample among other influencing parameters. The variation Y of the thermo-electrical response of the probe between the measurements performed far from the sample and measurements performed in contact with the sample is representative of this thermal conductance G s a m p l e and by the way of k of the sample. Two steps can then be followed in order to determine the measurement function h in (1) linking the intermediate measurand Y to the measurand k and the Q i .
The most challenging way is to build a complete theoretical physical model describing the heat transfer between the probe and the sample that remains complex and involves many influencing quantities [15]. The easiest way that is used in this paper is to build a calibration curve [23,24,25,26,27], based on the measurements of a set of bulk calibration materials with well-known thermal conductivities, and to fit the experimental points with a model deduced from phenomenological studies of the measurements [27,28]. The use of this model and of the calibration curve is based on a strong assumption: the heat transfer conditions are the same during calibration on bulk calibration materials as during measurement on the studied material. The advantage of this method is that it can ensure measurement traceability, provided that the thermal conductivity measurements of the calibration materials are themselves traceable. The International Vocabulary of Metrology (VIM) defines metrological traceability as the “property of a measurement result whereby the result can be related to a reference through a documented, unbroken chain of calibrations, each contributing to the measurement uncertainty” [22].
The evaluation of the uncertainty associated wih the estimated value of k requires the evaluation of uncertainties associated with ( Y , Q 1 , , Q n ) and the inversion procedure. Due to the absence of an explicit relationship between the measurand k and the uncertainty sources ( Y , Q 1 , , Q n ) , it is challenging to perform standard uncertainty evaluation using the traditional law of propagation of uncertainties [29]. One solution is to apply Bayesian statistical inversion in order to both predict the value k * of k corresponding to an observed value y * of Y and estimate the parameters of the SThM calibration curve and their associated uncertainties. The interest of this approach is the joint estimation of the calibration parameters and the thermal conductivities from fresh SThM measurements in addition to the calibration data. The marginal posterior distribution π ( k * | y , y * ) can be considered as the measurement result, from which can be extracted point estimates, standard uncertainties, and 95% credibility intervals.
Previous works performed at LNE (Laboratoire National de Métrologie et d’Essais) [30,31] present measurement and calibration procedures with the aim of performing quantitative and traceable measurements of k at nanoscale by the SThM technique. A measurement protocol named “dark mode”, where the laser of the optical force detection system of the SThM is switched off, has been proposed to avoid the bias induced by the overheating of the thermoresistive probe due to the laser beam [31]. An analysis of the measurement process and the uncertainty associated with the intermediate measure (corresponding to the electrical resistance variation Δ R of the probe) have been described in detail [30]. Then, a new intermediate measurand Y was defined, and the measurement and calibration protocols were improved in order to build the calibration curve (using a theoretical model and experimental measurements on bulk materials) to ensure the traceability to the International System of Units (SI) of measurements at nanoscale.
Based on Bayesian inversion, this paper presents the first complete uncertainty assessment for calibration of SThM apparatus for thermal conductivity k measurements. It details the measurement, calibration, and uncertainty assessment processes at the macro-scale and investigates the determination of thermal conductivity at the nano-scale, as well as the assessment of the corresponding measurement uncertainty as illustrated in Figure 1.

2. Materials and Methods

SThM is an atomic force microscope (AFM) with an instrumented probe acting as a thermal sensor. Therefore, SThM offers the possibility to characterise surfaces at the nanometric scale in terms of topography, such as standard AFM, as well as in terms of thermal properties. In order to perform thermal conductivity measurements, the probe acts both as a sensor and as a heater. When the probe maps the surface of a sample, its electrical resistance changes as a function of its own temperature, depending on, inter alia, heat transfer between the probe and the sample. In this section, the equipment, the used resistive probes, the SThM measurement protocol, the calibration protocol elaborated to perform SI traceable measurement, and the uncertainty assessment method are presented.

2.1. Measurement Equipment

2.1.1. SThM System

The LNE’s scanning thermal microscope is a commercial NTEGRA scanning probe microscope from the NT-MDT company. The dedicated Nova-Px software-v2013 enables accurate displacements of sample against the probe using x, y, and z directions of piezoelectric scanners. The contact force between the probe and the sample is controlled using a classical optical feedback system that includes a laser diode beam reflected on the cantilever and a quadrant of photodiodes. In addition to the standard deflection signal of the laser beam, the oscilloscopes included in the software enable the recording of signals from the thermal unit as current intensity or voltages.

2.1.2. Resistance Temperature Probe

The resistance temperature probes used for this study are constituted from a SiN-grooved cantilever with gold pads. A fine resistive ribbon of palladium (about 150 μ m length) is deposited on the tip of the cantilever, as illustrated in Figure 2. These probes, called KNT probes, are provided by Kelvin Nanotechnology. The new generation of KNT probes has a nominal electrical resistance R ( T a m b i e n t ) in the range of 340 Ω to 450 Ω (measurements performed on a set of 12 probes).

2.1.3. Thermal Unit

The probe is included in a homemade thermal unit that encloses a current generator whose amplitude I S T h M extends from 750 µA to 1350 µA and a Wheatstone bridge that is a well-established technique for accurate electrical resistance measurements. Figure 3 describes the homemade thermal unit with a scheme of the Wheatstone bridge adapted to the KNT resistance temperature probes. One leg of the bridge is comprised of the SThM resistance temperature probe in series with a fixed precision resistor R 1 . The other leg is comprised of an adjustable resistance R a in series with two precision resistors R 2 and R f . The precision resistor R 2 is fixed, and R f is a calibrated decade resistor (DB62-11K from IET Labs) with a Kelvin type 4-terminal configuration. The value of this resistor is adjusted to the lower decade of the nominal electrical resistance R ( T a m b i e n t ) for each probe (as an example, for a resistance temperature probe with a nominal electrical resistance of 412 Ω , value of R f is adjusted to 410 Ω ). The value of the adjustable resistance R a is manually monitored by a rotary knob arbitrarily scalable in 1000 graduations; the level of graduation for the knob is denoted as B B k in the following. The bridge balance voltage B B v is amplified by a factor A (110). Once R a is adjusted to meet the following equality:
R · R 2 = R 1 · [ R a + R f ] ,
the Wheatstone bridge is balanced, and the measured bridge balance value B B v , recorded by an external voltmeter, equals zero. An external voltmeter is used rather than the oscilloscope of the Nova-Px software-v2013 because an offset of 58.2 mV (amplified value by a factor 110) on the oscilloscope measurements when the bridge is balanced has been identified. This offset is due to parasitic resistance (the unequal resistance of the wires). In addition to the B B v voltage data recording, an external voltmeter enables the measurement of the probe voltage U. A complete description of the SThM setup can be found in paper [30]. Based on the knowledge of the electrical resistances ( R a , R f , R 1 , R 2 ) involved in the bridge, the amplificator factor A, the measurement of B B v voltage, and the U probe voltage, the electrical resistance R of the probe can be determined following Equation (3):
R = U R a + R f + R 2 ( R a + R f ) · R 1 R 2 R a + R f + R 2 U B B v A .

2.2. SThM Measurements

2.2.1. Active Mode Configuration

The SThM technique can be used in two different modes: the passive mode, where the probe is cold and comes in contact with the hot sample surface for temperature measurements, and the active mode, dedicated to thermal conductivity measurement, where the probe is hot, self-heated by the Joule effect, and comes in contact with the cold sample surface. In our configuration, we performed measurements in active mode, and our SThM measurement consists of estimating the variation of heat losses from the heated probe to its surrounding environment between the two configurations “out of contact” with the sample and “in contact” with the sample:
  • “out of contact” abbreviated in “oc” where the probe is placed far from the thermal influence of the sample. Furthermore, the electrical resistance R of the probe mostly depends on the convective and conductive heat losses between the probe and the ambient air, the current intensity in the probe, the conductive heat losses from the probe to the cantilever, and the heat induced by the laser diode beam illuminating the cantilever. The radiative heat losses can be neglected.
  • “in contact” abbreviated in “ic” where the probe is in contact with the sample surface. In this configuration, R depends on the same influencing parameters as in the “out of contact” configuration and on the heat transfers between the probe and the surface sample. These are functions of the thermal properties of the sample and the interface thermal resistance between the probe and the sample.
Bycomparing the probe signals in “oc” and “ic” configurations, a measured quantity value [22] Δ R can be defined associated with the probe used. This measured quantity value Δ R actually corresponds to a temperature drop of the probe between the “oc” and “ic” configurations:
Δ R = R o c R i c
and is, inter alia other parameters, function of the thermal conductivity k of the sample.

2.2.2. Definition of the Intermediate Measurand

As highlighted in the previous paragraph, the probe electrical resistance R ( T ) depends on its temperature T and follows a quite linear relationship with its temperature for low variation of temperature:
R ( T ) = R ( T a m b i e n t ) · [ 1 + α ( T T a m b i e n t ) ] ,
with α the temperature coefficient of the probe material and R ( T a m b i e n t ) the electrical resistance of the probe at room temperature. The measured quantity value Δ R is then expressed by:
Δ R = R ( T a m b i e n t ) · α · ( T o c T i c ) .
Unfortunately, experiments highlight the long-term drift of the value of R ( T a m b i e n t ) depending on the ageing of the probe. Indeed, wear of the tip affects the R ( T a m b i e n t ) value on many levels. Many landings cause wear of the apex that changes both the size of the solid-solid contact between the tip and the surface of the sample and the length of the resistive Pd ribbon. In addition, a high value of the current intensity can induce matter transport resulting from electron displacement. This can cause the Pd ribbon to be damaged. In metallic thin layers, the electro-migration is prevalent. One way to correct the influence of this drift is to work with a reference sample. By performing measurements on the reference sample before each measurement on a studied sample, we define an intermediate measurand Y m equal to the ratio between Δ R s a m p l e obtained on the studied sample and Δ R r e f obtained on the reference sample:
Y m = Δ R s a m p l e Δ R r e f = Δ R s a m p l e Δ R S i O 2 f = R o c s a m p l e R i c s a m p l e R o c S i O 2 f R i c S i O 2 f .
The reference sample must be chemically inert and must remain stable over time. Furthermore, SThM measurements require perfect surface conditions with the sample [14,27]. Finally, the material must have a thermal conductivity within the sensitivity range of the SThM technique, i.e., lower than 10 W   m 1   K 1 as we will discuss in Section 4.1. The reference sample selected by LNE is a piece of fused amorphous silica S i O 2 f with a roughness value Ra of 0.56 nm.
In order to evaluate the effect of repeatability and reproducibility conditions with regard to some uncertainty sources, our measurement process is repeated several times in different conditions described in Section 2.2.3. Each condition i = 1 , , 10 provides a corresponding measurement result Y m , i composed of a measured quantity value y m , i and its associated uncertainty u ( y m , i ) . Finally, the intermediate measurand Y is computed as the mean Y ¯ m of this set of ten individual measurement results Y m , i :
Y = Y ¯ m = i = 1 10 Y m , i 10 .

2.2.3. SThM Measurement Protocol

Sample Requirements

The sample dimensions must allow the landing of the SThM probe but be small enough to allow two samples to be inserted at the same time on the SThM sample holder. For this study, typical dimensions (diameter) were 0.5 and 1 cm. The sample must be protected from any effect likely to alter its properties (contamination, dust, oxidation, water contact, temperature changes). The contact between the tip and the sample is assumed to be identical from one sample to another, provided that the scanning surfaces are locally flat [27,32]. The sample surface should be mirror polished to avoid poor mechanical contacts between the tip and the surface and to assume that the contact resistance is the same for all samples.

Measuring Conditions

Temperature and relative humidity of the environment must be controlled as they influence directly the heat transfer between the tip and the environment. Humidity enhances the water meniscus formed at the tip/sample contact and, consequently, the thermal conductance of the contact [33].The SThM is located in an air-conditioned room where the temperature is set to 21 °C and the relative humidity is 50%. A station sensor located close to the scanning probe microscope monitors the environment. The variations are very slow over a day or a week, less than 0.5 °C and 3%, respectively. These conditions limit the influence of the hydrophilic nature of the different samples. The measurements were performed in air with thermal steady-state conditions. Some thermal drifts have been measured when performing measurements with the NTEGRA protection cover, which is quite airtight [30]. These thermal drifts are due to air conditioning and SThM heat sources such as step motors and piezo motors. A specific homemade enclosure viewed in Figure 4a has been developed to protect both the SThM and the thermal unit from air conditioning disruptions and to limit thermal drifts from internal heat sources.
Two samples are placed at the same time on the SThM sample holder, as illustrated in Figure 4b: the reference sample (fused silicon dioxide ( S i O 2 f ) and the studied sample. The two samples are close enough to assume that they are at the same temperature. The NTEGRA x, y directions of piezoelectric scanners provide displacement of the sample holder to place either of the samples under the probe.
Measurements require thermal steady-state conditions. Samples stay in the SThM for a stabilisation time of at least 2 h before starting measurements. To avoid any laser effects on bridge balance voltage B B v (probe overheating, thermal drift), measurements are performed in “dark mode” [31].
In dark mode, the laser of the probe guiding system is turned off to avoid continuous heating of the probe. The probe landing on the sample surface was performed step-by-step by the user. The contact of the probe with the sample surface is indicated by a strong discontinuity in the B B v voltage signal. The “oc” distance between probe and sample was 1000 steps (approximately 1.5 mm) during these experiments. This chosen distance of 1.5 mm is large enough to avoid the thermal influence of the sample in the dark mode.

Measurement Process

The measurements were performed with direct current (DC) heating of the resistive probe in a “constant current” configuration. At the beginning of the experiments, the value of the electrical current I S T h M is fixed arbitrarily to 900 µA. The value of the electrical current is maintained below the critical value specified by the supplier ( I S T h M < 2.5 mA ) in order to avoid damage or premature wear of the probe. In addition, a study of the influence of the value of the electrical current showed that working with an electrical current of 900 mA provided the best measurement repeatability (< 2 . 10 4 a . u . ) for the intermediate measurand Y [34]. The bridge balance is balanced by adjusting the value of R a with the rotary knob B B k , this adjustment of the B B k is strictly kept constant during all of the calibration experiments. For the preliminary adjustments of I S T h M and B B v , the probe is in an “out of contact” configuration.
After these preliminary settings, the experiment sequence for each measurement is:
  • after stabilisation (criteria of standard deviation < 10 4 V for B B v , o c r e f mean value calculated with measurements performed during a 100 s period), start recording of B B v , o c r e f and U o c r e f signals during a 100 s period with the probe in an “out of contact” configuration above the reference S i O 2 f sample;
  • land with “dark mode” on one position of the S i O 2 f reference sample, wait for stabilisation (with the same criteria as for the first step), record B B v , i c r e f and U i c r e f signals during a 100 s period with the probe in an “in contact” configuration;
  • remove the probe from contact and wait for stabilisation (with the same criteria as for the first step), record B B v , i c r e f and U i c r e f signals during a 100 s period with the probe in “out of contact” configuration; repeat these three operations for two other landings at the same position above and on the S i O 2 f reference sample (repeatability of measurements).
  • After 3 measurements at the same position on the S i O 2 f reference sample, from “out of contact” configuration, move to another position above the sample (reproducibility measurements). After stabilisation (criteria of standard deviation < 10 4 V for B B v , o c r e f mean value calculated with measurements performed during a 100 s period), start recording of B B v , o c r e f and U o c r e f signals during a 100 s period with the probe in “out of contact” configuration for the new position above the reference S i O 2 f sample;
  • land with “dark mode” on the new position of the S i O 2 f reference sample, wait for stabilisation (with the same criteria as for the first step), record B B v , i c r e f and U i c r e f signals during a 100 s period with the probe in “in contact” configuration;
  • remove the probe out of contact, waiting for stabilisation (with the same criteria as for the first step), record B B v , i c r e f and U i c r e f signals during a 100 s period with the probe in “out of contact” configuration;
  • from “out of contact” configuration, move to another position above the sample, and repeat the two steps described in the two last bullets for this third location on the sample.
  • After measurements on the S i O 2 f reference sample, perform measurements on the studied sample following the same protocol as for the S i O 2 f reference sample.
An example of the experimental U and B B v signal recording is presented in Figure 5. The five steps on the curve correspond to the measurements on the S i O 2 f reference sample, and the last five steps correspond to the measurements on the studied alumina calibration sample.
We obtain eleven B B v r e f and eleven U r e f measurement data (five in “in contact” configuration and six in “out of contact” configuration) on a S i O 2 f reference sample associated with, respectively, eleven B B v s a m p l e and eleven U s a m p l e measurement data on the studied sample. Using this measurement data, for each sample (reference and studied samples), we calculate eleven electrical resistance values for the probe following Equation (3) (five in “in contact” configuration and six in “out of contact” configuration). Then, we determine ten resistance differences per sample as defined in Equation (4). A regular residual drift is observed both on U and B B v measurements, mainly visible for “out of contact” data (represented by the red dashed line in Figure 5). This residual drift is induced by the electronics of the SThM (i.e., the motors of the piezoelectric stage). In order to reduce the influence of this drift of signals, we compute the y m , i (i from 1 to 10) measured quantity values from each individual data in a chronological way ( U r e f and B B v , r e f data from the first landing on the reference sample associated to U s a m p l e and B B v , s a m p l e data from the first landing on the studied sample and so one for next data. By this way, we obtain ten measurement results Y m , i , which are used to calculate the resulting Y intermediate measurand as the mean value following Equation (8).
To check the potential influence of landing or withdrawal conditions, we identified two types of resistance differences. The first one, Δ R l a n d i n g , identified as “landing”, corresponds to the resistance difference between measurements performed in the “out of contact” condition before the contact and the measurements performed during the “in contact” configuration. The second one, Δ R w i t h d r a w a l , identified as “withdrawal”, corresponds to the resistance difference between measurements performed in the “out of contact” condition after the contact and the measurements performed during the “in contact” configuration.
Δ R l a n d i n g = R o c , B e f o r e C o n t a c t R i c
Δ R w i t h d r a w a l = R o c , A f t e r C o n t a c t R i c
In conclusion, we compute ten measured quantity values y m , i for the intermediate measurand Y with repeatability and reproducibility conditions. The first six successive measurements. y m , i with i from one to six, are performed on the same location; the last four measurements, y m , i with i from seven to ten, are performed on two other locations in order to evaluate the potential influence of heterogeneity of the sample (reproducibility conditions). In addition, we also studied the potential influence of landing measurement conditions associated with odd i indexes and withdrawal measurement conditions associated with even i indexes. Measurement conditions and identification are summarized in Table 1.

2.3. SThM Calibration Protocol

Two ways are available in order to link the intermediate measurand Y to the thermal conductivity of the studied sample. The first solution is to establish a physical model describing heat transfers between the probe and the surface sample. Unfortunately, it is challenging to develop rigorous physical models at the nanoscale because the values of the influencing parameters are difficult to evaluate at these scales and not always reliable. Another solution is to establish a calibration curve based on materials with thermal conductivities measured at the macroscopic scale [23,24,27,28,35,36,37]. However, this method is based on a strong assumption, as highlighted in Section 1, that the heat transfer conditions are the same during both the calibration and sample measurement steps. Since calibration is performed at the macroscopic scale, heat transfers are diffusive. Particular care must be taken when measuring nanomaterials with a mean free path higher than the size of the studied structure (involving ballistic heat transfers).
The calibration function is mainly based on phenomenological study, and the calibration parameters are identified with Bayesian inversion from the experimental values and the mathematical model.

2.3.1. Definition of the Calibration Model

Based on phenomenological studies [24,28,38], a mathematical link has been established between the variation of the probe electrical resistance when it is put in contact with the material and the thermal conductivity of the material. These studies are based on the heat transfer from the hot probe (heated by the Joule effect) to the “cold” sample.
The thermal conductances involved in the probe-sample system during measurements are the thermal conductance relating to the probe heat losses through convection and radiation G e n v , the thermal conductance relating to the probe heat transfer through the cantilever G c a n t , the effective thermal conductance of the thermal contact between the probe and the sample (solid-solid thermal conduction, air and water meniscus thermal conduction G c o n t a c t and the thermal conductance of the sample G s a m p l e .
During “out of contact” measurement, the total thermal conductance G t o t , o c of the probe-sample system is:
G t o t , o c = G e n v , o c + G c a n t .
During “in contact” measurement, the total thermal conductance G t o t , i c of the probe-sample system is:
G t o t , i c = G e n v , i c + G c a n t + G c o n t a c t . G s a m p l e G c o n t a c t + G s a m p l e .
The variation of the thermal conductance of the probe-sample system between “oc” and “ic” configurations is consequently:
G t o t , i c G t o t , o c = G c o n t a c t . G s a m p l e G c o n t a c t + G s a m p l e + ( G e n v , i c G e n v , o c ) .
Assuming that heat flow in the sample is diffusive, that means that the mean-free path Λ of the energy carriers, defined as the average distance between two successive inelastic shocks of energy carriers, is lower than the radius b of the thermal contact between the probe and the sample; the thermal conductance G s a m p l e is given by:
G s a m p l e = 4 β k b
with β a specific coefficient corresponding to the influence of the shape of the thermal contact. As a result, the variation of thermal conductance Δ G t o t can be written as:
Δ G t o t = G c o n t a c t · k G c o n t a c t 4 β k b + k + Δ G e n v
Therefore, the variation of the electrical resistance of the probe between “oc” measurement and “ic” measurement can be described as:
R o c R i c = a R k b R + k + c R
where a R , b R and c R are parameters relative to the different thermal conductances involved during measurements. As R o c S i O 2 f R i c S i O 2 f does not depend on sample thermophysical properties, the k dependence of Δ R S a m p l e Δ R S i O 2 f remained unchanged when divided by a constant value. The calibration law of the intermediate measurand Y obeys thus the following law:
Y = Δ R S a m p l e Δ R S i O 2 f = a k b + k + c
Parameters a, b and c are estimated using Bayesian identification in order to estimate the uncertainty associated with these coefficients, which take into account the uncertainty associated with both the thermal conductivity of known materials and the intermediate measurand considered in this paper.
A representation of the calibration curve is provided in Figure 1, and the corresponding calibration model is identified as Equation (17).

2.3.2. Calibration Materials

Drawing on the results from the Quantiheat project [32], the choice of the calibration samples was based on their thermophysical properties, assuming that the thermal conductivity measured at the macroscale and at the nanoscale are comparable in terms of their mechanical stability and low roughness. To ensure the traceability of measurements, the materials have been chosen so that their thermal conductivity can be measured using traceable techniques. This requires them to be homogeneous and isotropic. Twelve samples of twelve different bulk materials (polymers, ceramics, and pure metals) with thermal conductivity k ranging from 0.1 W   m 1   K 1 to 100 W   m 1   K 1 have been selected for this study: poly(methyl methacrylate) (PMMA), poly-oxymethylene in copolymer (POM-C), borosilicate glass, two different grades of fused silicon dioxide ( S i O 2 f ) and ( S i O 2 N E G S 1 ), zircon oxide ( Z r O 2 ), titanium oxide ( T i O 2 ), alumina (poly crystal aluminium oxide A l 2 O 3 p ), sapphire ( α A l 2 O 3 ), germanium (Ge), p doped silicon ( S i p + + ), and zinc (Zn). All samples were mirror-polished to minimise roughness. The polymers were prepared by ultramicrotomy (cryogenic cutting). Table 2 presents materials, dimensions, structures, providers, measured thermal conductivity, and measured roughness of selected materials. Each sample is in the shape of a disc with a 10 mm diameter, except for the two polymer samples, which have a flat surface prepared by ultramicrotomy of around 1 mm 2 .
The thermal conductivity of each material was determined at 23 °C on bulk specimens (a few mm thick) using an indirect and traceable method. This method is based on the measurements of the thermal diffusivity by the laser flash method [39,40], the specific heat by differential scanning calorimetry, and the density by the Archimedean method. The expanded uncertainty associated with the measurement of thermal conductivity by this indirect method has been estimated at 5 % [41].

2.4. Method for the Evaluation of the Uncertainty Associated with the Estimation of the Intermediate Measurand

The intermediate measurand Y considered in this paper is the ratio of the resistance difference for the tested sample between the “ic” and “oc” configurations and the same difference obtained for the reference sample S i O 2 f (see Equation (8)). In order to evaluate dedicated uncertainty from repeatability and reproducibility conditions, Y is obtained as the combination of individual quantity measured values y m and their associated uncertainties. Based on the protocol described in Section 2.2.3, we obtained ten y m , i values for various conditions as described in Table 1. Some of these values are obtained from repeatability conditions (measurements at the same location), and others from reproducibility conditions (measurements at different locations or with two different configurations: landing or withdrawal conditions). In this way, we study the individual contribution of each uncertainty source (instrumental, repeatability, and reproducibility) in Section 3.2.
The estimated value of the intermediate measurand Y and its associated uncertainty include repeatability and reproducibility conditions influences. The different steps in the evaluation of the measurement uncertainty associated with Y are presented in Figure 6.
The following sections (Section 2.4.1, Section 2.4.2 and Section 2.4.3) describe the evaluation of measurement uncertainty on each individual measurand Y m , based on works (calculation, procedures, and conclusions) performed previously at LNE by Ramiandrisoa et al. [30], in which the evaluation of the measurement uncertainty associated with a single Δ R was considered. Then the next section (Section 2.4.4) describes the combination of y m , i measurements and uncertainty to report the intermediate measurand Y value for various conditions with associated standard uncertainty.

2.4.1. Modelling the Measurement Process for Individual Measurand

Four resistance measurements are involved in the mathematical model of the individual measurand Y m value as described in (7), each being determined thanks to the Wheatstone bridge in Figure 3. As a result, the expression for one of these four resistances (denoted as j from 1 to 4) is:
R j = U j R v + R 2 R v R 1 R 2 R v + R 2 U j B B j A
U j is the voltage (in V ), B B j is the bridge balance (in V ), R 1 is a fixed resistance of the bridge, given by the technical specifications (in Ω ), R 2 is a fixed resistance of the bridge, given by the technical specifications (in Ω ),
A = R 10 k R 1 k R 1 k + R 10 k R 1 k
is the amplification factor of the electrical chain (in Ω ) and
R v = R f + R v , m i n + B B k B B k m i n B B k m a x B B k m i n R v , m a x R v , m i n
is the variable resistance (in Ω ). B B k = 125 denotes the graduation of the know, with minimum and maximum graduations respectively denoted as B B k m i n = 0.5 and B B k m a x = 1003 . All three quantities have a negligible contribution to the uncertainty and are consequently considered as fixed.

2.4.2. Evaluating Input Quantities for Individual Measurand

In order to perform the Monte Carlo simulation, it is required to assign a Probability Density Function (PDF) to each input quantity in the mathematical model. The following subsections describe the choice made for every single input quantity.
  • Voltages:
The voltage measurements are denoted as U and B B v in Figure 3. Two 34401A multimeters measured one U voltage and the other one B B v voltage. Three errors are considered to be the most influential ones: the trueness, the quantification, and the repeatability.
  • Trueness of the multimeters: This error is the same for each measurement of a voltage, whether the sample is in or out of contact, and whether the unknown sample or the reference sample is measured, but is specific for each multimeter. Available information about the trueness error comes from the calibration certificate of each multimeter. These calibration certificates provide trueness corrections U t r u e and B B v , t r u e with an associated expanded uncertainty U ( U t r u e ) = U ( B B v , t r u e ) = 2.5 μV, using a coverage factor k = 2 . This correction is applied to the measurements, and a Gaussian probability distribution is assigned with a zero mean and
    u U t r u e = U U t r u e k = 1.25 × 10 6 V
    and
    u B B v , t r u e = U B B v , t r u e k = 1.25 × 10 6 V
    as standard deviations.
  • Quantification of the multimeters: The multimeters have the same quantification step q = 1 μV in the studied range. As a consequence, the quantification error lies in the interval 5 × 10 7 ; 5 × 10 7 . A rectangular probability distribution is assigned. However, this (unknown) error may be different for each voltage measurement. As a result, we define a different input quantity for each different voltage measurement.
  • Repeatability: In order to evaluate the repeatability of the voltage measurement, our measurement corresponds to the mean values U ¯ and B B v ¯ of the respective U voltage and B B v voltage for 100 measuring points (corresponding to a period of 100 s) associated with their respective standard deviations.
  • Measurement model for voltages: As a result, the measurement model used for each voltage measurement (in contact/out of contact) is:
    U ¯ i = U t r u e + U i q + U i R
    and
    B B ¯ v , i = B B v , t r u e + B B v , i q + B B v , i R
    where i is either equal to i c for the “in contact” voltage or to o c for the “out of contact” voltage, U t r u e and B B v , t r u e are the trueness corrections, U i q and B B v , i q are the quantification errors and, U i R and B B v , i R are the repeatability errors. An example of corresponding parameters are summarised in Table 3 for the PMMA sample.
  • Resistances involved in the Wheatstone bridge
In the mathematical model of the measurement process, 9 resistances are involved ( R 1 , R 2 , R f , R v , m a x , R 10 k , R 1 k , R 1 k , R 1 k and R 10 k ). In this section, we present the general mathematical model used to evaluate the uncertainty associated with a resistance, and we provide the summary for each considered resistance.
A resistance R is obtained as a reading with an associated trueness error, with a dilatation correction factor taking into account the difference between the temperature T L in the laboratory and the reference temperature T r e f = 25 °C:
R i = R 0 + T C R · ( T L T r e f ) · R 0 = R 0 + T C R · Δ T · R 0
The trueness of the resistance is evaluated thanks to the technical specifications: R 0 = R n o m ± a % , where R n o m is the nominal value of the resistance. As a consequence, R 0 is assigned a rectangular probability distribution over the interval [ R n o m a / 100 · R n o m ; R n o m a / 100 · R n o m ] . The Temperature Coefficient Ratio ( T C R ) is considered a fixed value. Regular measurements of the temperature in the laboratory show that it lies between 20.3 °C and 20.5 °C. As a consequence, Δ T is assigned a rectangular probability distribution over the interval [−4.5;−4.7].
In addition to the known resistances involved in the mathematical process, there is a parasitic resistance R p a r a s i t i c that is the sum of the parasitic resistances of the probe R p p r o b e and that of the bridge R p b r i d g e . This latter can be decomposed as the sum uncertainty sources: the reproducibility R p b R and the quantification error R p b q :
R p a r a s i t i c = R p p r o b e + R p b r i d g e = R p p r o b e + R p b R + R p b q
However, they are not considered in our mathematical model as they do not have any significant contribution to the uncertainty of the intermediate measurand when we consider the SThM indication to be Δ R S a m p l e Δ R S i O 2 f .
An example of assigning probability distributions to all input quantities involved in y m measurement process is given in Table 3 in Section 3.1.

2.4.3. Propagating Distributions for Individual Measurand

The evaluation of associated uncertainty y m measurement is performed using the Monte Carlo method (MCM), according to the principles of Supplement 1 to the GUM (Guide to the Expression of Uncertainty in Measurement) [22,42]. The previous quantification of the input quantities, described in Section 2.4.2, consisted in the choice of suitable Probability Distribution Functions (PDFs) based on the work of Ramiandrisoa [30]. Implementation of the propagation of these distributions is performed using LNE-MCM software-v2017 [43], with M = 10 6 simulations in order to get an overview of the PDFs for y m measured quantity value (with a best estimate, a standard uncertainty, and a coverage interval).
As a result, each individual measured quantity value y m is characterised, and the associated uncertainty of individual y m is evaluated.

2.4.4. Combining Reproducibility Measurements

In this study, the reproducibility measurements y m , i , u ( y m , i ) i = 1 , , 10 are combined using the random effects model defined as
y m , i = μ + λ i + ε i
where μ is the overall arithmetic mean response, the { λ i } i = 1 , , 10 are the effects of the reproducibility conditions (which are assumed to have a Gaussian distribution with mean 0 and standard deviation τ ), the { ε i } i = 1 , , 10 are random effects assumed Gaussian with mean 0 and standard deviation the reported standard uncertainties { u ( y m , i ) } i = 1 , , 10 , the λ i and the ε i are assumed to be independent.
Such a model is generally used to combine heterogeneous data, e.g., from interlaboratory studies or meta-analyses, where the λ i are referred to as “laboratory effects”. The parameter τ is often called “dark uncertainty” and is used to capture excess variability in the measurements with respect to the reported uncertainties.
In this application, we aim to identify all sources of uncertainty during measurements under both repeatability and reproducibility conditions. We suspect other influencing factors than material heterogeneity or landing/withdrawal conditions. The main suspected factor is the variation in force applied from one measurement to another, which could not be quantified within the scope of this study. In our case, “dark uncertainty” has a physical meaning related to underestimated u ( y m , i ) uncertainties.
In this study, the parameters μ are τ are estimated using the Bayesian approach of [44] with a non informative prior on μ , π ( μ ) 1 and the following prior for τ as recommended in the NIST Consensus Builder [44]: τ HalfCauchy ( mad ( { y m , i } i ) ) , where mad is the median absolute deviation.

2.5. Bayesian Approach to Estimate the Thermal Conductivity from SThM Measurements

2.5.1. Error-in-Variables Representation

In this section, we denote X = ( X 1 , , X N ) the vector of random variables modelling the thermal conductivities measured for the bulk materials and Y = ( Y 1 , , Y N ) the vector of random variables modelling the corresponding SThM measurements. Due to the normality arising from the GUM uncertainty propagation for both SThM and thermal conductivity measurements, the following error-in-variables representation (a regression model that accounts for measurement errors in the independent variables) is used:
x i N ( X i , u ( x i ) ) ,
y i N ( Y i , u ( y i ) ) ,
where x i and y i are respectively the reported estimates of the thermal conductivity and the SThM measurement for material i, u ( x i ) and u ( y i ) are their associated uncertainties respectively.
The relationship (17) between Y and X can be expressed as
Y i = h θ ( X i ) , i = 1 , , N ,
where h θ is the calibration curve and θ = ( a , b , c ) .
Denote X * the random variable representing the unknown thermal conductivity of a material under test and Y * the random variable modelling the reported SThM measurement result y * and its associated uncertainty u ( y * ) . Similarly to (29) and (30), we have
y * N ( Y * , u ( y * ) ) ,
Y * = h θ ( X * ) .
In the following, we denote μ , the vector containing all the reported SThM measurements
μ = ( Y * , Y 1 , , Y N ) = ( h θ ( X * ) , h θ ( X 1 ) , , h θ ( X N ) ) .
It is important to note that this approach can be readily extended to multiple (say M) new SThM measurements y * = ( y 1 * , , y M * ) in which case Y * and X * read as the vectors Y * = ( Y 1 * , , Y M * ) and X * = ( X 1 * , , X M * ) respectively.

2.5.2. Bayesian Paradigm

In the Bayesian paradigm [45], all quantities involved, namely X * , θ , X , Y , Y * , are viewed as parameters to be jointly updated in the inference process by the information contained in the measurements. It is important to note that, contrary to classical approaches to inversion, the SThM measurements performed on unknown materials y * are used to update knowledge on all parameters and are not only used for the prediction of X * . In other words, using a Bayesian approach allows us to simultaneously estimate the parameters θ of the calibration curve and make predictions from the curve.
According to (30), Y i explicitly depends on θ and X i , so that its posterior distribution can be obtained as a by-product of the Bayesian analysis from the samples of the joint posterior distribution of θ and X i . For this reason, we chose to remove Y and Y * from the Bayes formula.
Denoting x = ( x 1 , , x N ) and y = ( y * , y 1 , , y N ) , the Bayes formula gives the joint posterior distribution π X * , θ , X | x , y of all the quantities involved as
π X * , θ , X | x , y l x , y | X * , θ , X π X * , θ , X ,
where l x , y | X * , θ , X is the likelihood of the data and π X * , θ , X is the prior distribution of the parameters.
The overall objective of estimating the thermal conductivity from SThM measurements consists in estimating the so-called marginal posterior distribution of X * given all the measurements which reads π ( X * | x , y ) and is mathematically obtained from (34) as the integral
π ( X * | x , y ) = θ X π X * , θ , X | x , y d θ d X .

2.5.3. Likelihood

The likelihood can be factorized as
l x , y | X * , θ , X = l x | X × l y | X * , θ , X ,
where l x | X and l y | X * , θ , X are defined thereafter.
Denoting Σ x the covariance matrix of the Gaussian vector x, we have x N N X , Σ x and the associated part of the likelihood writes
l x | X 1 det Σ x exp 1 2 ( x X ) T Σ x 1 ( x X ) .
Denoting Σ y the covariance matrix of the Gaussian vector y, we have y N N + 1 ( Y , Σ y )
y * y 1 y N N N + 1 h θ ( X * ) h θ ( X 1 ) h θ ( X N ) , u 2 ( y * ) 0 0 0 0 u 2 ( y 1 ) 0 0 0 u 2 ( y N 1 ) 0 0 0 0 0 u 2 ( y N ) ,
and the associated part of the likelihood writes
l ( y | X * , θ , X ) 1 det Σ y exp 1 2 ( y μ ) T Σ y 1 ( y μ ) ,
where μ is defined in (33).
In this paper, we assume that the covariance matrices Σ x = diag ( u 2 ( x 1 ) , , u 2 ( x n ) ) and Σ y = diag ( u 2 ( y 1 ) , , u 2 ( y n ) ) are diagonal (meaning no covariance).

2.5.4. Prior Distribution

The joint prior distribution π X * , θ , X of all quantities involved can be expressed as the product
π X * , θ , X = π θ π X * π X ,
where π ( θ ) = π ( a ) π ( b ) π ( c ) expresses prior information on the parameters of the calibration curve a, b, c, π X * expresses prior belief on the sought thermal conductivity X * and π X represents the joint prior belief on the thermal conductivity measurements results stored in the vector X. If all measurements are assumed to be independent from each other, then π X can be factorised as the product of the individual densities of probability π ( X i ) .
A variety of distributions can be used to represent prior beliefs, ranging from poorly informative to informative [45]. In this study (see Section 3.3), we choose to use non-informative Jeffrey’s priors for π ( X ) namely π ( X i ) 1 , and mildly informative priors for π ( a ) , π ( b ) , π ( c ) , and π ( X * ) as Gaussian distributions with large variance.

2.5.5. Computing Posterior Distributions

Since analytical formulas for (34) and (35) are usually intractable, it is widely accepted to tackle the estimation of quantities of interest for these distributions (mean, standard deviation, coverage interval) using Markov Chain Monte Carlo (MCMC) simulations [46], from which estimations of the posterior distributions using histograms or kernel-based density estimates can also be obtained.
MCMC sampling is performed in software such as R or Python (among others), e.g., using RStan [47] as in this paper and requires statistical expertise for the tuning of the algorithms and the analysis of the results (correlation, convergence of Markov chains, etc.). A review of convergence diagnosis tools can be found in [48]. In this paper, we concentrate on the effective sample size and the Gelman-Rubin diagnostic (Rhat). In brief, the effective sample size gives the number of independent samples equivalent to a set of correlated Markov chain samples, and the output of the Gelman-Rubin diagnostic is the so-called potential scale reduction factor, which should be close to 1 and is computed from at least two chains running with over-dispersed starting points w.r.t. the posterior distribution. The great advantage of MCMC is that at each iteration, a sample from the joint posterior distribution is produced, while the collection of samples for each parameter is distributed according to its marginal posterior distribution.

3. Results

In this section, we implement the calibration methodology to obtain traceable estimates of the thermal conductivity using the SThM measurements and a Bayesian inversion procedure described in Section 2. We highlight that the calibration procedure is necessarily developed on bulk materials at macroscale, for which traceable thermal conductivity measurements can be obtained, and that the traceability at micrometric and nanometric scales is ensured by the SThM technique, provided that heat transfer regimes are the same between calibration and measurement, as will be discussed in Section 4.3. In Section 3.1, we display the measurement results obtained with the SThM on the calibration bulk materials presented in Table 2 according to the methodology developed in Section 2.2. In Section 3.2, we discuss the effect of the reproducibility conditions (landing and withdrawal conditions, heterogeneity of the sample) on the resulting SThM uncertainty. In Section 3.3, we illustrate the Bayesian methodology for both the identification of the parameters of the calibration curve and the prediction of traceable thermal conductivity. The approach is applied to the estimation of three thermal conductivities in the range [0.1–10] W   m 1   K 1 .

3.1. Experimental Measurements on Calibration Materials

Measurements have been performed on the twelve calibration samples presented in Table 2. As described in Section 2.2, a run of 5 measurements is performed on the S i O 2 f reference sample before each run of 5 measurements on the studied calibration sample.

Measurements of Input Quantities for Individual Measured Quantity and Their Associated PDFs

All input quantities ( B B V , U , B B k , ) have been measured as described in Section 2.4.2, and each measured quantity y m , i and the associated Probability Distributions Functions (PDFs) have been determined following the method described in Section 2.4.3. Table 3 summarises the PDFs assigned to the input quantities obtained for the y m , 1 measurement on the PMMA sample as an example.
Table 3. Summary of PDFs assigned to measured quantity value y m , 1 for the PMMA sample.
Table 3. Summary of PDFs assigned to measured quantity value y m , 1 for the PMMA sample.
Input QuantityUnitProbability DistributionMean ValueStandard DeviationLower BoundUpper Bound
U t r u e V Gaussian 0.004 1.25 × 10 6
U o c , q s V Rectangular 5 × 10 7 5 × 10 7
U o c , q r e f V Rectangular 5 × 10 7 5 × 10 7
U i c , q s V Rectangular 5 × 10 7 5 × 10 7
U i c , q r e f V Rectangular 5 × 10 7 5 × 10 7
U o c , R s V Gaussian 0.37547698 4.57 × 10 6
U o c , R r e f V Gaussian 0.37546387 4.08 × 10 6
U i c , R s V Gaussian 0.37529628 4.13 × 10 6
U i c , R r e f V Gaussian 0.37520403 4.03 × 10 6
B B v , t r u e V Gaussian 0.04 1.25 × 10 6
B B v , o c , q s V Rectangular 5 × 10 7 5 × 10 7
B B v , o c , q r e f V Rectangular 5 × 10 7 5 × 10 7
B B v , i c , q s V Rectangular 5 × 10 7 5 × 10 7
B B v , i c , q r e f V Rectangular 5 × 10 7 5 × 10 7
B B v , o c , R s V Gaussian 9.5697 × 10 3 5.66 × 10 5
B B v , o c , R r e f V Gaussian 8.2624 × 10 3 10.24 × 10 5
B B v , i c , R s V Gaussian 6.2327 × 10 3 5.07 × 10 5
B B v , i c , R r e f V Gaussian 14.4252 × 10 3 3.51 × 10 5
B B k a. u.Rectangular 124.5 125.5
B B k , m i n a. u.Fixed 0.5
B B k , m a x a. u.Fixed1003
R 1 Ω Rectangular9991001
R 2 Ω Rectangular9991001
R f Ω Gaussian 399.830 0.001
R v , m a x Ω Gaussian 198.119 0.001
R v , m i n Ω Gaussian 0.0938 0.001
R 1 k Ω Rectangular9991001
R 1 k Ω Rectangular9991001
R 10 k Ω Rectangular999910,001
R 10 k Ω Rectangular999910,001
R 1 k Ω Rectangular9991001
Based on these input quantities, we calculated the corresponding y m measured quantity value and determined the associated uncertainty by propagation of distributions as described in Section 2.4.3. As an example, the histogram for one value of y m , i is represented in Figure 7. It can be observed that the output quantity can be adequately described with a Gaussian behaviour. The Spearman study shows that the highest contribution to the variance on the individual measured quantity y m comes from the B B v measurements in “out of contact” conditions: B B v , w c s a m p l e and B B v , w c r e f contribute to 65% of instrumental variance on y m measurement. This highlights the need to properly manage the environmental conditions in order to maintain constant heat transfers between the probe and its surroundings during measurement.

3.2. Study of Influencing Factors Regarding Repeatability and Reproducibility Conditions of Measurement

In order to estimate the measurement precision, we performed replicate measurements on each sample under repeatability and reproducibility conditions. Usually, the measurement precision is expressed numerically as the standard deviation of the set of measurements performed under specified conditions. In the process of measuring the same sample, 5 stages of landing/withdrawal are performed, which yields 5 different measurements in each configuration (in contact/out of contact). It should be noted that, in order to take account of the potential heterogeneity of the sample, measurements are performed at three different locations on the sample. In addition, the potential influence of the landing or withdrawal of the probe has also been studied. As a result, the effects of the supposedly heterogeneity of the surface of the sample and the influence of probe movement are taken into account in the dispersion of the measurements.
As a result, we obtain ten measured values y m , i for each calibration sample with associated uncertainty. An example of ten measured values obtained for the PMMA sample is given in Table 4. From these results, we discuss in the following Section 3.2 the repeatability and reproducibility conditions, effects, and uncertainty associated with each Y.

3.2.1. Evaluation of Measurement Precision under Repeatability Conditions

For each calibration sample, we compute the standard deviation obtained with the repeatability condition for two sets of measurements: the landing condition measurement with the set composed of y m , 1 , y m , 3 , and y m , 5 identified measurements, and the withdrawal condition measurement with the set composed of y m , 2 , y m , 4 , and y m , 6 identified measurements y m , i as described in Table 1. Results are presented in Figure 8.
No specific trend is identified from Figure 8. Two materials ( A l 2 O 3 and Z n ) show higher standard deviation values in the repeatability condition (more than 1.5 10 2 [a.u.]) than the other materials. This is probably due to the influence of roughness. Indeed, these two materials present high roughness (respectively, 7.52 nm for A l 2 O 3 and 8.14 nm for Z n , see Table 2). Guen [27] has shown that surface roughness alters the mechanical contact between the tip and the surface sample, reducing the apparent contact radius. When the roughness of the sample increases, heat transfer through mechanical contact decreases by 30 % . As the signal decreases, the signal-to-noise ratio also decreases, and the dispersion of the measurements, evaluated by the standard deviation, increases with the roughness of the sample.
Nevertheless, it should be noted that the value of the standard deviation computed for the two sets of measurements (landing and withdrawal) is high compared with the instrumental uncertainty contribution. The values reported in Table 5 show that the relative uncertainty coming from the repeatability condition is at least double the relative uncertainty from instrumentation. This excessive dispersion of measurements suggests that the repeatability conditions are not fully respected. An influencing factor, not taken into account in this study, changes from one measurement to the next. This source of uncertainty will be discussed in Section 3.2.4.

3.2.2. Evaluation of Measurement Precision under Reproducibility Conditions: Study of Landing and Withdrawal Configurations

For each measurement, we compute the variation of resistance in two ways, as described in Section 2.2.3: landing conditions and withdrawal conditions. By this way, we check the potential influence of landing and withdrawal conditions on our measurements due to changes in heat transfers between the tip and its surroundings after contact with the sample (pollution or residual water film on the tip, variation of the surrounding temperature). We computed the mean value y ¯ m , l a n d i n g for the landing configuration with its associated expanded uncertainty and the mean value y ¯ m , w i t h d r a w a l for the withdrawal configuration with its associated expanded uncertainty. Results for the PMMA sample are presented in Figure 9. We see that there is no significant difference between the landing and withdrawal configurations for individual measurements. The mean values for y ¯ m , l a n d i n g measured quantity value in landing configuration and y ¯ m , w i t h d r a w a l in withdrawal configuration are comparable. The difference is in the same order of magnitude as the standard deviation, which means that there is no significant difference between landing and withdrawal configurations compared with the standard deviation computed from mean values.
The same analysis has been performed for all calibration samples in order to check if there is a significant trend in mean values or standard deviation between the landing configuration and the withdrawal configuration. The data are gathered in Figure 10. No significant trend has been observed either for the mean values or for the standard deviation values.

3.2.3. Evaluation of Measurement Precision under Reproducibility Conditions: Study of Heterogeneity of The Sample

Another influence factor that can induce an additional contribution to the uncertainty is the heterogeneity of the sample, which can be due to variations in roughness from one location to another, non-uniformity (variation of structure, size grain), oxide film, or pollution. In order to identify if the heterogeneity of our samples has a significant impact on our measurement spread, for each sample, we compare in Figure 11 the standard deviation obtained for the three measurements performed at the same location and standard deviation obtained for the measurements performed at three different locations.
No significant trend indicates that sets of points measured under reproducibility conditions have a higher contribution to the standard deviation than sets of points measured under repeatability conditions. To conclude, with our measurement protocol, no significant contribution to the standard deviation has been established either from the landing or withdrawal configuration or from the heterogeneity of the sample. That means other sources of uncertainty not clearly evaluated of the purpose of the previous experiments have to be considered.

3.2.4. Combination of Measurements in Repeatability and Reproducibility Conditions

For each studied sample, we obtained ten quantity-measured values Y m , i from the repeatability and reproducibility conditions. As discussed in previous sections, the variations in repeated observations of the measurand under apparently identical conditions in highlight underestimated the uncertainty associated with the measured quantity values. As a result, we computed all the ten values measured for each sample and determined the associated uncertainty by performing Bayesian consensus estimation of the mean value of each Y m , i as described in Section 2.4.4. In this way, we integrate all influencing factors, even non-identified uncertainty sources, as the variation of the force applied by the cantilever between each measurement or random variation of environmental conditions. As we perform landings without the laser, we have no feedback of the deflection of the cantilever, which means no fine control of the applied force. This can induce a strong variation in the interface thermal resistance between the tip and the surface of the sample and a strong variation in the thermal contact area. This could be the main influencing factor regarding our measurement protocol. Figure 12 presents the ten quantity measured values Y m , i obtained for the PMMA sample with their associated uncertainty as well as the mean value identified as the intermediate measurand y with its associated uncertainty.
The mean values and their associated uncertainty have been computed by performing Bayesian consensus estimation for each of the twelve calibration materials. Table 6 presents results, including intermediate measurand y, absolute uncertainty, and relative uncertainty.
As a result, our measurement protocol enables us to reach a relative standard uncertainty of 1.0 % at most. The highest values of uncertainty (≥0.7%) are obtained for materials with the highest roughness (Zn, Alumina, and POM-C). Regarding Z r O 2 , even if the roughness value of its surface is quite low (< 0.5 nm ), the relative standard uncertainty associated with the determination of Y Z r O 2 is estimated to be 0.7 % . A study of rough data shows that there was a variation in the thermal drift (discussed in Section 3.1 and illustrated on Figure 5) between the measurements performed on the reference S i O 2 sample and those performed on the Z r O 2 sample that increased the contribution of reproducibility to the uncertainty value. That confirms that special attention is required to the stability of environmental conditions during measurements in order to avoid increasing uncertainty.

3.3. Bayesian Identification of the Parameters

The Bayesian analysis proposed in Section 2.5 is carried out using the following poorly informative prior distributions for the parameters of the calibration curve: π ( a ) N ( 1 , 10 ) , π ( b ) N ( 1 , 10 ) , and π ( c ) N ( 1 , 10 ) .
The posterior distributions of the parameters a, b, and c of the calibration curve are presented in Figure 13. The corresponding calibration curve is illustrated in Figure 14 with its 95 % coverage intervals.
The calibration curve provided in Figure 14 suggests that the SThM technique is a promising technique for the determination of traceable thermal conductivities lower than 10 W   m 1   K 1 in the best case. Due to the shape of the calibration curve, the sensitivity of the technique highly decreases for thermal conductivity higher than 10 W   m 1   K 1 , which means that uncertainties on measurement for high thermally conductive materials should be significant values that will be discussed in Section 3.4.

3.4. Predictions and Associated Uncertainty Using the Calibration Curve

In this section, we study predictions of the thermal conductivity and their associated uncertainties X * using the calibration curve for arbitrary intermediate measurand Y * . In order to cover the range [ 0.1 , 10 ] W   m 1   K 1 of thermal conductivity, the following three artificial y * measurements are chosen to illustrate the prediction using the Bayesian methodology: y 1 * = 0.7 , y 2 * = 1.11 , and y 3 * = 1.12 . In this case, the vectors of SThM measurements Y * and predictions X * of the thermal conductivity write Y * = ( Y 1 * , Y 2 * , Y 3 * ) and X * = ( X 1 * , X 2 * , X 3 * ) respectively (see Section 2.5.1).
The following poorly informative prior distributions are chosen for the predictions π ( X 1 * ) N ( 0.5 , 10 ) , π ( X 2 * ) N ( 5 , 100 ) , π ( X 3 * ) N ( 10 , 100 ) and non informative Jeffrey’s prior distributions are chosen for the X i : π ( X i ) 1 for i 1 , , N .
In this study, we consider two uncertainty levels for y * , namely u ( y * ) = 0.005 , which corresponds to the median of observed uncertainties obtained with the SThM on the bulk materials in this study, and u ( y * ) = 0.002 , which corresponds to the lowest observed uncertainty. Since all parameters are jointly updated (see Section 2.5.2), we show that the uncertainty level of y * used for prediction has an effect on the estimations of all the parameters, in particular on those of the calibration curve, which is a desirable feature of the Bayesian inference. Summaries of the posterior distributions of all parameters ( a , b , c , X 1 , , X N , X 1 * , X 2 * , X 3 * , Y 1 , , Y N , Y 1 * , Y 2 * , Y 3 * ) with N = 12 are displayed in Appendix A Table A1 and Table A2 for u ( y * ) = 0.005 and u ( y * ) = 0.002 , respectively.
The effect of a lower uncertainty u ( y * ) = 0.002 on the predictions from the calibration curve is displayed in Section 4.1. The posterior distributions of the predictions and their associated 95 % coverage interval obtained for y 1 * = 0.7 , y 2 * = 1.11 , and y 3 * = 1.12 are presented in Figure 15 for the two uncertainty levels u ( y * ) = 0.002 and u ( y * ) = 0.005 . It should be noted that the posterior distributions of the predictions are not symmetric PDFs, in particular for the highest conductivity. The predicted values are given with their associated 95 % coverage interval. The highest density value for the predicted value is not located at the centre of the coverage interval. Therefore, it is not correct to attribute a standard uncertainty (absolute or relative) to the predicted value. Only the coverage interval provides a rigorous assessment of the uncertainty.

4. Discussions

4.1. Sensitivity of the Measurement Method

The calibration curve provided in Figure 14 suggests that the SThM technique is a promising technique for the determination of thermal conductivities lower than 10 W   m 1   K 1 . In this section, we propose to illustrate this finding with three different values of the Y measurements. Table 7 provides the results of the prediction obtained from simulation techniques based on MCMC for the three chosen values of the Y intermediate measurand. Two levels of standard uncertainty associated with Y measurements have been tested: 0.002 and 0.005 .
As expected, the uncertainty associated with the thermal conductivity increases as the Y intermediate measurand approaches the top of the rising part of the calibration curve in Figure 14. It appears that, currently, for values of Y higher than 1.2 , which corresponds to a thermal conductivity of 10   W   m 1   K 1 , the model does not enable the prediction of a thermal conductivity value with a sufficiently low uncertainty. The second important result presented in Table 7 is that the control of the uncertainty associated with the intermediate measurand Y is crucial. The simulation performed for a Y value of 1.2 with an associated uncertainty of 0.005 instead of 0.002 shows that it becomes impossible to determine the thermal conductivity with a sufficiently low uncertainty.

4.2. Improvement of Measurement Precision

As discussed in Section 3.2, the main contribution to the standard uncertainty is the dispersion of values in readability and reproducibility conditions. This statistical spread between measurements could be explained by heterogeneities at the surface of the sample (the contribution of reproducibility error to the variance seems to decrease a little bit when the roughness of materials decreases, see Table 2 and Table 6). Nevertheless, given the small decrease in the contribution of the reproducibility (at best, one point), it seems that the spread is also impacted by changes in the thermal probe-sample contact between each measurement (thermal interface resistance, thermal contact area, etc.). Since we follow the “dark mode” protocol to avoid the overheating induced by the laser, we do not manage accurately the applied tip force from one measurement to the next. Variations in applied tip force between measurements could induce changes in thermal interface resistance and thermal contact area. The next step of this work is to improve the reproducibility of our measurements by improving the management of the applied tip force.

4.3. Application to Nanomaterials

Estimating the thermal conductivity of an unknown material from the calibration curve requires that the heat transfer conditions are the same during the calibration and measurement. As the calibration curve has been established with bulk materials, the heat flows are mainly diffusive, which means that thermal conductivity can be directly identified only for nanostructures or materials whose dimensions are consistent with the Fourier regime with local thermodynamic equilibrium [49]. In the probe-sample-environment system, at least three dimensions have to be checked: the smallest dimension of the sample d m i n (such as the diameter of the nanowire or the thickness of the membrane), the mean-free path Λ of the studied material, and the effective solid-solid contact radius b e f f between the probe and the sample.
For systems where at least one of the dimensions d m i n or b e f f is lower than the mean-free path Λ , heat transport is completely different from that experienced in macroscopic systems. In that case, identification of the thermal conductance requires the development of specific model describing all heat transfers from the probe to the sample (including ballistic effects). To build the model, it is necessary to have a complete knowledge of the geometry of the probe (depending on the probe type, probe generation, or wear of the probe) [49,50], information about the thermal and electrical properties of the materials of the probe, and a perfect knowledge of the contact between the tip and the sample (surface roughness, water meniscus, boundary resistance, …) [15]. Unfortunately, it is really challenging to quantify all quantities involved in the measurement of the thermal conductivity of material and in the heat transfer model [50]. Indeed, measuring the thermal or electrical properties of the material constituents of the probe is tricky due to their small sizes. In addition, some parameters of the model and heat transfer mechanisms at a lower scale than the mean-free path are still the subject of fundamental research.

5. Conclusions

This work provides the first complete uncertainty assessment of thermal conductivity measurements by the SThM technique based on a calibration curve established with bulk calibration materials. This study shows that following the proposed protocols, it is possible to perform quantitative and traceable thermal conductivity measurements for materials with low thermal conductivity (under 10 W   m 1   K 1 ) with the SThM technique under specific conditions. Traceability is ensured by using calibration materials whose thermal conductivity measurements are themselves traceable. As traceability is established at the macroscale on bulk calibration materials, the measurement conditions must be the same between calibration and measurements (i.e., diffusive heat transfer regime). By ensuring strictly steady-state environmental conditions, minimising the roughness of studied material (ideally less than R a < 1 nm ), and minimising the dispersion of measurements to limit uncertainty on the measured quantity value to 0.2%, it is possible to reach an uncertainty value of less than 10% for the identified thermal conductivity value. For thermal conductivity greater than 10 W   m 1   K 1 , the current uncertainty values are too high from a metrological point of view to justify traceable measurement using the protocol developed for this study.
Regarding nanostructured materials, traceable measurement could be performed using the calibration and measurement protocols described in this paper when the heat transfer regimes stay the same between calibration and measurement, that is, in the diffusive regime. For a probe-sample system with one dimension lower than the mean-free path of the studied material, a dedicated model has to be developed, as highlighted in ref. [50].

Author Contributions

Conceptualization, N.F. and A.A.; methodology, N.F. and S.D. (Sarah Douri); experimental measurements, N.F. and S.D. (Sarah Douri); software, A.A and S.D. (Séverine Demeyer); validation, N.F., B.H. and S.D. (Séverine Demeyer); data analysis, N.F., A.A. and S.D. (Séverine Demeyer); data curation, N.F.; writing—original draft preparation, A.A., B.H. and N.F.; writing—review and editing, N.F., S.D. (Séverine Demeyer) and S.D. (Sarah Douri); visualization, N.F. and S.D. (Séverine Demeyer); supervision, B.H.; funding acquisition, N.F. and B.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the French metrology and the 19ENG05 NanoWires project. The 19ENG05 NanoWires project has received funding from the EMPIR programme co-financed by the Participating States and from the European Union’s Horizon 2020 research and innovation programme.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data that support the findings of this study are available on Zenedo repository: https://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.8188536; URL (accessed on 29 July 2023); Fleurence, Nolwenn; Douri, Sarah; Demeyer, Séverine; Allard, Alexandre and Hay, Bruno. (2023). Measurements on reference materials for SThM measurand versus thermal conductivity calibration curve (Version v1) [Data set].

Acknowledgments

The authors would like to thank the nanometrology team of LNE, in particular Alexandra Delvallée and Lydia Chibane for SEM and AFM characterisation of the samples and the probes, colleagues from the thermophysical team in particular Olivier Beaumont and Guillaume Davée for the thermal conductivity characterisation of samples. The authors are also grateful to Séverine Gomès from CETHIL-CNRS (France) for technical discussions and sharing experience from QUANTIHEAT project.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
AFMAtomic Force Microscopy
a. u.arbitrary unit
DCDirect current
GUMGuide to the expression of Uncertainty in Measurement
MCMMonte Carlo Method
MCMCMarkov Chain Monte Carlo
PDFProbability Distribution Functions
PMMApoly(methyl methacrylate)
POM-Cpoly-oxymethylene in copolymer
SEMScanning Electron Microscopy
SIInternational System of Units (SI for Système International)
SThMScanning Thermal Microscopy
TCRTemperature Coefficient Ratio

Nomenclature

Measurement ResultMeasured Quantity ValueUncertaintyDescription
Y m y m u ( y m ) individual measurand
Y m , i y m , i u ( y m , i ) individual measurand indexed by environmental measurement conditions
Yy u ( y ) mean value of measurand

Appendix A. Bayesian Estimates

Table A1. Posterior point estimates of all the quantities updated in the Bayesian inference for u ( Y ) = 0.005 .
Table A1. Posterior point estimates of all the quantities updated in the Bayesian inference for u ( Y ) = 0.005 .
MeanSD 2.50 % 25 % 50 % 75 % 97.50 % n_effRhat
a0.752480.021230.713320.737910.751540.766310.79741188911.000
b0.294790.016950.261920.283370.294560.306050.32837185341.000
c0.391780.021890.345780.377490.392740.406880.43207188321.000
X 1 116.982442.93226111.14404115.03346116.97302118.97162122.66783182631.000
X 2 93.511052.3366489.0063691.9187493.4917395.0935598.13226185881.000
X 3 52.080161.3006249.5182351.2120452.0804852.9490954.64027190551.000
X 4 36.811470.9283934.9985436.1856736.8127737.4411938.65114191781.000
X 5 29.775730.7358628.3321129.2815229.7734630.2694131.23389191631.000
X 6 9.192030.227848.746169.039599.189989.345949.63890188221.000
X 7 1.950290.046711.858111.918921.950291.981902.04136195031.000
X 8 1.416000.030441.356051.395741.415941.436361.47612188411.000
X 9 1.275430.028371.219521.256371.275411.294511.33108187961.000
X 10 1.060690.021771.018421.045921.060731.075401.10332180521.000
X 11 0.347490.006700.334560.342970.347520.351990.36062193471.000
X 12 0.182270.004510.173490.179180.182230.185340.19105193841.000
X 1 * 0.203870.007600.188870.198720.203960.208960.21867190611.000
X 2 * 6.692741.316024.806595.792996.462077.339919.85410172461.000
X 3 * 11.098525.613466.449798.356219.8387412.0934822.99046158981.000
Y 1 1.142370.002201.138011.140901.142361.143861.14668177981.000
Y 2 1.141900.002201.137561.140431.141901.143381.14620178171.000
Y 3 1.140030.002161.135751.138601.140041.141481.14427177471.000
Y 4 1.138290.002121.134091.136871.138301.139711.14244178641.000
Y 5 1.136890.002101.132721.135481.136891.138301.14100178591.000
Y 6 1.120900.001911.117141.119621.120911.122171.12462174631.000
Y 7 1.045580.002751.040101.043751.045591.047441.05089191081.000
Y 8 1.014790.002651.009581.013001.014771.016581.01997190421.000
Y 9 1.003190.003050.997131.001141.003211.005291.00906186281.000
Y 10 0.980860.002090.976740.979490.980870.982260.98490184851.000
Y 11 0.799420.003890.791710.796850.799400.802040.80706185761.000
Y 12 0.679980.002890.674320.678040.679990.681910.68565191861.000
Y 1 * 0.700050.005030.690180.696660.700090.703440.70986167201.000
Y 2 * 1.111550.005171.101411.108071.111541.115031.12171175681.000
Y 3 * 1.122560.005401.112091.118941.122401.126121.13366190171.000
Table A2. Posterior point estimates of all the quantities updated in the Bayesian inference for u ( Y ) = 0.002 .
Table A2. Posterior point estimates of all the quantities updated in the Bayesian inference for u ( Y ) = 0.002 .
MeanSD 2.50 % 25 % 50 % 75 % 97.50 % n_effRhat
a0.752570.021130.713900.737860.751450.766320.79607186631
b0.295000.016980.262600.283460.294710.306410.32894187101
c0.391770.021810.346780.377570.392880.406930.43158186331
X 1 116.981642.94698111.28035114.99642116.94231118.95543122.83754192531
X 2 93.503392.3220088.9802491.9429693.5264695.0489598.07742183171
X 3 52.063131.3019949.4933951.1753152.0701952.9461154.60207189961
X 4 36.814810.9259334.9873436.1840136.8139637.4341938.62662185051
X 5 29.758530.7406328.2938729.2588829.7603830.2595831.20828186351
X 6 9.191740.228028.742449.037499.193309.345119.63806184991
X 7 1.950600.046981.858511.919211.950351.982122.04297192101
X 8 1.416300.030241.357081.396101.416101.436171.47689189041
X 9 1.275000.028601.219191.255681.275211.294221.33215180751
X 10 1.060510.021941.017631.045621.060451.075251.10400195291
X 11 0.347590.006690.334400.343050.347610.352080.36065181541
X 12 0.182300.004480.173640.179240.182300.185310.19120185161
X 1 * 0.203900.005600.192960.200140.203900.207670.21493184981
X 2 * 6.254230.517525.360135.894226.210486.567987.39395187651
X 3 * 9.077031.130827.284838.281828.933039.7230011.65243187401
Y 1 1.142450.002181.138221.140961.142451.143931.14673181741
Y 2 1.141970.002171.137751.140491.141981.143441.14623181721
Y 3 1.140100.002131.135971.138651.140111.141551.14428180831
Y 4 1.138360.002101.134281.136931.138361.139791.14247181601
Y 5 1.136960.002071.132921.135541.136971.138361.14103181961
Y 6 1.120960.001881.117301.119691.120971.122221.12463184631
Y 7 1.045600.002761.040171.043751.045611.047481.05101191871
Y 8 1.014800.002671.009601.013011.014801.016591.02002187711
Y 9 1.003130.003060.997071.001061.003171.005231.00903184241
Y 10 0.980810.002090.976750.979390.980810.982230.98490190621
Y 11 0.799380.003870.791790.796770.799370.802020.80693179551
Y 12 0.679910.002900.674250.677920.679890.681830.68559191991
Y 1 * 0.700040.002020.696060.698670.700020.701410.70400190151
Y 2 * 1.110270.002001.106331.108941.110261.111601.11422190241
Y 3 * 1.120340.002041.116361.118961.120331.121701.12432187501

References

  1. El Sachat, A.; Alzina, F.; Sotomayor Torres, C.M.; Chavez-Angel, E. Heat Transport Control and Thermal Characterization of Low-Dimensional Materials: A Review. Nanomaterials 2021, 11, 175. [Google Scholar] [CrossRef]
  2. Majumdar, A. Scanning Thermal Microscopy. Ann. Rev. Mater. Sci. 1999, 29, 505–585. [Google Scholar] [CrossRef]
  3. Zhang, Y.; Zhu, W.; Hui, F.; Lanza, M.; Borca-Tasciuc, T.; Muñoz Rojo, M. A Review on Principles and Applications of Scanning Thermal Microscopy (SThM). Adv. Funct. Mater. 2020, 30, 1900892. [Google Scholar] [CrossRef]
  4. Mills, G.; Zhou, H.; Midha, A.; Donaldson, L.; Weaver, J.M.R. Scanning thermal microscopy using batch fabricated thermocouple probes. Appl. Phys. Lett. 1998, 72, 2900–2902. [Google Scholar] [CrossRef]
  5. Tovee, P.; Pumarol, M.; Zeze, D.; Kjoller, K.; Kolosov, O. Nanoscale spatial resolution probes for scanning thermal microscopy of solid state materials. J. Appl. Phys. 2012, 112, 114317. [Google Scholar] [CrossRef]
  6. Nordal, P.-E.; Kanstad, S.O. Photothermal Radiometry. Phys. Scr. 1979, 20, 659. [Google Scholar] [CrossRef]
  7. Phan, T.; Dilhaire, S.; Quintard, V.; Claeys, W.; Batsale, J.-C. Thermoreflectance measurements of transient temperature upon integrated circuits: Application to thermal conductivity identification. Microelectron. J. 1998, 29, 181–190. [Google Scholar] [CrossRef]
  8. Ocariz, A.; Sanchez-Lavega, A.; Salazar, A.; Fournier, D.; Boccara, A.C. Photothermal characterization of vertical and slanted thermal barriers: A quantitative comparison of mirage, thermoreflectance, and infrared radiometry. J. Appl. Phys. 1996, 80, 2968–2982. [Google Scholar] [CrossRef]
  9. Bontempi, A.; Thiery, L.; Teyssieux, D.; Briand, D.; Vairac, P. Quantitative thermal microscopy using thermoelectric probe in passive mode. Rev. Sci. Instrum. 2013, 84, 103703. [Google Scholar] [CrossRef]
  10. Guen, E.; Klapetek, P.; Puttock, R.; Hay, B.; Allard, A.; Maxwell, T.; Chapuis, P.-O.; Renahy, D.; Davee, G.; Valtr, M.; et al. SThM-based local thermomechanical analysis: Measurement intercomparison and uncertainty analysis. Int. J. Therm. Sci. 2020, 156, 106502. [Google Scholar] [CrossRef]
  11. Doumouro, J.; Perros, E.; Dodu, A.; Rahbany, N.; Leprat, D.; Krachmalnicoff, V.; Carminati, R.; Poirier, W.; De Wilde, Y. Quantitative Measurement of the Thermal Contact Resistance between a Glass Microsphere and a Plate. Phys. Rev. Appl. 2021, 15, 014063. [Google Scholar] [CrossRef]
  12. Gomès, S.; Assy, A.; Chapuis, P.-O. Scanning Thermal Microscopy: A Review. Phys. Status Solidi A 2015, 212, 477–494. [Google Scholar] [CrossRef]
  13. Zhang, Y.; Zhu, W.; Borca-Tasciuc, T. Sensitivity and spatial resolution for thermal conductivity measurements using noncontact scanning thermal microscopy with thermoresistive probes under ambient conditions. Oxf. Open Mater. Sci. 2021, 1, itab011. [Google Scholar] [CrossRef]
  14. Metzke, C.; Frammelsberger, W.; Weber, J.; Kühnel, F.; Zhu, K.; Lanza, M.; Benstetter, G. On the Limits of Scanning Thermal Microscopy of Ultrathin Films. Materials 2020, 13, 518. [Google Scholar] [CrossRef] [PubMed]
  15. Bodzenta, J.; Kaźmierczak-Bałata, A. Scanning thermal microscopy and its applications for quantitative thermal measurements. J. Appl. Phys. 2022, 132, 140902. [Google Scholar] [CrossRef]
  16. Amor, A.B.; Djomani, D.; Fakhfakh, M.; Dilhaire, S.; Vincent, L.; Grauby, S. Si and Ge allotrope heterostructured nanowires: Experimental evaluation of the thermal conductivity reduction. Nanotechnology 2019, 30, 375704. [Google Scholar] [CrossRef] [PubMed]
  17. Grauby, S.; Ben Amor, A.; Hallais, G.; Vincent, L.; Dilhaire, S. Imaging Thermoelectric Properties at the Nanoscale. Nanomaterials 2021, 11, 1199. [Google Scholar] [CrossRef] [PubMed]
  18. Li, Y.; Zhang, T.; Zhang, Y.; Zhao, C.; Zheng, N.; Yu, W. A comprehensive experimental study regarding size dependence on thermal conductivity of graphene oxide nanosheet. Int. Commun. Heat Mass Transf. 2022, 130, 105764. [Google Scholar] [CrossRef]
  19. Kaźmierczak-Bałata, A.; Grządziel, L.; Guziewicz, M.; Venkatachalapathy, V.; Kuznetsov, A.; Krzywiecki, M. Correlations of thermal properties with grain structure, morphology, and defect balance in nanoscale polycrystalline ZnO films. Appl. Surf. Sci. 2021, 546, 149095. [Google Scholar] [CrossRef]
  20. Sun, W.; Hamaoui, G.; Micusik, M.; Evgin, T.; Vykydalova, A.; Omastova, M.; Gomés, S. Investigation of the thermal conductivity enhancement mechanism of polymer composites with carbon-based fillers by scanning thermal microscopy. AIP Adv. 2022, 12, 105303. [Google Scholar] [CrossRef]
  21. Guide to the Expression of Uncertainty in Measurement—Part 6: Developing and Using Measurement Models; JCGM GUM-6: 2020; Joint Committee for Guides in Metrology (JCGM): Sèvres, France, 2020.
  22. BIPM. International Vocabulary of Metrology—Basic and General Concepts and Associated Terms (VIM), 3rd ed.; Joint Committee for Guides in Metrology (JCGM): Sèvres, France, 2012. [Google Scholar]
  23. Buzin, I.; Kamasa, P.; Pyda, M.; Wunderlich, B. Application of a Wollaston wire probe for quantitative thermal analysis. Thermochim. Acta 2002, 381, 9–18. [Google Scholar] [CrossRef]
  24. Lefevre, S.; Volz, S.; Saulnier, J.-B.; Fuentes, C.; Trannoy, N. Thermal conductivity calibration for hot wire based dc scanning thermal microscopy. Rev. Sci. Instrum. 2003, 74, 2418–2423. [Google Scholar] [CrossRef]
  25. Gomès, S.; Newby, P.; Canut, B.; Termentzidis, K.; Marty, O.; Fréchette, L.; Chantrenne, P.; Aimez, V.; Bluet, J.M.; Lysenko, V. Characterization of the thermal conductivity of insulating thin films by scanning thermal microscopy. Microelectron. J. 2013, 44, 1029–1034. [Google Scholar] [CrossRef]
  26. Dominika Trefon-Radziejewska, D.; Juszczyk, J.; Fleming, A.; Horny, N.; Antoniow, J.S.; Chirtoc, M.; Kaźmierczak-Bałata, A.; Bodzenta, J. Thermal characterization of metal phthalocyanine layers using photothermal radiometry and scanning thermal microscopy methods. Synth. Met. 2017, 232, 72–78. [Google Scholar] [CrossRef]
  27. Guen, E. Microscopie Thermique à Sonde Locale: étalonnages—Protocoles de Mesure et Applications Quantitatives sur des Matériaux Nanostructurés. Doctoral Thesis, Université de Lyon, Lyon, France, January 2020. [Google Scholar]
  28. Fischer, H. Quantitative determination of heat conductivities by scanning thermal microscopy. Thermochim. Acta 2005, 425, 69–74. [Google Scholar] [CrossRef]
  29. GUM 1995 with Minor Corrections Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement; JCGM Guide 100:2008; Joint Committee for Guides in Metrology (JCGM): Sèvre, France, 2008.
  30. Ramiandrisoa, L.; Allard, A.; Hay, B.; Gomès, S. Uncertainty assessment for measurements performed in the determination of thermal conductivity by scanning thermal microscopy. Meas. Sci. Technol. 2017, 28, 115010. [Google Scholar] [CrossRef]
  31. Ramiandrisoa, L.; Allard, A.; Joumani, Y.; Hay, B.; Gomés, S. A dark mode in scanning thermal microscopy,» Review of Scientific Instruments. Rev. Sci. Instrum. 2017, 188, 1125115. [Google Scholar]
  32. Gomès, S. QUANTItative Scanning Probe Microscopy Techniques for HEAT Transfer Management in Nanomaterials. [Final Research Report] QUANTIHEAT project—n° FP7-CP-IP – 604668. 2017. Available online: https://cordis.europa.eu/docs/results/604/604668/final1-quantiheat-final-report-1.pdf (accessed on 5 July 2023).
  33. Assy, A.; Lefèvre, S.; Chapuis, P.-O.; Gomès, S. Analysis of heat transfer in the water meniscus at the tip-sample contact in scanning thermal microscopy. J. Phys. D Appl. Phys. 2014, 47, 442001. [Google Scholar] [CrossRef]
  34. Douri, S.; Gomès, S.; Hameury, J.; Feltin, N.; Hay, B.; Fleurence, N. Study of influencing parameters for thermal conductivity measurements at the nanoscale by SThM technique. In Proceedings of the 20th International Metrology Congress, CIM 2021, Lyon, France, 7–9 September 2021. [Google Scholar]
  35. Saci, A.; Battaglia, J.-L.; De, I. Notice of Violation of IEEE Publication Principles: Accurate New Methodology in Scanning Thermal Microscopy. IEEE Trans. Nanotechnol. 2015, 14, 1035–1039. [Google Scholar] [CrossRef]
  36. Wilson, A.A.; Rojo, M.M.; Abad, B.; Perez, J.A.; Maiz, J.; Schomacker, J.; Martín-Gonzalez, M.; Borca-Tasciuc, D.-A.; Borca-Tasciuc, T. Thermal conductivity measurements of high and low thermal conductivity films using a scanning hot probe method in the 3ω mode and novel calibration strategies. Nanoscale 2015, 7, 15404–15412. [Google Scholar] [CrossRef]
  37. Bodzenta, J.; Juszczyk, J.; Kaźmierczak-Bałata, A.; Firek, P.; Fleming, A.; Chirtoc, M. Quantitative Thermal Microscopy Measurement with Thermal Probe Driven by dc+ac Current. Int. J. Thermophys. 2016, 37, 73. [Google Scholar] [CrossRef]
  38. Assy, A.; Gomès, S. Heat transfer at nanoscale contacts investigated with scanning thermal microscopy. Appl. Phys. Lett. 2015, 107, 043105. [Google Scholar] [CrossRef]
  39. Hay, B.; Filtz, J.-R.; Hameury, J.; Rongione, L. Uncertainty of thermal diffusivity measurements by laser flash method. Int. J. Thermophys. 2005, 26, 1883–1898. [Google Scholar] [CrossRef]
  40. Hay, B.; Beaumont, O.; Failleau, G.; Fleurence, N.; Grelard, M.; Razouk, R.; Davée, G.; Hameury, J. Uncertainty Assessment for Very High Temperature Thermal Diffusivity Measurements on Molybdenum, Tungsten and Isotropic Graphite. Int. J. Thermophys. 2022, 43, 2. [Google Scholar] [CrossRef]
  41. Hay, B.; Hameury, J.; Filtz, J.-R.; Haloua, F.; Morice, R. The metrological platform of LNE for measuring thermophysical properties of materials. High Temp.-High Press. 2010, 39, 181–208. [Google Scholar]
  42. Evaluation of Measurement Data—Supplement 1 to the “Guide to the Expression of Uncertainty in Measurement”—Propagation of Distributions Using a Monte Carlo Method; JCGM Guide 98-32008; Joint Committee for Guides in Metrology (JCGM): Sèvres, France, 2008.
  43. Laboratoire National de Métrologie et d’Essais. Available online: https://www.lne.fr/fr/logiciels/lne-mcm (accessed on 5 April 2023).
  44. Koepke, A.; Lafarge, T.; Possolo, A.; Toman, B. NIST Consensus Builder User Manual. 2020. Available online: https://consensus.nist.gov/app/nicob (accessed on 20 August 2023).
  45. Gelman, A.; Carlin, J.B.; Stern, H.S.; Dunson, D.B.; Vehtari, A.; Rubin, D.B. Bayesian Data Analysis, 3rd ed.; Chapman and Hall/CRC: Boca Raton, FL, USA, 2013. [Google Scholar]
  46. Robert, C.P.; Casella, G. Monte Carlo Statistical Methods, 2nd ed.; Springer Science+Business Media: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  47. Stan Development Team. RStan: The R Interface to Stan. R Package Version 2.26.13. 2023. Available online: https://mc-stan.org/ (accessed on 12 May 2023).
  48. Roy, V. Convergence Diagnostics for Markov Chain Monte Carlo. Annu. Rev. Stat. Its Appl. 2020, 7, 387–412. [Google Scholar] [CrossRef]
  49. Gomès, S. Microscopie thermique à balayage (SThM). Tech. L’ingén. 2023, R2770 V2:1–R2770 V2:30. [Google Scholar] [CrossRef]
  50. Douri, S.; Fleurence, N.; Hameury, J.; Gomès, S. Micro and Nanoscale Heat Transfer Investigation by 3D FEM for Scanning Thermal Microscopy. C’Nano 2023, Poitiers, France, 17 March 2023, Thematic Session: Nanoscale Heat Transfer—Measurement. Available online: https://cnano2023.sciencesconf.org/data/Book_of_Abstracts_NanoHeat.pdf (accessed on 27 July 2023).
Figure 1. Representation of the workflow for establishing traceability to the SI for thermal conductivity measurements by SThM technique for measurements in diffusive thermal regime.
Figure 1. Representation of the workflow for establishing traceability to the SI for thermal conductivity measurements by SThM technique for measurements in diffusive thermal regime.
Nanomaterials 13 02424 g001
Figure 2. SEM images of a KNT probe (2an type): (a) View of the SiN cantilever with the two gold pads (b) Zoom on the tip with the resistive Pt ribbon at deposited at the top end of the tip.
Figure 2. SEM images of a KNT probe (2an type): (a) View of the SiN cantilever with the two gold pads (b) Zoom on the tip with the resistive Pt ribbon at deposited at the top end of the tip.
Nanomaterials 13 02424 g002
Figure 3. Scheme of the thermal unit encloses an adjustable current generator and a Wheatstone bridge composed of two fixed resistances ( R 1 and R 2 ), two adjustable resistances: R f for coarse adjustments and R a for fine adjustments, the resistance temperature probe R and an amplifier setup A to amplify the bridge balance B B v voltage. The probe electrical resistance R is included in one leg of the Wheatstone bridge and the adjustable resistance R a in the opposite leg. R a can be set manually by rotary knob arbitrarily scalable in 1000 graduations. Value of the knob adjustment is denoted B B k .
Figure 3. Scheme of the thermal unit encloses an adjustable current generator and a Wheatstone bridge composed of two fixed resistances ( R 1 and R 2 ), two adjustable resistances: R f for coarse adjustments and R a for fine adjustments, the resistance temperature probe R and an amplifier setup A to amplify the bridge balance B B v voltage. The probe electrical resistance R is included in one leg of the Wheatstone bridge and the adjustable resistance R a in the opposite leg. R a can be set manually by rotary knob arbitrarily scalable in 1000 graduations. Value of the knob adjustment is denoted B B k .
Nanomaterials 13 02424 g003
Figure 4. Views of apparatus: (a) SThM and thermal unit inside the dedicated enclosure to decrease influence of thermal drifts and (b) view of the location of the probe above the S i O 2 f reference sample. Range of displacement of x, y direction piezoelectric scanners enables to switch from S i O 2 f reference sample to the studied sample ( Z r O 2 sample for the picture).
Figure 4. Views of apparatus: (a) SThM and thermal unit inside the dedicated enclosure to decrease influence of thermal drifts and (b) view of the location of the probe above the S i O 2 f reference sample. Range of displacement of x, y direction piezoelectric scanners enables to switch from S i O 2 f reference sample to the studied sample ( Z r O 2 sample for the picture).
Nanomaterials 13 02424 g004
Figure 5. Work flow of measurement sequence of U probe (cyan solid line) and B B v (orange dotted line) signals: five successive landings on the S i O 2 f reference sample following by five successive landings on the studied sample (alumina sample in this example). For each sample (reference and studied) the first three landings are performed on the same location and the last two on different locations. Measurement values used for R computation are identified in red for the S i O 2 f reference sample and in blue for the studied sample. The residual drift is highlighted with the red dashed line for U “out of contact” measurements.
Figure 5. Work flow of measurement sequence of U probe (cyan solid line) and B B v (orange dotted line) signals: five successive landings on the S i O 2 f reference sample following by five successive landings on the studied sample (alumina sample in this example). For each sample (reference and studied) the first three landings are performed on the same location and the last two on different locations. Measurement values used for R computation are identified in red for the S i O 2 f reference sample and in blue for the studied sample. The residual drift is highlighted with the red dashed line for U “out of contact” measurements.
Nanomaterials 13 02424 g005
Figure 6. Main steps for the evaluation of measurement uncertainty: In a first step: Modelling the measurement process for Y m from Section 2.4.1, Evaluating input quantities with associated Probability Distribution Functions for Y m in Section 2.4.2, Propagating with Monte Carlo method (MCM) in Section 2.4.3 to report single quantity measured values y m with its associated uncertainty u ( y m ) . In a second step: Use y m , i input quantities measured under various conditions indexed by i (i from 1 to 10) with measurement model described in Section 2.2.2 then combine the ten y m , i measurements and the associated uncertainties in Section 2.4.4 to report the y ¯ mean value with standard uncertainty u ( y ¯ ) .
Figure 6. Main steps for the evaluation of measurement uncertainty: In a first step: Modelling the measurement process for Y m from Section 2.4.1, Evaluating input quantities with associated Probability Distribution Functions for Y m in Section 2.4.2, Propagating with Monte Carlo method (MCM) in Section 2.4.3 to report single quantity measured values y m with its associated uncertainty u ( y m ) . In a second step: Use y m , i input quantities measured under various conditions indexed by i (i from 1 to 10) with measurement model described in Section 2.2.2 then combine the ten y m , i measurements and the associated uncertainties in Section 2.4.4 to report the y ¯ mean value with standard uncertainty u ( y ¯ ) .
Nanomaterials 13 02424 g006
Figure 7. Analyses of propagation of distributions on y m , 3 measured quantity value: (a) Probability density distribution (PDF) for the y m , 3 value (b) Spearman’s rank correlation coefficients.
Figure 7. Analyses of propagation of distributions on y m , 3 measured quantity value: (a) Probability density distribution (PDF) for the y m , 3 value (b) Spearman’s rank correlation coefficients.
Nanomaterials 13 02424 g007
Figure 8. Graphic representation of the standard deviation computed from a set of three measurements performed in repeatability condition for two types of measurement (landing and withdrawal) and for all the twelve calibration samples).
Figure 8. Graphic representation of the standard deviation computed from a set of three measurements performed in repeatability condition for two types of measurement (landing and withdrawal) and for all the twelve calibration samples).
Nanomaterials 13 02424 g008
Figure 9. Graphic representation of the three quantity measured values of y m , i , l a n d i n g ( i = 1 ; 3 and 5), measured at the same location, computed in landing configuration and the three quantity measured values of y m , i , w i t h d r a w a l ( i = 2 ; 4 and 6) computed in withdrawal configuration. Each data is indicated with its associated absolute uncertainty (coefficient k = 1 ) represented by black error bars. The blue circles correspond to landing measurement points and the orange square to withdrawal measurement points. The blue solid line represents the mean value for the landing configuration with its associated standard deviation (blue dashed lines) and the orange dash-dotted line represents the mean value for the withdrawal configuration with its associated standard deviation (orange dotted lines).
Figure 9. Graphic representation of the three quantity measured values of y m , i , l a n d i n g ( i = 1 ; 3 and 5), measured at the same location, computed in landing configuration and the three quantity measured values of y m , i , w i t h d r a w a l ( i = 2 ; 4 and 6) computed in withdrawal configuration. Each data is indicated with its associated absolute uncertainty (coefficient k = 1 ) represented by black error bars. The blue circles correspond to landing measurement points and the orange square to withdrawal measurement points. The blue solid line represents the mean value for the landing configuration with its associated standard deviation (blue dashed lines) and the orange dash-dotted line represents the mean value for the withdrawal configuration with its associated standard deviation (orange dotted lines).
Nanomaterials 13 02424 g009
Figure 10. Comparison of mean values and the standard deviation (black error bars) obtained in landing condition (blue rectangles) and withdrawal condition (orange rectangles) for all calibration samples.
Figure 10. Comparison of mean values and the standard deviation (black error bars) obtained in landing condition (blue rectangles) and withdrawal condition (orange rectangles) for all calibration samples.
Nanomaterials 13 02424 g010
Figure 11. Comparison of standard deviation obtained in repeatability conditions (same location) and reproducibility conditions (different locations).
Figure 11. Comparison of standard deviation obtained in repeatability conditions (same location) and reproducibility conditions (different locations).
Nanomaterials 13 02424 g011
Figure 12. Ten measured quantity values for PMMA sample.
Figure 12. Ten measured quantity values for PMMA sample.
Nanomaterials 13 02424 g012
Figure 13. Posterior distributions of the parameters a, b, c of the calibration curve.
Figure 13. Posterior distributions of the parameters a, b, c of the calibration curve.
Nanomaterials 13 02424 g013
Figure 14. Calibration curve (in black) obtained from the experimental data from Table 6 analysed with the Bayesian approach. Points are represented with their associated expanded uncertainty ( k = 2 ) for both axes. Red dashed lines represent the 95 % coverage intervals associated with the (estimated) calibration curve for each conductivity.
Figure 14. Calibration curve (in black) obtained from the experimental data from Table 6 analysed with the Bayesian approach. Points are represented with their associated expanded uncertainty ( k = 2 ) for both axes. Red dashed lines represent the 95 % coverage intervals associated with the (estimated) calibration curve for each conductivity.
Nanomaterials 13 02424 g014
Figure 15. Posterior distributions of the predictions and their associated 95 % coverage interval obtained for (a) Y = y 1 * , (b) Y = y 2 * and (c) Y = y 3 * for the two uncertainty levels u ( y * ) = 0.002 and u ( y * ) = 0.005 .
Figure 15. Posterior distributions of the predictions and their associated 95 % coverage interval obtained for (a) Y = y 1 * , (b) Y = y 2 * and (c) Y = y 3 * for the two uncertainty levels u ( y * ) = 0.002 and u ( y * ) = 0.005 .
Nanomaterials 13 02424 g015aNanomaterials 13 02424 g015b
Table 1. Identification of the ten measured values y m , i from various conditions of measurement indexed from i = 1 to i = 10 relative to the chronological acquisition, odd i indexes relative to landing measurement condition, even i indexes relative to withdrawal measurement condition.
Table 1. Identification of the ten measured values y m , i from various conditions of measurement indexed from i = 1 to i = 10 relative to the chronological acquisition, odd i indexes relative to landing measurement condition, even i indexes relative to withdrawal measurement condition.
LocationsLanding ConditionWithdrawal Condition
location n°1 y m , 1 y m , 2
y m , 3 y m , 4
y m , 5 y m , 6
location n°2 y m , 7 y m , 8
location n°3 y m , 9 y m , 10
Table 2. Thermal conductivity k (relative expanded uncertainty estimated to 5 % ) and roughness R a of calibrated materials measured at 23 °C. The sample thickness is identified as: 1 *, 2 ** and 5 *** mm.
Table 2. Thermal conductivity k (relative expanded uncertainty estimated to 5 % ) and roughness R a of calibrated materials measured at 23 °C. The sample thickness is identified as: 1 *, 2 ** and 5 *** mm.
SampleStructureProviderk [Wm 1 K 1 ] Ra [nm]
P M M A ***PolymerGoodfellow0.1875.04
P O M C ***PolymerRadiospare0.32911.7
B o r o s i l i c a t e g l a s s **AmorphousNeyco1.11< 0.5
S i O 2 f **AmorphousNeyco1.280.56
S i O 2 N E G S 1 **AmorphousNeyco1.40<1
Z r O 2 **Single crystalNeyco1.95< 0.5
T i O 2 **Single crystalNeyco9.15< 0.5
A l 2 O 3 p **Poly crystalNeyco29.87.52
S a p p h i r e *Single crystalCrystal GmbH36.9< 0.5
G e r m a n i u m **Single crystalCrystal GmbH52.0< 0.5
S i p + + **SemiconductorGoodfellow93.40.75
Z i n c **MetalNeyco1178.14
Table 4. Summary of the ten values obtain on the PMMA sample: the measured quantity value y m , i , the standard uncertainty (absolute and relative) and the 95% coverage interval are given.
Table 4. Summary of the ten values obtain on the PMMA sample: the measured quantity value y m , i , the standard uncertainty (absolute and relative) and the 95% coverage interval are given.
Identification
Measurement
y m , i Value
[a. u.]
Standard Uncertainty95% Coverage Interval
Abs Rel. (%) [a. u.] [a. u.]
y m , 1 0.6966 0.0047 0.68 0.6873 0.7058
y m , 2 0.6855 0.0033 0.48 0.6789 0.6919
y m , 3 0.6737 0.0028 0.42 0.6681 0.6792
y m , 4 0.6719 0.0025 0.38 0.6668 0.6768
y m , 5 0.6679 0.0026 0.39 0.6628 0.6731
y m , 6 0.6713 0.0026 0.38 0.6663 0.6764
y m , 7 0.6883 0.0032 0.46 0.6821 0.6944
y m , 8 0.6780 0.0034 0.50 0.6714 0.6848
y m , 9 0.6737 0.0033 0.49 0.6673 0.6801
y m , 10 0.6778 0.0030 0.44 0.6719 0.6836
Table 5. Comparison of instrumental standard uncertainty value associated to each measured quantity value to the standard deviation value computed for set of replicated measurements for each materials. The highest instrumental standard uncertainties (absolute and relative) and the computed standard deviation (absolute and relative) are given.
Table 5. Comparison of instrumental standard uncertainty value associated to each measured quantity value to the standard deviation value computed for set of replicated measurements for each materials. The highest instrumental standard uncertainties (absolute and relative) and the computed standard deviation (absolute and relative) are given.
Material
Sample
Measurement
Condition
Max. Standard UncertaintyStandard Deviation
Abs Rel. (%) Abs Rel. (%)
P M M A landing 0.0047 0.68 0.0152 2.23
withdrawal 0.0034 0.50 0.0080 1.19
P O M C landing 0.0031 0.38 0.0080 0.98
withdrawal 0.0034 0.41 0.0083 1.01
B o r o s i l i c a t e g l a s s landing 0.0045 0.46 0.0083 0.85
withdrawal 0.0045 0.47 0.0068 0.70
S i O 2 f landing 0.0049 0.49 0.0075 0.76
withdrawal 0.0043 0.43 0.0078 0.77
S i O 2 N E G S 1 landing 0.0051 0.50 0.0105 1.02
withdrawal 0.0044 0.44 0.0074 0.70
Z r O 2 landing 0.0052 0.49 0.074 0.70
withdrawal 0.0047 0.46 0.0074 0.72
T i O 2 landing 0.0056 0.49 0.0101 0.88
withdrawal 0.0057 0.50 0.0100 0.88
A l u m i n a landing 0.0051 0.45 0.0338 3.00
withdrawal 0.0050 0.45 0.024 2.14
S a p p h i r e landing 0.0038 0.34 0.0061 0.54
withdrawal 0.0039 0.35 0.0096 0.86
G e r m a n i u m landing 0.0045 0.39 0.0118 1.03
withdrawal 0.0051 0.45 0.0082 0.72
S i p + + landing 0.0042 0.36 0.0082 0.71
withdrawal 0.0041 0.35 0.0094 0.81
Z i n c landing 0.0055 0.48 0.0181 1.59
withdrawal 0.0055 0.49 0.0264 2.34
Table 6. Summary of the experimental measurements on calibration samples: the intermediate measurand mean value Y, the standard uncertainty u ( Y ) (absolute and relative).
Table 6. Summary of the experimental measurements on calibration samples: the intermediate measurand mean value Y, the standard uncertainty u ( Y ) (absolute and relative).
SampleThermal ConductivityY Intermediate Measurand u ( Y ) Standard Uncertainty
(Wm 1 K 1 ) Mean Value (a.u.)AbsRel.(%)
P M M A 0.187 0.6780 0.0029 0.4
P O M C 0.329 0.8145 0.0055 0.7
B o r o s i l i c a t e g l a s s 1.11 0.9780 0.0023 0.2
S i O 2 f 1.28 1.0019 0.0048 0.5
S i O 2 N E G S 1 1.40 1.0173 0.0038 0.4
Z r O 2 1.95 1.0457 0.0072 0.7
T i O 2 9.15 1.1316 0.0057 0.5
A l u m i n a 29.8 1.1140 0.0091 0.8
S a p p h i r e 36.9 1.1241 0.0045 0.4
G e r m a n i u m 52.0 1.1460 0.0035 0.3
S i p + + 93.4 1.1548 0.0044 0.4
Z i n c 117 1.1158 0.0111 1.0
Table 7. Prediction of the thermal conductivity for different values of the direct measurement. The indications of the standard uncertainty (absolute and relative) are given only for the order of magnitude, only the coverage interval give the rigorous estimation of the uncertainty level as discussed in previous section (Section 3.4).
Table 7. Prediction of the thermal conductivity for different values of the direct measurement. The indications of the standard uncertainty (absolute and relative) are given only for the order of magnitude, only the coverage interval give the rigorous estimation of the uncertainty level as discussed in previous section (Section 3.4).
y 0 u y 0 k 0
Wm 1 K 1
u k 0
Wm 1 K 1
u k 0
(%)
95 % Coverage Interval
Wm 1 K 1
0.7 0.005 0.20407 0.00753 3.7 0.18923 ; 0.21897
0.7 0.002 0.20390 0.00560 2.7 0.19296 ; 0.21493
1.1 0.005 6.69603 1.29795 19.4 4.83485 ; 9.77990
1.1 0.002 6.25423 0.51752 8.3 5.36013 ; 7.39395
1.2 0.005 11.01500 5.21046 47.3 6.48843 ; 22.57146
1.2 0.002 9.07703 1.13082 12.5 7.28483 ; 11.65243
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fleurence, N.; Demeyer, S.; Allard, A.; Douri, S.; Hay, B. Quantitative Measurement of Thermal Conductivity by SThM Technique: Measurements, Calibration Protocols and Uncertainty Evaluation. Nanomaterials 2023, 13, 2424. https://0-doi-org.brum.beds.ac.uk/10.3390/nano13172424

AMA Style

Fleurence N, Demeyer S, Allard A, Douri S, Hay B. Quantitative Measurement of Thermal Conductivity by SThM Technique: Measurements, Calibration Protocols and Uncertainty Evaluation. Nanomaterials. 2023; 13(17):2424. https://0-doi-org.brum.beds.ac.uk/10.3390/nano13172424

Chicago/Turabian Style

Fleurence, Nolwenn, Séverine Demeyer, Alexandre Allard, Sarah Douri, and Bruno Hay. 2023. "Quantitative Measurement of Thermal Conductivity by SThM Technique: Measurements, Calibration Protocols and Uncertainty Evaluation" Nanomaterials 13, no. 17: 2424. https://0-doi-org.brum.beds.ac.uk/10.3390/nano13172424

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop