Vapor–Liquid Equilibrium Simulations of Hydrocarbons Using Molecular Dynamics with Long-Range Lennard-Jones Interactions

Abstract

Because of exposure of injected fuel to extreme temperatures and pressures in modern diesel and jet engines, there is much interest in understanding vapor−liquid equilibrium (VLE), interfacial behavior, and critical properties of hydrocarbon mixtures. Molecular dynamics (MD) simulations with long-range dispersion interactions can be used to study liquid−vapor interfacial properties. In this work, the accuracy of three force fields (CHARMM, OPLS-AA, and TraPPE-UA) is assessed by using MD simulations with the Lennard-Jones particle-mesh Ewald method to examine the VLE behavior of 12 pure hydrocarbons. While the TraPPE-UA potential yields significantly better critical temperature predictions, all three force fields predict critical densities and pressures with about the same accuracy. Surface tensions calculated using the TraPPE-UA force field are too high compared to experiment; values using the CHARMM and OPLS-AA potentials tend to be more accurate at lower temperatures and then decrease too quickly as temperature increases. Simulations of binary mixtures show that both the CHARMM and TraPPE-UA potentials predict the same liquid-and vapor-phase compositions as a function of overall mole fraction. Not surprisingly, the composition of the interfacial region is correlated with the surface tension values of the pure components. Pressure versus mole fraction phase diagrams were quantitatively reproduced by TraPPE-UA and only qualitatively reproduced by the CHARMM force field. Given the lack of phase equilibrium data used in the parameterization of the CHARMM potential, its ability to reproduce many aspects of VLE behavior is impressive. This work demonstrates the power of MD simulations for examining the phase behavior of hydrocarbons under conditions that are difficult to achieve experimentally.

Correction to “Molecular Dynamics Simulations of the Vapor–Liquid Equilibria in CO2/n-Pentane, Propane/n-Pentane, and Propane/n-Hexane Binary Mixtures”

“…Our conclusions that the comparison is affected by the computed pressure are valid, but the MD overestimates the pressure, rather than underestimating it. The tendency of the CHARMM force field to give high vapor pressures has been noted previously …”

Relative Resolution: A Computationally Efficient Implementation in LAMMPS

“…32 The LJ potential is frequently used by the LAMMPS community in various applications, [33][34][35][36][37][38] especially in studying elementary alkanes, as well as oligomeric and polymeric hydrocarbons. [39][40][41][42][43][44] The main benefit of LAMMPS for RelRes is that it inherently supports multiple neighbor lists (based on various cutting distances), and thus, the computational efficiency for RelRes is successfully achieved. Ultimately, we exemplify our LAMMPS algorithm on several alkane liquids (propane, neopentane, alternate cooligomer, block copolymer, etc.).…”

Solubility of Nitrogen in Methane, Ethane, and Mixtures of Methane and Ethane at Titan-Like Conditions: A Molecular Dynamics Study

“…Standard methods for calculating solubility in the limit of infinite dilution is the Widom insertion method . Because the solubility of nitrogen is very high in both methane and ethane, we have performed molecular dynamics vapor–liquid equilibrium (VLE) simulations − of binary mixtures of methane and nitrogen and ethane and nitrogen as well as ternary mixtures of methane, ethane, and nitrogen to investigate the solubility of nitrogen in methane, ethane, and mixtures of methane and ethane at Titan-like conditions. The paper is organized as follows: in the Methods, we describe our simulation approaches for all of the thermodynamic conditions studied here.…”