Preferential adsorption from liquid water-ethanol mixtures in alumina pores

Adsorption Kinetics of Poly(benzyl acrylate) Chains onto Alumina Interface during the Flow-Driven Translocation through Cylindrical Nanochannels

In this work, the adsorption kinetics of the PBAN/AAO system under flushing condition has been investigated, where PBAN and AAO represent poly(benzyl acrylate) and anodic alumina oxide (AAO, average pore radius R0 ≈ 10 nm) nanochannel, respectively. Our specially designed double-pump flushing system is proved to eliminate the overshoot phenomenon and in situ monitor transmembrane pressure (ΔP) as a function of flushing time (t) and flow rate (Q), which gives the effective pore radius (R), cross-sectional coverage factor (χ = [1 - (R/R0)2]), and characteristic ratio (rc) of the increments of χ during each adsorption/desorption cycle at a given bulk solution concentration (Cbulk). Our findings include: (1) by gradient increasing Cbulk from 10 to 200 mg/L at Q = 10 mL/h, the shortest PBA40 displays a saturation adsorption behavior when Cbulk ≥ 80 mg/L and t ≥ 2000 s, which agrees well with the prediction of blob model, whereas for the longer PBAN chains, the chain length (N) and concentration-dependent adsorption tendency get stronger as N increases from 40 to 620 at t ≥ 2000 s, in particular, R/R0 ∼ N-0.20 is observed at Cbulk = 140 mg/L; (2) by focusing on the platform χ in the saturation adsorption regime (χsat), the longer PBAN displays a stronger adsorption trend with partially reversible feature at Q = 5.0 mL/h, namely, as N increases from 40 to 620, χsat increases from 0.15 to 0.83 at Cbulk = 100 mg/L, where rc changes from 0.25 ± 0.10 to 0.80 ± 0.10 as the adsorption/desorption flushing cycle increases from 1 to 8 at Cbulk = 100 mg/L; (3) by further assuming a solvent nonpenetrating and nondraining adsorption layer, χsat determined in the case of curved surface can be comparable to the physical meaning of adsorption thickness (Δad) in the case of flat-surface adsorption, and the fitting result indicates χsat ∼ Δad ∼ N0.58, falling between Δad ∼ N1/2 and Δad ∼ N1.0 predicted by the mean-field and scaling theories for real multichain adsorption, respectively. Overall, the present work not only clarifies some controversies but also provides unambiguous evidence supporting the existence of tightly adsorbed internal and loosely adsorbed external layers.

Long term phase separation dynamics in liquid crystal-enriched microdroplets obtained from binary fluid mixtures

The dynamics of long term phase separation in binary liquid mixtures remains a subject of fundamental interest. Here, we study a binary liquid mixture, where the minority phase is confined to a liquid crystal (LC)-rich droplet, by investigating the evolution of size, defect and mesogen alignment over time. We track the binary liquid mixture evolving towards equilibrium by visualising the configuration of the liquid crystal droplet through polarisation microscopy. We compare our experimental findings with computational simulations and elucidate differences between bulk phases and confined droplets based on the respective thermodynamics of phase separation. Our work provides insights on how phase transitions on the microscale can deviate from bulk phase diagrams with relevance to other material systems, such as the liquid-liquid phase separation of polymer and protein solutions.

Ethanol-Induced Aggregation of Nonpolar Nanoparticles in Water/Ethanol Mixed Solvents

The dispersity of nonpolar nanoparticles (NPs) in water/ethanol mixed solvents was studied using molecular dynamics simulations. Based on the rule of "like dissolves like," nonpolar NPs should be dispersed better in a solvent with a lower polarity. As the mole fraction of ethanol in a mixed solvent (R) increases from 0% (pure water) to 100% (pure ethanol), the polarity of the mixed solvent is indicated to decrease monotonically. However, the dispersity of nonpolar NP does not increase monotonically: it first decreases after the addition of a small fraction of ethanol (R < 8.0%) and then markedly increases as R further grows. When there is a small amount of ethanol, the ethanol molecules around aggregated NPs tend to simultaneously make contact with multiple NPs, which can increase the tendency of NP aggregation. Furthermore, with a considerable ethanol ratio, the interaction of the solvent with NPs becomes notably strong, which facilitates the dissolution of NPs. Our findings may help to better understand the mechanism of dispersion of NPs in mixed solvents and may provide a useful precipitation technology for NP production.

