Molecular Dynamics Simulations of CO2/Water/Quartz Interfacial Properties: Impact of CO2 Dissolution in Water

Water and Carbon Dioxide Capillary Bridges in Nanoscale Slit Pores: Effects of Temperature, Pressure, and Salt Concentration on the Water Contact Angle

We perform molecular dynamics (MD) simulations of a nanoscale water capillary bridge (WCB) surrounded by carbon dioxide over a wide range of temperatures and pressures (T = 280-400 K and carbon dioxide pressures P CO 2 ≈ 0-80 MPa). The water-carbon dioxide system is confined by two parallel silica-based surfaces (hydroxylated β-cristobalite) separated by h = 5 nm. The aim of this work is to study the WCB contact angle (θc) as a function of T and P CO 2 . Our simulations indicate that θc varies weakly with temperature and pressure: Δθc ≈ 10-20° for P CO 2 increasing from ≈0 to 80 MPa (T = 320 K); Δθc ≈ -10° for T increasing from 320 to 360 K (with a fixed amount of carbon dioxide). Interestingly, at all conditions studied, a thin film of water (1-2 water layers-thick) forms under the carbon dioxide volume. Our MD simulations suggest that this is due to the enhanced ability of water, relative to carbon dioxide, to form hydrogen-bonds with the walls. We also study the effects of adding salt (NaCl) to the WCB and corresponding θc. It is found that at the salt concentrations studied (mole fractions xNa = xCl = 3.50, 9.81%), the NaCl forms a large crystallite within the WCB with the ions avoiding the water-carbon dioxide interface and the walls surface. This results in θc being insensitive to the presence of NaCl.

Nanoscale Patterning of Surface Nanobubbles by Focused Ion Beam

Surface nanobubbles forming on hydrophobic surfaces in water present an exciting opportunity as potential agents of top-down and bottom-up nanopatterning. The formation and characteristics of surface nanobubbles are strongly influenced by the physical and chemical properties of the substrate. In this study, focused ion beam (FIB) milling is used for the first time to spatially control the nucleation of surface nanobubbles with 75 nm precision. The spontaneous formation of nanobubbles on alternating lines of a self-assembled monolayer (octadecyltrichlorosilane) patterned by FIB is detected by atomic force microscopy. The effect of chemical vs topographical surface heterogeneity on the formation of nanobubbles is investigated by comparing samples with OTS coating applied pre- vs post-FIB patterning. The results confirm that nanoscale FIB-based patterning can effectively control surface nanobubble position by means of chemical heterogeneity. The effect of FIB milling on nanobubble morphology and properties, including contact angle and gas oversaturation, is also reported. Molecular dynamics simulations provide further insight into the effects of FIB amorphization on surface nanobubble formation. Combined experimental and simulation investigations offer insights to guide future nanobubble-based patterning using FIB milling.

Contrasting Views of the Electric Double Layer in Electrochemical CO2 Reduction: Continuum Models vs Molecular Dynamics

In the field of electrochemical CO2 reduction, both continuum models and molecular dynamics (MD) models have been used to understand the electric double layer (EDL). MD often focuses on the region within a few nm of the electrode, while continuum models can span up to the device level (cm). Still, both methods model the EDL, and for a cohesive picture of the CO2 electrolysis system, the two methods should agree in the regions where they overlap length scales. To this end, we make a direct comparison between state-of-the-art continuum models and classical MD simulations under the conditions of CO2 reduction on a Ag electrode. For continuum modeling, this includes the Poisson-Nernst-Planck formulation with steric (finite ion size) effects, and in MD the electrode is modeled with the constant potential method. The comparison yields numerous differences between the two modeling methods. MD shows cations forming two adsorbed layers, including a fully hydrated outer layer and a partial hydration layer closer to the electrode surface. The strength of the inner adsorbed layer increases with cation size (Li+ < Na+ < K+ < Cs+) and with more negative applied potentials. Continuum models that include steric effects predict CO2 to be mostly excluded within 1 nm of the cathode due to tightly packed cations, yet we find little evidence to support these predictions from the MD results. In fact, MD shows that the concentration of CO2 increases within a few Å of the cathode surface due to interactions with the Ag electrode, a factor not included in continuum models. The EDL capacitance is computed from the MD results, showing values in the range of 7-9 μF cm-2, irrespective of the electrolyte concentration, cation identity, or applied potential. The direct comparison between the two modeling methods is meant to show the areas of agreement and disagreement between the two views of the EDL, so as to improve and better align these models.

