Molecular structure and dynamics of an aqueous sodium chloride solution in nano-pores between portlandite surfaces: a molecular dynamics study

Giant Effect of CO2 Injection on Multiphase Fluid Adsorption and Shale Gas Production: Evidence from Molecular Dynamics

Carbon dioxide (CO2) injection in unconventional gas-bearing shale reservoirs is a promising method for enhancing methane recovery efficiency and mitigating greenhouse gas emissions. The majority of methane is adsorbed within the micropores and nanopores (≤50 nm) of shale, which possess extensive surface areas and abundant adsorption sites for the sequestration system. To comprehensively discover the underlying mechanism of enhanced gas recovery (EGR) through CO2 injection, molecular dynamics (MD) provides a promising way for establishing the shale models to address the multiphase, multicomponent fluid flow behaviors in shale nanopores. This study proposes an innovative method for building a more practical shale matrix model that approaches natural underground environments. The grand canonical Monte Carlo (GCMC) method elucidates gas adsorption and sequestration processes in shale gas reservoirs under various subsurface conditions. The findings reveal that previously overlooked pore slits have a significant impact on both gas adsorption and recovery efficiency. Based on the simulation comparisons of absolute and excess uptakes inside the kerogen matrix and the shale slits, it demonstrates that nanopores within the kerogen matrix dominate the gas adsorption while slits dominate the gas storage. Regarding multiphase, multicomponent fluid flow in shale nanopores, moisture negatively influences gas adsorption and carbon storage while promoting methane recovery efficiency by CO2 injection. Additionally, saline solution and ethane further impede gas adsorption while facilitating displacement. Overall, this work elucidates the substantial effect of CO2 injection on fluid transport in shale formations and advances the comprehensive understanding of microscopic gas flow and recovery mechanisms with atomic precision for low-carbon energy development.

Cadmium pollution exacerbated by drought: Insights from the nanoscale interaction at the clay mineral surface

Drought is a global environmental problem, while the effect of drought-induced unsaturation on the fate of heavy metal ions is still poorly understood, particularly the lack of mechanistic information at the molecular level. This study used molecular dynamics simulations to investigate nanoscale interactions at the montmorillonite surface under different moisture conditions. Compared to the saturated condition, drought increased the amounts and strength of Cd2+ ions adsorbed on the montmorillonite (MMT) surface while decreased the diffusivity, which was especially obvious in extreme drought conditions (θv=21%-7%). This is closely related to the compressed electric double layer, overcompensation of surface charge, and increased ion pair interactions, resulting from the confinement of water films under drought stress. Further analysis showed that the decrease of hydration effect was responsible for the exacerbated cadmium pollution. Therefore, this study may break the stereotypes about the interactions between heavy metal ions and soil minerals. The results suggest that water management (e.g., irrigation) may be prioritized before beginning heavy metal remediation.

Impact of Hydrostatic Pressure on Molecular Structure and Dynamics of the Sodium and Chloride Ions in Portlandite Nanopores

In order to address the issues of energy depletion, more resources are being searched for in the deep sea. Therefore, research into how the deep-sea environment affects cement-based materials for underwater infrastructure is required. This paper examines the impact of ocean depth (0, 500, 1000, and 1500 m) on the ion interaction processes in concrete nanopores using molecular dynamics simulations. At the portlandite interface, the local structural and kinetic characteristics of ions and water molecules are examined. The findings show that the portlandite surface hydrophilicity is unaffected by increasing depth. The density profile and coordination number of ions alter as depth increases, and the diffusion speed noticeably decreases. The main cause of the ions' reduced diffusion velocity is expected to be the low temperature. This work offers a thorough understanding of the cement hydration products' microstructure in deep sea, which may help explain why cement-based underwater infrastructure deteriorates over time.

Salinity-Driven Structural and Viscosity Modulation of Confined Polar Oil Phases by Carbonated Brine Films: Novel Insights from Molecular Dynamics

