Non-equilibrium molecular dynamics calculation of heat conduction in liquid and through liquid-gas interface

Significantly Enhanced Thermal Conductivity of hBN/PTFE Composites: A Comprehensive Study of Filler Size and Dispersion

High-temperature polymers are attractive for applications in extreme temperatures, where they maintain their mechanical flexibility and electrical insulating properties. However, their heat dissipation capability is limited due to their intrinsically low thermal conductivities. Hexagonal boron nitride (hBN) is a chemically inert, thermally stable, and electrically insulative compound with a high thermal conductivity, making it an ideal candidate as a filler within a high-temperature polymer matrix to increase the thermal conductivity. This study evaluates the effect of filler size and dispersion on thermal conductivity by producing homogeneous composite samples using a combination of solvent mixing and resonant acoustic mixing (RAM). We carefully characterized our samples, including the spread of the size distribution, and observed that the smaller sized hBN centered around 5 μm was able to integrate more seamlessly into the polytetrafluoroethylene (PTFE) matrix with particle size in the 15 μm range and hence outperformed 30 μm, in contrast to the conventional wisdom, which asserts that larger fillers universally perform better than smaller ones. Our thermal conductivity of hBN/PTFE composites at 30 wt % is 2× higher than the literature values. Notably, we reached the record-high value of 3.5 W/m K at 40 wt % with an onset of percolation at 20 wt %, attributed to optimized hBN dispersion that facilitates the formation of thermal percolation. Our findings provide general guidelines to enhance the thermal conductivity of polymer composites for thermal management, ranging from power transmission to microelectronics cooling.

A supervised learning-assisted multi-scale study for thermal and mechanical behavior of porous Silica

This paper presents a comprehensive investigation of mesoporous Silica utilizing a multi-scale modeling approach under periodic boundary conditions integrated with machine learning algorithms. The study begins with Molecular Dynamics (MD) simulations to extract Silica's elastic properties and thermal conductivity at the nano-scale, employing the Tersoff potential. Subsequently, the derived material characteristics are applied to a series of generated porous Representative Volume Elements (RVEs) at the microscale. This phase involves the exploration of porosity and void shape effects on Silica's thermal and mechanical properties, considering inhomogeneities' distributions along the X-axis and random dispersion of pore cells within a three-dimensional space. Furthermore, the influence of pore shape is examined by defining open and closed-cell models, encompassing spherical and ellipsoidal voids with aspect ratios of 2 and 4. To predict the properties of porous Silica, a shallow Artificial Neural Network (ANN) is deployed, utilizing geometric parameters of the RVEs and porosity. Subsequently, it is revealed that Silica's thermal and mechanical behavior is linked to pore geometry, distribution, and porosity model. Finally, to classify the behavior of porous Silica into three categories, quasi-isotropic, orthotropic, and transversely-isotropic, three methodologies of decision tree approach, K-Nearest Neighbors (KNN) algorithm, and Support Vector Machines (SVMs) are employed. Among these, SVMs employing a quadratic kernel function demonstrate robust performance in categorizing the thermal and mechanical behavior of porous Silica.

Thermal Conductivity Calculation in Organic Liquids: Application to Poly-α-Olefin

In this work, we aim to understand and predict the thermal properties of automotive lubricants using non-equilibrium molecular dynamics. After a previous study on model materials for the mechanical parts of a car engine, we now focus on the thermal conductivity κ of the poly-α-olefin base oil (PAO4) using the well-known sink and source method to study the response of the system to an imposed heat flux. We present a detailed methodology for the calculation of κ, taking into account specific constraints related to the system under study, such as large steady-state fluctuations and rapidly growing stationarization times. We provide thermal conductivity results using four different force fields, including OPLS-AA, PCFF and COMPASS, in a temperature range of 300 to 500 K, which corresponds to the typical operating range of a car engine. The results are compared to experimental measurements performed on the commercial compound using the laser flash method. Agreement at room temperature is shown to be excellent for our in-house force field.

Significant thermal rectification induced by phonon mismatch of functional groups in a single-molecule junction

