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

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.

Atom Types Independent Molecular Mechanics Method for Predicting the Conformational Energy of Small Molecules.

We previously implemented a well-known qualitative chemical principle into an accurate quantitative model computing relative potential energies of conformers. According to this principle, hyperconjugation strength correlates with electronegativity of donors and acceptors. While this earlier version of our model applies to σ bonds, lone pairs, disregarded in this earlier version, also have a major impact on the conformational preferences of molecules. Among the well-established principles used by organic chemists to rationalize some organic chemical behaviors are the anomeric effect, the alpha effect, basicity, and nucleophilicity. These effects are directly related to the presence of lone pairs. We report herein our effort to incorporate lone pairs into our model to extend its applicability domain to any saturated small molecules. The developed model H-TEQ 2 has been validated on a wide variety of molecules from polyaromatic molecules to carbohydrates and molecules with high heteroatoms/carbon ratios. Interestingly, this method, in contrast to common force field-based methods, does not rely on atom types and is virtually applicable to any organic molecules.

Elucidating Hyperconjugation from Electronegativity to Predict Drug Conformational Energy in a High Throughput Manner.

Computational chemists use structure-based drug design and molecular dynamics of drug/protein complexes which require an accurate description of the conformational space of drugs. Organic chemists use qualitative chemical principles such as the effect of electronegativity on hyperconjugation, the impact of steric clashes on stereochemical outcome of reactions, and the consequence of resonance on the shape of molecules to rationalize experimental observations. While computational chemists speak about electron densities and molecular orbitals, organic chemists speak about partial charges and localized molecular orbitals. Attempts to reconcile these two parallel approaches such as programs for natural bond orbitals and intrinsic atomic orbitals computing Lewis structures-like orbitals and reaction mechanism have appeared. In the past, we have shown that encoding and quantifying chemistry knowledge and qualitative principles can lead to predictive methods. In the same vein, we thought to understand the conformational behaviors of molecules and to encode this knowledge back into a molecular mechanics tool computing conformational potential energy and to develop an alternative to atom types and training of force fields on large sets of molecules. Herein, we describe a conceptually new approach to model torsion energies based on fundamental chemistry principles. To demonstrate our approach, torsional energy parameters were derived on-the-fly from atomic properties. When the torsional energy terms implemented in GAFF, Parm@Frosst, and MMFF94 were substituted by our method, the accuracy of these force fields to reproduce MP2-derived torsional energy profiles and their transferability to a variety of functional groups and drug fragments were overall improved. In addition, our method did not rely on atom types and consequently did not suffer from poor automated atom type assignments.