Interfacial tension and wettability alteration during hydrogen and carbon dioxide storage in depleted gas reservoirs

The storage of CO2 and hydrogen within depleted gas and oil reservoirs holds immense potential for mitigating greenhouse gas emissions and advancing renewable energy initiatives. However, achieving effective storage necessitates a thorough comprehension of the dynamic interplay between interfacial tension and wettability alteration under varying conditions. This comprehensive review investigates the multifaceted influence of several critical parameters on the alterations of IFT and wettability during the injection and storage of CO2 and hydrogen. Through a meticulous analysis of pressure, temperature, treatment duration, pH levels, the presence of nanoparticles, organic acids, anionic surfactants, and rock characteristics, this review elucidates the intricate mechanisms governing the changes in IFT and wettability within reservoir environments. By synthesizing recent experimental and theoretical advancements, this review aims to provide a holistic understanding of the processes underlying IFT and wettability alteration, thereby facilitating the optimization of storage efficiency and the long-term viability of depleted reservoirs as carbon capture and storage or hydrogen storage solutions. The insights gleaned from this analysis offer invaluable guidance for researchers, engineers, and policymakers engaged in harnessing the potential of depleted reservoirs for sustainable energy solutions and environmental conservation. This synthesis of knowledge serves as a foundational resource for future research endeavors aimed at enhancing the efficacy and reliability of CO2 and hydrogen storage in depleted reservoirs.

Atomic-Scale Insights into the Phase Behavior of Carbon Dioxide and Water from 313 to 573 K and 8 to 30 MPa

We performed molecular dynamics (MD) simulations of CO2 + H2O systems by employing widely used force fields (EPM2, TraPPE, and PPL models for CO2; SPC/E and TIP4P/2005 models for H2O). The phase behavior observed in our MD simulations is consistent with the coexistence lines obtained from previous experiments and SAFT-based theoretical models for the equations of state. Our structural analysis reveals a pronounced correlation between phase transitions and the structural orderliness. Specifically, the coordination number of Ow (oxygen in H2O) around other Ow significantly correlates with phase changes. In contrast, coordination numbers pertaining to the CO2 molecules show less sensitivity to the thermodynamic state of the system. Furthermore, our data indicate that a predominant number of H2O molecules exist as monomers without forming hydrogen bonds, particularly in a CO2-rich mixture, signaling a breakdown in the hydrogen bond network's orderliness, as evidenced by a marked decrease in tetrahedrality. These insights are crucial for a deeper atomic-level understanding of phase behaviors, contributing to the well-grounded design of CO2 injection under high-pressure and high-temperature conditions, where an atomic-scale perspective of the phase behavior is still lacking.

The impact of salinity on the cohesion process of quartz substracts: A molecular dynamics study

The moisture with salt ions adsorbed on the mineral soil surface is crucial to the cohesion process when the media is exposed to marine or coastal environments. However, the impact of salinity on the cohesion of soils is not well studied at the nanoscale. In this study, the salinity effect was investigated by studying the wettability and capillary force of NaCl solutions on quartz substrates via a molecular dynamics-based approach. Besides, a new visualization method was proposed to measure the contact angle of liquid droplets from the aspect of nanoscale. The results indicated that salt ions can weaken the wettability of the liquid on the quartz surface and inhibit the capillary force. Compared with water, the liquid with a 10% NaCl solution can achieve a capillary force reduction of around 70%, resulting in a detrimental effect on the cohesion of soils. Overall, this study enhanced the understanding of the nanoscale salinity effect on the cohesion process and provided insights into the modification of the mechanical properties of soils from the aspect of nanoscale.

Interfacial properties of the brine + carbon dioxide + oil + silica system