The structural and dynamic properties of fluids under confinement in a porous medium differ from their bulk properties. This study delves into the surface structuring and hydrodynamic characteristics of oil/thin film carbonated brine two-phase within a calcite channel upon salinity variation. To this end, both equilibrium and non-equilibrium molecular dynamics simulations are utilized to unveil the effect of the carboxylic acid component (benzoic acid) in a simple model oil (decane) confined between two thin films of carbonated brine on the oil-brine-calcite characteristics. The salinity effect was scrutinized under four saline carbonated waters, deionized carbonated water (DCW), carbonated low-salinity brine (CLSB, 30,000 ppm), carbonated seawater (CSW, 60,000 ppm), and carbonated high-salinity brine (CHSB, 180,000 ppm). An electrical double layer (EDL) is observed at varying salinities, comprising a Stern-like positive layer (formed by Na+ ions) followed by a negative one (formed by Cl- ions primarily residing on top of the adsorbed sodium cations). By lowering the salinity, the Na+ ions cover the interface regions (brine-calcite and brine-oil), depleting within the brine bulk region. The lowest positive surface charge on the rock surface was found in salinity corresponding to seawater. Two distinct Na+ peaks at the oleic phase interface have been observed in the carbonated high-salinity brine system, enhancing the adsorption of polar molecules at the thin brine film interfaces. There is a pronounced EDL formation at the oleic phase interface in the case of CSW, resulting in a strong interface region containing ions and functional fractions. Likewise, the oil region confined by CSW exhibited the lowest apparent viscosity, attributed to the optimized salinity distribution and inclination of benzoic acid fractions uniformly at the brine-oil interface, acting as a slippery surface. Moreover, the results reveal that the presence of polar fractions could increase the oil phase's apparent viscosity, and introducing ions to this system reduces the polar molecules' destructive effect on the apparent viscosity of the oil region. Therefore, the fluidity of confined systems is modulated by both composition of the brine and oil phases.

Computational Simulations of Adsorption Behavior of Anionic Surfactants at the Portlandite-Water Interface under Sulfate and Calcium Ions

The adsorption behaviors of two kinds of anionic surfactants (called HSO4 and HPO4, respectively) with different negatively charged hydrophilic head groups (sulfate and phosphate groups) under different concentrations of sulfate and calcium ions at the portlandite-water interface are investigated by molecular dynamics simulations. Although the adsorption strength of HPO4 is much greater than that of HSO4, the desorption energy of HSO4 is slightly greater at an early stage of desorption due to a more perpendicular orientation and denser packing of hydrophobic tail chains. After adding ions, the sulfate ion has a significant weakening effect due to competitive adsorption, and the negative influence of the calcium ion is weaker, and it even slightly promotes the adsorption at low concentration. Due to the stronger electrostatic interaction of phosphate head groups with the portlandite surface, adsorption strength and adsorption stability for HPO4 are always greater than that of HSO4 under the interference of sulfate ions. The competitive adsorption of the sulfate ion significantly weakens the interaction of hydrophilic head groups with portlandite and the dense packing of two surfactants. The calcium ion with low concentration approaches the portlandite surface and acts as an ion bridge to slightly enhance the adsorption of the surfactant. The ion bridging effect is stronger in the HPO4 system than in the HSO4 system.

Interaction of Nitrite Ions with Hydrated Portlandite Surfaces: Atomistic Computer Simulation Study

The nitrite admixtures in cement and concrete are used as corrosion inhibitors for steel reinforcement and also as anti-freezing agents. The characterization of the protective properties should account for the decrease in the concentration of free NO₂^- ions in the pores of cement concretes due to their adsorption. Here we applied the classical molecular dynamics computer simulation approach to quantitatively study the molecular scale mechanisms of nitrite adsorption from NaNO₂ aqueous solution on a portlandite surface. We used a new parameterization to model the hydrated NO₂^- ions in combination with the recently upgraded ClayFF force field (ClayFF-MOH) for the structure of portlandite. The new NO₂^- parameterization makes it possible to reproduce the properties of hydrated NO₂^- ions in good agreement with experimental data. In addition, the ClayFF-MOH model improves the description of the portlandite structure by explicitly taking into account the bending of Ca-O-H angles in the crystal and on its surface. The simulations showed that despite the formation of a well-structured water layer on the portlandite (001) crystal surface, NO₂^- ions can be strongly adsorbed. The nitrite adsorption is primarily due to the formation of hydrogen bonds between the structural hydroxyls on the portlandite surface and both the nitrogen and oxygen atoms of the NO₂^- ions. Due to that, the ions do not form surface adsorption complexes with a single well-defined structure but can assume various local coordinations. However, in all cases, the adsorbed ions did not show significant surface diffusional mobility. Moreover, we demonstrated that the nitrite ions can be adsorbed both near the previously-adsorbed hydrated Na⁺ ions as surface ion pairs, but also separately from the cations.

