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.     

Transport Coefficients Of Ar—Kr Mixtures by Molecular Dynamics Computer Simulation

Abstract

Equilibrium molecular dynamics computer simulations have been used to determine the transport coefficients of model Ar—Kr mixtures, which are represented by Lennard-Jones pair potentials with Lorentz—Berthelot rules for the cross-species interactions. The component self-diffusion and mutual-diffusion coefficients are calculated from time correlation functions and mean square displacements. Time correlation functions are used to evaluate the shear and bulk viscosity, thermal conductivity and the thermal diffusion coefficient (Soret/Dufour coefficient). In the case of the thermal transport coefficients, the partial enthalpy of the two species is required at each state point to define the heat flux rigorously. We obtain this and the partial volume (and species resolved chemical potential) using particle-exchange (and particle insertion) techniques implemented in separate [NPT] simulations at the same state point.

The viscoelasticity of the fluids is characterised by the relaxation times for bulk and shear stress relaxation. The results are for dense liquids close to the triple point temperature and density. Agreement with experiment and previous simulation is particularly good for the density of the mixtures, the shear modulus, shear viscosity, shear stress relaxation time and thermal conductivity. As for the single component noble gas fluids (simulated and experiment) there is a significant qualitative difference in the temperature and, for mixtures, composition dependence of the bulk viscosity.     

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.     

Thermal conductivity of methane hydrate from experiment and molecular simulation

A single-sided transient plane source technique has been used to determine the thermal conductivity and thermal diffusivity of a compacted methane hydrate sample over the temperature range of 261.5-277.4 K and at gas-phase pressures ranging from 3.8 to 14.2 MPa. The average thermal conductivity, 0.68 +/- 0.01 W/(m K), and thermal diffusivity, 2.04 x 10(-7) +/- 0.04 x 10(-7) m2/s, values are, respectively, higher and lower than previously reported values. Equilibrium molecular dynamics (MD) simulations of methane hydrate have also been performed in the NPT ensemble to estimate the thermal conductivity for methane compositions ranging from 80 to 100% of the maximum theoretical occupation, at 276 K and at pressures ranging from 0.1 to 100 MPa. Calculations were performed with three rigid potential models for water, namely, SPC/E, TIP4P-Ew, and TIP4P-FQ, the last of which includes the effects of polarizability. The thermal conductivities predicted from MD simulations were in reasonable agreement with experimental results, ranging from about 0.52 to 0.77 W/(m K) for the different potential models with the polarizable water model giving the best agreement with experiments. The MD simulation method was validated by comparing calculated and experimental thermal conductivity values for ice and liquid water. The simulations were in reasonable agreement with experimental data. The simulations predict a slight increase in the thermal conductivity with decreasing methane occupation of the hydrate cages. The thermal conductivity was found to be essentially independent of pressure in both simulations and experiments. Our experimental and simulation thermal conductivity results provide data to help predict gas hydrate stability in sediments for the purposes of production or estimating methane release into the environment due to gradual warming.     

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.     

Dependence of thermal diffusion factor of binary mixtures to the thermodynamic state by NEMD simulation

In this paper, we investigate the dependence of thermal diffusion factor and thermal conductivity to the temperature, density and mole fraction in Lennard–Jones binary mixtures of isotopes, noble gases and SF₆–noble gases by non-equilibrium molecular dynamics simulations.The results for the isotopic mixtures indicated that the density has a crucial effect on the dependence of thermal diffusion factor on the temperature. For isotope system at low density, thermal diffusion factor increased with temperature then remains constant at higher temperatures and the slope of thermal diffusion factor vs. temperature is positive while at higher density, thermal diffusion factor decreased with temperature and then fluctuate. For noble gas mixtures, thermal diffusion factor reduces with increasing of temperature and remain constant at high temperatures. For SF₆–Ar system, thermal diffusion factor has a negative slope and reduced with increasing of temperature, but remain nearly constant at high temperatures. For Xe–SF₆ thermal diffusion factor changed sign and the slope of thermal diffusion factor vs. temperature was negative. The results also show that thermal conductivity increases with temperature for all systems.The dependence of thermal diffusion factor to mole fraction of heavier component also investigated. The inverse of thermal diffusion factor versus mole fraction of heavier component is linear for isotope mixtures at thermodynamic conditions: (a) Low temperature, large mass ratio and all densities. (b) High temperature, large mass ratio and low densities. For Ne–Kr mixture, the inverse of thermal diffusion factor shows a linear dependence to the mole fraction of heavier component in moderate temperatures and all densities. For SF₆–Ar and Xe–SF₆ mixtures, the inverse of thermal diffusion factor has linear behaviour at moderate temperatures and low density and high temperature and low density, respectively.     

