Nonequilibrium molecular dynamics calculation of the thermal conductivity of amorphous polyamide-6,6

The thermal conductivity of the amorphous phase of polyamide-6,6 is investigated by nonequilibrium molecular dynamics simulations. Two different algorithms are used, reverse nonequilibrium molecular dynamics and the dual-thermostat method. Particular attention is paid to the force field used. Four different models are tested, flexible and rigid bonds and all-atom and united-atom descriptions. They mainly differ in the number of high-frequency degrees of freedom retained. The calculated thermal conductivity depends systematically on the number of degrees of freedom of the model. This dependence is traced to the quantum nature of the fast molecular vibrations, which are incorrectly described by classical mechanics. The best agreement with experiment is achieved for a united-atom model with all bonds kept rigid. It could be shown that both the thermal conductivity and the heat capacity of a model show a similar (but not equal) dependence on its number of degrees of freedom. Hence, the computationally more convenient heat capacity can be used for force field optimization. A side result is the anisotropy of the thermal conductivity in stretched polyamide; heat conduction is faster parallel to the drawing direction than perpendicular to it.