Synergistic effects between a non-ionic and an anionic surfactant on the micellization process and the adsorption at liquid/air surfaces.

Predicting the behaviour of solutions with surfactants of significantly different critical micelle concentration (CMC) values remains a challenge. The study of the molecular interactions within micelles and interfaces in surfactant combinations used in everyday products is essential to understand these complex systems. In this work, the equilibrium and dynamic surface tension in the presence of mixed non-ionic (tristyrylphenol ethoxylates) and anionic (sodium benzene sulfonate with alkyl chain lengths of C10-C13) surfactants, commonly encountered as delivery systems in agrochemicals, were studied and their CMC values were determined. For the surfactant mixtures, four molar ratios were examined: nEOT/nNaDDBS = 0.01, 0.1, 1, 4 and two different cases were analysed, the premixed and the add one by one surfactant. The surface tension for single surfactants stabilised quickly, while the mixtures needed a long time to reach equilibrium; up to 15 h for the premixed mixtures and 40 min when surfactants were added one by one. The CMC values for the nEOT/nNaDDBS = 0.01, 0.1 premixed surfactant mixtures were found to be in between the CMC values of the single surfactants, but those for the nEOT/nNaDDBS = 1 and 4 mixtures were lower than the CMCs of both single surfactants. Calculations based on the regular solution theory suggested that there are attractive forces in the mixed micelles and at the interface layers, while the supramolecular assemblies in the bulk (i.e., micelles) and at interfaces (surfactant films) are preferentially enriched in EOT.

Water Wettability Coupled with Film Growth on Realistic Cyclopentane Hydrate Surfaces.

Although the wettability of hydrate surfaces and hydrate film growth are key to understanding hydrate agglomeration and pipeline plugging, a quantitative understanding of the coupled behavior between both phenomena is lacking. In situ measurements of wettability coupled with film growth were performed for cyclopentane hydrate surfaces in cyclopentane at atmospheric pressure and temperatures between 1.5-6.8 °C. Results were obtained as a function of annealing (conversion) time and subcooling. Hydrate surface wettability decreased as annealing time increased, while hydrate film growth rate was unaffected by annealing time at any subcooling. The results are interpreted as a manifestation of the hydrate surface porosity, which depends on annealing time and controls water spreading on the hydrate surface. The wettability generally decreased as the subcooling increased because higher subcooling yields rougher hydrate surfaces, making it harder for water to spread. However, this effect is balanced by hydrate growth rates, which increase with subcooling. Also affecting the results, surface heating from heat release (from exothermic crystallization) allows excess surface water to promote spreading. The hydrate film growth rate on water droplets increased with subcooling, as expected from a higher driving force. At any subcooling, the instantaneous hydrate growth rate decreased over time, likely from heat transfer limitations. A new phenomenon was observed, where the angle at the three-phase point increases from the initial contact angle upon hydrate film growth, named the crystallization angle. This is attributed to the water droplet trying to spread while the thin film is weak enough to be redirected. Once the hydrate film grows and forms a "wall" around the droplet, it cannot be moved, and further growth yields a crater on the droplet surface, attributed to water penetrating the hydrate surface pore structures. This fundamental behavior has many flow assurance implications since it affects the interactions between the agglomerating hydrate particles and water droplets.

Structural and dynamic properties of some aqueous salt solutions.

Aqueous salt solutions are utilized and encountered in wide-ranging technological applications and natural settings. Towards improved understanding of the effect of salts on the dynamic properties of such systems, dilute aqueous salt solutions (up to 1 molar concentration) are investigated here, via experiments and molecular simulations. Four salts are considered: sodium chloride, for which published results are readily available for comparison, ammonium acetate, barium acetate and barium nitrate, for which published data are scarce. In the present work, molecular dynamics (MD) simulations are conducted to quantify viscosity and water self-diffusion coefficients, together with rheometry and Pulsed Field Gradient Spin Echo (PFGSE)-NMR experiments for validation. Simulation predictions are consistent with experimental observations in terms of trend and magnitude of salt-specific effects. Combining insights from the approaches considered, an interpretation of the results is proposed whereby the capacity of salts to influence bulk dynamics arises from their molecular interfacial area and strength of interaction with first hydration-shell water molecules. For the concentration range investigated, the interpretation could be useful in formulating aqueous systems for applications including the manufacturing of advanced catalysts.

Effects of surface contamination on the interfacial properties of CO2/water/calcite systems.