Single-molecule junctions (SMJs) may bring exotic physical effects. In this work, a significant thermal rectification effect is observed in a cross-dimensional system, comprising a diamond, a single-molecule junction, and a carbon nanotube (CNT). The molecular dynamics simulations indicate that the interfacial thermal resistance varies with the direction of heat flow, the orientation of the crystal planes of the diamond, and the length of the CNT. We find that the thermal rectification ratio escalates with the length of the CNT, achieving a peak value of 730% with the CNT length of 200 nm. A detailed analysis of phonon vibrations suggests that the primary cause of thermal rectification is the mismatched vibrations between the biphenyl and carbonyl groups. This discovery may offer theoretical insights for both the experimental exploration and practical application of SMJs in efficient thermal management strategy for high power and highly integrated chips.

The Effect of Copper-Graphene Composite Architecture on Thermal Transport Efficiency

This paper presents the results of molecular dynamic modeling, revealing that inserting confined graphene layers into copper crystal reduces the thermal conductivity of the whole composite, and the coefficient of thermal conductivity κ decreases upon an increase in the number of graphene layers. The injection of one, two, and three layers of 15 nm graphene leads to a change in the coefficient of thermal conductivity from 380 W/(m·K) down to 205.9, 179.1, and 163.6 W/(m·K), respectively. Decreasing the length of graphene layers leads to a decrease in the density of defects on which heat is dissipated. With one, two, and three layers of 8 nm graphene, the coefficient of thermal conductivity of the composite is equal to 272.6, 246.8, and 240.8 W/(m·K), appropriately. Meanwhile the introduction of an infinite graphene layer results in the growth of κ to 414.2-803.3 W/(m·K).

Analysis of the mechanism of enhanced heat transfer by nanofluids

CONTEXT

Industrial production and humans cannot exist without energy, but the low efficiency of the heat transfer in the excessive use of energy is the most significant aspect of energy saving and emission reduction. Molecular dynamics simulation methods are devoted to simulate the heat transfer efficiency of a nanofluid system with different particle sizes, and the heat transfer enhancement mechanism of the nanofluid is simulated and studied from a microscopic perspective. The analysis showed that as nanoparticle size increases, the thermal conductivity of the Al-Ar nanofluid tends to decrease, but all of them are still higher than the thermal conductivity of the liquid argon system. According to the findings of the density and radial distribution function analyses, it can be seen that the microstructure of the system changes after putting solid nanoparticles to the original fluid. This alteration in the system's microstructure is the primary component responsible for the increased heat transfer efficiency of nanofluids.

METHODS

In this paper, based on the theory of molecular dynamics, the simulation calculations were mainly performed using LAMMPS software, which is a commonly used open source computational program in the field of MD simulation research. VMD is used for visualization and analysis. The Lennard-Jones potential function was used in the simulation to accurately describe the forces acting between the atoms.

Development of Molecular Dynamics and Research Progress in the Study of Slag

Molecular dynamics is a method of studying microstructure and properties by calculating and simulating the movement and interaction of molecules. The molecular dynamics simulation method has become an important method for studying the structural and dynamic characteristics of slag systems and can make up for the shortcomings of existing detection methods and experiments. Firstly, this paper analyzes the development process and application fields of molecular dynamics, summarizes the general simulation steps and software algorithms of molecular dynamics simulation methods, and discusses the advantages and disadvantages of the algorithms and the common functions of the software. Secondly, the research status and application progress of molecular dynamics simulation methods in the study of phosphate, silicate, aluminate and aluminosilicate are introduced. On this basis, a method of combining molecular dynamics simulation with laboratory experiments is proposed, which will help obtain more accurate simulation results. This review provides theoretical guidance and a technical framework for the effective analysis of the microstructure of different slag systems via molecular dynamics, so as to finally meet the needs of iron and steel enterprises in producing high-quality steel grades.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

Mechanical and thermal properties of graphyne-coated carbon nanotubes: a molecular dynamics simulation on one-dimensional all-carbon van der Waals heterostructures

