Nadine
Gächter
*a,
Fabian
Könemann
a,
Masiar
Sistani
b,
Maximilian G.
Bartmann
b,
Marilyne
Sousa
a,
Philipp
Staudinger
a,
Alois
Lugstein
b and
Bernd
Gotsmann
*a
aIBM Research - Zurich, 8803 Rüschlikon, Switzerland. E-mail: nadine.gaechter@unibas.ch; bgo@zurich.ibm.com
bInstitute of Solid State Electronics - TU Wien, 1040 Vienna, Austria
First published on 8th October 2020
The thermoelectric properties of a nanoscale germanium segment connected by aluminium nanowires are studied using scanning thermal microscopy. The germanium segment of 168 nm length features atomically sharp interfaces to the aluminium wires and is surrounded by an Al2O3 shell. The temperature distribution along the self-heated nanowire is measured as a function of the applied electrical current, for both Joule and Peltier effects. An analysis is developed that is able to extract the thermal and thermoelectric properties including thermal conductivity, the thermal boundary resistance to the substrate and the Peltier coefficient from a single measurement. Our investigations demonstrate the potential of quantitative measurements of temperature around self-heated devices and structures down to the scattering length of heat carriers.
The research in this area is marked by the difficulty in performing nanoscopic heat transport measurements.1 Thermal contact resistances dominate at small length scales translate into systematic challenges for measuring the local temperature or heat flux in a structure. The thermal contacts to nanowire samples, for example, have turned out to be a major issue.5,8 It has been proposed that this can be mitigated using extended measurement series for example of thermal transport as a function of length, sometimes called transmission line method9 or related.10 Apart from being time-intensive, these solutions rely on other assumptions such as the ability to fabricate reproducible test structures with the same contact resistances to electrodes or thermometers. For semiconducting nanowires this is, oftentimes, a major difficulty.11
Scanning thermal microscopy (SThM) has been used to measure thermal properties of thermoelectric materials.12 Motivated by the recent success of using SThM,13 we here describe the development of a method to systematically extract thermal and thermoelectric transport properties from the spatial temperature information. In contrast to transport measurements using external heaters and thermometers, here we analyse the temperature distribution in self-heated nanostructures to extract information through fitting with appropriate models. While thermal conductivity measurements with SThM are often hampered by the thermal resistance between probing tip and sample, we apply periodic self-heating to circumvent this problem and are able to measure the sample temperatures.8 In this way, we can not only mitigate but also quantify the effect of thermal contact resistances inside and around the device. Thereby, we can address length scales not readily accessible using other methods (such as a segment and its interfaces within a nanowire). Additionally, the temperature distribution is measured in its topographic context by atomic force microscopy (AFM). It is by that possible to align and compare scans under different measurement conditions. Also, the method is independent on the geometry and material of the device. This means, that there is no need for altering the sample architecture for the thermal measurements.
There are examples of interpreting the temperature field of nanowire structures using SThM to extract thermal conductivity.8,14–16 Here, we extend the method to extract the thermal conductivity of a thermoelectric material and the metal leads, the thermal interfaces to a substrate, and the Peltier coefficient from the observation of the temperature fields. The nanowire device consists of a germanium (Ge) segment which is monolithically integrated in a single crystalline aluminium (c-Al) nanowire with atomically sharp interfaces17 and a diameter of 35 nm. A dark field scanning transmission electron microscope (DF-STEM) image of the nanowire cross-section, a SEM image and a sketch of the device are shown in Fig. 1(b–d). Further information on the fabrication, electrical and structural characterization of the device can be found in the ESI† and in ref. 18–20. The thermoelectric figure of merit can thus be calculated from a single measurement collecting the complementary information from the Joule and the Peltier signals. The analysis concentrates on the temperature profile in the lateral direction of the wire for the different operating currents. The Ge segment length of 168 nm is interesting to study as it is long compared to typical thermionic emission barriers and shorter than thermoelectric materials typically studied. Our study reveals thermalization effects of carriers in the metal electrodes, which is the typical picture used to explain thermionic emission, while less often considered in the context of conventional thermoelectric effects.
A major difficulty in the determination of the sample temperature Tsample lies in quantifying the thermal resistance of the tip, the tip–sample contact and the spreading into the substrate. They are all summarized in Rts, which then depends on the surface material, the precise shape of the tip and largely on the size of the touching point that changes with the surface granularity and the topography of the sample. To overcome this problem, the nanowire is operated with an electrical AC-bias of frequency fmod = 1234 Hz that creates a continuously modulated temperature field. Rts is then eliminated in the equations by simultaneous measurement of the time-averaged sensor signal ΔVDC and the de-modulated sensor signal amplitude ΔVAC. An illustration of the set-up is depicted in Fig. 1a. The sample temperature difference to ambient temperature, ΔTsample, for linear devices is then given by8
(1) |
An additional feature of the lock-in detection is the possibility to separate the temperature change that is caused by Joule heating and the Peltier effect, respectively. Joule heating, on the one hand, is proportional to the dissipated electrical power and the signal appears at 2fmod. On the other hand, the Peltier effect is directly proportional to the current and produces a thermal signal at 1fmod. Therefore, a detection of SThM-signals at different frequencies allows the separation of Joule and Peltier effects for linear resistors (see below and eqn (4) and (5) for details).
Fig. 2 Scanning thermal microscopy data of the wire operating with an AC-current of amplitude I = 26 μA. (a) AFM topography. (b) Thermal signal amplitude taken at 2 fmod in response to Joule heating. (c) Phase of the 2 fmod signal. (d) Joule thermal map calculated using eqn (1) and the scans in (b) and (e). (e) DC-thermal measurement proportional to the thermal resistance of the tip–sample contact. (f) Thermal signal amplitude taken at 1 fmod in response to Peltier heating/cooling. (g) Phase of the 1 fmod AC signal, where ±5 V equals ±90°. (h) Peltier thermal map calculated using eqn (1) and the scans in (e), (f) and (g). |
The following analysis concentrates on the line-profiles extracted from the two-dimensional temperature maps of Fig. 2(d) and (h) and according ones for other currents between 13 and 26 μA. The temperature profiles along the long axis of the nanowire are shown in Fig. 3 for both Joule and Peltier signals.
(2) |
Firstly, the radial temperature distribution can be neglected. This assumption works well for wires with high aspect ratio as in this case. For this composite system we further require that the thermal interface resistance between the core and Al2O3-shell is smaller than that between wire and substrate. While this is a plausible assumption (see below), our study is limited in studying lateral intra-wire effects.
Secondly, the heat spreading into the substrate from a line heat source has no analytical solution. To be able to simplify this into a single thermal conductance term per unit length between a nanowire and the substrate (gΔT) requires, that the thermal interface resistance to the substrate is larger than the spreading resistances within the wire and the substrate. As a first indication, we observe a discontinuity in the horizontal temperature profile between wire and substrate. Fig. 4a shows that the temperature rise along a line section perpendicular to the wire is approximately constant across the nanowire surface and also five to ten times larger than on the substrate. Note, that the wire appears to have a width of 150 nm (instead of the actual 65 nm) due to the well-known convolution effects with the tip shape. In the convolution area and due to the finite diameter of the tip–surface the temperature distribution appears smeared out over some tens of nanometers in this line scan, in contrast to a higher resolution along the nanowire axis.
Thirdly, eqn (2) implies that the transport is essentially diffusive. This is a well established assumption for thermal transport in nanowires of these dimensions with one notable exception where reduced dimensions lead to a ballistic regime.9 However, in this study the observed temperature profiles within the segments are typical for diffusive transport and cannot be explained using ballistic effects. They can be observed only around the interfaces and boundaries, where the spatial resolution of the scan is higher than the carrier scattering length. Near the Ge–Al interfaces, the thermalization lengths of non-equilibrium charge carriers lead to a spatial distribution of the thermoelectric effect.7,15,23 The effective Peltier coefficient is then the integration over the spatial distribution of the local Peltier coefficient
(3) |
Finally, the AC-modulation period of the driving current is orders of magnitudes shorter than any thermal time constant of the system. Therefore, the scan is taken at steady state and the last term of eqn (2) is negligible.
For a system whose electric response is sufficiently linear, the total temperature distribution is disentangled to T(x,t) = Tambient + ΔTPeltier(x)sin(2πft) + ΔTJoule(x)sin(4πfmodt) for the current I(t) = I0 sin(2πfmodt). The respective thermal profiles are, then, described independently by
(4) |
(5) |
An initial qualitative observation of the line profiles, shown in Fig. 3, already allows for important conclusions. The Ge and Al sections of the wire exhibit very different characteristics. In the Ge segment, the temperature profile is approximately parabolic, as expected for a diffusive 1D system with uniform Joule dissipation. A heat source term in eqn (2) generally leads to a concave temperature profile. In contrast, the profile in the Al sections is convex, indicative of a situation which is dominated by heat dissipation and loss to the substrate. Here, the approximately exponential decay is expected for a uniform wire without any local Joule dissipation. The Al–Ge interface marks the Joule profiles by a sudden change of slope. When extracting the points of maximum temperature gradient for the Joule signals, their relative distances coincides for all measurements exactly with the length of the Ge segment. Also, after the deflection points there is a transition length of about 40–50 nm in the Al-segments. This area corresponds to the phonon–electron thermalization length that is blurred by a combination of spatial resolution and parallel heat transport in the oxide shell (the figure is in the ESI†). In the Peltier profile we observe the maximum and the minimum of the temperatures some tens of nanometers next to the Al–Ge interface due the spatial distribution of the thermionic effect at the order of the carrier scattering length.
(6) |
The indices j and p denote Joule and Peltier terms, respectively. As the temperature profile reaches room temperature before reaching the electrodes, the analytic solution to this differential equation is given by
(7) |
The analysis is pushed further by deriving κAl and g independently. It is observed, that the wire temperature of the Joule profiles decays to room temperature before reaching the electrodes. In other words, the entire heat that is generated over electric power dissipation in the wire (Pdis = I2ρLGe/Acore) is equal to the heat lost to the substrate. By making use of a finite element integral over all the measuring points on the nanowire, the substrate coupling is then the only unknown quantity. The conservation of heat requires
(8) |
We note that one reaches the same result by considering continuity of heat flux at the Ge–Al interface in addition to the 1D heat diffusion equation. The heat transported into the Al segment is given by Al = T(x0)·κAlAwire. The temperature gradient is then calculated from the exponential fit and extrapolated to the Al–Ge interface. Energy conservation implies:
Pdis = 2Al + g·ΔTavgGe·LGe, | (9) |
Next, we turn to the thermal conductivity of the Al segment. Using the extracted values of g and the ratio κ/g, we arrive at a thermal conductivity for the Al segments. To interpret, we need to consider that the value from the fit is an average over the cross section of both the Al core and the Al2O3 shell which are weighed by the respective thermal conductivities and obtain for the core κAl ≈ 150 W (m K)−1. We use the measured cross-sectional areas from the STEM image and the tabulated values for Al2O3. To compare with expected values we apply the Wiedemann Franz law, which results in an electrical resistivity of ρAl ≈ 50 × 10−9 Ω m about twice as high than theoretical predictions for bulk c-Al25,26 which is expected due to increased surface scattering in nanostructures. Indeed, this value is three times lower than previously determined for pure c-Al wires.19 However, in an SThM scan of such an operating c-Al wire a dominant local heat source is identified in the center (see the temperature profiles in the ESI†). This may be expected according to the wire fabrication and indicates a grain boundary. Therefore, the expected values for the electrical conductivity of the c-Al segments without any grain boundary are significantly higher and in good agreement with our experimental value. This result is shows that the influence of grain boundary scattering on the electrical resistivity exceeds the impact of surface scattering in this device. This is typical for polycrystalline nanowires of these dimensions.27 Finally, we would like to mention the importance of taking into account the contribution of an oxide shell to the heat transport in a nanowire.
(10) |
Here, κGe is an average value over core and shell in the Ge segment. This second order differential equation is analytically solved by
(11) |
It is interesting to note, that in this case of relatively small values of m, the temperature profile looks parabolic. For g = 0, supposing that the temperature profile is dominated by Joule heating and heat evacuation towards the electrodes, the solution to eqn (10) would indeed be a parabola. Such a fit seems to match the data well. However, the extracted value is given by κGe = 23.7 ± 3.3 W m−1 K−1 higher by 40% and we conclude that the hyperbolic model should not be simplified.
To interpret the result, we consider that the extracted value for κGe comprises the contributions of both oxide shell and Ge core. The area of the shell is double the area of the core, however, the exact values are measured on the STEM image. Both Al2O3 and Ge crystals have a larger conductivity in bulk, of about 30 and 50 to 60 W m−1 K−1, respectively.28,29 Both values, however, are expected to be significantly reduced, due to the boundary scattering of the phonons at reduced dimensions30 and the measured value is consistent with expectations.
The extraction of the Peltier coefficient is based on the conservation of heat flux at steady state. In the Al wires and near the Al–Ge interfaces, there is either a heat sink or a heat source caused by the Peltier effect GeAl = Π·I or AlGe = −Π·I, respectively, where AlGe > 0 and GeAl < 0, and Π is the Peltier coefficient. An illustration is seen in Fig. 5b. The injected heat from these regions is then equal to the sum of the heat flux towards the Al-electrodes and the heat flux over the Ge-segment to the other heat source/sink. (We can again neglect on the short length scale λ, as discussed above.) By use of Fourier's law the relation is then written as
= (−κAl∇TAl − κGe∇TGe)·Awire. | (12) |
In order to take into account the above described microscopic extension of the Peltier coefficient we consider the temperature profile beyond the length scale away from the Al–Ge interface. Then we can treat the temperature profiles as solutions to eqn (5) in sections without source terms. The temperature gradient in the Al-segment is the derivative with respect to the position x of the exponential fit to the Peltier profile evaluated at the junction as , see eqn (7). By extrapolating the exponential fit we conceptually concentrate the extended heat source in a point at the AlGe interface. The temperature gradient in the Ge-segment is obtained from a linear fit to the data between heat source and heat sink. The linearity of the temperature profile within the Ge section is a result from the strong temperature gradient and the κ/g ratio which justifies neglecting substrate loss g in this region. The effective Peltier coefficient is then calculated as Π = 267 ± 25 mW A−1, resulting in a Seebeck coefficient of S = 790 ± 95 μV K−1 for the temperature at the Al–Ge interface.
It is interesting to relate the effective Peltier coefficient to the Schottky barrier height to test the validity of the derived value. The cooling power by thermionic emission over a barrier is given by the total current times the average energy of the carriers over the barrier
(13) |
Hence, the barrier height may be calculated from our measurements. We find φc = 325 ± 45 meV, which is in agreement with previously determined barrier heights from gating experiments of similar devices.20 Independent of metal type and doping concentrations, metal–germanium junctions form Schottky contacts and exhibit very strong Fermi-level pinning close to the valence band.34,35 Finally, the thermoelectric figure of merit is calculated for the scan taken at a current of I = 26 μA to be ZT = 0.020 ± 0.005. This finalizes the thermoelectric characterization of the Al–Ge heterostructure nanowire with a segment length of 168 nm from a single measurement.
Lastly, we have a closer look at the Ge–Al interface regions. We acknowledge that extracting quantitative results is hampered by the oxide shell, which does not comprise an interface. However, important observations can nevertheless be made. The conjecture that the Peltier source term has an exponential spatial distribution has testable consequences. As a result of the distributed source term (eqn (3)), the point of maximum temperature should be shifted with respect to the Al–Ge interface. Using the differential equation with this source term and the material parameters extracted above, we calculate that the maximum temperature should be about 46 nm separated from the interface using an equilibration length of 22 nm. This shift is a good estimation with the position of maximum temperature extracted from the experimental temperature profile of 58 ± 10 nm.
Continuous support from Anna Fontcuberta i Morral, Ilaria Zardo, Marta De Luca, Kirsten Moselund, Heike Riel and Steffen Reidt is gratefully acknowledged.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/D0NR05504B |
This journal is © The Royal Society of Chemistry 2020 |