Interfacial properties of hydrocarbon/water systems predicted by molecular dynamic simulations

J Chem Phys. 2019 Mar 21;150(11):114703. doi: 10.1063/1.5078739.

Abstract

The presence of small hydrocarbons is known to reduce the interfacial tension of the gas-water interface, and this phenomenon can affect the formation of the clathrate hydrates of these gases. In this work, the interfacial behavior of the pure methane-, ethane-, and propane-water, and the ternary 90:7:3 mol. % gas mixture of (methane + ethane + propane)-water were studied with molecular dynamics simulations. The interfacial tension, γ, and z-density profiles for the gases and water from simulations of the gas-water systems were determined at the temperatures of 275.15 and 298.15 K, and pressures up to 10 MPa for methane and up to near the experimental saturation pressures of ethane and propane. The goal is to accurately calculate the interfacial tension for the hydrocarbon/water systems and to analyze the molecular behaviors at the interfaces which lead to the observed trends. At the same hydrostatic gas phase pressure, propane, ethane, and methane reduce the gas-water interfacial tension in that order. The local density of the gas molecules at the interface is enhanced relative to the bulk gas, and it was determined that about 13%-20%, 33%-40%, and 54%-59% of the gas molecules in the simulation congregated at the interfaces for the CH4-, C2H6-, and C3H8-water systems, respectively, at the different simulated hydrostatic pressure ranges. For all gases in the pressure range studied, a complete monolayer of gas had not formed at the water interface. Furthermore, a dynamic equilibrium with fast exchange between molecules at the interface and in the gas phase was observed. For the gas mixture, deviations were observed between total calculated interfacial tension, γmix, and the "ideal mixture" value, ∑xiγi,pure, calculated from the interfacial tensions of the pure gases, where xi is the mole fraction of each substance in the simulation. Some possible implications of the results on the mechanism of clathrate hydrate formation are discussed.