The mechanical and thermal properties of a hybrid nanotube consisting of a coaxial carbon nanotube (CNT) inside a graphyne nanotube (GNT), i.e., CNT@GNT, are investigated in this paper by using molecular dynamics simulations. The results show that the mechanical properties of CNT@GNT under uniaxial tension depend on the nanotube chirality of its components. Specifically, the Young's modulus of the CNT@GNT structure with an inner zigzag CNT is larger than that of its counterpart with an armchair CNT, while CNT@GNT with an armchair CNT and a zigzag GNT is found to possess the largest tensile strength and fracture strain. In addition, a unique fracture behavior of the successive rupture of its two components is observed in CNT@GNT. The thermal conductivity of CNT@GNT is found to be almost independent of the nanotube chirality of its components but increases as the length and diameter of the CNT@GNT increase. Moreover, strain engineering is shown as an effective avenue to modulate the thermal conductivity of CNT@GNT, which can be enhanced by tension but reduced by compression. The analysis of the phonon spectrum and spectral energy density demonstrates that this strain effect originates from changes of the phonon group velocity and phonon scattering in the strained CNT@GNT.

Molecular dynamics simulation of thermal transport across a solid/liquid interface created by a meniscus

Understandings heat transfer across a solid/liquid interface is crucial for establishing novel thermal control pathways in a range of energy applications. One of the major problems raised in this context is the impact of the three-phase contact line between solid, liquid, and gas on heat flux perturbations at the nanoscale. The focus of this research is the thermal transport via nanosized meniscus restricted between two solid walls. The molecular dynamics approach was used to consider different wetting states of the meniscus by varying the interaction potential between atoms of the substrate and the liquid. The influence of the meniscus size on the energy exchange between two solid walls was also studied. It was discovered that possessing a three-phase contact line reduces the interfacial boundary resistance between solid and liquid. Furthermore, the finite element method was employed to connect atomistic simulations with continuum mechanics. We show that the wetting angle and interfacial boundary resistance are essential important parameters for multiscale analysis of thermal engineering issues with precise microscale parametrization.

Wide-Temperature Tunable Phonon Thermal Switch Based on Ferroelectric Domain Walls of Tetragonal KTN Single Crystal

Ferroelectric domain walls (DWs) of perovskite oxide materials, which can be written and erased by an external electric field, offer the possibility to dynamically manipulate phonon scattering and thermal flux behavior. Different from previous ferroelectric materials, such as BaTiO₃, PbTiO₃, etc., with an immutable and low Curie temperature. The Curie temperature of perovskite oxide KTa_1-xNb_xO₃ (KTN) crystal can be tuned by altering the Ta/Nb ratio. In this work, the ferroelectric KTa_0.6Nb_0.4O₃ (KTN) single crystal is obtained by the Czochralski method. To understand the role of ferroelectric domains in thermal transport behavior, we perform a nonequilibrium molecular dynamics (NEMD) calculation on monodomain and 90° DWs of KTN at room temperature. The calculated thermal conductivity of monodomain KTN is 9.84 W/(m·k), consistent with experimental results of 8.96 W/(m·k), and distinctly decreased with the number of DWs indicating the outstanding performance of the thermal switch. We further evaluate the thermal boundary resistance (TBR) of KTN DWs. An interfacial thermal resistance value of 2.29 × 10^-9 K·m²/W and a large thermal switch ratio of 4.76 was obtained for a single DW of KTN. Our study shows that the ferroelectric KTN can provide great potential for the application of thermal switch at room temperature and over a broad temperature range.

Graphene and 2D Hexagonal Boron Nitride Heterostructure for Thermal Management in Actively Tunable Manner

Thermal management is a critical task for highly integrated or high-power semiconductor devices. Low dimensional materials including graphene and single-layer hexagonal boron nitride (BN) are attractive candidates for this task because of their high thermal conductivity, semi-conductivity and other excellent physical properties. The similarities in crystal structure and chemistry between graphene and boron nitride provide the possibility of constructing graphene/BN heterostructures bearing unique functions. In this paper, we investigated the interfacial thermal transport properties of graphene/BN nanosheets via non-equilibrium molecular dynamics (NEMD) simulations. We observed a significant thermal rectification behavior of these graphene/BN nanosheets, and the rectification ratio increased with the system length increases up to 117%. This phenomenon is attributed to the mismatch of out-of-plane phonon vibration modes in two directions at the interface. In addition, we explored the underlying mechanism of the length dependence of the thermal transport properties. The results show promise for the thermal management of this two-dimensional heterostructure in an actively tunable manner.

