A new molecular dynamics simulation method, with coupling to external baths, is used to perform equilibrium simulations on polyamide-6,6 trimers nanoconfined between graphene surfaces, in equilibrium with the bulk polymer. The method is coupled with the reverse nonequilibrium molecular dynamics simulation technique to exchange heat in the direction normal to the surfaces. To be able to study the effect of confinement on the heat conductance in nanoconfined pores, in this work a number of simulations on systems with different pore sizes are done. It is concluded that the coefficient of heat conductivity depends on the degree of polymer layering between the surfaces and on the pore width. Our results further indicate a considerable temperature drop at the interface between the surfaces and polymer. The calculated Kapitza lengths depend on the intersurface distance and on the layering of the polymer nanoconfined between the surfaces.

While polymers are known as thermal insulators, recent studies show that stretched single chains of polymers have a very high thermal conductivity. In this work, our new simulation scheme for simulation of heat flow in nanoconfined fluids [H. Eslami, L. Mohammadzadeh, and N. Mehdipour, J. Chem. Phys. 135, 064703 (2011)] is employed to study the effect of chain ordering (stretching) on the rate of heat transfer in polyamide-6,6 nanoconfined between graphene surfaces. Our results for the heat flow in the parallel direction (the plane of surfaces) show that the coefficient of thermal conductivity depends on the intersurface distance and is much higher than that of the bulk polymer. A comparison of results in this work with our former findings on the heat flow in the perpendicular direction, with the coefficient of heat conductivity less than the bulk sample, reveal that well-organized polymer layers between the confining surfaces show an anisotropic heat conduction; the heat conduction in the direction parallel to the surfaces is much higher than that in the perpendicular direction. The origin of such anisotropy in nanometric heat flow is shown to be the dramatic anisotropy in chain conformations (chain stretching) beside the confining surfaces. The results indicate that the coefficients of heat conductivity in both directions, normal and parallel to the surfaces, depend on the degree of polymer layering between the surfaces and the pore width.

Reverse nonequilibrium molecular dynamics simulations were done to quantify the effect of the inclusion of carbon nanotubes (CNTs) in the Polyamide-6,6 matrix on the enhancement in the thermal conductivity of polymer. Two types of systems were simulated; systems in which polymer chains were in contact with a single CNT, and those in which polymer chains were in contact with four CNTs, linked together via polymer linkers at different linkage fractions. In both cases, heat transfer in both perpendicular and parallel (to the CNT axis) directions were studied. To examine the effect of surface curvature (area) on the heat transfer between CNT and polymer, systems containing CNTs of various diameters were simulated. We found a large interfacial thermal resistance at the CNT-polymer boundary. The interfacial thermal resistance depends on the surface area of the CNT (lower resistances were seen at the interface of flatter CNTs) and is reduced by linking CNTs together via polymer chains, with the magnitude of the reduction depending on the linkage fraction. The thermal conductivity of polymer in the perpendicular direction depends on the surface proximity; it is lower at closer distances to the CNT surface and converges to the bulk value at distances as large as 2 nm. The chains at the interface of CNT conduct heat more in the parallel than in the perpendicular directions. The magnitude of this thermal conductivity anisotropy reduces with decreasing the CNT diameter and increasing the linkage fraction. Finally, microscopic parameters obtained from simulations were used to investigate macroscopic thermal conductivities of polymer nanocomposites within the framework of effective medium approximation.

Grand canonical ensemble molecular dynamics simulation is employed to calculate the solubility of water in polyamide-6,6. It is shown that performing two separate simulations, one in the polymeric phase and one in the gaseous phase, is sufficient to find the phase coexistence point. In this method, the chemical potential of water in the polymer phase is expanded as a first-order Taylor series in terms of pressure. Knowing the chemical potential of water in the polymer phase in terms of pressure, another simulation for water in the gaseous phase, in the grand canonical ensemble, is done in which the target chemical potential is set in terms of pressure in the gas phase. The phase coexistence point can easily be calculated from the results of these two independent simulations. Our calculated sorption isotherms and solubility coefficients of water in polyamide-6,6, over a wide range of temperatures and pressures, agree with experimental data.

No abstract