Thermal expansion and impurity effects on lattice thermal conductivity of solid argon

Thermal expansion and impurity effects on the lattice thermal conductivity of solid argon have been investigated with equilibrium molecular dynamics simulation. Thermal conductivity is simulated over the temperature range of 20-80 K. Thermal expansion effects, which strongly reduce thermal conductivity, are incorporated into the simulations using experimentally measured lattice constants of solid argon at different temperatures. It is found that the experimentally measured deviations from a T(-1) high-temperature dependence in thermal conductivity can be quantitatively attributed to thermal expansion effects. Phonon scattering on defects also contributes to the deviations. Comparison of simulation results on argon lattices with vacancy and impurity defects to those predicted from the theoretical models of Klemens and Ashegi et al. demonstrates that phonon scattering on impurities due to lattice strain is stronger than that due to differences in mass between the defect and the surrounding matrix. In addition, the results indicate the utility of molecular dynamics simulation for determining parameters in theoretical impurity scattering models under a wide range of conditions. It is also confirmed from the simulation results that thermal conductivity is not sensitive to the impurity concentration at high temperatures.     

Modifying thermal transport in electrically conducting polymers: effects of stretching and combining polymer chains

If their thermal conductivity can be lowered, polyacetylene (PA) and polyaniline (PANI) offer examples of electrically conducting polymers that can have potential use as thermoelectrics. Thermal transport in such polymers is primarily influenced by bonded interactions and chain orientations relative to the direction of heat transfer. We employ molecular dynamics simulations to investigate two mechanisms to control the phonon thermal transport in PANI and PA, namely, (1) mechanical strain and (2) polymer combinations. The molecular configurations of PA and PANI have a significant influence on their thermal transport characteristics. The axial thermal conductivity increases when a polymer is axially stretched but decreases under transverse tension. Since the strain dependence of the thermal conductivity is related to the phonon scattering among neighboring polymer chains, this behavior is examined through Herman's orientation factor that quantifies the degree of chain alignment in a given direction. The conductivity is enhanced as adjacent chains become more aligned along the direction of heat conduction but diminishes when they are orthogonally oriented to it. Physically combining these polymers reduces the thermal conductivity, which reaches a minimum value for a 2:3 PANI/PA chain ratio.     

A gentler approach to RNEMD: nonisotropic velocity scaling for computing thermal conductivity and shear viscosity

We present a new method for introducing stable nonequilibrium velocity and temperature gradients in molecular dynamics simulations of heterogeneous systems. This method extends earlier reverse nonequilibrium molecular dynamics (RNEMD) methods which use momentum exchange swapping moves. The standard swapping moves can create nonthermal velocity distributions and are difficult to use for interfacial calculations. By using nonisotropic velocity scaling (NIVS) on the molecules in specific regions of a system, it is possible to impose momentum or thermal flux between regions of a simulation while conserving the linear momentum and total energy of the system. To test the method, we have computed the thermal conductivity of model liquid and solid systems as well as the interfacial thermal conductivity of a metal-water interface. We find that the NIVS-RNEMD improves the problematic velocity distributions that develop in other RNEMD methods.     

Thermal Conductance of Graphene-Titanium Interface: A Molecular Simulation

Titanium is a commonly used material in aviation, aerospace, and military applications, due to the outstanding mechanical properties of titanium and its alloys. However, its relatively low thermal conductivity restricts its extended usage. The use of graphene as a filler shows great potential for the enhancement of thermal conductivity in titanium-based metal-matrix composites (MMCs). We used classical molecular dynamics (MD) simulation methods to explore the thermal conductance at the titanium–graphene (Ti/Gr) interface for its thermal boundary conductance, which plays an important role in the thermal properties of Ti-based MMCs. The effects of system size, layer number, temperature, and strain were considered. The results show that the thermal boundary conductance (TBC) decreases with an increasing layer number and reaches a plateau at n = 5. TBC falls under tensile strain and, in turn, it grows with compressive strain. The variation of TBC is explained qualitatively by the interfacial atomic vibration coupling factor. Our findings also provide insights into ways to optimize future thermal management based on Ti-based MMCs materials.     

Molecular theory of thermal conductivity of the Lennard-Jones fluid

In this paper the thermal conductivity of the Lennard-Jones fluid is calculated by applying the combination of the density-fluctuation theory, the modified free volume theory of diffusion, and the generic van der Waals equation of state. A Monte Carlo simulation method is used to compute the equilibrium pair-correlation function necessary for computing the mean free volume and the coefficient in the potential-energy and virial contributions to the thermal conductivity. The theoretical results are compared with our own molecular dynamics simulation results and with those reported in the literature. They agree in good accuracy over wide ranges of density and temperature examined in molecular dynamics simulations. Thus the combined theory represents a molecular theory of thermal conductivity of the Lennard-Jones fluid and by extension simple fluids, which enables us to compute the nonequilibrium quantity by means of the Monte Carlo simulations for the equilibrium pair-correlation function.     