Transport properties of 2F <==> F2 in a temperature gradient as studied by molecular dynamics simulations

We calculate transport properties of a reacting mixture of F and F(2) from results of non-equilibrium molecular dynamics simulations. The reaction investigated is controlled by thermal diffusion and is close to local chemical equilibrium. The simulations show that a formulation of the transport problem in terms of classical non-equilibrium thermodynamics theory is sound. The chemical reaction has a large effect on the magnitude and temperature dependence of the thermal conductivity and the interdiffusion coefficient. The increase in the thermal conductivity in the presence of the chemical reaction, can be understood as a response to an imposed temperature gradient, which reduces the entropy production. The heat of transfer for the Soret stationary state was more than 100 kJ mol(-1), meaning that the Dufour and Soret effects are non-negligible in reacting mixtures. This sheds new light on the transport properties of reacting mixtures.     

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.     

Thermal conductivity of solid argon at high pressure and high temperature: a molecular dynamics study

The thermal conductivity of solid argon at high-pressure (up to 50 GPa) and high-temperature (up to 2000 K) has been calculated by equilibrium molecular dynamics simulations using the Green-Kubo formalism and an exponential-6 interatomic potential. A simple empirical expression is given for its pressure and temperature dependence. The results are compared with predictions based on kinetic theory. The relative change of the thermal conductivity lambda with density rho is found to be consistent with a partial differential ln lambda/ partial differential ln rho slope of approximately 6 in a wide range of pressures and temperatures, in good agreement with predictions based on kinetic theory.     

Transport properties of Mie(14,7) fluids: molecular dynamics simulation and theory

An extensive computer simulation study is presented for the self-diffusion coefficient, the shear viscosity, and the thermal conductivity of Mie(14,7) fluids. The time-correlation function formalism of Green-Kubo is utilized in conjunction with molecular dynamics (MD) simulations. In addition to molecular simulations, the results of a recent study [A. Eskandari Nasrabad, J. Chem. Phys. 128, 154514 (2008)] for the mean free volume are applied to calculate the self-diffusion coefficients within a free volume theory framework. A detailed comparison between the MD simulation and free volume theory results for the diffusion coefficient is given. The density fluctuation theory of shear viscosity is used to compute the shear viscosity and the results are compared to those from MD simulations. The density and temperature dependences of different time-correlation functions and transport coefficients are studied and discussed.     

Mechanisms for thermal conduction in hydrogen hydrate

Extensive equilibrium molecular dynamics simulations have been performed to investigate thermal conduction mechanisms via the Green-Kubo approach for (type II) hydrogen hydrate, at 0.05 kbar and between 30 and 250 K, for both lightly filled H(2) hydrates (1s4l) and for more densely filled H(2) systems (2s4l), in which four H(2) molecules are present in the large cavities, with respective single- and double-occupation of the small cages. The TIP4P water model was used in conjunction with a fully atomistic hydrogen potential along with long-range Ewald electrostatics. It was found that substantially less damping in guest-host energy transfer is present in hydrogen hydrate as is observed in common type I clathrates (e.g., methane hydrate), but more akin in to previous results for type II and H methane hydrate polymorphs. This gives rise to larger thermal conductivities relative to common type I hydrates, and also larger than type II and H methane hydrate polymorphs, and a more crystal-like temperature dependence of the thermal conductivity.     

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.     

Toward the correction of effective electrostatic forces in explicit-solvent molecular dynamics simulations: restraints on solvent-generated electrostatic potential and solvent polarization