Molecular dynamics simulations of the H2O + CO2 + aromatic hydrocarbon and H2O + CO2 + benzene + silica (hydrophilic) systems are performed to gain insights into CO2-enhanced oil recovery (EOR) processes. For comparison purposes, an overview of the previous simulation studies of the interfacial properties of the brine + CO2 + alkane + silica system is also presented. In general, the water contact angle (CA) of the H2O + CO2 + silica (hydrophilic) system increased with pressure and decreased with temperature. The CAs of the H2O + hydrocarbon + silica (hydrophilic) system are not significantly affected by temperature and pressure. The simulated CAs were in the ranges of about 58°-77° and 81°-93° for the H2O + hexane + silica (hydrophilic) and the H2O + aromatic hydrocarbon + silica (hydrophilic) systems, respectively. In general, these CAs were not significantly influenced by the addition of CO2. The simulated CAs were in the ranges of about 51.4°-95.0°, 69.1°-86.0°, and 72.0°-87.9° for the brine + CO2 + silica (hydrophilic), brine + hexane + silica (hydrophilic), and brine + CO2 + hexane + silica (hydrophilic) systems, respectively. All these CAs increased with increasing NaCl concentration. The adhesion tension of the brine + silica (hydrophilic) system in the presence of CO2 and/or hexane decreased with increasing salt concentration. The simulated CAs were in the range of about 117°-139° for the H2O + alkane + silica (hydrophobic) system. These CAs are increased by the addition of CO2. At high pressures, the distributions of H2O normal to the silica (hydrophobic) surface in the droplet region of the H2O + silica system were found to be strongly affected by the presence of CO2. These insights might be key for optimizing the performance of the miscible CO2 water-alternating-gas injection schemes widely used for EOR.

Water Film Structure and Wettability of Different Quartz Surfaces: Hydrogen Bonding Across Various Cutting Planes

Quartz is ubiquitous in subsurface formations. The crystal faces have different atomic arrangements. Knowledge of the molecular structures on the surface of quartz is key in many processes. Molecular dynamics simulations are conducted to investigate the atomic arrangement effect on the water film structure, ion adsorption, and wettability at three different α-quartz surfaces. The interfacial structures depend on the quartz surface. Intrasurface hydrogen bonding between surface silanols differs in the α-quartz surface. At the (0001) surface, the OH density is 9.58 nm-2, consisting of Q2 units with two hydroxyl groups per silicone atom. At the (101̅0)-β surface, the OH density is 7.54 nm-2, consisting of Q3 units with one hydroxyl group per silicone atom; there is significant intrasurface hydrogen bonding. At the (101̅0)-α surface, the OH density is 7.54 nm-2, consisting of Q2 units; however, there is little intrasurface hydrogen bonding. The intrasurface hydrogen bonding results in the exposure of hydrogen-bond acceptors to the aqueous phase, causing water molecules to have an H-up (hydrogen toward surface) orientation. This orientation can be found at the (0001) and (101̅0)-β surfaces; it is related to the degree of ordering at the surface. The ordering at the (0001) and (101̅0)-β surfaces is higher than that at the (101̅0)-α surface. In aqueous systems with ions, cation adsorption is the most dominant at the (0001) surface due to the largest surface density of the intrasurface hydrogen bonding, providing interaction sites for cations to be adsorbed. We observe a pronounced decrease in water film thickness from the ions at the (0001) surface only, likely due to significant cation adsorption. In this work, we demonstrate that the hydrogen-bond network, which varies from the plane cut, affects the water film structure and ion adsorption. The contact is nearly zero despite the changes in the film thickness and molecular structure at the temperature of 318 K.

Effect of microstructure in mesoporous adsorbents on the adsorption of low concentrations of VOCs: An experimental and simulation study

This study evaluated the adsorption of five volatile organic compounds (VOCs) on Opoka, precipitated silica, and palygorskite, to elucidate the effect of their pore size on VOCs adsorption. The adsorption capacity of these adsorbents is not only highly correlated with their surface area and pore volume, but also notably improved by the presence of micropores. The variation in adsorption capacity for different VOCs was primarily influenced by their boiling point and polarity. Palygorskite, which had the smallest total pore volume (0.357 cm3/g) but the largest micropore volume (0.043 cm3/g) among the three adsorbents, exhibited the highest adsorption capacity for all tested VOCs. Additionally, the study constructed slit pore models of palygorskite with micropores (0.5 and 1.5 nm) and mesopores (3.0 and 6.0 nm), calculated and discussed the heat of adsorption, concentration distribution, and interaction energy of VOCs adsorbed on different pore models. The results revealed that the adsorption heat, concentration distribution, total interaction energy, and van der Waals energy decrease with increasing pore size. The concentration of VOCs in 0.5 nm pore was nearly three times that in 6.0 nm pore. This work can also provide guidance for further research on using adsorbents with mixed microporous and mesoporous structures to control VOCs.

