An Integrated Experimental and Computational Platform to Explore Gas Hydrate Promotion, Inhibition, Rheology, and Mechanical Properties at McGill University: A Review

(1) Background: Gas hydrates are historically notable due to their prevalence and influence on operational difficulties in the oil and gas industry. Recently, new technologies involving the formation of gas hydrates to accomplish various applications have been proposed. This has created new motivation for the characterization of rheological and mechanical properties and the study of molecular phenomena in gas hydrates systems, particularly in the absence of oil and under pre-nucleation conditions. (2) Methodology: This work reviews advances in research on the promotion, inhibition, rheology, and mechanical properties of gas hydrates obtained through an integrated material synthesis-property characterization-multi-scale theoretical and computational platform at McGill University. (3) Discussion: This work highlights the findings from previous experimental work by our group and identifies some of their inherent physical limitations. The role of computational research methods in extending experimental results and observations in the context of mechanical properties of gas hydrates is presented. (4) Summary and Future perspective: Experimental limitations due to the length and time scales of physical phenomena associated with gas hydrates were identified, and future steps implementing the integrated experimental-computational platform to address the limitations presented here were outlined.

Molecular dynamics study of wetting of alkanes on water: from high temperature to the supercooled region and the influence of second inflection points of interfacial tensions

To explore the wetting behavior of alkanes on bulk water interfaces, molecular dynamics (MD) simulations were carried out for united-atom PYS alkane models, and for SPC/E and TIP4P/2005 water models over a wide temperature range. The MD results at each temperature were used to find (1) the surface tension of the alkanes (octane, nonane) and water, and (2) the interfacial tensions of the alkane-water systems. These quantities were then used to calculate the spreading coefficient (S) and contact angle (θc) for each alkane on water. At higher temperatures, the contact angle of octane and nonane on water is found to behave in accord with conventional expectations, i.e., it decreases with increasing temperature for both water models as each system approaches the usual high-temperature transition to perfect wetting. At lower temperatures, we found an unusual temperature dependence of S and θc for each PYS alkane on SPC/E water. In contrast to conventional expectations, θc decreases with a decrease in the temperature. For octane-SPC/E water, this unusual behavior of θc occurs due to the presence of second inflection points (SIP) in the vapor-water and the octane-water interfacial tensions, whereas the SIP effect is much less important for the nonane-water system. The unusual temperature dependence of θc observed for nonane on SPC/E water is also found for nonane on TIP4P/2005 water. On the other hand, such unusual wetting behavior has not been observed in the PYS octane-TIP4P/2005 water system, except possibly for the two lowest temperatures studied.

Molecular dynamics simulations of interfacial properties of the CO2-water and CO2-CH4-water systems

Molecular dynamics simulations were performed to study the interfacial behavior of the pure carbon dioxide-water system and a binary 40:60 mol. % gas mixture of (carbon dioxide + methane)-water at the temperatures of 275.15 K and 298.15 K and pressures near 4 MPa for CO₂ and up to 10 MPa for methane. The simulations are used to study the dynamic equilibrium of the gases at the water-gas interface, to determine the z-density profiles for the gases and water, and calculate the interfacial tension γ under the different temperature/pressure conditions close to those of the formation of clathrate hydrates of these gases. At the same hydrostatic gas phase pressure, the CO₂-water interface has a lower interfacial tension than the CH₄-water interface. A greater number of CO₂ molecules, as much as three times more than methane at the same pressure, were adsorbed at the interfacial layer, which reflects the stronger electrostatic quadrupolar and van der Waals interactions between CO₂ and water molecules at the interface. The water surfaces are covered by less than a monolayer of gas even when the pressure of the system goes near the saturation pressure of CO₂. The surface adsorbed molecules are in dynamic equilibrium with the bulk gas and with exchange between the gas and interface regions occurring repeatedly within the timescale of the simulations. The effects of the changes in the CO₂-water interfacial tension with external temperature and pressure conditions on the formation of the clathrate hydrates and other CO₂ capture and sequestration processes are discussed.

Association of alkanes with the aqueous liquid-vapor interface: a reference system for interpreting hydrophobicity generally through interfacial fluctuations

We report free energy calculations and fluctuation profiles of single alkanes (from methane to pentane) along the direction normal to the air-water interface. The induced fluctuations and the interfacial stabilities of alkanes are found to be correlated and similar to the results of inorganic monovalent ions (Ou et al., J. Phys. Chem. B, 2013, 117, 11732). This suggests that hydrophobic solvation of solutes and ions is important in determining the adsorption behavior.

Bubble bursting as an aerosol generation mechanism during an oil spill in the deep-sea environment: molecular dynamics simulations of oil alkanes and dispersants in atmospheric air/salt water interfaces

Potential of mean force (PMF) calculations and molecular dynamics (MD) simulations were performed to investigate the properties of oil n-alkanes [i.e., n-pentadecane (C15), n-icosane (C20) and n-triacontane (C30)], as well as several surfactant species [i.e., the standard anionic surfactant sodium dodecyl sulfate (SDS), and three model dispersants similar to the Tween and Span species present in Corexit 9500A] at air/salt water interfaces. This study was motivated by the 2010 Deepwater Horizon (DWH) oil spill, and our simulation results show that, from the thermodynamic point of view, the n-alkanes and the model dispersants have a strong preference to remain at the air/salt water interface, as indicated by the presence of deep free energy minima at these interfaces. The free energy minimum of these n-alkanes becomes deeper as their chain length increases, and as the concentration of surfactant species at the interface increases. The n-alkanes tend to adopt a flat orientation and form aggregates at the bare air/salt water interface. When this interface is coated with surfactants, the n-alkanes tend to adopt more tilted orientations with respect to the vector normal to the interface. These simulation results are consistent with the experimental findings reported in the accompanying paper [Ehrenhauser et al., Environ. Sci.: Processes Impacts 2013, in press, (DOI: 10.1039/c3em00390f)]. The fact that these long-chain n-alkanes show a strong thermodynamic preference to remain at the air/salt water interfaces, especially if these interfaces are coated with surfactants, makes these species very likely to adsorb at the surface of bubbles or droplets and be ejected to the atmosphere by sea surface processes such as whitecaps (breaking waves) and bubble bursting. Finally, the experimental finding that more oil hydrocarbons are ejected when Corexit 9500A is present in the system is consistent with the deeper free energy minima observed for the n-alkanes at the air/salt water interface at increasing concentrations of surfactant species.

Recent advances in Many Body Dissipative Particles Dynamics simulations of liquid-vapor interfaces

Many Body Dissipative Particles Dynamics (MDPD) simulation is a novel promising mesoscopic method to model the liquid-vapor interfaces. Based upon works of Paganobarraga and Frenkel (J. Chem. Phys. 15, 5015 (2001)) and Trofimov (J. Chem. Phys. 117, 9383 (2002)) and of Warren (Phys. Rev. E 68, 066702 (2003)) this method has been critically reviewed during this last decade. We propose here to give an overview of the Many Body Dissipative Particles Dynamic simulation within the framework of the liquid-vapor interfaces. We recall the theoretical background of MDPD and we present some recent results of systems of interest such as water liquid-vapor interfaces and salt effect on water surface tension. Additionally we discuss the ability of MDPD to capture the mechanisms at the mesoscopic scale through the formation of micelles and the coalescence of a nanodroplet water on water surface.