Despite considerable advances in computing power, atomistic simulations under nonperiodic boundary conditions, with Coulombic electrostatic interactions and in systems large enough to reduce finite-size associated errors in thermodynamic quantities to within the thermal energy, are still not affordable. As a result, periodic boundary conditions, systems of microscopic size and effective electrostatic interaction functions are frequently resorted to. Ensuing artifacts in thermodynamic quantities are nowadays routinely corrected a posteriori, but the underlying configurational sampling still descends from spurious forces. The present study addresses this problem through the introduction of on-the-fly corrections to the physical forces during an atomistic molecular dynamics simulation. Two different approaches are suggested, where the force corrections are derived from special potential energy terms. In the first approach, the solvent-generated electrostatic potential sampled at a given atom site is restrained to a target value involving corrections for electrostatic artifacts. In the second approach, the long-range regime of the solvent polarization around a given atom site is restrained to the Born polarization, i.e., the solvent polarization corresponding to the ideal situation of a macroscopic system under nonperiodic boundary conditions and governed by Coulombic electrostatic interactions. The restraints are applied to the explicit-water simulation of a hydrated sodium ion, and the effect of the restraints on the structural and energetic properties of the solvent is illustrated. Furthermore, by means of the calculation of the charging free energy of a hydrated sodium ion, it is shown how the electrostatic potential restraint translates into the on-the-fly consideration of the corresponding free-energy correction terms. It is discussed how the restraints can be generalized to situations involving several solute particles. Although the present study considers a very simple system only, it is an important step toward the on-the-fly elimination of finite-size and approximate-electrostatic artifacts during atomistic molecular dynamics simulations.     

Molecular dynamics simulations of the thermal conductivity of methane hydrate