Interfacial properties of the hexane + carbon dioxide + water system in the presence of hydrophilic silica

Molecular dynamics simulations were conducted to study the interfacial behavior of the CO₂ + H₂O and hexane + CO₂ + H₂O systems in the presence of hydrophilic silica at geological conditions. Simulation results for the CO₂ + H₂O and hexane + CO₂ + H₂O systems are in reasonable agreement with the theoretical predictions based on the density functional theory. In general, the interfacial tension (IFT) of the CO₂ + H₂O system exponentially (linearly) decreased with increasing pressure (temperature). The IFTs of the hexane + CO₂ + H₂O (two-phase) system decreased with the increasing mole fraction of CO₂ in the hexane/CO₂-rich phase x_CO2 . Here, the negative surface excesses of hexane lead to a general increase in the IFTs with increasing pressure. The effect of pressure on these IFTs decreased with increasing x_CO2 due to the positive surface excesses of carbon dioxide. The simulated water contact angles of the CO₂ + H₂O + silica system fall in the range from 43.8° to 76.0°, which is in reasonable agreement with the experimental results. These contact angles increased with pressure and decreased with temperature. Here, the adhesion tensions are influenced by the variations in fluid-fluid IFT and contact angle. The simulated water contact angles of the hexane + H₂O + silica system fall in the range from 58.0° to 77.0° and are not much affected by the addition of CO₂. These contact angles increased with pressure, and the pressure effect was less pronounced at lower temperatures. Here, the adhesion tensions are mostly influenced by variations in the fluid-fluid IFTs. In all studied cases, CO₂ molecules could penetrate into the interfacial region between the water droplet and the silica surface.

Hydrogeochemical characterization and suitability of groundwater for drinking and irrigation in Menzel Bourguiba aquifers (Northeastern Tunisia)

The aquifer system of Menzel Bourguiba, located in northeastern Tunisia, is an important groundwater resource that, for the last decades, has been a subject of excessive use by agricultural and industrial activities leading to the degradation of its water quality. To better understand the hydrochemistry of these aquifers, identify the main factors controlling their mineralization, and assess their suitability for drinking and agricultural purposes, 35 groundwater samples collected during the campaign of October 2019 were analyzed by different physicochemical methods. The results showed that the hydrochemical facies of the study area is characterized by two major types SO4-Cl-Ca and Na-Cl. The total dissolved solids (TDS) values range from 0.5 to 1.9 g/L in the shallow aquifer and from 0.6 to 2 g/L in the deep aquifer. The high rates are recognized downstream of the studied basin. The interpretation of multivariate statistical analysis and geochemical approaches revealed that rock-water interaction, dissolution and precipitation process, ionic exchange mechanisms, water irrigation return, and marine intrusion are the main factors controlling the mineralization of the study area groundwater. The comparison of the analytical data to the Tunisian and World Health Organization (WHO) standards unveiled that the majority of groundwater samples are unsuitable for drinking because of high nitrate contents, elevated electrical conductivity (EC), TDS values, and high concentrations of some chemical parameters. In addition, the sodium absorption ratio and the sodium percent parameters indicated that most of the analyzed samples are inappropriate for irrigation uses.

Capillary bridges between unsaturated nano-mineral particles: a molecular dynamics study

Capillary bridges play an important role in the process of cohesion, which is crucial for wet granular media, and engineering of pharmaceuticals and food processing. However, the understanding of capillary bridges at the nanoscale remains unclear because the mechanical performance of nanoscale capillary bridges cannot be fully captured and explained by classical capillary theory. We applied a novel molecular dynamic simulation to investigate the dynamic formation process of nanoscale capillary bridges between quartz asperities. In comparison with classical capillary theory, our results suggested that the application of the toroidal approximation and gorge method will break down at the scale of 1 nm. Below this threshold, a pronounced oscillation in the adhesive force was observed due to inconsistent distribution of water molecules in the capillary bridges. Moreover, we found a non-linear correlation between the adhesive force and the saturation degree. Different from the cohesive stress of sandy soil as a function of saturation degree, we identified an optimal saturation range of 0.5-0.7 instead of 0.2-0.9 for the sandy soil. Our findings enhance the understanding of capillary bridges and provide new insights into the capillary force between particles in the fields of geotechnical engineering, food-process engineering, the pharmaceutical industry and nanotechnology.