Mechanical and thermal properties of carbon nanotubes in carbon nanotube fibers under tension-torsion loading

In carbon nanotube fibers (CNFs) fabricated by spinning methods, it is well-known that the mechanical and thermal performances of CNFs are highly dependent on the mechanical and thermal properties of the inherent CNTs. Furthermore, long CNTs are usually preferred to assemble CNFs because the interaction and entanglement between long CNTs are effectively stronger than between short CNTs. However, in CNFs fabricated using long CNTs, the interior carbon nanotubes (CNTs) inevitably undergo both tension and torsion loading when they are stretched, which would influence the mechanical and thermal performances of CNFs. Here, molecular dynamics (MD) simulations were carried out to study the mechanical and thermal properties of individual CNTs under tension-torsion loading. As for mechanical properties, it was found that both the fracture strength and Young's modulus of CNTs decreased as the twist angle α increased. Besides, step-wise fracture happened due to stress concentration when the twisted CNTs are stretched. On the other hand, it could be seen that the thermal conductivity of CNTs decreased as α increased. This work presents the systematic investigation of the mechanical and thermal properties of CNTs under tension-torsion loading and provides a theoretical guideline for the design and fabrication of CNFs.

Local Thermodynamic Description of Isothermal Single-Phase Flow in Simple Porous Media

Darcy’s law for porous media transport is given a new local thermodynamic basis in terms of the grand potential of confined fluids. The local effective pressure gradient is determined using non-equilibrium molecular dynamics, and the hydraulic conductivity and permeability are investigated. The transport coefficients are determined for single-phase flow in face-centered cubic lattices of solid spheres. The porosity changed from that in the closest packing of spheres to near unity in a pure fluid, while the fluid mass density varied from that of a dilute gas to a dense liquid. The permeability varied between $$5.7 \times {10^{-20}} \hbox {m}^2$$ 5.7 × 10 - 20 m 2 and $$5.5 \times {10^{-17}} \hbox {m}^2$$ 5.5 × 10 - 17 m 2 , showing a porosity-dependent Klinkenberg effect. Both transport coefficients depended on the average fluid mass density and porosity but in different ways. These results set the stage for a non-equilibrium thermodynamic investigation of coupled transport of multi-phase fluids in complex media.

Prediction of Thermal Conductivities of Rubbers by MD Simulations-New Insights

In this article, two main approaches to the prediction of thermal conductivities by molecular dynamics (MD) simulations are discussed, namely non-equilibrium molecular dynamics simulations (NEMD) and the application of the Green–Kubo formula, i.e., EMD. NEMD methods are more affected by size effects than EMD methods. The thermal conductivities of silicone rubbers in special were found as a function of the degree of crosslinking. Moreover, the thermal conductivities of thermoplastic polyurethane as function of the mass fraction of soft segments were obtained by those MD simulations. All results are in good agreement with data from the experimental literature. After the analysis of normalized heat flux autocorrelation functions, it has been revealed that heat in the polymers is mainly transferred by low-frequency phonons. Simulation details as well as advantages and disadvantages of the single methods are discussed in the article.

Tuning the interfacial friction force and thermal conductance by altering phonon properties at contact interface

Controlling friction force and thermal conductance at solid/solid interface is of great importance but remains a significant challenge. In this work, we propose a method to control the matching degree of phonon spectra at the interface through modifying the atomic mass of contact materials, thereby regulating the interfacial friction force and thermal conductance. Results of Debye theory and molecular dynamics simulations show that the cutoff frequency of phonon spectrum decreases with increasing atomic mass. Thus, two contact surfaces with equal atomic mass have same vibrational characteristics, so that more phonons could pass through the interface. In these regards, the coupling strength of phonon modes on contact surfaces makes it possible to gain insight into the nonmonotonic variation of interfacial friction force and thermal conductance. Our investigations suggest that the overlap of phonon modes increases energy scattering channels and therefore phonon transmission at the interface, and finally, an enhanced energy dissipation in friction and heat transfer ability at interface.