Nonequilibrium molecular dynamics simulations with the nonpolarizable SPC/E (Berendsen et al., J. Phys. Chem. 1987, 91, 6269) and the polarizable COS/G2 (Yu and van Gunsteren, J. Chem. Phys. 2004, 121, 9549) force fields have been employed to calculate the thermal conductivity and other associated properties of methane hydrate over a temperature range from 30 to 260 K. The calculated results are compared to experimental data over this same range. The values of the thermal conductivity calculated with the COS/G2 model are closer to the experimental values than are those calculated with the nonpolarizable SPC/E model. The calculations match the temperature trend in the experimental data at temperatures below 50 K; however, they exhibit a slight decrease in thermal conductivity at higher temperatures in comparison to an opposite trend in the experimental data. The calculated thermal conductivity values are found to be relatively insensitive to the occupancy of the cages except at low (T相似文献   

Influence of solvent on temperature and thermal peak broadening in capillary zone electrophoresis

The present paper deals with the role of the solvent on thermal peak broadening. One main solvent property that determines the magnitude of the temperature gradient due to the generation of Joule heat in capillary zone electrophoresis is the thermal conductivity. As organic solvents have lower thermal conductivity than water (methanol and acetonitrile, e.g., nearly by a factor of 3) it can be hypothesized that the temperature gradient inside the capillary is more pronounced in organic solvents compared to an aqueous solution. On the other hand, the temperature dependence of the ion mobility (which is responsible for the velocity profile and thus for thermal peak broadening) is smaller in organic solvents. To get insight into the thermal effect of the solvent, first the temperature of a solution in a cylindrical tube was calculated utilizing the heat balance equation. It was shown that the two theoretical models most common in the literature (based on the analytical solution or on an assumption of the parabolic temperature profile in the tube, respectively) give the same results. The latter model was chosen for the further calculations, adding a quadratic term to express the electric conductivity as a function of the temperature. The temperature at the inner capillary wall and center as function of the capillary dimensions and the electric power was computed for electrolytes with a given conductivity at 25.0 degrees C with water, methanol, and acetonitrile as solvents. Capillary cooling systems used were circulating liquid cooling, enforced air-cooling, and natural convection in still air. The mean temperature (averaged over the cross section) resulting from Joule heating was compared with experimentally determined temperatures established upon application of an electric field; the latter temperature was derived from the measurement of the electric conductance of the background electrolyte solution and its (measured) temperature dependence. All investigations were carried out with solutions of the same initial electric conductivity (about 0.5 S.m(-1) at 25.0 degrees C). Agreement is found for natural convection conditions, and the deviation between theoretical and experimental results for the forced air and circulated liquid cooling systems can be related to the poorly defined thermal conditions of the capillaries in commercial instrumentation (with a part in a thermostated cassette and a part outside). For given conditions the temperature gradients in the organic solvents exceed largely those in water, independent of the type of cooling. As a consequence, the thermal plate height is significantly larger in organic solvents, at least under conditions where the deviation from the Nernst-Einstein limiting case is not too high. However, even for the maximum applicable field strengths the thermal plate height contributions are negligible compared to longitudinal diffusion in all solvents.     

Thermal conductivity of nanotubes revisited: effects of chirality, isotope impurity, tube length, and temperature

We study the dependence of the thermal conductivity of single-walled nanotubes on chirality, isotope impurity, tube length, and temperature by nonequilibrium molecular-dynamics method with accurate potentials. It is found that, contrary to electronic conductivity, the thermal conductivity is insensitive to the chirality. The isotope impurity, however, can reduce the thermal conductivity up to 60% and change the temperature dependence behavior. We also found that the tube length dependence of thermal conductivity is different for nanotubes of different radii at different temperatures.     

Molecular simulations of the pressure, temperature, and chemical potential dependencies of clay swelling

A new method for the determination of clay swelling thermodynamics from computer simulation is discussed and evaluated. This method allows for the determination of temperature, pressure, and water chemical potential dependence of clay swelling from simulations at a single thermodynamic state point. The temperature dependence and pressure dependence of clay swelling are shown to be directly related to the composite system entropy and volume change, respectively, that accompany swelling. Expressions for the chemical potential dependence of clay swelling are used to determine constant pressure layer spacing and adsorption isotherms, quantities that are well suited for comparison with experimental measurements. This method is evaluated through grand isoshear ensemble simulations of Na-montmorillonite, a prototypical swelling clay. Approximations associated with all expressions are discussed with explicit calculations used to demonstrate their regimes of validity.     

Validation of intermolecular pair potential model of SiH4: molecular-dynamics simulation for saturated liquid density and thermal transport properties

We demonstrate a validation of the intermolecular pair potential model of SiH(4), which is constructed from ab initio molecular-orbital calculations and expressed as the sum of the exponential and the London dispersion terms. The saturated liquid densities of SiH(4) are calculated for temperatures from 100 to 225 K by molecular-dynamics (MD) simulation. The average deviation between the experiment and the MD simulation using the present potential model is 3.9%, while the deviations exceed 10% for other well-known potential models such as the five-center Lennard-Jones (LJ) model. Subsequently, the shear viscosity, the thermal conductivity, and the self-diffusion coefficient of liquid SiH(4) are calculated by an equilibrium MD simulation with the Green-Kubo formula from 100 to 225 K. The average deviations from experiment are 11.8% and 13.7% for the shear viscosity and the thermal conductivity, respectively. Comparing the present model with an empirical one-center LJ model, it turns out that the rotational energy transfer through the intermolecular potential energy, which comes from the anisotropic potential energy, plays an important role in the thermal conductivity of liquid SiH(4). These results indicate that the present intermolecular potential model has an ability to give realistic pictures for liquid SiH(4) through molecular simulations.     

The origin of the conductivity maximum in molten salts. I. Bismuth chloride

A new theory is presented to explain the conductivity maxima of molten salts (versus temperature and pressure). In the new theory, conductivity is due to ions hopping from counterion to counterion, and its temperature dependence can be explained with an ordinary Arrhenius equation in which the frequency prefactor A (for hopping opportunities) and activation energy E(a) (for hopping) are density dependent. The conductivity maximum is due to competing effects: as density decreases, the frequency of opportunities for hopping increases, but the probability that an opportunity is successfully hopped decreases due to rising E(a) caused by the increased hopping distance. The theory is successfully applied to molten bismuth (III) chloride, and supported by density-functional based molecular dynamics simulations which not only reproduce the conductivity maximum, but disprove the long-standing conjecture that this liquid features an equilibrium between BiCl(3) molecules, and BiCl(2)(+) and BiCl(4)(-) ions that shifts to the left with increasing temperature.