Microscopic properties and stabilization mechanism of a supercritical carbon dioxide microemulsion with extremely high water content

HYPOTHESIS

Developing the supercritical carbon dioxide microemulsion with a broad water content (W₀) window can provide more possibility for designing highly efficient chemical processes, which is challenging due to the lack of comprehension about its formation mechanism. Molecular dynamics simulation method is expected to reveal the microscopic stabilization mechanism of high-W₀ microemulsions.

SIMULATIONS

All-atom molecular dynamics simulations of the ternary systems with varied W₀ stabilized by 4FG(EO)₂ surfactant were designed according to phase behavior experiments. A systematic investigation was performed concerning the self-assembling, equilibrium morphology and detailed microstructure of the microemulsion droplet. An in-depth comparative study about the distribution of both H₂O and CO₂, the interfacial behaviors of 4FG(EO)₂, as well as the microscopic interactions was conducted.

FINDINGS

For the first time, direct evidence was provided for the formation of water-in-carbon dioxide microemulsion with extremely high W₀ (80) under the effect of 4FG(EO)₂. Furthermore, a unique interfacial phenomenon, i. e. CO₂ accumulating at the interface, was revealed to be responsible for the formation and enhanced stability of the nanosized droplet with high W₀. This should set a new guiding star for synthesizing and selecting effective interfacial modifiers to create high-W₀ microemulsions.

Bulk and Interfacial Properties of the Decane + Water System in the Presence of Methane, Carbon Dioxide, and Their Mixture

Molecular dynamics simulations are carried out to study the two-phase behavior of the n-decane + water system in the presence of methane, carbon dioxide, and their mixture under reservoir conditions. The simulation studies were complemented by theoretical modeling using the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state (EoS) and density gradient theory. Our results show that the presence of methane and carbon dioxide decreases the interfacial tension (IFT) of the decane + water system. In general, the IFT increases with increasing pressure and decreasing temperature for the methane + decane + water and carbon dioxide + decane + water systems, similar to what has been found for the corresponding decane + water system. The most important finding of this study is that the presence of carbon dioxide decreases the IFT of the methane + decane + water system. The atomic density profiles provide evidence of the local accumulation of methane and carbon dioxide at the interface, in most of the studied systems. The results of this study show the preferential dissolution in the water-rich phase and enrichment at the interface for carbon dioxide in the methane + carbon dioxide + decane + water system. This indicates the preferential interaction of water with carbon dioxide relative to methane and decane. Notably, there is an enrichment of the interface by decane at high mole fractions of methane in the methane/decane-rich or methane/carbon dioxide/decane-rich phase. Overall, the solubility of methane and carbon dioxide in the water-rich phase increases with increasing pressure and temperature. Additionally, we find that the overall performance of the PC-SAFT EoS and the cubic-plus-association EoS is similar with respect to the calculation of bulk and interfacial properties of these systems.

Interaction of low salinity surfactant nanofluids with carbonate surfaces and molecular level dynamics at fluid-fluid interface at ScCO2 loading

HYPOTHESIS

The advanced low salinity aqueous formulations are yet to be validated as an injection fluid for enhanced oil recovery (EOR) from the carbonate reservoirs and CO₂ geosequestration. Interaction of various ionic species present in the novel low salinity surfactant nanofluids with scCO₂/CO₂ saturated aqueous phase interface and at the interface of CO₂ saturated aqueous phase/mixed wet (with CO₂ and Decane) limestone surface at the conditions of low salinity at reservoir conditions are to yet to be understood.

EXPERIMENTS

This study, carried out for the first time in low salinity at scCO₂ loading conditions at 20 MPa pressure and 343 K temperature, comprises of wettability study of the limestone surface by aqueous phase contact angle measurements using ZrO₂ nanoparticles (in the concentration range of 100-2000 mg/L) and 0.82 mM Hexadecyltrimethylammonium bromide (CTAB) surfactant. Molecular dynamics simulations results were used to understand the underlying mechanism of wettability alteration and interfacial tension (IFT) change.

FINDINGS

