Two-molecule theory of polyethylene liquids

Two-molecule theory refers to a class of microscopic, self-consistent field theories for the radial distribution function in classical molecular liquids. The version examined here can be considered as one of the very few formally derived closures to the reference interaction site model (RISM) equation. The theory is applied to polyethylene liquids, computing their equilibrium structural and thermodynamic properties at melt densities. The equation for the radial distribution function, which is represented as an average over the accessible states of two molecules in an external field that mimics the effects of the other molecules in the liquid, is computed by Monte Carlo simulation along with the intramolecular structure function. An improved direct sampling algorithm is utilized to speed the equilibration. Polyethylene chains of 24 and 66 united atom CH2 units are studied. The results are compared to full, many-chain molecular dynamics (MD) simulations and self-consistent polymer-RISM (PRISM) theory with the atomic Percus-Yevick (PY) closure under the same conditions. It is shown that the two-molecule theory produces results that are close to those of MD and is thus able to overcome defects of PRISM-PY theory and predict more accurate liquid structure at both short and long ranges. Predictions for the equation of state are also discussed.