Understanding the wetting properties of reservoir rocks can be of great benefit for advanced applications such as the effective trapping and geological storage of CO2. Despite their importance, not all mechanisms responsible for wetting mineral surfaces in subsurface environments are well understood. Factors such as temperature, pressure and salinity are often studied, achieving results with little unanimity; other possible factors are left somewhat unexplored. One such factor is the effect of contamination. In the present study, the effects of adding a non-aqueous organic contaminant, ethanol, on the CO2-water interfacial tension (IFT) and the CO2/water/calcite contact angle were investigated using molecular dynamics simulations. Within the conditions studied, relatively small amounts of ethanol cause a significant decrease in the CO2-water IFTs, as well as a pronounced increase in the water-calcite-CO2 three phase contact angle. The latter result is due to the decrease of the IFT between CO2 and water and the strong adsorption of ethanol on the solid substrate. These findings could be helpful for explaining how impurities can affect experimental data and could lead to effective carbon sequestration strategies.

Water/oil interfacial tension reduction - an interfacial entropy driven process.

The interfacial tension (IFT) of a fluid-fluid interface plays an important role in a wide range of applications and processes. When low IFT is desired, surface active compounds (e.g. surfactants) can be added to the system. Numerous attempts have been made to relate changes in IFT arising from such compounds to the specific nature of the interface. However, the IFT is controlled by an interplay of factors such as temperature and molecular structure of surface-active compounds, which make it difficult to predict IFT as those conditions change. In this study, we present the results from molecular dynamics simulations revealing the specific role surfactants play in IFT. We find that, in addition to reducing direct contact between the two fluids, surfactants serve to increase the disorder at the interface (related to interfacial entropy) and consequently reduce the water/oil IFT, especially when surfactants are present at high surface density. Our results suggest that surfactants that yield more disordered interfacial films (e.g. with flexible and/or unsaturated tails) reduce the water/oil IFT more effectively than surfactants which yield highly ordered interfacial films. Our results shed light on some of the factors that control IFT and could have important practical implications in industrial applications such as the design of cosmetics, food products, and detergents.

Aqueous films on pore surfaces mediate adsorption and transport of gases through crowded nanopores.

Interactions of trapped reservoir gases within organic-rich and brine-bearing sedimentary rocks have direct relevance to many geoenergy applications. Extracting generalizable information from experimental campaigns is hindered by the fact that geological systems are extremely complex. However, modern computational tools offer the opportunity of studying systems with controlled complexity, in an effort to better understand the mechanisms at play. Employing molecular dynamics, we examine here adsorption and transport of gases containing CH4 and either CO2 or H2S within amorphous silica nanopores filled with benzene. We explicitly quantify the effect of small amounts of water/brines at geological temperature and pressure conditions. Because of wetting, the presence of brines lessens the adsorption capacity of the aromatic-filled pore. The simulation results show salt-specific effects on the transport properties of the gases when either KCl or CaCl2 brines are considered, although adsorption was not affected. The acid gases considered either facilitate or hinder CH4 transport depending on whether they are more or less preferentially adsorbed within the pore as compared to benzene, and this effect is mediated by the presence of water/brines. Our simulation results could be used to extract thermodynamic quantities that in the future will help to optimize transport of various gases through organic-rich and brine-bearing sedimentary rocks, which is likely to have a positive impact on both hydrocarbon production and carbon sequestration applications. As a first step, a phenomenological model is presented here, which allows one to predict permeability based on interatomic energies.

Correction: Effects of droplet size and surfactants on anchoring in liquid crystal nanodroplets.

Correction for 'Effects of droplet size and surfactants on anchoring in liquid crystal nanodroplets' by Zeynep Sumer et al., Soft Matter, 2019, 15, 3914-3922, DOI: .

Self-assembly of mono- and poly-dispersed nanoparticles on emulsion droplets: antagonistic vs. synergistic effects as a function of particle size.

In this work, using Dissipative Particle Dynamics simulations, we provide fundamental insights into the self-assembly of nanoparticles (NPs) on droplet surfaces in an oil-in-water emulsion. We highlight the effect of particle size on the arrangement of NPs for different interparticle interactions. NPs of two different sizes were considered. In general, when the NP-NP interaction is changed from repulsive to attractive, a transition in the NP arrangement occurs from weekly-connected networks to clusters of NPs separated by particle-free domains. When NP-NP interactions are strongly attractive, NPs yield small 3D aggregates on the droplet surface. These arrangements seem to agree with experimental observations reported in the literature. In addition, our simulations suggest that small NPs are able to diffuse more easily on the droplet surface, which leads to prompt self-organisation, while large NPs are more likely to form metastable structures, perhaps because of slow mobility and strong adsorption to the interface. Our analysis suggests that thermal fluctuations could provide the activation energy for the small NPs to escape local minima in the free energy landscape. The results obtained for systems containing NPs of two sizes provide evidence of size segregation on the droplet surface, which could be useful when NP self-assemblies are used, for example, to template supra-molecular materials. However, analysis of the simulated trajectories suggests that the results depend strongly on the initial configuration, as the larger NPs seem to impose barriers for the small NPs to adsorb and diffuse on the droplet surface.