Nonequilibrium molecular dynamics method based on coarse-graining formalism: Application to a nonuniform temperature field system

A nonequilibrium molecular dynamics method is proposed to produce nonequilibrium states flexibly. In this method, virtual points are set in a simulation box, and coarse-grained physical quantities at these points are constrained using Gauss's principle of least constraint. The coarse-grained physical quantities are evaluated by averaging microscopic quantities with an appropriate weight. To obtain the weight to evaluate the coarse-grained physical quantities, a shape function matrix is initially constructed from the particle configuration. This matrix expresses an interpolation of the physical quantities at particle positions from the coarse-grained quantities at the virtual points. Then, a matrix form of the weight is calculated as the Moore-Penrose pseudoinverse matrix for the shape function matrix. This method is applied to constrain the coarse-grained kinetic energy and produce a nonuniform temperature field in the system. The temperature profile at a nonequilibrium steady state depends on the method for constructing the shape function matrix. In particular, a local temperature coincides with the coarse-grained temperature when the shape function matrix is constructed based on a higher-order interpolation.

Anisotropic thermal conductivity and corrugated patterns in single-layer black phosphorus nanoribbon subjected to shear loading: a molecular dynamics study

Single-layer black phosphorus (SLBP) also known as phosphorene is a recently introduced two-dimensional material with unique structure and promising physical properties that has drawn considerable attention in the field of nanodevices. This structure demonstrates a high anisotropy in mechanical and thermal behavior along zigzag (ZZ) and armchair (AC) principal in-plane directions. Here in this study, it is shown that implementing shear strain on 10 nm × 50 nm SLBP nanoribbons (SLBPNRs) along ZZ and AC directions, the anisotropy leads to different corrugated patterns on the pristine structure. Applying non-equilibrium molecular dynamics under a parameterized Stillinger-Weber potential for modelling SLBP, thermal conductivity (TC) behavior of the sheared SLBPNRs with corrugated patterns are examined. The results show a higher amplitude and wavelength of the corregations on the ZZ-aligned SLBPNRs, which is around two times higher than that of AC-aligned counterparts. Although, it is also shown that unlike some other 2D materials, such as graphene, the wrinkling does not have such a significant effect on TC of SLBP. The phonon density of states results obtained in this work as well as phonon dispersion curves by first-principle calculations in other works concrete this finding. The results show small frequency shifts in both high- and low-frequency phonons, which are not strong enough to affect TC in SLBPNRs. This interesting thermal property of SLBP under shear strain suggests the great potential application of these corrugated structures in nanodevices without any loss of TC abilities.

In silico all-atom approach to thermodiffusion in dilute aqueous solutions

Thermodiffusion (or thermophoresis) is the phenomenon by which the spatial distributions of constituents of liquid or gas phases become inhomogeneous in response to a temperature gradient. It has been evidenced in a variety of systems and has many practical applications as well as implications in the context of the origins of life. A complete molecular picture of thermophoresis is still missing, and phenomenological approaches are often employed to account for the experimental observations. In particular, the amplitude of the resulting concentration-gradients (quantified by the Soret coefficient) depends on many factors that are not straightforwardly rationalized. All-atom molecular dynamics simulations appear as an exquisite tool to shed light on the molecular origins for this phenomenon in molecular systems, but the practical implementation of thermophoretic settings in silico poses significant challenges. Here, we propose a robust approach to tackle thermophoresis in dilute realistic solutions at the molecular level. We rely on a recent enhanced heat-exchange algorithm to generate temperature-gradients. We carefully assess the convergence of thermophoretic simulations in dilute aqueous solutions. We show that simulations typically need to be propagated on long timescales (hundreds of nanoseconds). We find that the magnitude of the temperature gradient and the box sizes have little effect on the measured Soret coefficients. Practical guidelines are derived from such observations. Provided with this reliable setup, we discuss the results of thermophoretic simulations on several examples of molecular, neutral solutes, which we find in very good agreement with experimental measurements regarding the concentration-, mass-, and temperature-dependence of the Soret coefficient.

