Search

In this study, we employ molecular simulations to investigate the enhancement in thermal conductance at the solid/liquid interface in the presence of a meniscus reported previously (Klochko et al., Phys. Chem. Chem. Phys. 25(4):3298-3308, 2023). We vary the solid/liquid interaction strength at Lennard-Jones interfaces for both confined liquid and meniscus systems, finding that the presence of a meniscus yields an enhancement in the interfacial thermal conductance across all wettabilities. However, the magnitude of the enhancement is found to depend on the surface wettability, initially rising monotonously for low to moderate wettabilities, followed by a sharp rise between moderate and high wettabilities. The spectral decomposition of heat flux formalism was applied to understand the nature of this phenomenon further. By computing the in-plane and out-of-plane components of the heat fluxes within both the interfacial solid and liquid, we show that the initial monotonous rise in conductance enhancement predominantly stems from a rise in the coupling of out-of-plane vibrations within both the solid and the liquid. In contrast, the subsequent sharp rise at more wetting interfaces is linked to sharp increases in the utilization of the in-plane modes of the solid and liquid. These observations result from the interplay between the solid/liquid adhesive forces and the liquid/vapor interfacial tension. Our results can aid engineers in optimizing thermal transport at realistic interfaces, which is critical to designing effective cooling solutions for electronics, among other applications.

Nano-engineering crystalline materials can be used to tailor their thermal properties. By adding new nanoscale phonon scattering centers and controlling their size, one can effectively decrease the phonon mean free path and hence the thermal conductivity of a fully crystalline material. In this letter, we use the 3$\omega$ method in the temperature range of 100-300 K to experimentally report on the more than threefold reduction of the thermal conductivity of an epitaxially-grown crystalline germanium thin film with embedded polydispersed crystalline \ch{Ge3Mn5} nano-inclusions with diameters ranging from 5 to 25~nm. A detailed analysis of the structure of the thin film coupled with Monte Carlo simulations of phonon transport highlight the role of the nano-inclusions volume fraction in the reduction of the phononic contribution to the thermal conductivity, in particular its temperature dependence, leading to a phonon mean free path that is set by geometrical constraints., Comment: Applied Physics Letters, In press

There is still a gap in understanding phonons scattering by geometrical defects at the nanoscale, and it remains a significant challenge for heat transfer management in nanoscale devices and systems. In this study, we aim to explore the characteristics of phonon scattering by a single pore to gain insights into thermal transport in nanostructures. The paper outlines a methodology for assessing the spatial distribution of the magnitude of the radial, azimuthal, and polar components of the velocity of scattered phonons by a spherical pore. We demonstrated that the size parameter, commonly employed in electromagnetic wave scattering theory, is vital in determining the scattering regime. Specifically, we show that calculated scattering efficiency has the same pattern as one commonly obtained in classical wave scattering theory. However, we found that crystallographic directions are pivotal in shaping the scattering patterns, especially in the regions where scattering patterns are defined by the Mie resonances. This observation holds significance in understanding the influence of phonon coherence on thermal transport in nanostructured materials.

Understanding the physics of a three-phase contact line between gas, liquid, and solid is important for numerous applications. At the macroscale, the three-phase contact line response to an external force action is often characterized by a contact angle hysteresis, and several models are presented in the literature for its description. Yet, there is still a need for more information about such model applications at the nanoscale. In this study, a molecular dynamics approach was used to investigate the shape of a liquid droplet under an external force for different wetting regimes. In addition, an analytic model for describing the droplet shape was developed. It gives us the possibility to evaluate the receding and advancing wetting angle accurately. With our modeling, we found that the interplay between capillary forces and viscous forces is crucial to characterize the droplet shape at the nanoscale. In this frame, the importance of the rolling movement of the interface between liquid and vapor was pointed out. We also demonstrate that in the range of the external forces when capillary forces are most significant compared to others, hysteresis is well described by the macroscale Cox-Voinov model.

This paper studies thermal transport in nanoporous silicon with a significant specific surface area. First, the equilibrium molecular dynamics approach was used to obtain the dependence of thermal conductivity on a specific surface area. Then, a modified phonon transport kinetic theory-based approach was developed to analyze thermal conductivity. Two models were used to evaluate the phonon mean free path in the porous materials. The first model approximates the dependence of the mean free path only with the specific surface area, and the second one investigates the dependence of the mean free path variation with the porosity in the peculiar case of a highly porous matrix. Both models approximate molecular dynamics data well for the smaller porosity values, while the first model fails for large porosities. The second model matches well with molecular dynamics simulations for all considered ranges of the porosities. This work illustrates that the phonon mean free path dependence with the porosity/volume fraction of composite materials is essential for describing thermal transport in systems with significant surface-to-volume fractions.

