Supercritical Carbon Dioxide at Smectite Mineral–Water Interfaces: Molecular Dynamics and Adaptive Biasing Force Investigation of CO₂/H₂O Mixtures Nanoconfined in Na-Montmorillonite

The carbon dioxide (CO 2 ) retention capacity and adsorption/desorption energetics of layered nanoporous oxide materials depend critically on the hydration level and the nature of molecular interactions among H 2 O, CO 2, charge-balancing cations, and the oxide/hydroxide layers. Molecular-scale understanding of the structure, dynamics, and interfacial energetics of H 2 O/CO 2 binary mixtures confined in the interlayer nanopores is paramount to geological CO 2 storage efforts in clay-rich materials. This Article investigates the effects of supercritical CO 2 (scCO 2 ) in the hydrated interlayer galleries of the hydrophilic smectite mineral (Na-montmorillonite) under geochemically relevant conditions using classical molecular dynamics simulations and enhanced sampling free energy methods. For the compositions investigated, the interactions among the cations, intercalated fluid species, and the basal surfaces result in structures with H 2 O and CO 2 coexisting in a single layer at the center of the interlayer. The water molecules in this central H 2 O/CO 2 layer cluster around and hydrate Na + ions desorbed from the basal surfaces, whereas CO 2 –CO 2 hydrophobic interactions favor mutual clustering of CO 2 molecules. This arrangement results in dynamic percolation paths that facilitate single file-like anisotropic lateral diffusion of CO 2 . The water clusters around the Na + ions act as two-dimensional nanopores for the diffusion of Na + between the basal surfaces and across the central H 2 O/CO 2 layer, whereas the CO 2 -rich regions are not permeable to Na + . The near-surface Na + ions occur in two distinct types of coordination environments with distinct NMR spectral fingerprints. Type-I near-surface Na + ions are coordinated by two basal oxygen atoms and four water molecules, whereas for type-II one of the coordinating water molecules is replaced by a CO 2 molecule. The activation energies for a H 2 O and a CO 2 molecule to move out of the first coordination shell of a near-surface Na + are ∼2.75 and ∼0.5 kcal/mol, respectively. The activation barriers for site-hopping of a H 2 O molecule within the first coordination shell of near-surface and displaced Na + ions are ∼1.6 kcal/mol whereas those for site-hopping of CO 2 around the near-surface and displaced Na + ions are ∼1.8 and ∼3.5 kcal/mol, respectively. The results provide a detailed picture of the interlayer structure and energetics of diffusional motion of cations and intercalates.