Quasi-Isotropically Thermal Conductive, Highly Transparent, Insulating and Super-Flexible Polymer Films Achieved by Cross Linked 2D Hexagonal Boron Nitride Nanosheets

Polymer-based thermal management materials (TIMs) show great potentials as TIMs due to their excellent properties, such as high insulation, easy processing, and good flexibility. However, the limited thermal conductivity seriously hinders their practical applications in high heat generation devices. Herein, highly transparent, insulating, and super-flexible cellulose reinforced polyvinyl alcohol/nylon12 modified hexagonal boron nitride nanosheet (PVA/(CNC/PA-BNNS)) films with quasi-isotropic thermal conductivity are successfully fabricated through a vacuum filtration and subsequent self-assembly process. A special structure composed of horizontal stacked hexagonal boron nitride nanosheets (h-BNNSs) connected by their warping edges in longitudinal direction, which is strengthened by cellulose nanocrystals, is formed in PVA matrix during self-assembly process. This special structure makes the PVA/(CNC/PA-BNNS) films show excellent thermal conductivity with an in-plane thermal conductivity of 14.21 W m^-1 K^-1 and a through-plane thermal conductivity of 7.29 W m^-1 K^-1 . Additionally, the thermal conductive anisotropic constants of the as-obtained PVA/(CNC/PA-BNNS) films are in the range of 1 to 4 when the h-BNNS contents change from 0 to 60 wt%, exhibiting quasi-isotropic thermal conductivity. More importantly, the PVA/(CNC/PA-BNNS) films exhibit excellent transparency, super flexibility, outstanding mechanical strength, and electric insulation, making them very promising as TIMs for highly efficient heat dissipation of diverse electronic devices.

Modulation of Interfacial Thermal Transport between Fumed Silica Nanoparticles by Surface Chemical Functionalization for Advanced Thermal Insulation

Since solid-state heat transport in a highly porous nanocomposite strongly depends on the thermal boundary conductance (TBC) between constituent nanomaterials, further suppression of the TBC is important for improving performance of thermal insulators. Here, targeting a nanocomposite fabricated by stamping fumed silica nanoparticles, we perform a wide variety of surface functionalizations on fumed silica nanoparticles by a silane coupling method and investigate the impact on the thermal conductivity (K_m). The K_m of the silica nanocomposite is approximately 20 and 9 mW/m/K under atmospheric and vacuum conditions at the material density of 0.2 g/cm³ without surface functionalization, respectively, and the experimental results indicate that the K_m can be modulated depending on the chemical structure of molecules. The surface modification with a linear alkyl chain of optimal length significantly suppresses K_m by approximately 30%, and the suppression can be further enhanced to approximately 50% with an infrared opacifier. The magnitude of suppression was found to sensitively depend on the length of the terminal chain. The magnitude is also related to the number of reactive silanol groups in the chemical structure, where the surface modification with fluorocarbon gives the largest suppression. The surface hydrophobization merits thermal insulation through significant suppression of the TBC, presumably by reducing the water molecules that otherwise would serve as heat conduction channels at the interface. On the other hand, when the chain length is long, the suppression is counteracted by the enhanced phonon transmission through the silane coupling molecules that grow with the chain length. This is supported by the analytical model and present simulation results, leading to prediction of the optimal chemical structure for better thermal insulation.

Heat conduction in polymer chains with controlled end-to-end distance

The low thermal conductance of polymers is one of the major drawbacks for many polymer-based products. However, a single polymer chain when stretched can have high thermal conductivities. We use non-equilibrium molecular dynamics simulations to study the steady-state thermal conductance along finite macromolecules under mechanical control of the end-to-end distance. We find that the nature of heat transport along such chains strongly depends on mechanical tuning, leading to significantly different heat conductions and temperature profiles along the chain in the compressed-chain and stretched-chain limits. This transition between modes of behaviors appears to be a threshold phenomenon: at relatively small end-to-end distances, the thermal conductance remains almost constant as one stretches the polymer chain. At given critical end-to-end distances, thermal conductances start to increase, reaching the fully extended chain values. Correlated with this behavior are two observations: first, the temperature bias falls mostly at contacts in the fully stretched chain, while part of it falls along the molecule in the compressed limit. Second, the heat conduction does not change significantly with the chain length in the stretched-chain limit but decreases dramatically when this length increases in the compressed molecule. This suggests that heat transfer along stretched chains is mostly ballistic, while in the compressed chain, heat is transferred by diffusive mechanisms. Significantly, these trends persist also for a large range of molecular structures and force fields, and the changing behavior correlates well with mode localization properties. Similar studies conducted with disordered chains and bundles of several chains show remnants of the same behavior.

