A Replica Exchange Molecular Dynamics Simulation of a Single Polyethylene Chain: Temperature Dependence of Structural Properties and Chain Conformational Study at the Equilibrium Melting Temperature

The conformational properties of a finite length polyethylene chain were explored over a wide range of temperatures using a replica exchange molecular dynamics simulation providing high quality simulation data representative for the equilibrium behavior of the chain molecule. The radial distribution function (RDF) and the structure factor S(q) of the chain as a function of temperature are analyzed in detail. The different characteristic peaks in the RDF and S(q) were assigned to specific distances in the chain and structural changes occurring with the temperature. In S(q), a peak characteristic for the order in the solid state was found and used to determine the equilibrium melting temperature. A detailed scaling analysis of the structure factor covering the full q range was performed according to the work of Hammouda. In the Theta region, a quantitative analysis of the full structure factor was done using the equivalent Kuhn chain, which enabled us to assign the Theta region of our chain and to demonstrate, in our particular case, the failure of the Gaussian chain approach. The chain conformational properties at the equilibrium melting temperature are discussed using conformational distribution functions, using the largest principal component of the radius of gyration and shape parameters as order parameters. We demonstrate that for the system studied here, the Landau free energy expression based on this conformational distribution information leads to erroneous conclusions concerning the thermodynamic transition behavior. Finally, we focus on the instantaneous conformational properties at the equilibrium melting temperature and give a detailed analysis of the conformational shapes using different shape parameters and a simulation snapshot. We show that the chain does not only take the lamellar rod-like and globular conformational shapes, typical of the solid and liquid states, but can also explore many other conformational states, including the toroidal conformational state. It is the first demonstration that a flexible molecule like PE can also take a toroidal conformational state, which is normally linked to stiffer chains.