Hydration Structures on γ-Alumina Surfaces With and Without Electrolytes Probed by Atomistic Molecular Dynamics Simulations

A wide range of systems, both engineered and natural, feature aqueous electrolyte solutions at interfaces. In this study, the structure and dynamics of water at the two prevalent crystallographic terminations of gamma-alumina, [110] and [100], and the influence of salts─sodium chloride, ammonium acetate, barium acetate, and barium nitrate on such properties─were investigated using equilibrium molecular dynamics simulations. The resulting interfacial phenomena were quantified from simulation trajectories via atomic density profiles, angle probability distributions, residence times, 2-D density distributions within the hydration layers, and hydrogen bond density profiles. Analysis and interpretation of the results are supported by simulation snapshots. Taken together, our results show stronger interaction and closer association of water with the [110] surface, compared to [100], while ion-induced disruption of interfacial water structure was more prevalent at the [100] surface. For the latter, a stronger association of cations is observed, namely sodium and ammonium, and ion adsorption appears determined by their size. The differences in surface-water interactions between the two terminations are linked to their respective surface features and distributions of surface groups, with atomistic-scale roughness of the [110] surface promoting closer association of interfacial water. The results highlight the fundamental role of surface characteristics in determining surface-water interactions, and the resulting effects on ion-surface and ion-water interactions. Since the two terminations of gamma-alumina considered represent interfaces of significance to numerous industrial applications, the results provide insights relevant for catalyst preparation and adsorption-based water treatment, among other applications.

Nanoconfinement facilitates reactions of carbon dioxide in supercritical water

The reactions of CO2 in water under extreme pressure-temperature conditions are of great importance to the carbon storage and transport below Earth's surface, which substantially affect the carbon budget in the atmosphere. Previous studies focus on the CO2(aq) solutions in the bulk phase, but underground aqueous solutions are often confined to the nanoscale, and nanoconfinement and solid-liquid interfaces may substantially affect chemical speciation and reaction mechanisms, which are poorly known on the molecular scale. Here, we apply extensive ab initio molecular dynamics simulations to study aqueous carbon solutions nanoconfined by graphene and stishovite (SiO2) at 10 GPa and 1000 ~ 1400 K. We find that CO2(aq) reacts more in nanoconfinement than in bulk. The stishovite-water interface makes the solutions more acidic, which shifts the chemical equilibria, and the interface chemistry also significantly affects the reaction mechanisms. Our findings suggest that CO2(aq) in deep Earth is more active than previously thought, and confining CO2 and water in nanopores may enhance the efficiency of mineral carbonation.

Insight into the Structure and Dynamics of Ethanol-Water Binary Mixture Confined in Nanochannel by Mica and Graphene

Investigation of the structural properties and dynamics of fluid mixture confined in nanochannels has become an essential topic in many fields due to potential applications in nanofluidic devices and biological systems. Here, we study the ethanol-water blend confined between the mica and single or multilayer graphene for different slit pore widths, ethanol content, and temperatures. Our molecular dynamics simulation indicates that water molecules are adsorbed at the mica surface, while ethanol molecules prefer to be adsorbed near the graphene surface. We find that distinct layers of ethanol molecules form as the channel width and ethanol content in the mixture are increased. The diffusion of confined ethanol and water molecules depends on the nanopore widths, concentrations, and temperatures. Interestingly, at nanopore widths of 1.0 and 1.3 nm, the mobility of confined ethanol molecules is greater than that of water molecules for all ethanol concentrations. In contrast, at pore width of 0.7 nm, the opposite behavior is observed at lower concentrations of ethanol (x_EtOH = 0.1 and 0.3) in the mixture. Furthermore, the diffusivity of ethanol and water in the mixtures increases with increasing the temperatures. The hydrogen bond and cluster analysis imply the segregation of water molecules near the mica surface, while ethanol molecules are near the opposite pore wall (graphene).

Partition Constant of Binary Mixtures for the Equilibrium between a Bulk and a Confined Phase

