Unseeded, spontaneous nucleation of spherulitic magnesium calcite

Interface dissolution kinetics and porosity formation of calcite and dolomite (110) and (104) planes: An implication to the stability of geologic carbon sequestration

Geologic carbon sequestration (GCS) via injecting CO2 into deep carbonate reservoirs (mainly calcite and dolomite) is a promising strategy to reduce CO2 level. However, the dissolution or precipitation of calcite/dolomite planes on minerals/solution interface during long-term GCS process develops intergranular porosity and thus affects the permeability and stability of reservoirs. To investigate this process, both calcite and dolomite were dissolved in acetic and carbonic acids. A diffusion-controlled process was identified, with greater diffusion rates in acetic acid than that in carbonic acid. Quantified planes activity of both minerals follows (110) > (116) > (101) > (113) > (018) > (104) through density functional theory. Accomplished with preferential dissolution of calcite (110) planes in carbonic acid, calcite crystals precipitated with (104) planes at 423.15 K, under which, more calcite crystals were observed on dolomite surface, producing Ca-deplete surface. Molecular dynamic calculations showed higher dissolution rates of calcite/dolomite (110) planes than (104). In addition, the dissolution coefficients of Ca2+ were approximately triple of that Mg2+. Therefore, this study reveals the interface dissolution mechanisms of calcite and dolomite, especially on (110) and (104) planes at an atomic level, for the first time, providing better understanding for the stability of long-term GCS process.

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

HYPOTHESIS

The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). We showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers.

EXPERIMENTS

First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances.

FINDINGS

Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.