Capillary suction under unsaturated condition drives strong specific affinity of ions at the surface of clay mineral

Mineral-solution interface is of great importance in many soil and geochemical processes as well as industrial applications. Most relevant studies were based on saturated condition and given the corresponding theory, model, and mechanism. However, soils are usually in the non-saturation with different capillary suction. Our study presents substantially different scenery for ions interacting with mineral surface under unsaturated condition using molecular dynamics method. Under partially hydrated state, both cations (Ca²⁺) and anions (Cl^-) can be adsorbed as outer-sphere complexes at the montmorillonite surface, and the number significantly increased with the increase of unsaturated degree. Ions preferred to interact with clay mineral instead of water molecules under unsaturated state, and the mobility of both cations and anions substantially decreased with the increase of capillary suction as reflected by the diffusion coefficient analysis. Potential of mean force calculations further clearly revealed that the adsorption strength of both Ca²⁺ and Cl^- increased with capillary suction. Such an increase was more obvious for Cl^- compared to Ca²⁺, despite the adsorption strength of Cl^- was much weaker than Ca²⁺ at a certain capillary suction. Therefore, it is the capillary suction under unsaturated condition that drives the strong specific affinity of ions at the surface of clay mineral, which was tightly related to the steric effect of confined water film, the destruction of EDL structure, and the cation-anion pair interaction. This suggests that our common understanding of mineral-solution interaction should be largely improved.

Hydration Structure at the Calcite-Water (10.4) Interface in the Presence of Rubidium Chloride

Solid-liquid interfaces are of significant importance in a multitude of geochemical and technological fields. More specifically, the solvation structure plays a decisive role in the properties of the interfaces. Atomic force microscopy (AFM) has been used to resolve the interfacial hydration structure in the presence and absence of ions. Despite many studies investigating the calcite-water interface, the impact of ions on the hydration structure at this interface has rarely been studied. Here, we investigate the calcite-water interface at various concentrations (ranging from 0 to 5 M) of rubidium chloride (RbCl) using three-dimensional atomic force microscopy (3D AFM). We present molecularly resolved images of the hydration structure at the interface. Interestingly, the characteristic pattern of the hydration structure appears similar regardless of the RbCl concentration. The presence of the ions is detected in an indirect manner by more frequent contrast changes and slice displacements.

Adsorption mechanism of Pb2+ in montmorillonite nanopore under various temperatures and concentrations

Adsorption of lead (Pb²⁺) onto the montmorillonite (Mt) surface is one of the key approaches to remove Pb²⁺ in geological and environmental engineering. Temperature and initial Pb²⁺ concentration are two essential factors that influence the adsorption capacity of Mt on absorbing Pb²⁺. However, the nanoscale governing mechanism of temperature and initial concentration on Pb²⁺ adsorbing of Mt is still unclear. This research performed comprehensively molecular dynamics (MD) simulations to investigate how temperature and initial concentration affect the dynamic Pb²⁺ adsorption of Mt nanopore. The Pb²⁺ removal ratio shows a two-stage variation with the increase of initial Pb²⁺ concentration. Temperature controls the maximum initial Pb²⁺ concentration for complete Pb²⁺ removal by changing the maximum adsorption energy of Mt. Temperature also influences the maximum adsorption capacity and Pb²⁺ removal ratio of Mt nanopore indirectly by changing diffusion and hydration state of Pb²⁺. The initial Pb²⁺ concentration corresponding to the maximum adsorption energy coincides with the maximum initial Pb²⁺ concentration determined by the Pb²⁺ removal ratio. Lower adsorption energy and higher level of hydration and diffusion make Pb²⁺ absorbing on Mt surface become more difficult, reducing the Pb²⁺ adsorbing capacity of Mt. The initial Pb²⁺ concentration influences adsorption capacity and Pb²⁺ removal ratio not only via altering the quantity of Pb²⁺ but also through controlling the adsorption energy of Mt, as well as the diffusion and hydration state of Pb²⁺. With the increase of initial Pb²⁺ concentration, the hydration of Pb²⁺ is weakened while the adsorption energy of Mt and diffusion of Pb²⁺ are enhanced.

