A ReaxFF Potential for Modeling Organic Matter Degradation with Oxybromine Oxidants.

Oxidation of organic matter with oxybromine oxidants is ushering in a new era of enhanced hydrocarbon recovery. While these potent reagents are being tested in laboratory and field experiments, there is a pressing demand to delineate the molecular processes governing oxidation reactions at geological depth. Here, we parameterize a ReaxFF potential to model the oxidative decompositions of aliphatic and aromatic hydrocarbons in the presence of water-NaBr solutions that contain oxybromine (BrOn)- oxidizers. Our parameterization results in a reliable empirical bond-order potential that accurately calculates bond energies, exhibiting an RMSE of ∼1.18 eV, corresponding to 1.36 % average error. Reproducing bond dissociation and binding energies from Density Functional Theory (DFT), our parameterization proves transferable to aqueous environments. This H/C/O/Na/Br ReaxFF potential accurately reproduces the oxidation pathways of small hydrocarbons with oxybromine oxidizers. This force field captures proton and oxygen transfer, C-C bond tautomerization, and cleavage, leading to ring-opening and chain fragmentation. Molecular dynamic simulations demonstrate the oxidative degradation of aromatic and aliphatic kerogen-like moieties in bulk solutions. We envision that such reactive force fields will be useful to understand better the oxidation reactions of organic matter formed in geological reservoirs for enhanced shale gas recovery and improved carbon dioxide treatments.

Reactive force fields for aqueous and interfacial magnesium carbonate formation.

We develop Mg/C/O/H ReaxFF parameter sets for two environments: an aqueous force field for magnesium ions in solution and an interfacial force field for minerals and mineral-water interfaces. Since magnesium is highly ionic, we choose to fix the magnesium charge and model its interaction with C/O/H through Coulomb, Lennard-Jones, and Buckingham potentials. We parameterize the forcefields against several crystal structures, including brucite, magnesite, magnesia, magnesium hydride, and magnesium carbide, as well as Mg2+ water binding energies for the aqueous forcefield. Then, we test the forcefield for other magnesium-containing crystals, solvent separated and contact ion-pairs and single-molecule/multilayer water adsorption energies on mineral surfaces. We also apply the forcefield to the forsterite-water and brucite-water interface that contains a bicarbonate ion. We observe that a long-range proton transfer mechanism deprotonates the bicarbonate ion to carbonate at the interface. Free energy calculations show that carbonate can attach to the magnesium surface with an energy barrier of about 0.22 eV, consistent with the free energy required for aqueous Mg-CO3 ion pairing. Also, the diffusion constant of the hydroxide ions in the water layers formed on the forsterite surface are shown to be anisotropic and heterogeneous. These findings can help explain the experimentally observed fast nucleation and growth of magnesite at low temperature at the mineral-water-CO2 interface in water-poor conditions.

Densification of the interlayer spacing governs the nanomechanical properties of calcium-silicate-hydrate.

Calciuam-silicate-hydrate (C-S-H) is the principal binding phase in modern concrete. Molecular simulations imply that its nanoscale stiffness is 'defect-driven', i.e., dominated by crystallographic defects such as bridging site vacancies in its silicate chains. However, experimental validation of this result is difficult due to the hierarchically porous nature of C-S-H down to nanometers. Here, we integrate high pressure X-ray diffraction and atomistic simulations to correlate the anisotropic deformation of nanocrystalline C-S-H to its atomic-scale structure, which is changed by varying the Ca-to-Si molar ratio. Contrary to the 'defect-driven' hypothesis, we clearly observe stiffening of C-S-H with increasing Ca/Si in the range 0.8 ≤ Ca/Si ≤ 1.3, despite increasing numbers of vacancies in its silicate chains. The deformation of these chains along the b-axis occurs mainly through tilting of the Si-O-Si dihedral angle rather than shortening of the Si-O bond, and consequently there is no correlation between the incompressibilities of the a- and b-axes and the Ca/Si. On the contrary, the intrinsic stiffness of C-S-H solid is inversely correlated with the thickness of its interlayer space. This work provides direct experimental evidence to conduct more realistic modelling of C-S-H-based cementitious material.

Rigidity transition in materials: hardness is driven by weak atomic constraints.

Understanding the composition dependence of the hardness in materials is of primary importance for infrastructures and handled devices. Stimulated by the need for stronger protective screens, topological constraint theory has recently been used to predict the hardness in glasses. Herein, we report that the concept of rigidity transition can be extended to a broader range of materials than just glass. We show that hardness depends linearly on the number of angular constraints, which, compared to radial interactions, constitute the weaker ones acting between the atoms. This leads to a predictive model for hardness, generally applicable to any crystalline or glassy material.

Anomalous composition-dependent dynamics of nanoconfined water in the interlayer of disordered calcium-silicates.

With shear interest in nanoporous materials, the ultraconfining interlayer spacing of calcium-silicate-hydrate (C-S-H) provides an excellent medium to study reactivity, structure, and dynamic properties of water. In this paper, we present how substrate composition affects chemo-physical properties of water in ultraconfined hydrophilic media. This is achieved by performing molecular dynamics simulation on a set of 150 realistic models with different compositions of calcium and silicon contents. It is demonstrated that the substrate chemistry directly affects the structural properties of water molecules. The motion of confined water shows a multi-stage dynamics which is characteristic of supercooled liquids and glassy phases. Inhomogeneity in that dynamics is used to differentiate between mobile and immobile water molecules. Furthermore, it is shown that the mobility of water molecules is composition-dependent. Similar to the pressure-driven self-diffusivity anomaly observed in bulk water, we report the first study on composition-driven diffusion anomaly, the self diffusivity increases with increasing confined water density in C-S-H. Such anomalous behavior is explained by the decrease in the typical activation energy required for a water molecule to escape its dynamical cage.