Effects of droplet size and surfactants on anchoring in liquid crystal nanodroplets.

Liquid crystal (LC) droplets attract scientific attention for many advanced applications, including, but not limited to optical and sensing devices. To aid experimental advancements, theoretical calculations have been conducted to quantify molecular driving forces responsible for the collective behaviour of LC molecules within micrometer-size spherical droplets. To quantify the LC molecular anchoring within spherical physical constraints, molecular simulations at atomistic resolution would be useful. In an attempt to bridge the gap between computational capabilities and experimental interest, coarse-grained simulations are used here to study nematic LC nanodroplets dispersed in water. A LC phase diagram is generated as a function of droplet size and temperature. The effect of adding surfactants on LC anchoring was quantified, considering surfactants of different molecular features. When few surfactants are present, they self-assemble at the droplet boojums regardless of their molecular features. All surfactants tested shifted LC orientation from bipolar to uniaxial. When the surfactants have a hydrophobic tail of sufficient length, they cause deviations from the spherical symmetry of LC droplets. Increasing the concentration of these surfactants enhances such phenomenon. Simulations were also conducted to assess the ability of the surfactants to prevent the agglomeration between two LC droplets. The results showed that coalescence was inevitable at all conditions and suggested that large enough surfactant concentrations can delay the phenomenon. The results presented could be helpful for designing novel surface-active compounds to develop optical and/or sensing devices at conditions in which mutual solubility between water and LCs is low.

Effects of water on the stochastic motions of propane confined in MCM-41-S pores.

Hydrocarbons confined in porous media find applications in a wide variety of industries and therefore their diffusive behavior is widely studied. Most of the porous media found in natural environments are laden with water, which might affect the confined hydrocarbons. To quantify the effect of hydration, we report here a combined quasielastic neutron scattering (QENS) and molecular dynamics (MD) simulation study on the dynamics of propane confined in the 1.5 nm-wide micropores of MCM-41-S in the presence of water at 230 and 250 K. To eliminate the strong incoherent signal from water and emphasize the propane signal we have used heavy water (D2O). QENS data show two dynamically different populations of propane in MCM-41-S and suggest that the presence of water hinders the diffusion of propane. Weak elastic contributions to the QENS spectra suggest that only long-range translational motion of propane molecules contributes to the quasielastic broadening. MD simulations carried out using a model cylindrical silica pore of 1.6 nm diameter filled with water and propane agree with the experimental finding of water hindering the diffusion of propane. Further, the simulation results suggest that the slowing down of propane motions is a function of the water content within the pore and is stronger at higher water contents. At high water content, the structure and the dynamics, both translational and rotational, of propane are severely impacted. Simulation data suggest that the rotational motion of the propane molecule occurs on time scales much faster than those accessible with the QENS instrument used, and thus explain the weak elastic contribution to the QENS spectra measured in the experiments. This study shows the effects of hydration on the structure and dynamics of volatiles in porous media, which are of interest for fundamental understanding and applied studies of confined fluids.

Atomic Force Microscopy Force Mapping Analysis of an Adsorbed Surfactant above and below the Critical Micelle Concentration.