Computing the Heat Conductivity of Fluids from Density Fluctuations

Equilibrium molecular dynamics simulations, in combination with the Green-Kubo (GK) method, have been extensively used to compute the thermal conductivity of liquids. However, the GK method relies on an ambiguous definition of the microscopic heat flux, which depends on how one chooses to distribute energies over atoms. This ambiguity makes it problematic to employ the GK method for systems with nonpairwise interactions. In this work, we show that the hydrodynamic description of thermally driven density fluctuations can be used to obtain the thermal conductivity of a bulk fluid unambiguously, thereby bypassing the need to define the heat flux. We verify that, for a model fluid with only pairwise interactions, our method yields estimates of thermal conductivity consistent with the GK approach. We apply our approach to compute the thermal conductivity of a nonpairwise additive water model at supercritical conditions, and of a liquid hydrogen system described by a machine-learning interatomic potential, at 33 GPa and 2000 K.

Response of glassy liquids to thermal gradients

The Soret effect, i.e., the flow of matter caused by a temperature gradient, is studied in a glass-forming binary Lennard-Jones (LJ) mixture, using nonequilibrium molecular dynamics computer simulation. The transport processes associated with this effect are thermal diffusion and interdiffusion. While interdiffusion processes exhibit a drastic slowing down when approaching the glass transition, thermal diffusion appears to be a fast process even in the glass. We show that the Soret effect becomes more pronounced in the vicinity of the glass transition, due to the decoupling between thermal diffusion and interdiffusion as well as the chemical ordering in the considered LJ mixture. This is reflected in the occurrence of large concentration gradients, nonlinear concentration profiles, and long-lived nonstationary structures.

Influence of thermostatting on nonequilibrium molecular dynamics simulations of heat conduction in solids

Nonequilibrium molecular dynamics (NEMD) has been extensively used to study thermal transport at various length scales in many materials. In this method, two local thermostats at different temperatures are used to generate a nonequilibrium steady state with a constant heat flux. Conventionally, the thermal conductivity of a finite system is calculated as the ratio between the heat flux and the temperature gradient extracted from the linear part of the temperature profile away from the local thermostats. Here, we show that, with a proper choice of the thermostat, the nonlinear part of the temperature profile should actually not be excluded in thermal transport calculations. We compare NEMD results against those from the atomistic Green's function method in the ballistic regime and those from the homogeneous nonequilibrium molecular dynamics method in the ballistic-to-diffusive regime. These comparisons suggest that in all the transport regimes, one should directly calculate the thermal conductance from the temperature difference between the heat source and sink and, if needed, convert it into the thermal conductivity by multiplying it with the system length. Furthermore, we find that the Langevin thermostat outperforms the Nosé-Hoover (chain) thermostat in NEMD simulations because of its stochastic and local nature. We show that this is particularly important for studying asymmetric carbon-based nanostructures, for which the Nosé-Hoover thermostat can produce artifacts leading to unphysical thermal rectification.

Interfacial Thermal Conductivity and Its Anisotropy

There is a significant effort in miniaturizing nanodevices, such as semi-conductors, currently underway. However, a major challenge that is a significant bottleneck is dissipating heat generated in these energy-intensive nanodevices. In addition to being a serious operational concern (high temperatures can interfere with their efficient operation), it is a serious safety concern, as has been documented in recent reports of explosions resulting from many such overheated devices. A significant barrier to heat dissipation is the interfacial films present in these nanodevices. These interfacial films generally are not an issue in macro-devices. The research presented in this paper was an attempt to understand these interfacial resistances at the molecular level, and present possibilities for enhancing the heat dissipation rates in interfaces. We demonstrated that the thermal resistances of these interfaces were strongly anisotropic; i.e., the resistance parallel to the interface was significantly smaller than the resistance perpendicular to the interface. While the latter is well-known—usually referred to as Kapitza resistance—the anisotropy and the parallel component have previously been investigated only for solid-solid interfaces. We used molecular dynamics simulations to investigate the density profiles at the interface as a function of temperature and temperature gradient, to reveal the underlying physics of the anisotropy of thermal conductivity at solid-liquid, liquid-liquid, and solid-solid interfaces.