Understanding specific ion effects and the Hofmeister series

Specific ion effects (SIE), encompassing the Hofmeister Series, have been known for more than 130 years since Hofmeister and Lewith's foundational work. SIEs are ubiquitous and are observed across the medical, biological, chemical and industrial sciences. Nevertheless, no general predictive theory has yet been able to explain ion specificity across these fields; it remains impossible to predict when, how, and to what magnitude, a SIE will be observed. In part, this is due to the complexity of real systems in which ions, counterions, solvents and cosolutes all play varying roles, which give rise to anomalies and reversals in anticipated SIEs. Herein we review the historical explanations for SIE in water and the key ion properties that have been attributed to them. Systems where the Hofmeister series is perturbed or reversed are explored, as is the behaviour of ions at the liquid-vapour interface. We discuss SIEs in mixed electrolytes, nonaqueous solvents, and in highly concentrated electrolyte solutions - exciting frontiers in this field with particular relevance to biological and electrochemical applications. We conclude the perspective by summarising the challenges and opportunities facing this SIE research that highlight potential pathways towards a general predictive theory of SIE.

A Deep Look into the Dynamics of Saltwater Imbibition in a Calcite Nanochannel: Temperature Impacts Capillarity Regimes

This research concerns fundamentals of spontaneous transport of saltwater (1 mol·dm^-3 NaCl solution) in nanopores of calcium carbonates. A fully atomistic model was adopted to scrutinize the temperature dependence of flow regimes during solution transport under CaCO₃ nanoconfinement. The early time of capillary filling is inertia-dominated, and the solution penetrates with a near-planar meniscus at constant velocity. Following a transition period, the meniscus angle falls to a stabilized value, characterizing the capillary-viscous advancement in the calcite channel. At this stage, brine displacement follows a parabolic relationship consistent with the classical Lucas-Washburn (LW) theory. Approaching the slit outlet, the meniscus contact lines spread widely on the solid substrate and brine leaves the channel at a constant rate, in oppose to the LW law. The brine imbibition rate considerably increases at higher temperatures as a result of lower viscosity and greater tendency to form wetting layers on slit walls. We also pointed out a longer primary inertial regime and delayed onset of the viscous-capillary regime at higher temperatures. Throughout the whole span of capillary displacement, transport of sodium and chloride ions is tied to dynamics and diffusion of the water phase, even at the mineral interface. The results presented in this study are of broad implications in diverse science and technological applications.

The Impact of Salinity on the Interfacial Structuring of an Aromatic Acid at the Calcite/Brine Interface: An Atomistic View on Low Salinity Effect

This study aims to elucidate the impact of salinity on the interactions governing the adsorption of polar aromatic oil compounds onto calcite. To this end, molecular dynamics simulations were employed to assess adsorption of a model polar organic molecule (deprotonated benzoic acid, benzoate) on the calcite surface in NaCl brines of different concentration levels, namely, deionized water (DW), low-salinity water (LS, 5000 ppm), and sea water (SW; 45,000 ppm). Calcite was found to be completely covered by several well-ordered water layers. The top hydration layer is very compact and prevents direct adsorption of benzoates onto the substrate. Instead, Na⁺ ions form a distinct positively charged layer by adhering on the calcite substrate through inner-sphere complexion mode. Cl^- ions mostly lodge on top of the adsorbed sodium cations, forming a negatively charged layer. The distribution of ions at the calcite/brine interface thus exhibits the features of an electrical double layer, composed of a Stern-like positive layer followed by a negative one. The positive charged layer attracts benzoates toward the surface. As such, the sodium ions attached onto the calcite can act as adsorption sites to connect benzoates to the surface. By increasing the salinity, more Na⁺ ions adsorb onto the calcite surface, and the density of benzoate molecules at the interface is enhanced as a result of more Na⁺ bridging ions. The monotonic salinity-dependent adsorption of benzoate molecules on calcite offers a mechanism driving additional oil recovery upon injection of diluted brine into subsurface carbonate reservoirs.