Phonon thermal conductivity in nanolaminated composite metals via molecular dynamics

We use nonequilibrium molecular dynamics to characterize the phonon contribution to thermal conduction of Al nanostructures and the role of interfaces in metallic nanocomposites. We characterize the lattice thermal conductivity of pure Al samples as a function of size and temperature from which we obtain, using kinetic theory, the temperature dependence of the phonon mean free path. We also calculated the thermal conductivity of AlAl* and AlNi nanolaminate composites (where Al* differs from Al only in its mass) for various periodic sizes and compositions as well as the associated interfacial thermal resistivities (ITRs). We find that simple, additive models provide good estimates of the thermal conductivities of the nanocomposites in terms of those of the individual components and interfaces if size effects on the behavior of the individual components are considered. The additive models provide important insight to the decrease in thermal conductivity of the nanolaminates as their periodicity (thickness of a bilayer) is reduced to a size comparable with the phonon mean free path and break down when this characteristic size is reduced further. At this point the system can be regarded as homogeneous and the conductivity increases with decreasing periodicity of the laminates. We also observe that the ITR depends on the direction of the heat flux; this is the first molecular level characterization of such thermal diode behavior in a realistic three dimensional material.     

Molecular dynamics simulations of self-diffusion coefficient and thermal conductivity of methane at low and moderate densities

This article demonstrates a highly accurate molecular dynamics (MD) simulation of thermal conductivity of methane using an ab initio intermolecular potential. The quantum effects of the vibrational contribution to thermal conductivity are more efficiently accounted for in the present MD model by an analytical correction term as compared to by the Monte Carlo method. The average deviations between the calculated thermal conductivity and the experimental data are 0.92% for dilute methane and 1.29% for methane at moderate densities, as compared to approximately 20% or more in existing MD calculations. The results demonstrate the importance of considering vibrational contribution to the thermal conductivity which is mainly through the self-diffusion process.     

Scaling behavior in on‐chip field‐amplified sample stacking

Field amplified sample stacking (FASS) uses differential electrophoretic velocity of analyte ions in the high‐conductivity background electrolyte zone and low conductivity sample zone for increasing the analyte concentration. The stacking rate of analyte ions in FASS is limited by molecular diffusion and convective dispersion due to nonuniform electroosmotic flow (EOF). We present a theoretical scaling analysis of stacking dynamics in FASS and its validation with a large set of on‐chip sample stacking experiments and numerical simulations. Through scaling analysis, we have identified two stacking regimes that are relevant for on‐chip FASS, depending upon whether the broadening of the stacked peak is dominated by axial diffusion or convective dispersion. We show that these two regimes are characterized by distinct length and time scales, based on which we obtain simplified nondimensional relations for the temporal growth of peak concentration and width in FASS. We first verify the theoretical scaling behavior in diffusion‐ and convection‐dominated regimes using numerical simulations. Thereafter, we show that the experimental data of temporal growth of peak concentration and width at varying electric fields, conductivity gradients, and EOF exhibit the theoretically predicted scaling behavior. The scaling behavior described in this work provides insights into the effect of varying experimental parameters, such as electric field, conductivity gradient, electroosmotic mobility, and electrophoretic mobility of the analyte on the dynamics of on‐chip FASS.     

Nonequilibrium molecular dynamics simulations of the thermal conductivity of water: A systematic investigation of the SPC/E and TIP4P/2005 models

We report an extensive nonequilibrium molecular dynamics investigation of the thermal conductivity of water using two of the most accurate rigid nonpolarizable empirical models available, SPC/E and TIP4P/2005. Our study covers liquid and supercritical states. Both models predict the anomalous increase of the thermal conductivity with temperature and the thermal conductivity maximum, hence confirming their ability to reproduce the complex anomalous behaviour of water. The performance of the models strongly depends on the thermodynamic state investigated, and best agreement with experiment is obtained for states close to the liquid coexistence line and at high densities and temperatures. Considering the simplicity of these two models the overall agreement with experiments is remarkable. Our results show that explicit polarizability and molecular flexibility are not needed to reproduce the anomalous heat conduction of water.     

Prediction of the transport properties of a polyatomic gas

