Search

Crystallisation of alkaline earth metal carbonates from water is important for biomineralisation and environmental geochemistry. Here large-scale computer simulations are a useful approach to complement experimental studies by providing atomistic insights and even by quantitatively determining the thermodynamics of individual steps. However, this is dependent on the existence of force field models that are sufficiently accurate while being computationally efficient enough to sample complex systems. Here we introduce a revised force field for aqueous alkaline earth metal carbonates that reproduces both the solubilities of the crystalline anhydrous minerals, as well as the hydration free energies of the ions. The model is also designed to run efficiently on GPUs thereby reducing the cost of such simulations. The performance of the revised force field is compared against previous results for important properties relevant to crystallisation, including ion-pairing and mineral-water interfacial structure and dynamics.

In this study, the adaption of the recently published molecular GFN-FF for periodic boundary conditions (pGFN-FF) is described through the use of neighbor lists combined with appropriate charge sums to handle any dimensionality from 1D polymers to 2D surfaces and 3D solids. Numerical integration over the Brillouin zone for the calculation of π bond orders of periodic fragments is also included. Aside from adapting the GFN-FF method to handle periodicity, improvements to the method are proposed in regard to the calculation of topological charges through the inclusion of a screened Coulomb term that leads to more physical charges and avoids a number of pathological cases. Short-range damping of three-body dispersion is also included to avoid collapse of some structures. Analytic second derivatives are also formulated with respect to both Cartesian and strain variables, including prescreening of terms to accelerate the dispersion/coordination number contribution to the Hessian. The modified pGFN-FF scheme is then applied to a wide range of different materials in order to examine how well this universal model performs.

Density functional theory (DFT) with semiempirical dispersion corrections (DFT-D2) has been used to calculate the binding energy of a CO2 molecule on the calcite {10.4} surface for different positions and orientations. This generated potential energy landscape was then used to parametrize a classical force field. From this, we used metadynamics (MTD) to derive free energy profiles at 300 and 350 K for CO2 binding to calcite, CO2 binding with Ca2+, and pairing of two CO2 molecules, all for aqueous conditions. We subsequently performed classical molecular dynamics (MD) simulations of CO2 and water on the {10.4} surface at pressures and temperatures relevant for CO2 geological storage. Density profiles show characteristic structured water layering at the calcite surface and two distinct phases of water and CO2. We have also calculated the densities of the CO2-rich and water-rich phases and thereby determined the mutual solubilities. For all the pressures and temperatures in the studied range, CO2 was unable ...

Assessment of the risks and environmental impacts of carbon geosequestration requires knowledge about the wetting behavior of mineral surfaces in the presence of CO2 and the pore fluids. In this context, the interfacial tension (IFT) between CO2 and the aqueous fluid and the contact angle, theta, 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 CO2–water IFT and the CO2–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 CO2 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 water 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 CO2 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 CO2. Our results suggest that the contact angle of a cylindrical water droplet on calcite {10.4}, in the presence of CO2, 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 degrees, indicating that calcite is strongly wetted by water.

Knowledge about the interaction between fluids and solids and the interfacial tension (IFT) that results is important for predicting behavior and properties in industrial systems and in nature, such as in rock formations before, during, and after CO2 injection for long-term storage. Many authors have studied the effect of the environmental variables on the IFT in the CO2–H2O system. However, experimental measurements above CO2 supercritical conditions are scarce and sometimes contradictory. Molecular modeling is a valuable tool for complementing experimental IFT determination, and it can help us interpret results and gain insight under conditions where experiments are difficult or impossible. Here, we report predictions for CO2–water interfacial tension performed using density functional theory (DFT) combined with the COSMO-RS implicit solvent model. We predicted the IFT dependence as a function of pressure (0–50 MPa), temperature (273–383 K), and salinity (0–5 M NaCl). The results agree well with literat...

Hetero-nanoparticles represent an important family of composite nanomaterials that in the past years are attracting ever-growing interest. Here, we report a new strategy for the synthesis of water dispersible cobalt ferrite nano particles (CoxFe3-xO4 NPs) decorated with ultrasmall (2-3 nm). gold nanoparticles (Au NPs). The synthetic procedure is based on the use of 2,3-meso-dimercaptosuccinic acid (DMSA), which plays a double role. First, it transfers cobalt ferrite NPs from the organic phase to aqueous media. Second, the DMSA reductive power promotes the in situ nucleation of gold NPs in proximity of the magnetic NP surface. Following this procedure, we achieved a water dispersible nanosystem (CoxFe3-xO4-DMSA-Au NPs) which combines the cobalt ferrite magnetic properties with the catalytic features of ultrasmall Au NPs. We showed that CoxFe3-xO4-DMSA-Au NPs act as an efficient nanocatalyst to reduce 4-nitrophenol to 4-aminophenol and that they can be magnetically recovered and recycled. It is noteworthy that such nanosystem is more catalytically active than Au NPs with equal size. Finally, a complete structural and chemical characterization of the hetero-NPs is provided.