Force curves collected using an atomic force microscope (AFM) in the presence of adsorbed surfactants are often used to draw conclusions about adsorbed film packing, rigidity, and thickness. However, some noteworthy features of such force curve characteristics have yet to be thoroughly investigated and explained. In this work, we collected force curves from tetradecyltrimethylammonium bromide films adsorbed on highly oriented pyrolytic graphite (HOPG), silica, and silica that had been hydrophobized by functionalization with dichlorodimethyl silane. Breakthrough events in the force curves from several different trials were compared to show that the breakthrough distance, often reported as the adsorbed film thickness, increased with concentration below the critical micelle concentration (CMC) but was approximately 3.5 nm on all surfaces between 2× and 10× CMC; an unexpected result because of the different surface chemistries for the three surfaces. We employed an AFM probe with a different force constant ( k) value as well as a colloidal probe and the breakthrough distance remained approximately 3.5 nm in all cases. Gradient mapping, a variant of force mapping, was also implemented on the three surfaces and resulted in a new technique for visualizing adsorbed surfactant in situ. The resulting maps showed patches of adsorbed surfactant below the CMC and revealed that with increasing concentration, the size of the patches increased resulting in full coverage near and above the CMC. These results are, to our knowledge, the first time force mapping has been used to spatially track patches of adsorbed surfactant. Finally, layers of surfactants on an AFM tip were investigated by collecting a force map on a single AFM tip using the tip of a separate AFM probe. A breakthrough event was observed between the tips, indicating that a layer of surfactant was present on at least one, if not both tips.

Emergent Properties of Antiagglomerant Films Control Methane Transport: Implications for Hydrate Management.

The relationship between collective properties and performance of antiagglomerants (AAs) used in hydrate management is handled using molecular dynamics simulations and enhanced sampling techniques. A thin film of AAs adsorbed at the interface between one flat sII methane hydrate substrate and a fluid hydrocarbon mixture containing methane and n-dodecane is studied. The AA considered is a surface-active compound with a complex hydrophilic head that contains both amide and tertiary ammonium cation groups and hydrophobic tails. At a sufficiently high AA density, the interplay between the surfactant layer and the liquid hydrocarbon excludes methane from the interfacial region. In this scenario, we combine metadynamics and umbrella sampling frameworks to study accurately the free-energy landscape and the equilibrium rates associated with the transport of one methane molecule across the AA film. We observe that the local configurational changes of the liquid hydrocarbon packed within the AA film are associated with high free-energy barriers for methane transport. The time scales estimated for the transport of methane across the AA film can be, in some cases, comparable to those reported in the literature for the growth of hydrates, suggesting that one possible mechanism by which AAs delay the formation of hydrate plugs could be providing a barrier to methane transport. Considering the interplay between the structural design and collective properties of AAs might be of relevance to improve their performance in flow assurance.

Manipulating molecular order in nematic liquid crystal capillary bridges via surfactant adsorption: guiding principles from dissipative particle dynamics simulations.

The ability of liquid crystals (LCs) to change orientational order is used in applications, ranging from sensors to displays. The aim of this work is to computationally investigate how surfactant adsorption on cylindrical LC bridges can be used to control such orientational order. Building from classical fundamental lessons, understanding the ordering of mesogens along a preferred axis with the help of molecular modelling contributes to investigations of systems that could be a platform for LC-based sensing applications. The coarse-grained dissipative particle dynamics (DPD) simulation method is implemented here, because it allows us to quantify the effect of molecular features on the properties of meso-scopic systems containing LC bridges, an aqueous solvent, and surfactants at various concentrations. Three surfactant types are modelled with short, medium, and long tail lengths, respectively. All surfactants adsorb at the LC-water interface. It is found that the length of the surfactant hydrophobic tail determines the effectiveness by which the LC order is affected. Short tails are not as effective as long ones. Surfactants with long tails affect the LC order, but, in agreement with experiments, predominantly only within a short distance from the LC-water interface. For these surfactants, the surface density at the LC-water interface is an important knob that can be used to control the order of the LCs. As the effective LC-surfactant interactions change, so does the distribution of the surfactants at the interface. Consistent with theoretical expectations, the results presented here elucidate the effect of molecular features on the anchoring mechanism between surfactants and mesogens within cylindrical bridges dispersed in aqueous systems and could be helpful for designing novel surface-active compounds in the development of advanced sensing devices based on LCs.

Surfactant Assemblies on Selected Nanostructured Surfaces: Evidence, Driving Forces, and Applications.

Surfactant adsorption at solid-liquid interfaces is critical for a number of applications of vast industrial interest and can also be used to seed surface-modification processes. Many of the surfaces of interest are nanostructured, as they might present surface roughness at the molecular scale, chemical heterogeneity, as well as a combination of both surface roughness and chemical heterogeneity. These effects provide lateral confinement on the surfactant aggregates. It is of interest to quantify how much surfactant adsorbs on such nanostructured surfaces and how the surfactant aggregates vary as the degree of lateral confinement changes. This review focuses on experimental evidence on selected substrates, including gold- and carbon-based substrates, suggesting that lateral confinement can have pronounced effects both on the amount adsorbed and on the morphology of the aggregates as well as on a systematic study, via diverse simulation approaches, on the effect of lateral confinement on the structure of the surfactant aggregates. Atomistic and coarse-grained simulations conducted for surfactants on graphene sheets and carbon nanotubes are reviewed, as well as coarse-grained simulations for surfactant adsorption on nanostructured surfaces. Finally, we suggest a few possible extensions of these studies that could positively impact a few practical applications. In particular, the simultaneous effect of lateral confinement and of the coadsorption of molecular compounds within the surface aggregates is expected to yield interesting fundamental results with long-lasting consequences in applications ranging from drug delivery to the design of advanced materials.