All- and one-particle distribution functions at nonequilibrium steady state under thermal gradient

We provide a concrete expression for the phase-space distribution function at nonequilibrium steady state under a constant thermal gradient, which is a typical system of the nonequilibrium molecular dynamics simulation of heat conduction. First, the phase-space distribution function of all particles in a local volume is formulated. Our formulation explicitly takes into account the entropy production due to the change in equilibrium thermodynamic variables in addition to the traditional entropy production described by the spatial gradients and fluxes of equilibrium thermodynamic variables. This treatment is necessary to explain the nonequilibrium response of a quantity that has no equilibrium correlation with mass and heat fluxes and is essential to correctly deduce one-particle distribution functions from the all-particle one. From the all-particle distribution function, we derive the Green-Kubo relations that express the one-particle distribution functions of density and velocity in terms of equilibrium correlation functions and verify these expressions using the molecular dynamics simulation of a Lennard-Jones liquid. These nonequilibrium one-particle distribution functions are sufficiently tractable for practical use, such as for the analytical evaluation of the nonequilibrium average of physical quantities.

The Effect of Thermal Contact Number on the Tube⁻Tube Contact Conductance of Single-Walled Carbon Nanotubes

The contact conductance of single, double, and triple thermal contacts of single-walled carbon nanotubes (SWCNTs) was investigated using molecular dynamics simulations. Our results showed that the effect of the thermal contact number on the contact conductance was not as strong as previously reported. The percentages of contact conductance of double and triple thermal contacts were about 72% and 67%, respectively, compared to that of a single thermal contact. Moreover, we found that the contact conductance of the double and triple thermal contacts was associated with the SWCNT length and the positional relationship of the thermal contacts.

Molecular dynamics simulations of the thermal conductivity of graphene for application in wearable devices

Graphene has been highlighted as a great potential material in wearable devices, owing to its extraordinary properties such as mechanical softness, high electrical conductivity and ultra-thin thickness. However, there are remaining challenges in thermal dissipation of graphene under such complicated conditions, which significantly affect the performance of portable electronics. Using molecular dynamics simulations, we have performed systematic analysis of thermal performance for the application in wearable devices in terms of graphene with defects, under uniaxial tensile strain, and vertical stress. Three kinds morphology of defects (horizontal line defect, circular defect, and vertical line defect) are constructed by deleting atoms on the pristine graphene plane. The thermal conductivity is related to the projected defected area perpendicular to the direction of the heat current. The relative thermal conductivity displays a cubic decreasing trend with the increase of uniaxial tensile strain. Besides, the thermal conductivity of graphene is not only related to the deformation quantity, but also related to the type of compression region. Our results show that the thermal conductivity decreases a lot under local stress with larger vertical deformation, while no obvious decline is observed under the global stress. This study aims to provide guidelines and ballpark estimates for experimentalists fabricating flexible devices from graphene.

Thermo-orientation in fluids of arbitrarily shaped particles

Recent nonequilibrium Molecular Dynamics (NEMD) simulations revealed preferential orientation, induced by a temperature gradient, in fluids of uncharged dumbbell-like particles. The magnitude of this phenomenon, called thermo-orientation, was found to be linear in the applied temperature gradient and to increase with the difference in shape or mass between the two beads of the particles. The underlying mechanism and the microscopic determinants of the phenomenon are not obvious. Here, after examination of the general symmetry requirements for thermo-orientation, we have extended the NEMD simulations to uncharged particles of various shapes and mass distribution, including chiral cases. The numerical results are rationalized by a microscopic model, based on the assumption of local equilibrium. This allows us to correlate the thermo-orientation response of arbitrarily shaped particles to quantities that characterize their shape and mass distribution.