An ab initio molecular potential model is employed in this paper to show its excellent predictability for the transport properties of a polyatomic gas from molecular dynamics simulations. A quantum mechanical treatment of molecular vibrational energies is included in the Green and Kubo integral formulas for the calculation of the thermal conductivity by the Metropolis Monte Carlo method. Using CO₂ gas as an example, the fluid transport properties in the temperature range of 300–1000 K are calculated without using any experimental data. The accuracy of the calculated transport properties is significantly improved by the present model, especially for the thermal conductivity. The average deviations of the calculated results from the experimental data for self-diffusion coefficient, shear viscosity, thermal conductivity are, respectively, 2.32%, 0.71% and 2.30%.     

A universal gauge for thermal conductivity of silicon nanowires with different cross sectional geometries

By using molecular dynamics simulations, we study thermal conductivity of silicon nanowires (SiNWs) with different cross sectional geometries. It is found that thermal conductivity decreases monotonically with the increase of surface-to-volume ratio (SVR). More interestingly, a simple universal linear dependence of thermal conductivity on SVR is observed for SiNWs with modest cross sectional area (larger than 20 nm(2)), regardless of the cross sectional geometry. As a result, among different shaped SiNWs with the same cross sectional area, the one with triangular cross section has the lowest thermal conductivity. Our study provides not only a universal gauge for thermal conductivity among different cross sectional geometries, but also a designing guidance to tune thermal conductivity by geometry.     

Phonon coherent resonance and its effect on thermal transport in core-shell nanowires

We study heat current autocorrelation function and thermal conductivity in core-shell nanowires by using molecular dynamics simulations. Interestingly, a remarkable oscillation effect in heat current autocorrelation function is observed in core-shell NWs, while the same effect is absent in pure silicon nanowires, nanotube structures, and random doped nanowires. Detailed characterizations of the oscillation signal reveal that this intriguing oscillation is caused by the coherent resonance effect of the transverse and longitudinal phonon modes. This phonon resonance results in the localization of the longitudinal modes, which leads to the reduction of thermal conductivity in core-shell nanowires. Our study reveals that a coherent mechanism can be used to tune thermal conductivity in core-shell nanowires.     

Equilibrium and nonequilibrium molecular dynamics simulations of the thermal conductivity of molten alkali halides

The thermal conductivity of molten NaCl and KCl was calculated through the Evans-Gillan nonequilibrium molecular dynamics (NEMD) algorithm and Green-Kubo equilibrium molecular dynamics (EMD) simulations. The EMD simulations were performed for a "binary" ionic mixture and the NEMD simulations assumed a pure system for reasons discussed in this work. The cross thermoelectric coefficient obtained from Green-Kubo EMD simulations is discussed in terms of the homogeneous thermoelectric power or Seebeck coefficient of these materials. The thermal conductivity obtained from NEMD simulations is found to be in very good agreement with that obtained through Green-Kubo EMD simulations for a binary ionic mixture. This result points to a possible cancellation between the neglected "partial enthalpy" contribution to the heat flux associated with the interdiffusion of one species through the other and that part of the thermal conductivity related to the coupled fluxes of charge and heat in "binary" ionic mixtures.     

Heat conduction in chain polymer liquids: molecular dynamics study on the contributions of inter- and intramolecular energy transfer

In this paper, the molecular mechanisms which determine the thermal conductivity of long chain polymer liquids are discussed, based on the results observed in molecular dynamics simulations. Linear n-alkanes, which are typical polymer molecules, were chosen as the target of our studies. Non-equilibrium molecular dynamics simulations of bulk liquid n-alkanes under a constant temperature gradient were performed. Saturated liquids of n-alkanes with six different chain lengths were examined at the same reduced temperature (0.7T(c)), and the contributions of inter- and intramolecular energy transfer to heat conduction flux, which were identified as components of heat flux by the authors' previous study [J. Chem. Phys. 128, 044504 (2008)], were observed. The present study compared n-alkane liquids with various molecular lengths at the same reduced temperature and corresponding saturated densities, and found that the contribution of intramolecular energy transfer to the total heat flux, relative to that of intermolecular energy transfer, increased with the molecular length. The study revealed that in long chain polymer liquids, thermal energy is mainly transferred in the space along the stiff intramolecular bonds. This finding implies a connection between anisotropic thermal conductivity and the orientation of molecules in various organized structures with long polymer molecules aligned in a certain direction, which includes confined polymer liquids and self-organized structures such as membranes of amphiphilic molecules in water.     

Thermal conductivity of solid argon from molecular dynamics simulations

The thermal conductivity of solid argon in the classical limit has been calculated by equilibrium molecular dynamic simulations using the Green-Kubo formalism and a Lennard-Jones interatomic potential. Contrary to previous theoretical reports, we find that the computed thermal conductivities are in good agreement with experimental data. The computed values are also in agreement with the high-temperature limit of the three-phonon scattering contribution to the thermal conductivity. We find that finite-size effects are negligible and that phonon lifetimes have two characteristic time scales, so that agreement with kinetic theory is obtained only after appropriate averaging of the calculated phonon lifetimes.     