This study reveals that a low dosage (100 mg/L) of ZrO₂ nanoparticles forming ZrO₂-CTAB nano-complexes helps in wettability alteration of the rock surface to more water-wetting state; certain ionic species augment this effect when used in appropriate concentration. Also, these nano-complexes helps in scCO₂/CO₂ saturated aqueous phase IFT reduction. This study can be used to design advanced low salinity injection fluids for water alternating gas injection for EOR and CO₂ geosequestration projects.

CO2 wetting on pillar-nanostructured substrates

CO₂ capture by dropwise CO₂ condensation on cold solid surfaces is a promising technology. Understanding the role of the nanoscale surface and topographical features of CO₂ droplet wetting characteristics is of importance for CO₂ capture by this technology, but this remains unexplored as of yet. Here, using large-scale molecular dynamics (MD) simulations, the contact angle and wetting behaviors of CO₂ droplets on pillar-structured Cu-like surfaces are investigated for the first time. Dynamic wetting simulations show that, by changing the strength of the solid-liquid attraction [Formula: see text] a smooth Cu-like surface offers a transition from CO₂-philic to CO₂-phobic. By periodically pillared roughening of the Cu-like surfaces, however, a higher contact angle and a smaller spreading exponent of a liquid CO₂ droplet are realized. Particularly, a critical crossover of CO₂-philic to CO₂-phobic can appear. The wetting of the pillared surfaces by a liquid CO₂ droplet proceeds non-uniformly. A liquid CO₂ droplet is capable of exhibiting a transition from the Cassie state to the Wenzel state with increasing [Formula: see text] increasing inter-pillar distance, and increasing pillar width. The wetting morphologies of the metastable Wenzel state of a CO₂ droplet are very different from each other. The findings will inform the ongoing design of CO₂-phobic solid surfaces for practical dropwise condensation-based CO₂ capture applications.

High-performance carbon dioxide capture and storage by multi-functional sphingosine kinase inhibitors through a CO2-philic membrane

Carbon dioxide (CO₂) capture using environmentally friendly sphingosine-based materials was theoretically studied.

Transitions between nanomechanical and continuum mechanical contacts: new insights from liquid structure

The use of continuum mechanics to describe contacts involving nanoscale and atomic interactions has been one of the key controversies in nanoscience, tribology, and petrophysical and geological studies. By applying a novel nonequilibrium molecular dynamics scheme to wet quartz contacts, this study revealed the key transitions between continuum electrostatic, nanomechanical and Hertzian contact behaviors at around one nm of surface separation, which results in critical contact pressure fluctuations between -30 and 100 MPa. Using a novel liquid-structure analysis scheme based on the spatial distribution of water molecules, the nanomechanical behavior was found to originate from the collapse and localization of layers of water molecules. Moreover, the role of surface curvature on this effect was also quantified and explained based on a new topological descriptor. The findings of this study enrich our understanding of wet contacts and have a wide range of applications from the nanoscale to macroscale.

Wetting Properties of the CO2-Water-Calcite System via Molecular Simulations: Shape and Size Effects

Assessment of the risks and environmental impacts of carbon geosequestration requires knowledge about the wetting behavior of mineral surfaces in the presence of CO₂ and the pore fluids. In this context, the interfacial tension (IFT) between CO₂ and the aqueous fluid and the contact angle, θ, with the pore mineral surfaces are the two key parameters that control the capillary pressure in the pores of the candidate host rock. Knowledge of these two parameters and their dependence on the local conditions of pressure, temperature, and salinity is essential for the correct prediction of structural and residual trapping. We have performed classical molecular dynamics simulations to predict the CO₂-water IFT and the CO₂-water-calcite contact angle. The IFT results are consistent with previous simulations, where simple point charge water models have been shown to underestimate the water surface tension, thus affecting the simulated IFT values. When combined with the EPM2 CO₂ model, the SPC/Fw water model indeed underestimates the IFT in the low-pressure region at all temperatures studied. On the other hand, at high pressure and low temperature, the IFT is overestimated by ∼5 mN/m. Literature data regarding the CO₂/water/calcite contact angle on calcite are contradictory. Using our new set of force field parameters, we performed NVT simulations at 323 K and 20 MPa to calculate the contact angle of a water droplet on the calcite {10.4} surface in a CO₂ atmosphere. We performed simulations for both spherical and cylindrical droplet configurations for different initial radii to study the size dependence of the water contact angle on calcite in the presence of CO₂. Our results suggest that the contact angle of a cylindrical droplet, is independent of droplet size, for droplets with a radius of 50 Å or more. On the contrary, spherical droplets make a contact angle that is strongly influenced by their size. At the largest size explored in this study, both spherical and cylindrical droplets converge to the same contact angle, 38°, indicating that calcite is strongly wetted by water.