It is well-known that the thermodynamic, kinetic and structural properties of fluids, and in particular of water and its solutions, can be drastically affected in nanospaces. A possible consequence of nanoscale confinement of a solution is the partial segregation of its components. Thereby, confinement in nanoporous materials (NPM) has been proposed as a means for the separation of mixtures. In fact, separation science can take great advantage of NPM due to the tunability of their properties as a function of nanostructure, morphology, pore size, and surface chemistry. Alcohol-water mixtures are in this context among the most relevant systems. However, a quantitative thermodynamic description allowing for the prediction of the segregation capabilities as a function of the material-solution characteristics is missing. In the present study we attempt to fill this vacancy, by contributing a thermodynamic treatment for the calculation of the partition coefficient in confinement. Combining the multilayer adsorption model for binary mixtures with the Young equation, we conclude that the liquid-vapor surface tension and the contact angle of the pure substances can be used to predict the separation ability of a particular material for a given mixture to a semiquantitative extent. Moreover, we develop a Kelvin-type equation that relates the partition coefficient to the radius of the pore, the contact angle, and the liquid-vapor surface tensions of the constituents. To assess the validity of our thermodynamic formulation, coarse grained molecular dynamics simulations were performed on models of alcohol-water mixtures confined in cylindrical pores. To this end, a coarse-grained amphiphilic molecule was parametrized to be used in conjunction with the mW potential for water. This amphiphilic model reproduces some of the properties of methanol such as enthalpy of vaporization and liquid-vapor surface tension, and the minimum of the excess enthalpy for the aqueous solution. The partition coefficient turns out to be highly dependent on the molar fraction, on the interaction between the components and the confining matrix, and on the radius of the pore. A remarkable agreement between the theory and the simulations is found for pores of radius larger than 15 Å.

Transferable Gaussian Attractive Potentials for Organic/Oxide Interfaces

Organic/oxide interfaces play an important role in many areas of chemistry and in particular for lubrication and corrosion. Molecular dynamics simulations are the method of choice for providing complementary insight to experiments. However, the force fields used to simulate the interaction between molecules and oxide surfaces tend to capture only weak physisorption interactions, discarding the stabilizing Lewis acid/base interactions. We here propose a simple complement to the straightforward molecular mechanics description based on "out-of-the-box" Lennard-Jones potentials and electrostatic interactions: the addition of an attractive Gaussian potential between reactive sites of the surface and heteroatoms of adsorbed organic molecules, leading to the Gaussian Lennard-Jones (GLJ) potential. The interactions of four oxygenated and four amine molecules with the typical and widespread hematite and γ-alumina surfaces are investigated. The root mean square deviation (RMSD) for all probed molecules is only 5.7 kcal/mol, which corresponds to an error of 23% over hematite. On γ-alumina, the RMSD is 11.2 kcal/mol using a single parameter for all five chemically inequivalent surface aluminum atoms. Applying GLJ to the simulation of organic films on oxide surfaces demonstrates that the mobility of the surfactants is overestimated by the simplistic LJ potential, while GLJ and other qualitatively correct potentials show a strong structuration and slow dynamics of the surface films, as could be expected from the first-principles adsorption energies for model head groups.

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 CH₄ and either CO₂ or H₂S 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 CaCl₂ brines are considered, although adsorption was not affected. The acid gases considered either facilitate or hinder CH₄ 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.

Simulation study on the structural and dynamic properties of ethanol confined in nanochannels

The structural and dynamic properties of ethanol molecules in nanochannels of various diameters have been studied using molecular dynamics simulations.

Quantifying Pore Width Effects on Diffusivity via a Novel 3D Stochastic Approach with Input from Atomistic Molecular Dynamics Simulations

The increased production of unconventional hydrocarbons emphasizes the need to understand the transport of fluids through narrow pores. Although it is well-known that confinement affects fluids structure and transport, it is not yet possible to quantitatively predict properties such as diffusivity as a function of pore width in the range of 1-50 nm. Such pores are commonly found not only in shale rocks but also in a wide range of engineering materials, including catalysts. We propose here a novel and computationally efficient methodology to obtain accurate diffusion coefficient predictions as a function of pore width for pores carved out of common materials, such as silica, alumina, magnesium oxide, calcite, and muscovite. We implement atomistic molecular dynamics (MD) simulations to quantify fluid structure and transport within 5 nm-wide pores, with particular focus on the diffusion coefficient within different pore regions. We then use these data as input to a bespoke stochastic kinetic Monte Carlo (KMC) model, developed to predict fluid transport in mesopores. The KMC model is used to extrapolate the fluid diffusivity for pores of increasing width. We validate the approach against atomistic MD simulation results obtained for wider pores. When applied to supercritical methane in slit-shaped pores, our methodology yields data within 10% of the atomistic simulation results, with significant savings in computational time. The proposed methodology, which combines the advantages of MD and KMC simulations, is used to generate a digital library for the diffusivity of gases as a function of pore chemistry and pore width and could be relevant for a number of applications, from the prediction of hydrocarbon transport in shale rocks to the optimization of catalysts, when surface-fluid interactions impact transport.

