Seiji Yamashita

Effects of molecular structure on microscopic heat transport in chain polymer liquids

In this paper, we discuss the molecular mechanism of the heat conduction in a liquid, based on nonequilibrium molecular dynamics simulations of a systematic series of linear- and branched alkane liquids, as a continuation of our previous study on linear alkane [T. Ohara et al., J. Chem. Phys. 135, 034507 (2011)]. The thermal conductivities for these alkanes in a saturated liquid state at the same reduced temperature (0.7Tc) obtained from the simulations are compared in relation to the structural difference of the liquids. In order to connect the thermal energy transport characteristics with molecular structures, we introduce the new concept of the interatomic path of heat transfer (atomistic heat path, AHP), which is defined for each type of inter- and intramolecular interaction. It is found that the efficiency of intermolecular AHP is sensitive to the structure of the first neighbor shell, whereas that of intramolecular AHP is similar for different alkane species. The dependence of thermal conductivity on different lengths of the main and side chain can be understood from the natures of these inter- and intramolecular AHPs.

Equivalence of the EMD- and NEMD-based decomposition of thermal conductivity into microscopic building blocks

Thermal conductivity of a material can be comprehended as being composed of microscopic building blocks relevant to the energy transfer due to a specific microscopic process or structure. The building block is called the partial thermal conductivity (PTC). The concept of PTC is essential to evaluate the contributions of various molecular mechanisms to heat conduction and has been providing detailed knowledge of the contribution. The PTC can be evaluated by equilibrium molecular dynamics (EMD) and non-equilibrium molecular dynamics (NEMD) in different manners: the EMD evaluation utilizes the autocorrelation of spontaneous heat fluxes in an equilibrium state whereas the NEMD one is based on stationary heat fluxes in a non-equilibrium state. However, it has not been fully discussed whether the two methods give the same PTC or not. In the present study, we formulate a Green-Kubo relation, which is necessary for EMD to calculate the PTCs equivalent to those by NEMD. Unlike the existing theories, our formulation is based on the local equilibrium hypothesis to describe a clear connection between EMD and NEMD simulations. The equivalence of the two derivations of PTCs is confirmed by the numerical results for liquid methane and butane. The present establishment of the EMD-NEMD correspondence makes the MD analysis of PTCs a robust way to clarify the microscopic origins of thermal conductivity.