Molecular dynamics study of the swelling patterns of Na/Cs-, Na/Mg-montmorillonites and hydration of interlayer cations

Montmorillonite Membranes with Tunable Ion Transport by Controlling Interlayer Spacing

Membranes incorporating two-dimensional (2D) materials have shown great potential for water purification and energy storage and conversion applications. Their ordered interlayer galleries can be modified for their tunable chemical and structural properties. Montmorillonite (MMT) is an earth-abundant phyllosilicate mineral that can be exfoliated into 2D flakes and reassembled into membranes. However, the poor water stability and random interlayer spacing of MMT caused by weak interlamellar interactions pose challenges for practical membrane applications. Herein, we demonstrate a facile approach to fabricating 2D MMT membranes with alkanediamines as cross-linkers. The incorporation of diamine molecules of different lengths enables controllable interlayer spacing and strengthens interlamellar connections, leading to tunable ion transport properties and boosted membrane stability in aqueous environments.

Interaction between Hydrated Smectite Clay Particles as a Function of Salinity (0-1 M) and Counterion Type (Na, K, Ca)

Swelling clay minerals control the hydrologic and mechanical properties of many soils, sediments, and sedimentary rocks. This important and well-known phenomenon remains challenging to predict because it emerges from complex multiscale couplings between aqueous chemistry and colloidal interaction mechanics in nanoporous clay assemblages, for which predictive models remain elusive. In particular, the predominant theory of colloidal interactions across fluid films, the widely used Derjaguin-Landau-Verwey-Overbeek model, fails to predict the ubiquitous existence of stable swelling states at interparticle distances below 3 nm that are stabilized by specific inter-atomic interactions in overlapping electrical double layers between the charged clay surfaces. Atomistic simulations have the potential to generate detailed insights into the mechanisms of these interactions. Recently, we developed a metadynamics-based molecular dynamics simulation methodology that can predict the free energy of interaction between parallel smectite clay particles in a wide range of interparticle distances (from 0.3 to 3 nm) and salinities (from 0.0 to 1.0 M NaCl). Here, we extend this work by characterizing the sensitivity of interparticle interactions to counterion type (Na, K, Ca). We establish a detailed picture of the free energy of interaction of parallel clay particles across water films as the sum of five interaction mechanisms with different sensitivities to salinity, counterion type, and interparticle distance.

Preparation of Montmorillonite Nanosheets with a High Aspect Ratio through Heating/Rehydrating and Gas-Pushing Exfoliation

Montmorillonite (MMT) is an abundant silicate mineral with ultrahigh stability. The exfoliation of stacked MMT into high-aspect-ratio nanosheets is of crucial importance for various applications such as toxic gas suppression, barrier property enhancement, flame retardancy, and ion conduction. In this work, we develop a new heating/rehydrating and gas-pushing method that can successfully exfoliate MMT into nanosheets with aspect ratios (600-5000) far higher than the currently reported values (1-120). The MMT first goes through a "starvation pretreatment" under different heating temperatures to improve its hydrophilicity and is then rehydrated in a hydrogen peroxide solution. The hydrogen peroxide in the MMT interlayer space can decompose into water and oxygen bubbles, thus finally leading to the exfoliation via gas-pushing while preserving the large lateral size (mainly in the range of 1-6 μm) of the nanosheets. By changing the pretreatment temperature and pH value of the hydrogen peroxide solution, the exfoliation performance can be tuned. This simple and low-cost exfoliation method is promising to achieve the mass production of MMT nanosheets with a high aspect ratio and may promote its application in various fields such as energy conversion, drug delivery, and photocatalysis.

Effect of Temperature on the Diffusion and Sorption of Cations in Clay Vermiculite

The MD method for modeling vermiculite containing Na⁺, Rb⁺, Cs⁺, Mg²⁺, and Ba²⁺ cations shows the following: With a weak swelling of clay, the temperature has no significant effect on the diffusion of water and cations through vermiculite. With a high content of water in vermiculite, the effect of temperature on the diffusion coefficient of water is greater than that of cations. We studied the structure of RDF ions in Na⁺-vermiculite, in which some of the cations are replaced by Rb⁺, Cs⁺, Mg²⁺, and Ba²⁺. Cations of alkali and alkaline earth metals compete with Na⁺ ions for adsorption sites on the surface of the clay layer. The alkaline earth metal cations are in the middle between the clay layers due to their higher charge and stronger hydration. In this case, Na⁺ is localized at the surface of the clay layer. Thus, cations of alkaline earth metals have little effect on the temperature dependence of the diffusion coefficient Na⁺.

