Equation of state modeling and force field-based molecular dynamics simulations of supercritical polyethylene +...

The understanding of polymer solution thermodynamics and characterization of pressure effects on fundamental polymer physics of macromolecular systems is significant in the manufacturing of polyolefins. Consequently, numerous experimental and theoretical efforts have been made towards understanding phase behavior of polymer solutions at elevated pressures. Despite this progress, only limited efforts are directed towards understanding the underlying phenomena behind the influence of high pressure upon the thermophysical properties of ternary polymer solutions at a molecular level. The present paper, therefore, reports on the influence of supercritical ethylene on the density of PE + hydrocarbon solvent system by exploring ternary mixtures of PE + hexane + ethylene for ethylene concentrations up to 10 wt% at varied temperatures and in a pressure range from 100 to 1000 bar via fully-atomistic molecular dynamics (MD) simulations. Additionally, the modified Sanchez-Lacombe equation of state (EOS) model is iteratively solved to capture the pressure, concentration, and temperature dependence of ternary PE solution density. It is shown that the small amounts of ethylene dissolved in the liquid mixtures of PE + hexane significantly decreases the polymer solution density. The presence of unreacted monomer in the solution polymerization process utilized in PE manufacturing was found to substantially lower the PE solution density particularly at the lower end of the investigated pressure range. This noteworthy reduction in mixture density as a consequence impacts design and operation of the liquid-liquid phase separator in manufacturing of PE via solution polymerization. Another key point to bear in mind is that the mixture density exhibits fairly less sensitivity to ethylene amount as external pressure raises. Nevertheless, pressure, solvent composition, and temperature dependence of density display less sensitivity as pressure increases. In relation to the characterization of the impact of addition of ethylene an atomistic-level insight is provided, which proves to be of great value in revealing intermolecular interactions in the binary subsystems of polymer/solvent/monomer. The MD computations are shown to be in excellent agreement with the theoretical EOS model, confirming the validity of the proposed methodology. Furthermore, the adopted OPLS-AA has been found a reliable atomistic force field, which provides detailed molecular information on the thermophysical properties of polyolefin in hydrocarbon solutions. Ultimately, it is demonstrated that the MD simulations complement parametric EOS predictions and costly experimental approaches.