Wettability of rock/CO2/brine and rock/oil/CO2-enriched-brine systems:Critical parametric analysis and future outlook

CO₂ geo-sequestration is a promising technology to permanently store CO₂ in geological formations to control the atmospheric carbon footprint. In addition, CO₂ is frequently utilized in enhanced oil recovery operations to accelerate oil production. Both, CO₂ geo-storage and EOR, are significantly influenced by the wettability of the associated rock/CO₂/brine systems. Wettability drives the multiphase flow dynamics, and microscopic fluid distribution in the reservoir. Furthermore, while wettability is known to be influenced by varying in-situ conditions and surface chemistry of the rock/mineral, the current state-of-the-art indicates wider variabilities of the wetting states. This article, therefore, critically reviews the published datasets on CO₂ wettability of geological formations. Essentially, the rock/CO₂/brine and rock/crude-oil/CO₂-enriched-brine contact angle datasets for the important reservoir rocks (i.e. sandstone and carbonate rocks), as well as for the key minerals quartz and calcite are considered. Also, the parameters that influence wettability are critically analyzed, and the associated parametric trends are discussed and summarized. Finally, we identify pertinent research gaps and define the outlook of future research. The review, therefore, establishes a repository of the recent contact angle data, which thus assists to enhance our current understanding of the subject.

Application of the CLAYFF and the DREIDING Force Fields for Modeling of Alkylated Quartz Surfaces

To extend applicability and to overcome limitations of combining rules for nonbond potential parameters, in this study, CLAYFF and DREIDING force fields are coupled at the level of atomic site charges to model quartz surfaces with chemisorpt hydrocarbons. Density functional theory and Bader charge analysis are applied to calculate charges of atoms of the OC bond connecting a quartz crystal and an alkyl group. The study demonstrates that the hydrogen atom of the quartz surface hydroxyl group can be removed and its charge can be redistributed among the oxygen and carbon atoms of the OC bond in a manner consistent with the results calculated at the density functional level of theory. Augmented with modified charges of the OC bond, force fields can then be applied to a practical problem of evaluation of the contact angle of a water droplet on alkylated quartz surfaces in a carbon dioxide environment, which is relevant for carbon geo-sequestration and in a broader context of oil and gas recovery. Alkylated quartz surfaces have been shown to be extremely hydrophobic even when the surface density of hydroxyl groups is close to the highest naturally observed density of 6.2 OH groups per square nanometer.

Atomistic Study of Dynamic Contact Angles in CO2-Water-Silica System

The dynamic wetting for the CO₂-water-silica system occurring in deep reservoirs is complex because of the interactions among multiple phases. This work aims to quantify the contact angle of CO₂-water flow in the silica channel at six different flow velocities using molecular dynamics. The dynamic contact angle values at different contact line velocities are obtained for the CO₂-water-silica system. By calculating the rates of the adsorption-desorption process of CO₂ and water molecules on the silica surface using molecular dynamics simulations, it has been found that the results of the dynamic contact angle can be explained by the molecular kinetic theory and predicted from the equilibrium molecular simulations. Moreover, the capillary pressure at different contact line velocities is predicted according to the Young-Laplace equation. The change in contact angles at different velocities is compared with empirical equations in terms of capillary number. The results of this study can help us better understand the dynamic process of the multiphase flow at the nanoscale under realistic reservoir conditions.

Characterizing surface wetting and interfacial properties using enhanced sampling (SWIPES)

We introduce an accurate and efficient method for characterizing surface wetting and interfacial properties, such as the contact angle made by a liquid droplet on a solid surface, and the vapor-liquid surface tension of a fluid. The method makes use of molecular simulations in conjunction with the indirect umbrella sampling technique to systematically wet the surface and estimate the corresponding free energy. To illustrate the method, we study the wetting of a family of Lennard-Jones surfaces by water. For surfaces with a wide range of attractions for water, we estimate contact angles using our method, and compare them with contact angles obtained using droplet shapes. Notably, our method is able to capture the transition from partial to complete wetting as surface-water attractions are increased. Moreover, the method is straightforward to implement and is computationally efficient, providing accurate contact angle estimates in roughly 5 nanoseconds of simulation time.