Dissociation mechanism of methane hydrate by CaCl2: an experimental and molecular dynamics study

The formation of gas hydrate is a serious threat to the safe and effective completion of deepwater drilling and transportation operations, although it is considered as a potential energy resource. The inorganic salts are generally used as thermodynamic inhibitors; CaCl₂ as a common additive in drilling fluids exhibits unique properties. In this study, we explored the dissociation mechanism of CH₄ hydrate in CaCl₂ solutions at the macroscopic and microscopic scale using experiment and molecular dynamics (MD) simulation. The experimental results showed that CaCl₂ accelerated the dissociation rate of CH₄ hydrate. The dissociation rate of CH₄ hydrate increased with the increase of CaCl₂ concentration at large depressurization pressure and was mainly affected by pressure when the depressurization pressure was lower. MD simulations were used to give an atomic scale interpretation of the macroscopic results obtained from the experiment. The results showed that the addition of CaCl₂ destroyed the resistance liquid film formed during CH₄ hydrate dissociation, thus accelerating the dissociation process, in good agreement with experimental results. HIGHLIGHTS: • The amount of CaCl₂ affects CH₄ hydrate dissociation at large depressurization pressure. • The dissociation of CH₄ hydrate at low depressurization pressure is dependent on pressure. • Ca²⁺ destroys effectively the resistance liquid film produced during hydrate dissociation. • MD simulation results are in agreement with those of the experiment.

Molecular Dynamics Simulation of Clay Hydration Inhibition of Deep Shale

In the process of the exploitation of deep oil and gas resources, shale wellbore stability control faces great challenges under complex temperature and pressure conditions. It is difficult to reflect the micro mechanism and process of the action of inorganic salt on shale hydration with the traditional experimental evaluation technology on the macro effect of restraining shale hydration. Aiming at the characteristics of clay minerals of deep shale, the molecular dynamics models of four typical cations (K+, NH4+, Cs+ and Ca2+) inhibiting the hydration of clay minerals have been established by the use of the molecular dynamics simulation method. Moreover, the micro dynamics mechanism of typical inorganic cations inhibiting the hydration of clay minerals has been systematically evaluated, as has the law of cation hydration inhibition performance in response to temperature, pressure and ion type. The research indicates that the cations can promote the contraction of interlayer spacing, compress fluid intrusion channels, reduce the intrusion ability of water molecules, increase the negative charge balance ability and reduce the interlayer electrostatic repulsion force. With the increase in temperature, the inhibition of the cations on montmorillonite hydration is weakened, while the effect of pressure is opposite. Through the molecular dynamics simulation under different temperatures and pressures, we can systematically understand the microcosmic dynamics mechanism of restraining the hydration of clay in deep shale and provide theoretical guidance for the microcosmic control of clay hydration.

Molecular Dynamics Simulation of Ion Adsorption and Ligand Exchange on an Orthoclase Surface

Orthoclase (K-feldspar) is one of the natural inorganic materials, which shows remarkable potential toward removing heavy metal ions from aqueous solutions. Understanding the interactions of the orthoclase and metal ions is important in the treatment of saline wastewater. In this paper, molecular dynamics simulations were used to prove the adsorption of different ions onto orthoclase. The adsorption isotherms show that orthoclase has remarkable efficiency in the removal of cations at low ion concentrations. Aluminol groups are the preferential adsorption sites of cations due to higher negative charges. The adsorption types and adsorption sites are influenced by the valence, radius, and hydration stability of ions. Monovalent cations can be adsorbed in the cavities, whereas divalent cations cannot. The hydrated cation may form an outer-sphere complex or an inner-sphere complex in association with the loss of hydration water. Na+, K+, and Ca2+ ions mainly undergo inner-sphere adsorption and Mg2+ ions prefer outer-sphere adsorption. On the basis of simulation results, the mechanism of ion removal in the presence of orthoclase is demonstrated at a molecular level.

Hydration and swelling: a theoretical investigation on the cooperativity effect of H-bonding interactions between p-hydroxy hydroxymethyl calix[4]/[5]arene and H2O by many-body interaction and density functional reactivity theory