Propane-Water Mixtures Confined within Cylindrical Silica Nanopores: Structural and Dynamical Properties Probed by Molecular Dynamics.

Despite the multiple length and time scales over which fluid-mineral interactions occur, interfacial phenomena control the exchange of matter and impact the nature of multiphase flow, as well as the reactivity of C-O-H fluids in geologic systems. In general, the properties of confined fluids, and their influence on porous geologic phenomena are much less well understood compared to those of bulk fluids. We used equilibrium molecular dynamics simulations to study fluid systems composed of propane and water, at different compositions, confined within cylindrical pores of diameter ∼16 Šcarved out of amorphous silica. The simulations are conducted within a single cylindrical pore. In the simulated system all the dangling silicon and oxygen atoms were saturated with hydroxyl groups and hydrogen atoms, respectively, yielding a total surface density of 3.8 -OH/nm2. Simulations were performed at 300 K, at different bulk propane pressures, and varying the composition of the system. The structure of the confined fluids was quantified in terms of the molecular distribution of the various molecules within the pore as well as their orientation. This allowed us to quantify the hydrogen bond network and to observe the segregation of propane near the pore center. Transport properties were quantified in terms of the mean square displacement in the direction parallel to the pore axis, which allows us to extract self-diffusion coefficients. The diffusivity of propane in the cylindrical pore was found to depend on pressure, as well as on the amount of water present. It was found that the propane self-diffusion coefficient decreases with increasing water loading because of the formation of water bridges across the silica pores, at sufficiently high water content, which hinder propane transport. The rotational diffusion, the lifespan of hydrogen bonds, and the residence time of water molecules at contact with the silica substrate were quantified from the simulated trajectories using the appropriate autocorrelation functions. The simulations contribute to a better understanding of the molecular phenomena relevant to the behavior of fluids in the subsurface.

Evidence of Structure-Performance Relation for Surfactants Used as Antiagglomerants for Hydrate Management.

Molecular dynamics simulations were employed to study the structure of molecularly thin films of antiagglomerants adsorbed at the interface between sII methane hydrates and a liquid hydrocarbon. The liquid hydrocarbon was composed of dissolved methane and higher-molecular-weight alkane such as n-hexane, n-octane, and n-dodecane. The antiagglomerants considered were surface-active compounds with three hydrophobic tails and a complex hydrophilic head that contains both amide and tertiary ammonium cation groups. The length of the hydrophobic tails and the surface density of the compounds were changed systematically. The results were analyzed in terms of the preferential orientation of the antiagglomerants, density distributions of various molecular compounds, and other molecular-level properties. At low surface densities, the hydrophobic tails do not show preferred orientation, irrespectively of the tail length. At sufficiently high surface densities, our simulations show pronounced differences in the structure of the interfacial film depending on the molecular features and on the type of hydrocarbons present in the system. Some antiagglomerants are found to pack densely at the interface and exclude methane from the interfacial region. Under these conditions, the antiagglomerant film resembles a frozen interface. The hydrophobic tails of the antiagglomerants that show this feature has a length comparable to that of the n-dodecane in the liquid phase. It is possible that the structured interfacial layer is in part responsible for determining the performance of antiagglomerants in flow-assurance applications. The simulation results are compared against experimental data obtained with the rocking cell apparatus. It was found that the antiagglomerants for which our simulations suggest evidence of a frozen interface at sufficiently high surface densities are those that show better performance in rocking cell experiments.

Surfactants adsorption on crossing stripes and steps.

