Min-Seung
Jo
ab,
Beom-Jun
Kim
a,
Myung-Kun
Chung
a,
Se-Yoon
Jung
a,
Min-Ho
Seo
cd,
Jae-Young
Yoo
e,
Jae-Soon
Yang
a,
Sung-Ho
Kim
a and
Jun-Bo
Yoon
*a
aSchool of Electrical Engineering, Korea Advanced Institute of Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea. E-mail: jbyoon@kaist.ac.kr
bCenter for Bio-Integrated Electronics, Northwestern University, Evanston, IL 60208, USA
cDepartment of Information Convergence Engineering, College of Information and Biomedical Engineering, Pusan National University, Republic of Korea
dSchool of Biomedical Convergence Engineering, College of Information and Biomedical Engineering, Pusan National University, Republic of Korea
eDepartment of Semiconductor Convergence Engineering, Sungkyunkwan University, Suwon 16419, Republic of Korea
First published on 11th February 2025
Thermal management at the nanoscale offers potential advancements in power-efficient geometrical design; however, both conduction- and convection-based structural optimization have yet to be fully investigated due to dimensional limitations. In this study, we analytically designed a structured configuration within a regime comparable to the mean free path of heat-transferring carriers. The optimally designed nanowire configuration, featuring aligned nanowires with narrow gaps (∼22 nm), was designed based on the analytic calculation. A tailored nanofabrication technique enabled a reliable geometrical parametric study, experimentally validating the proposed theoretical design. Finally, by engineering both conduction heat loss and air molecular interactions in convective heat loss at the nanoscale, we suggest an optimized heater for atmospheric conditions based on scaling factor-induced thermophysical properties, surpassing the efficiency of traditional film-based heaters.
In terms of thermal conduction, when materials are scaled down to the nanoscale, comparable to the mean free path of heat-transferring carriers such as electrons and phonons, the resulting surface scattering effectively suppresses heat loss. Additionally, the inherently small cross-sectional area of nanomaterials further minimizes heat loss, making them highly effective when considered as low-power heating elements. Indeed, the use of nanomaterials in thermal management has facilitated the development of power-efficient heating elements, which successfully demonstrated low-power gas sensing electronics.11–13 Nevertheless, convection has often been overlooked in these designs, primarily due to concerns about its potential to enlarge the air conductive heat loss.
With respect to thermal convection, as the size decreases, the increased surface-to-volume ratio induces the thinner thermal boundary and greater heat loss path around the exposed surface areas, negatively affecting power-efficient heating. Although heat transfer based on molecular collisions decreases due to the reduced collision frequency in smaller characteristic length regimes,14,15 these complex molecular interactions have not been effectively applied to heating components in the transition regime (0.1 < Knudsen number (Kn) < 10). In this context, convective heat transfer based on geometrical scaling has been primarily explored for individual nanostructures,16,17 with experimental studies limited to the sub-micro scale and conditions with an extended air mean free path through pressure control.10,18–20 Meanwhile, advanced geometrical designs that modify arrangements have been proposed to further suppress molecular collisions for more efficient thermal management.21 Specifically, structures arranged with narrow gaps comparable to the air mean free path can significantly suppress molecular diffusion, resulting in reduced effective heat loss in the assembly of structure. However, this arrangement of structures to restrict convective heat transfer has only been explored at the microscale under vacuum conditions, which artificially extend the air mean free path at that scale.19 These partially explored heating configurations, aimed at suppressing thermal conduction at the nanoscale and convective heat transfer through geometric arrangements, have yet to be fully realized, indicating that engineered design and demonstration are required for optimization.
Here, the optimized efficient heating architecture at the nanoscale is designed to address both thermal conduction and convective heat loss for low-power heating under atmospheric conditions. In pursuit of this goal, a highly aligned arrangement of thin nanowire (<100 nm) with extremely narrow nanogaps (<30 nm), comparable to both electron and air mean free path of heat-transferring carriers, is proposed. An analytic model was established to predict the geometry-dependent conductive and convective heat loss in heating components, and experimental validation was conducted using a reliable fabrication method of platinum (Pt) nanowire array. Ultimately, the optimized heating configuration was successfully demonstrated by parametric study, exhibiting strong agreement with the developed analytic model. The proposed design presents more efficient heating than traditional thin film heater even at the same joule-heating current density.
First, conductive heat (Qcond), or beam conduction, is transferred to the substrate along the axial direction of the beam through the electrodes and is proportional to the thermal conductivity (k) and cross-sectional area of the beam (Ac) (Fig. 1b(i)). As mentioned earlier, in materials smaller than 100 nm in thickness, surface scattering of heat-transferring carriers, electrons in metals, becomes more intense. This increased surface scattering raises electrical resistivity, which is directly proportional to thermal resistivity, thereby reducing beam conduction.
Second, convective heat (Qconv) is transferred from the material surface to the surrounding air (Fig. 1b(ii)). As the material size decreases, convective heat transfer driven by bulk fluidic motion decreases, while heat transfer due to air molecular interactions increases.10,18,23 Particularly in the high Knudsen number regime (Kn > 0.1), thinner effective thermal boundaries caused by surface interactions lead to sharp temperature gradients at the surface, enhancing overall heat transfer in nanosized material surfaces.18,24 However, when nanosized materials are placed close to each other, comparable to the air mean free path (λair), molecular interactions within each other's thermal boundaries can be suppressed, reducing the increased convective heat transfer coefficient.21 In this regard, designing geometrical configurations with narrow gaps between nanostructures can effectively mitigate beam and air conduction heat loss, resulting in improved thermal management for the optimization of heating efficiency.
In the optimization of joule-heating configuration, consideration of heat loss alongside the heating length is necessary, as both beam conduction and air conduction are significantly influenced by the axial length. Under the steady state condition of fixed room temperature at both ends of the electrodes, the joule-heated beam exhibits a parabolic temperature profile.23 The temperature along the nanowire's longitudinal direction, depending on the operating power, can be calculated using the steady-state thermal equilibrium equation (eqn (1)).25
![]() | (1) |
We further established an analytic model to account for the variation in the heat transfer coefficient at the nanoscale. Based on initial experimental data across different dimensions and the two-layer model in high Knudsen number regime,19,26 the conductive and convective heat transfer coefficients are calculated. First, in terms of conductive heat transfer, thermal conductivity was calculated based on electrical resistivity measurements using the Wiedemann–Franz law27 (Fig. 2b). The reduced thermal conductivity, observed when the thickness of the film and the characteristic length (D) of the nanowire decreased below 100 nm, was implemented. Here, for analytic calculation, the defined characteristic length and the structural parameter are shown in ESI Fig. S2.†
Next, the convective heat transfer coefficient for nanosized materials in the transition regime at atmospheric pressure is expressed as h = Nu·kair/D. The Nusselt number (Nu) for Pt is calculated in the transition regime, considering both ballistic transport and diffusion, based on an experimentally determined slip length.18,20 Air conductivity (kair) is used as 0.0259 W m−1 K−1.10,28 As a result, for nanosized materials at room temperature, the convective heat transfer coefficient increases as the characteristic length decreases, due to diffusion-dominant molecular interactions within the thinner thermal boundary layer (Fig. 2c). Then, as the gap between the material arrays decreased, the convective heat transfer coefficient was calculated as a ratio compared to arrays spaced far enough apart to avoid mutual influence (Fig. 2d). Based on the two-layer model,18,19 in the continuum layer outside the material surface, air molecular collisions are dominant. When the continuum layers of adjacent two wires become close enough to overlap, we assumed that heat loss path is blocked due to molecular collisions between them, and the total convective heat transfer coefficient for the array was decreased accordingly (ESI Fig. S3†). Consequently, the total convective heat transfer coefficient decreases sharply as the gap narrows, particularly when it falls below 2λair.
Assuming that reduced thermal conductivity of the metal and air conductivity is applied for heating film, power consumption according to the heating length is expected to follow as shown in Fig. 2e and f. As the conductive heat loss through the nanowires axial direction decreases, the dominant beam conduction at shorter lengths is expected to decrease, leading to reduced power consumption (Fig. 2e). On the other hand, as air conduction decreases, resulting in lower convective heat transfer, the heat loss at longer lengths is anticipated to diminish (Fig. 2f). Note that the conductive heat transfer coefficient, influenced by the characteristic length, also affects the exposed heat loss path due to the changed surface area. To achieve an optimized heating component, an analytic model was developed that takes these factors into comprehensive consideration.
To experimentally validate the developed analytic predictions, we employed a reliable method to fabricate nanowire arrays with precisely controlled gap and size. The fabrication process is based on physical vapor deposition (PVD) on a nanograting film with a predefined period, previously developed to realize a suspended nanowire array, nanolene,29 as schematically shown in Fig. 3a. First, Si nanograting substrates with different periods (200 nm, 400 nm, and 600 nm) were prepared using krypton fluoride (KrF) lithography and a subsequent reactive ion etching (RIE) process. For the 200 nm nanograting substrate, we additionally applied spacer lithography30 and size-reduction method31 to achieve a linewidth below 100 nm, as the resolution of KrF lithography is limited to 150 nm. By starting with different periods of Si substrate, the gap between the nanowires was precisely controlled. Next, UV lithography is used for pattering the area where the nanowires would be, and Pt was obliquely deposited onto the substrate resulting in a C-shaped cross-section. Following deposition, a lift-off process was implemented to remove the patterned photoresist, leaving Pt nanowires only in the areas defined by the photoresist. The characteristic length of the nanowires was controlled by adjusting the deposition thickness and angle. Subsequently, using the conventional lift-off process, gold (Au) electrodes with a chromium (Cr) adhesion layer were formed on the patterned areas on both ends of the Pt nanowires. Here, the surface roughness of Pt as-deposited was measured as Sq = 0.431 nm using atomic force microscopy (ESI Fig. S4†), and we conducted experiment under consistent conditions to verify the analytical design without the influence of surface roughness. Finally, the Si substrate was dry-etched using xenon difluoride (XeF2) gas to make the nanowire array suspended in the air.
The fabricated nanowire array was visually inspected using an optical microscope (Fig. 3b) and scanning electron microscope (SEM) (Fig. 3c). Four electrodes were configured for each nanowire array to perform 4-point probing, allowing for accurate measurement of the induced electrical power on the nanowires while eliminating the effects of contact resistance. The fabricated nanowire arrays were reliably formed with consistent dimensions over a large area of the Si substrate (inset, Fig. 3b). A cross-sectional image was obtained by cutting the nanowire array cross-section after sputtering an additional Pt layer inside the SEM as a passivation layer to protect against the ion beam (Fig. 3c, upper left). The nanowires, featuring high mechanical bending stiffness due to their inherent C-shaped cross-section, remained perfectly separated, even at spacings of less than 100 nm (Fig. 3c, upper right). Note that this C-shaped configuration of nanowires offers advantages in mechanical stability, ensuring reliable device applications through thermal and mechanical durability tests.29,32
Using this reliable fabrication method, we realized nanowires with specific sizes and gaps to independently compare beam conduction and air conduction, ultimately optimizing heat loss. The actual cross-sections of the fabricated nanowire arrays, which enable a parametric study, are shown in Fig. 3d. By forming nanowires with different thicknesses on a fixed period of 600 nm, the size-dependent beam conduction, due to surface scattering, can be verified. For gaps larger than 200 nm, the effect of air conduction on the gap size is expected to be negligible. Subsequently, by controlling the gap between nanowires on Si substrates with different periods of Si nanograting substrate, while maintaining the same thickness of Pt, the gap-dependent air conduction effect can be studied under the same current density. By achieving gaps between nanowires smaller than twice the air mean free path (∼136 nm),15 suppressed air conduction is expected under atmospheric pressure (760 Torr).
The fabricated nanowires were electrically examined, as shown in Fig. 4a and b. For the thinnest nanowires with a thickness of 25 nm, the I–V curve demonstrates reliable electrical performance, even with a length of 100 μm, as resistance increases with length (Fig. 4a). We then compared the electrical resistivity with respect to the nanowire and film thickness (Fig. 4b). The resistivity increased proportionally to the reciprocal of the thickness, following the Fuch–Sondheimer equations,33 which account for the increased electrical resistivity in thin films due to surface scattering. Even with the same thickness, nanowires exhibited higher electrical resistivity compared to films because of their confinement along multiple sides. We electrically measured the temperature coefficient of resistance by controlling the ambient temperature conditions and used those values to normalize the measured resistance to temperature. To clearly compare the heating efficiency with respect to the dimensions, we compared the nanowire array using the same current density, taking a 50 nm film with a width of 24 μm as the reference. Before making this comparison, we measured the cross-sectional area of the individual nanowire from the SEM images and then normalized to the same cross-sectional area (Fig. 4c). For further calculations based on the nanowire dimensional parameters, the perimeter and cross-sectional area were directly calculated from the SEM images (ESI, Fig. S5†).
![]() | ||
Fig. 4 (a) I–V measurement data for different lengths of nanowires with a thickness of 25 nm, as shown in Fig. 3. (b) Electrical resistivity as a function of nanowire and film thickness. (c) Measured cross-sectional areas of nanowires with thicknesses of 25 nm, 50 nm, and 75 nm, along with their corresponding SEM images (measurements from three different nanowires, scale bar = 200 nm). (d) Calculated (dashed line) and measured (dots) operating power required to reach 100 °C for heaters of different lengths, using a 50 nm film with a width of 24 μm as the reference. (e) Operating power extracted from (d) at a length of 10 μm and at the length with minimum operating power. |
The experimental results and analytic calculations of the operating power required to reach an average beam temperature of 100 °C as a function of length are shown in Fig. 4d and ESI Fig. S6,† corresponding to the nanowires shown in Fig. 3d(i)–d(iii). As expected in Fig. 2e, due to the increased thermal resistivity of the nanowires, beam conduction decreased, as evident from the power consumption at shorter lengths. However, as the heater length increases, the power consumption of the nanowire array with a thickness of 25 nm surpasses that of the film. This is attributed to the smaller characteristic length of the thin nanowires, which leads to increased air conduction of individual nanowire. The extracted power consumption at a length of 10 μm and at the optimal length is shown in Fig. 4e. At 10 μm-heater length, power consumption decreased in the order of the film, followed by nanowires with smaller thicknesses. However, the minimum power consumption at the optimal length was observed with a thickness of 75 nm. This suggests that simply reducing the size of nanowires to suppress beam conduction does not necessarily lead to lower power consumption, as there is an optimal thickness, compared to films, which balances the increased air conduction.
With the 75 nm thickness nanowires, we further compared the operating power required to reach 100 °C at different gaps, corresponding to the nanowires shown in Fig. 3d(iii)–3d(v). According to our analytic calculations, the effective convective heat transfer coefficient of the nanowire array decreases as the gap becomes smaller, compared to that of a single nanowire (Fig. 5a). The fabricated nanowires with smaller gaps resulted in nanowires having smaller characteristic lengths. Hence, the convective coefficient increases as the individual nanowires are made denser, since it is inversely proportional to the characteristic length. However, when the gap is reduced to below 100 nm, the effective heat transfer coefficient can become even lower. The operating power of these nanowires with a reduced convective heat transfer coefficient was compared across different lengths in both experimental data and analytic calculations (Fig. 5b). The proposed calculations, which account for the reduced convective heat transfer coefficient with respect to characteristic length and gap, exhibited strong agreement with the measurements. As predicted in Fig. 2f, the reduction in power consumption due to decreased air conduction is evidenced in longer heater lengths. This reduction is clearly confirmed in the comparison of power consumption at exceedingly long lengths, such as 500 μm, where air conduction becomes the dominant heat loss path over beam conduction (Fig. 5c and ESI Fig. S7†). The optimal power consumption for the nanowires with the same thickness occurred at a longer optimal length. This extension of the optimal length indicates that the balance point between beam conduction and air conduction has shifted due to the reduced convection. Finally, for the optimally designed nanowire thickness, gap, and length, significant agreement of the proportional increase in temperature with operating power is confirmed in both analytic calculations and experimental measurements (Fig. 5d).
Footnote |
† Electronic supplementary information (ESI) available: Analytic calculation of power consumption in Pt film heater, geometrical parameters and symbols for structural design, schematic of the air-conduction model, calculated and measured power consumption of the fabricated nanowire heaters. See DOI: https://doi.org/10.1039/d4na00894d |
This journal is © The Royal Society of Chemistry 2025 |