In order to explore the nature of the hydration and swelling of superabsorbent resin, a theoretical investigation into the cooperativity effect of the H-bonding interactions in the hydrates of four model compounds that can be regarded as the units of hydroquinone formaldehyde resin (HFR) (i.e., O-hydroxymethyl-1,4-dihydroxybenzene, methylene di-O-hydroxymethyl-1,4-dihydroxybenzene, p-hydroxy hydroxymethyl calix[4]arene and p-hydroxy hydroxymethyl calix[5]arene) was carried out by many-body interaction and density functional reactivity theory. The HFR···H₂O···H₂O complexes, in which the H₂O···H₂O moieties are bound with both the hydroxyl groups of HFR, are the most stable. For the HFR(H₂O)_n clusters, the interaction energy per building block is increased as the number of the size n increases, indicating the cooperativity effect. Therefore, a deduction is given that the cooperativity effects of the H-bonding interactions play an important role in the process of the hydration and swelling of HFR, and the swelling behavior is mainly attributed to the cooperativity effects which arised from the interactions between the H₂O molecules. The origin of the cooperativity effect was examined employing several information-theoretic quantities in the density functional reactivity theory. The degree of swelling of HFR was quantitated using a measure of volume. Graphical abstract.

Effect of Layer Charge Density on Hydration Properties of Montmorillonite: Molecular Dynamics Simulation and Experimental Study

Four kinds of Ca-montmorillonite with different layer charge density were used to study the effect of charge density on their hydration properties by molecular dynamics simulation and experiments. The research results of Z-density distribution of water molecules, H_w (hydrogen in water molecules), and Ca in the interlayer of montmorillonite show that the hydration properties of montmorillonite are closely related to its layer charge density. If the charge density is low, the water molecules in the interlayers are mainly concentrated on the sides of the central axis about –1.3 Å and 1.5 Å. As the charge density increases from 0.38_semi-cell to 0.69_semi-cell, the water molecules are distributed −2.5 Å and 2.4 Å away from the siloxane surface (Si-O), the concentration of water molecules near the central axis decreases, and at the same time, Ca²⁺ appears to gradually shift from the vicinity of the central axis to the Si-O surface on both sides in the montmorillonite layer. The simulation results of the radial distribution function (RDF) of the Ca-H_w, Ca-O_w (oxygen in water molecules), and Ca-O_t (the oxygen in the tetrahedron) show that the Ca²⁺ and O_w are more tightly packed together than that of H_w; with the increase of the charge density, due to the fact that the negative charge sites on the Si-O surface increase, under the action of electrostatic attraction, some of the Ca²⁺ are pulled towards the Si-O surface, which is more obvious when the layer charge density of the montmorillonite is higher. The results of the RDF of the O_t-H_w show that with the increase of charge density, the number of hydrogen bonds formed by O_t and H_w in the interlayers increase, and under the action of hydrogen bonding force, the water molecules near the central axis are pulled towards the two sides of Si-O surface. As a result, the arrangement of water molecules is more compact, and the structure is obvious. Correspondingly, the self-diffusion coefficient shows that the higher the layer charge density, the lower the self-diffusion coefficient of water molecules in interlayers is and the worse the hydration performance of montmorillonite. The experimental results of the experiments fit well with the above simulation results.

Effect of Layer Charge Characteristics on the Distribution Characteristics of H2O and Ca2+ in Ca-Montmorillonites Interlayer Space: Molecular Dynamics Simulation

The charge characteristics of montmorillonite have significant effects on its hydration and application performances. In this study, a molecular dynamics simulation method was used to study the influence of the charge position and charge density of montmorillonite on the distribution of H2O and Ca2+ in layers. The results showed that when the layer charge is mainly derived from the substitution among ions in the tetrahedron, a large number of Hw and Ot are combined into a hydrogen bond in the interlayer, thus the water molecules are more compactly arranged and the diffusion of water molecules among the layers is reduced. In addition, the Ca2+ are diffused to the sides by a concentrated distribution in the central axis of the layer. As the charge density of the montmorillonite increases, the polarity of the Si-O surface increases, which lesds to the deterioration of the diffusibility of the water molecules and the structure of the water molecules in the interlayers is more stable. The increase in the layer charge density lesds to the expansion of the isomorphic substitution range of the crystal structure, which results in a more dispersed distribution of Ca2+ among the layers under the action of electrostatic attraction between the substituted negative sites and the Ca2+.