Using coarse-grained dissipative particle dynamics (DPD) simulations, we systematically study the effect of surface heterogeneity on surfactant adsorption. Here we investigate the adsorption and aggregation of surfactants on hydrophobic stripes crossing each other perpendicularly (i.e., crossing stripes) and on hydrophobic steps. The results are compared with those obtained for isolated stripes. We find that on crossing stripes of moderate stripe widths (e.g., L = 0.61LS, 1.22LS and 1.83LS, where LS is the length of one surfactant molecule) the crossing region hinders the formation of defect-free adsorbed surfactant structures. By increasing the stripe width and/or by increasing the length of one of the two perpendicularly crossing stripes (i.e., lowering the surface density of defects/intersections), the crossing region is found to have a weaker effect on the features of the adsorbed structures. Regarding surfactant adsorption on steps, our simulation results show that the self-assembled aggregates can be stretched along the step corner, and the resultant elastic deformation can hinder adsorption. This qualitative observation can facilitate a description of surfactant adsorption that takes into consideration also the deformation of the self-assembled film. As suggested by such a general model, increasing the convex angle of the step, increasing the size of the surfactant head groups, and changing other physical parameters can reduce the elastic energy penalty, and yield larger amounts of surfactants adsorbed. The results presented could assist in understanding and sometimes predicting surfactant adsorption on heterogeneous surfaces, suggest methods to formulate surfactant mixtures to control surface coverage on heterogeneous surfaces, and perhaps facilitate new methods for the fabrication of nano-structured surfaces.

Molecular dynamics simulations of propane in slit shaped silica nano-pores: direct comparison with quasielastic neutron scattering experiments.

Molecular motion under confinement has important implications for a variety of applications including gas recovery and catalysis. Propane confined in mesoporous silica aerogel as studied using quasielastic neutron scattering (QENS) showed anomalous pressure dependence in its diffusion coefficient (J. Phys. Chem. C, 2015, 119, 18188). Molecular dynamics (MD) simulations are often employed to complement the information obtained from QENS experiments. Here, we report an MD simulation study to probe the anomalous pressure dependence of propane diffusion in silica aerogel. Comparison is attempted based on the self-diffusion coefficients and on the time scales of the decay of the simulated intermediate scattering functions. While the self-diffusion coefficients obtained from the simulated mean squared displacement profiles do not exhibit the anomalous pressure dependence observed in the experiments, the time scales of the decay of the intermediate scattering functions calculated from the simulation data match the corresponding quantities obtained in the QENS experiment and thus confirm the anomalous pressure dependence of the diffusion coefficient. The origin of the anomaly in pressure dependence lies in the presence of an adsorbed layer of propane molecules that seems to dominate the confined propane dynamics at low pressure, thereby lowering the diffusion coefficient. Further, time scales for rotational motion obtained from the simulations explain the absence of rotational contribution to the QENS spectra in the experiments. In particular, the rotational motion of the simulated propane molecules is found to exhibit large angular jumps at lower pressure. The present MD simulation work thus reveals important new insights into the origin of anomalous pressure dependence of propane diffusivity in silica mesopores and supplements the information obtained experimentally by QENS data.

Numerical analysis of Pickering emulsion stability: insights from ABMD simulations.

The issue of the stability of Pickering emulsions is tackled at a mesoscopic level using dissipative particle dynamics simulations within the Adiabatic Biased Molecular Dynamics framework. We consider the early stage of the coalescence process between two spherical water droplets in a decane solvent. The droplets are stabilized by Janus nanoparticles of different shapes (spherical and ellipsoidal) with different three-phase contact angles. Given a sufficiently dense layer of particles on the droplets, we show that the stabilization mechanism strongly depends on the collision speed. This is consistent with a coalescence mechanism governed by the rheology of the interfacial region. When the system is forced to coalesce sufficiently slowly, we investigate at a mesoscopic level how the ability of the nanoparticles to stabilize Pickering emulsions is discriminated by nanoparticle mobility and the associated caging effect. These properties are both related to the interparticle interaction and the hydrodynamic resistance in the liquid film between the approaching interfaces.

Factors governing the behaviour of aqueous methane in narrow pores.

All-atom equilibrium molecular dynamics simulations were employed to investigate the behaviour of aqueous methane confined in 1-nm-wide pores obtained from different materials. Models for silica, alumina and magnesium oxide were used to construct the slit-shaped pores. The results show that methane solubility in confined water strongly depends on the confining material, with silica yielding the highest solubility in the systems considered here. The molecular structure of confined water differs within the three pores, and density fluctuations reveal that the silica pore is effectively less 'hydrophilic' than the other two pores considered. Comparing the water fluctuation autocorrelation function with local diffusion coefficients of methane across the hydrated pores we observed a direct proportional coupling between methane and water dynamics. These simulation results help to understand the behaviour of gas in water confined within narrow subsurface formations, with possible implications for fluid transport.