Water properties under nano-scale confinement

Water is the universal solvent and plays a critical role in all known geological and biological processes. Confining water in nano-scale domains, as encountered in sedimentary rocks, in biological, and in engineered systems, leads to the deviations in water’s physicochemical properties relative to those measured for the non-confined phase. In our comprehensive analysis, we demonstrate that nano-scale confinement leads to the decrease in the melting/freezing point temperature, density, and surface tension of confined water. With increasing degree of spatial confinement the population of networked water, as evidenced by alterations in the O-H stretching modes, increases. These analyses were performed on two groups of mesoporous silica materials, which allows to separate pore size effects from surface chemistry effects. The observed systematic effects of nano-scale confinement on the physical properties of water are driven by alterations to water’s hydrogen-bonding network—influenced by water interactions with the silica surface — and has implications for how we understand the chemical and physical properties of liquids confined in porous materials.

Atomistic insights into cesium chloride solution transport through the ultra-confined calcium–silicate–hydrate channel

A new capillary transport model is proposed by modifying the original Lucas–Washburn function.

Structure, orientation, and dynamics of water-soluble ions adsorbed to basal surfaces of calcium monosulfoaluminate hydrates

Transport of water molecules and chloride ions in nanopores of hydrated cement paste (HCP) is proven to adversely affect the long-term durability of reinforced concrete structures exposed to seawater or deicing salts. The resistance against chloride attack is primarily associated with the chloride binding capacity of the main HCP constituents. Experimental tests revealed that AFm phases of HCP play a central role in binding the chloride ions. However, many aspects of AFm-solution interactions were largely unknown, especially at their interfaces. This was the motivation of the current study, in which the atomistic processes underlying the transport of water-soluble ions are investigated in detail using the classical molecular dynamics (MD) method. To this end, an aqueous layer containing various concentrations of sodium chloride solution is sandwiched between two basal surfaces of calcium monosulfoaluminate hydrate, which is the most abundant phase of AFm. The adsorption mechanisms of water molecules and diffusing ions are then characterized for inner- and outer-sphere distance ranges from the basal surfaces of monosulfoaluminate. It is found that the self-diffusion coefficient of the chloride and sodium ions present in the outer-sphere range are 83% and 47% larger than those residing in the inner-sphere range. With increasing the distance from the solid surface, an increase in the self-diffusion coefficient is captured. This increase in mobility is larger for chloride ions than sodium ions. This can be understood based on the observation that the inner- and outer-sphere complex formation are the governing adsorption mechanisms for the chloride and sodium ions, respectively.

Atomic-scale investigation of physical adsorption of water molecules and aggressive ions to ettringite's surfaces

The strength and durability of cementitious composite materials are adversely affected by the ingress of water molecules and aggressive ions into their intrinsic meso- and nano-pore spaces. Among various phases of hydrated cement paste (HCP), aluminum-rich phases play an important role in controlling the diffusivity of aqueous solutions, which can contain aggressive ions. To this date, however, there has been no systematic study to understand the adsorption mechanisms and chloride binding capacity of the aluminum-rich phases of HCP. This research gap has been the motivation of the current study to investigate the physical adsorption characteristics of ettringite as the main aluminum-rich phase of HCP and the primary hydrated product of calcium sulfoaluminate cement. Through a set of Molecular Dynamics simulations supported by macro-scale experimental tests, a fundamental insight into the molecular origins of the diffusion of water molecules, as well as sodium and chloride ions, in contact with ettringite is provided. As the primary objective of this study is to evaluate the transport properties at and near solution/solid interfaces, the molecular adsorption mechanisms are characterized for inner- and outer-sphere distances from the solid substrate. With an in-depth understanding of the structure and dynamics of water molecules and aggressive ions in contact with ettringite's surfaces, the outcome of this study provides reliable measures of physical adsorption, binding capacity, and self-diffusion coefficient, which can be further employed to introduce strategies to avoid the degradation of a wide variety of cementitious materials exposed to harsh environmental conditions.

Influence of aluminates on the structure and dynamics of water and ions in the nanometer channel of calcium silicate hydrate (C–S–H) gel

Al species incorporated in silicate chains enhance hydrophilicity and cation immobilization ability of the C–S–H gel.