Modeling CO2-Water-Mineral Wettability and Mineralization for Carbon Geosequestration

Carbon dioxide (CO₂) capture and storage (CCS) is an important climate change mitigation option along with improved energy efficiency, renewable energy, and nuclear energy. CO₂ geosequestration, that is, to store CO₂ under the subsurface of Earth, is feasible because the world's sedimentary basins have high capacity and are often located in the same region of the world as emission sources. How CO₂ interacts with the connate water and minerals is the focus of this Account. There are four trapping mechanisms that keep CO₂ in the pores of subsurface rocks: (1) structural trapping, (2) residual trapping, (3) dissolution trapping, and (4) mineral trapping. The first two are dominated by capillary action, where wettability controls CO₂ and water two-phase flow in porous media. We review state-of-the-art studies on CO₂/water/mineral wettability, which was found to depend on pressure and temperature conditions, salt concentration in aqueous solutions, mineral surface chemistry, and geometry. We then review some recent advances in mineral trapping. First, we show that it is possible to reproduce the CO₂/water/mineral wettability at a wide range of pressures using molecular dynamics (MD) simulations. As the pressure increases, CO₂ gas transforms into a supercritical fluid or liquid at ∼7.4 MPa depending on the environmental temperature. This transition leads to a substantial decrease of the interfacial tension between CO₂ and reservoir brine (or pure water). However, the wettability of CO₂/water/rock systems depends on the type of rock surface. Recently, we investigated the contact angle of CO₂/water/silica systems with two different silica surfaces using MD simulations. We found that contact angle increased with pressure for the hydrophobic (siloxane) surface while it was almost constant for the hydrophilic (silanol) surface, in excellent agreement with experimental observations. Furthermore, we found that the CO₂ thin films at the CO₂-hydrophilic silica and CO₂-H₂O interfaces displayed a linear correlation, which can in turn explain the constant contact angle on the hydrophilic silica surface. In view of the literature and our study results, a few recommendations seem necessary to construct a molecular system suitable to study wettability with MD simulations. Future work should be conducted to determine the influence of brine salinity on the wettability of minerals with high cation exchange capacity. Mineral trapping is believed to be an extremely slow process, likely taking thousands of years. However, a recent pilot study demonstrated that CO₂ mineralization occurs within 2 years in highly reactive basalt reservoirs. A first-principles MD study has also shown that carbonation reactions occur rapidly at the surface oxygen sites of a reactive mineral. We observed carbonate ions on both a newly cleaved quartz surface (without hydrolysis), and a basalt andesine surface after hydrolysis in a CO₂-rich environment. Future work should consider the influence of water, gas impurities, and mineral cation type on carbonation.

CO2 activating hydrocarbon transport across nanopore throat: insights from molecular dynamics simulation

In tight oil reservoirs, nanopore throat acting as the narrowest section of fluidic channel determines the oil transport performance; injecting CO₂ is found to significantly promote the oil flow.

Water Contact Angle Dependence with Hydroxyl Functional Groups on Silica Surfaces under CO2 Sequestration Conditions

Functional groups on silica surfaces under CO2 sequestration conditions are complex due to reactions among supercritical CO2, brine and silica. Molecular dynamics simulations have been performed to investigate the effects of hydroxyl functional groups on wettability. It has been found that wettability shows a strong dependence on functional groups on silica surfaces: silanol number density, space distribution, and deprotonation/protonation degree. For neutral silica surfaces with crystalline structure (Q(3), Q(3)/Q(4), Q(4)), as silanol number density decreases, contact angle increases from 33.5° to 146.7° at 10.5 MPa and 318 K. When Q(3) surface changes to an amorphous structure, water contact angle increases 20°. Water contact angle decreases about 12° when 9% of silanol groups on Q(3) surface are deprotonated. When the deprotonation degree increases to 50%, water contact angle decreases to 0. The dependence of wettability on silica surface functional groups was used to analyze contact angle measurement ambiguity in literature. The composition of silica surfaces is complicated under CO2 sequestration conditions, the results found in this study may help to better understand wettability of CO2/brine/silica system.