Understanding methane/carbon dioxide partitioning in clay nano- and meso-pores with constant reservoir composition molecular dynamics modeling

CRC-MD simulations show that nanopores in shales bounded by clay minerals have a strong preference for CO₂ relative to CH₄.

Separation of water-ethanol solutions with carbon nanotubes and electric fields

Bioethanol has been used as an alternative energy source for transportation vehicles to reduce the use of fossil fuels. The separation of water-ethanol solutions from fermentation processes is still an important issue in the production of anhydrous ethanol. Using molecular dynamics simulations, we investigate the effect of axial electric fields on the separation of water-ethanol solutions with carbon nanotubes (CNTs). In the absence of an electric field, CNT-ethanol van der Waals interactions allow ethanol to fill the CNTs in preference to water, i.e., a separation effect for ethanol. However, as the CNT diameter increases, this ethanol separation effect significantly decreases owing to a decrease in the strength of the van der Waals interactions. In contrast, under an electric field, the energy of the electrostatic interactions within the water molecule structure induces water molecules to fill the CNTs in preference to ethanol, i.e., a separation effect for water. More importantly, the electrostatic interactions are dependent on the water molecule structure in the CNT instead of the CNT diameter. As a result, the separation effect observed under an electric field does not diminish over a wide CNT diameter range. Moreover, CNTs and electric fields can be used to separate methanol-ethanol solutions too. Under an electric field, methanol preferentially fills CNTs over ethanol in a wide CNT diameter range.

Evidence for confinement induced phase separation in ethanol-water mixture: a positron annihilation study

We report an experimental evidence for the phase separation of ethanol-water mixture confined in mesoporous silica with different pore size using positron annihilation lifetime spectroscopy (PALS). A bulk-like liquid in the core of the pore and a distinct interfacial region near the pore surface have been identified based on ortho-positronium lifetime components. The lifetime corresponding to the core liquid shows similar behavior to the bulk liquid mixture while the interfacial lifetime shows an abrupt rise within a particular range of ethanol concentration depending on the pore size. This abrupt increase is attributed to the appearance of excess free-volume near the interfacial region. The excess free-volume is originated due to microphase separation of confined ethanol-water primarily at the vicinity of the pore wall. We envisage that probing free-volume changes at the interface using PALS is a sensitive way to investigate microphase separation under nanoconfinement.

Kinetic-contact-driven gigantic energy transfer in a two-dimensional Lennard-Jones fluid confined to a rotating pore

A two-dimensional Lennard-Jones system in a circular and rotating container has been studied by means of molecular dynamics technique. A nonequilibrium transition to the rotating stage has been detected in a delayed time since an instant switching of the frame rotation. This transition is attributed to the increase of the density at the wall because of the centrifugal force. At the same time the phase transition occurs, the inner system changes its configuration of the solid-state type into the liquid type. Impact of angular frequency and molecular roughness on the transport properties of the nonrotating and rotating systems is analyzed.

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.

Confined Water Determines Transport Properties of Guest Molecules in Narrow Pores

We computed the transport of methane through 1 nm wide slit-shaped pores carved out of solid substrates. Models for silica, magnesium oxide, and alumina were used as solid substrates. The pores were filled with water. The results show that the methane permeability through the hydrated pores is strongly dependent on the solid substrate. Detailed analysis of the simulated systems reveals that local properties of confined water, including its structure, and more importantly, evolution of solvation free energy and hydrogen bond structure are responsible for the pronounced differences observed. The simulations are extended to multicomponent systems representative of natural gas, containing methane, ethane, and H2S. The results show that all pores considered have high affinity for H2S, moderate affinity for methane, and low affinity for ethane. The H2S/methane transport selectivity through the hydrated alumina pore is comparable, or superior, to that reported for existing commercial membranes. A multiscale approach was then implemented to demonstrate that a Smoluchowski one-dimensional model is able to reproduce the molecular-level results for short pores when appropriate values for the local self-diffusion coefficients are used as input parameters. We propose that the model can be extended to predict methane transport through uniform hydrated pores of macroscopic length. When verified by experiments, our simulation results could have important implications in applications such as natural gas sweetening and predictions of methane migration through hydraulically fractured shale formations.

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.