This paper investigates thermal transport in a nanocomposite system "porous silicon matrix filled with ionic liquid". First, the thermal conductivity and heat capacity of two imidazolium and one ammonium ionic liquids were evaluated using the photoacoustic approach in piezoelectric configuration and differential scanning calorimetry, respectively. Then, the thermal transport properties of the composite system "ionic liquid confined inside porous silicon matrix" were investigated with the photoacoustic approach in gas-microphone configuration. It was found that a significant enhancement of the thermal conductivity of the composite when compared to the individual components, i.e. (i) more than two times for pristine porous silicon and (ii) more than eight times for ionic liquids. These results provide new paths for innovative solutions in the field of thermal management in highly efficient energy storage devices.

Understandings heat transfer across a solid/liquid interface is important to develop new pathways to improve thermal management in various energy applications. One of the important questions that arises in this context is the impact of three-phase contact line between solid, liquid and gas on the perturbations of the heat fluxes at the nanoscale. Therefore, this paper is devoted to the investigations of features of thermal transport across nanosized meniscus constrained between two solid walls. Different wetting states of the meniscus were considered with molecular dynamics approach by the variation of the interactional potential between atoms of the substrate and the liquid. The effect of the size of the meniscus on the exchange of energy between two solid walls was also investigated. It was shown that the presence of a three phase contact line leads to a decrease of the interfacial boundary resistance between solid and liquid. Further, investigations with the finite element method were used to link atomistic simulations with the continuum mechanics. We demonstrate that the wetting angle and the interfacial boundary resistance are the required key-parameters to perform multiscale simulations of such engineering problems with an accurate microscale parametrization.

Material thermal properties characterization at nanoscales remains a challenge even if progresses were done in developing specific characterization techniques like the Scanning Thermal Microscopy (SThM). In the present work, we propose a detailed procedure based on the combined use of a SThM probe characterization and its Finite Element Modeling (FEM) to recover in-operando 3omega measurements achieved under high vacuum. This approach is based on a two-step methodology: (i) a fine description of the probe s electrical and frequency behaviors in out of contact mode to determine intrinsic parameters of the SThM tip, (b) a minimization of the free parameter of our model, i.e. the contact thermal resistance, by comparing 3omega measurements to our simulations of the probe operating in contact mode. Such an approach allows us to accurately measure thermal interface resistances of the probe as a function of the strength applied between the tip and the surface for three different materials (silicon, silica and gold). In addition, FEM modeling provides insights about the 3omega-SThM technique sensitivity, as a function of probe / sample interface resistance to measure material thermal conductivity, paving the way to quantitative SThM measurements., Comment: 29 pages, 10 figures

Germanium manganese compounds exhibit a variety of stable and metastable phases with different stoichiometry. These materials entail interesting electronic, magnetic and thermal properties both in their bulk form and as heterostructures. Here we develop and validate a transferable machine learning potential, based on the high-dimensional neural network formalism, to enable the study of Mn$_x$Ge$_y$ materials over a wide range of compositions. We show that a neural network potential fitted on a minimal training set reproduces successfully the structural and vibrational properties and the thermal conductivity of systems with different local chemical environments, and it can be used to predict phononic effects in nanoscale heterostructures., Comment: 29 pages, 9 figures

The thermal transport properties of porous silicon and nano-hybrid "porous silicon/water" systems are presented here. The thermal conductivity was evaluated with equilibrium molecular dynamics technique for porous systems made of spherical voids or water-filled cavities. We revealed large thermal conductivity enhancement in the nano-hybrid systems as compared to their dry porous counterparts, which cannot be captured by effective media theory. This rise of thermal conductivity is related to the increases of the specific surface of the liquid/solid interface. We demonstrated that significant difference for more than two folds of thermal conductivity of pristine porous silicon and "porous silicon liquid/composite" is due to the liquid density fluctuation close to "solid/liquid interface" (layering effect). This effect is getting more important for the high specific surface of the interfacial area. Specifically, the enhancement of the effective thermal conductivity is 50 % for specific surface area of 0.3 (1/nm), and it increases further upon the increase of the surface to volume ratio. Our study provides valuable insights into the thermal properties of hybrid liquid/solid nanocomposites and about the importance of confined liquids within nanoporous materials.

In this work, we investigated tensile and compression forces effect on the thermal conductivity of silicon. We used equilibrium molecular dynamics approach for the evaluation of thermal conductivity considering different interatomic potentials. More specifically, we tested Stillinger-Weber, Tersoff, Environment-Dependent Interatomic Potential and Modified Embedded Atom Method potentials for the description of silicon atom motion under different strain and temperature conditions. Additionally, we extracted phonon density of states and dispersion curves from molecular dynamics simulations. These data were used for direct calculations of thermal conductivity considering the kinetic theory approach. Comparison of molecular dynamics and kinetic theory simulations results as a function of strain and temperature allowed us to investigate the different factors affecting the thermal conductivity of strained silicon.