Centroid molecular dynamics approach to the transport properties of liquid para-hydrogen over the wide temperature range

Fundamental transport properties of liquid para-hydrogen (p-H(2)), i.e., diffusion coefficients, thermal conductivity, shear viscosity, and bulk viscosity, have been evaluated by means of the path integral centroid molecular dynamics (CMD) calculations. These transport properties have been obtained over the wide temperature range, 14-32 K. Calculated values of the diffusion coefficients and the shear viscosity are in good agreement with the experimental values at all the investigated temperatures. Although a relatively large deviation is found for the thermal conductivity, the calculated values are less than three times the amount of the experimental values at any temperature. On the other hand, the classical molecular dynamics has led all the transport properties to much larger deviation. For the bulk viscosity of liquid p-H(2), which was never known from experiments, the present CMD has given a clear temperature dependence. In addition, from the comparison based on the principle of corresponding states, it has been shown that the marked deviation of the transport properties of liquid p-H(2) from the feature which is expected from the molecular parameters is due to the quantum effect.     

Thermal conductivity of ordered molecular water

The authors use molecular dynamics simulation to investigate the thermal transport characteristics of water with various degree of orientational and translational orders induced by the application of an electric field. The authors observe that the orientational ordering of the water dipole moments has a minor effect on the thermal conductivity. However, electric-field-induced crystallization and associated translational order result in approximately a three fold increase of thermal conductivity with respect to the base water, i.e., to values comparable with those characterizing ice crystal structures.     

Investigation of nanocomposite PNIPAm-g-graphene with pre-lower critical solution temperature rapid thermal response function for chip adaptive cooling: A molecular dynamics and computational fluid dynamics simulations

Temperature-sensitive hydrogels have been widely used for rapid adaptive cooling in electronic device thermal management with promising applications. However, existing temperature-sensitive hydrogels can only regulate the flow in the chip cooling system after the ambient temperature reaches their lower critical solution temperature (LCST). Before reaching LCST, effective rapid heat dissipation for electronic chips is not achievable. This study aims to develop a temperature-sensitive hydrogel that can provide assisted adaptive cooling for electronic chips before reaching its LCST. This requires the hydrogel to have a thermal conductivity far surpassing existing hydrogel materials. Using the temperature-sensitive hydrogel PNIPAm and graphene molecules as base materials, this research utilized molecular dynamics simulations to graft graphene molecules onto PNIPAm molecules in different ways, resulting in the temperature-sensitive hydrogel material PNIPAm-g-graphene. Non-equilibrium molecular dynamics (NEMD) was employed to calculate the thermal conductivity of this material under different temperature conditions. The results indicate that the thermal conductivity of PNIPAm-g-graphene can reach up to 1.95474 W/m K (graphene grafted at  CH₃ functional group, temperature at 280 K). Compared to the thermal conductivity of PNIPAm under the same conditions (0.45 W/m K), the increase in thermal conductivity is significant, demonstrating excellent thermal conductivity compared to PNIPAm. Subsequently, this study analyzed the underlying mechanisms of different thermal conductivities in materials obtained by grafting graphene molecules at different points using the method of overlap in Phonon Density of States Curves (PDOS) from the perspective of interfacial thermal conduction. Finally, through computational fluid dynamics (CFD) simulations, this study investigates the chip's adaptive cooling performance with PNIPAm-g-graphene material. The results show that, compared to traditional temperature-sensitive hydrogels, PNIPAm-g-graphene can achieve efficient adaptive cooling of chip hotspots before the cooling fluid temperature reaches its LCST value. This finding is significant for the field of chip cooling. The study proposes a new method for rapid, adaptive cooling of chip hotspots and explores its feasibility from the perspectives of molecular dynamics and CFD simulation. It holds importance in the thermal management of electronic devices and the rapid adaptive cooling of electronic chips.     

Equilibrium calculations of viscosity and thermal conductivity across a solid-liquid interface using boundary fluctuations

We calculate viscosity and thermal conductivity in systems of Lennard-Jones particles consisting of coexisting solid and liquid with different interface wetting properties using the recently developed equilibrium boundary fluctuation theory. We compare the slip length and equivalent liquid length obtained from these calculations with those obtained from nonequilibrium molecular dynamics. The equilibrium and nonequilibrium calculations of the slip length and the sum of the thermal equivalent lengths are in good agreement. We conclude that for both interfacial properties, the nonequilibrium simulations were probing the linear response. The significant dependence of the intrinsic equivalence length on the interfacial temperature difference used to generate the thermal gradient is explained as a consequence of the different thermodynamic states of the two interfaces.