In this work we present a molecular dynamics investigation of thermal transport in a silica-gallium nitride nanocomposite. A surprising enhancement of the thermal conductivity for crystalline volume fractions larger than 5% is found, which cannot be predicted by an effective medium approach, not even including percolation effects, the model systematically leading to an underestimation of the effective thermal conductivity. The behavior can instead be reproduced if an effective volume fraction twice larger than the real one is assumed, which translates in a percolation effect surprisingly stronger than the usual one. Such scenario can be understood in terms of a phonon tunneling between inclusions, enhanced by the iso-orientation of all particles. Indeed, if a misorientation is introduced, the thermal conductivity strongly decreases. We also show that a percolating nanocomposite clearly stand in a different position than other nanocomopsites, where thermal transport is domimnated by the interface scattering, and where parameters such as the interface density play a major role, differently from our case.

We report measurements and Monte Carlo simulations of thermal conductivity of porous 100nm- thick silicon membranes, in which size, shape and position of the pores were varied randomly. Measurements using 2-laser Raman thermometry on both plain membrane and porous membranes revealed more than 10-fold reduction of thermal conductivity compared to bulk silicon and six-fold reduction compared to non-patterned membrane for the sample with 37% filling fraction. Using Monte Carlo solution of the Boltzmann transport equation for phonons we compared different possibilities of pore organization and its influence on the thermal conductivity of the samples. The simulations confirmed that the strongest reduction of thermal conductivity is achieved for a distribution of pores with arbitrary shapes that partly overlap. Up to 15% reduction of thermal conductivity with respect to the purely circular pores was predicted for a porous membrane with 37% filling fraction. The effect of pore shape, distribution and surface roughness is further discussed.

We report on thermal transport properties in 2 and 3 dimensions interconnected nanowire networks (strings and nodes). The thermal conductivity of these nanostructures decreases in increasing the distance of the nodes, reaching ultra-low values. This effect is much more pronounced in 3D networks due to increased porosity, surface to volume ratio and the enhanced backscattering at 3D nodes compared to 2D nodes. We propose a model to estimate the thermal resistance related to the 2D and 3D interconnections in order to provide an analytic description of thermal conductivity of such nanowire networks; the latter is in good agreement with Molecular Dynamic results.

We report on thermal transport properties of wurtzite GaN in the presence of dislocations, by using molecular dynamics simulations. A variety of isolated dislocations in a nanowire configuration were analyzed and found to reduce considerably the thermal conductivity while impacting its temperature dependence in a different manner. We demonstrate that isolated screw dislocations reduce the thermal conductivity by a factor of two, while the influence of edge dislocations is less pronounced. The relative reduction of thermal conductivity is correlated with the strain energy of each of the five studied types of dislocations and the nature of the bonds around the dislocation core. The temperature dependence of the thermal conductivity follows a physical law described by a T$^{-1}$ variation in combination with an exponent factor which depends on the material's nature, the type and the structural characteristics of the dislocation's core. Furthermore, the impact of the dislocations density on the thermal conductivity of bulk GaN is examined. The variation and even the absolute values of the total thermal conductivity as a function of the dislocation density is similar for both types of dislocations. The thermal conductivity tensors along the parallel and perpendicular directions to the dislocation lines are analyzed. The discrepancy of the anisotropy of the thermal conductivity grows in increasing the density of dislocations and it is more pronounced for the systems with edge dislocations.

We study heat transport in semiconductor nanostructures by solving the Boltzmann Transport Equation (BTE) by means of the Discrete Ordinate Method (DOM). Relaxation time and phase and group velocitiy spectral dependencies are taken into account. The Holland model of phonon relaxation time is revisited and recalculated from dispersion relations (taken in litterature) in order to match bulk silicon and germanium values. This improved model is then used to predict silicon nanowire and nanofilm thermal properties in both ballistic and mesoscopic regimes.

Heat transport at nanoscales in semiconductors is investigated with a statistical method. The Boltzmann Transport Equation (BTE) which characterize phonons motion and interaction within the crystal lattice has been simulated with a Monte Carlo technique. Our model takes into account media frequency properties through the dispersion curves for longitudinal and transverse acoustic branches. The BTE collisional term involving phonons scattering processes is simulated with the Relaxation Times Approximation theory. A new distribution function accounting for the collisional processes has been developed in order to respect energy conservation during phonons scattering events. This non deterministic approach provides satisfactory results in what concerns phonons transport in both ballistic and diffusion regimes. The simulation code has been tested with silicon and germanium thin films; temperature propagation within samples is presented and compared to analytical solutions (in the diffusion regime). The two materials bulk thermal conductivity is retrieved for temperature ranging between 100 K and 500 K. Heat transfer within a plane wall with a large thermal gradient (250 K-500 K) is proposed in order to expose the model ability to simulate conductivity thermal dependence on heat exchange at nanoscales. Finally, size effects and validity of heat conduction law are investigated for several slab thicknesses.