Search

Models predicting wettability alteration of mineral-brine-oil interfaces during low-salinity-waterflooding (LSW) should account for the elevated temperatures typically found in oil reservoirs. For the first time, high temperature ζ-potential (zeta potential) data for silica are collected and used to interpret surface chemistries and interactions at reservoir-like conditions to predict temperature’s effect on wettability alteration. Mobility data for amorphous silica in varying NaCl(aq) concentrations at 25, 100, and 150 °C and neutral pH were obtained through microelectrophoresis experiments. Calculated ζ-potentials were fit with surface complexation model (SCM) parameters to predict electrical double layer (EDL) parameters based upon the Gouy-Chapman-Stern-Grahame (GCSG) model. ζ-potentials increased with increasing temperature (around 50% increase from 25 to 150 °C) and decreasing NaCl concentrations (10−1–10−4 mol kg−1). These trends, along with Derjaguin-Verwey-Landau-Overbeek (DLVO) theory, suggests that overall repulsive forces extend farther from the surface at low salinity and higher temperatures, implying greater wetting thickness/surface wettability in these environments. The resulting surface concentration calculations suggest that LSW is most impactful up to 10−2 mol kg−1 of salt, and that additional dilution below 10−3 mol kg−1 will negligibly impact oil recovery, particularly at reservoir temperatures above 100 °C. The analysis provides a framework for treating more complex reservoir systems, such as carbonates in multivalent brines.

Detailed analysis of the adsorption of oxalic acid ions, that is, oxalate and hydrogenoxalate, on the rutile (110) surface was carried out using molecular dynamics augmented by free energy calculations and supported by ab initio calculations. The predicted adsorption on perfect nonhydroxylated and hydroxylated surfaces with surface charge density from neutral to +0.208 C/m2 corresponding to pH values of about 6 and 3.7, respectively, agrees with experimental adsorption data and charge-distribution multisite ion complexation model predictions obtained using the most favorable surface complexes identified in our simulations. We found that outer-sphere complexes are the most favorable, owing to strong hydrogen binding of oxalic acid ions with surface hydroxyls and physisorbed water. The monodentate complex, the most stable among inner-sphere complexes, was about 15 kJ/mol higher in energy, but separated by a large energy barrier. Other inner-sphere complexes, including some previously suggested in the literature as likely adsorption structures such as bidentate and chelate complexes, were found to be unstable both by classical and by ab initio modeling. Both the surfaces and (hydrogen)oxalate ions were modeled using charges scaled to 75% of the nominal values in accord with the electronic continuum theory and our earlier parameterization of (hydrogen)oxalate ions, which showed that nominal charges exaggerate ion-water interactions.

A multistart optimization algorithm for surface complexation equilibrium parameters (MUSE) was applied to a large and diverse data set for chromate adsorption on iron (oxy)hydroxides (ferrihydrite and goethite). Within the Basic Stern and the charge-distribution multisite complexation (CD-MUSIC) framework, chromate binding constants and the Stern Layer capacitance were optimized simultaneously to develop a consistent parameter set for surface complexation models. This analysis resulted in three main conclusions regarding the model parameters: (a) There is no single set of parameter values that describes such diverse data sets when modeled independently. (b) Parameter differences among the data sets are mainly due to different amounts of total sites, i.e., surface area and surface coverages, rather than structural differences between the iron (oxy)hydroxides. (c) Unified equilibrium constants can be extracted if total site dependencies are taken into account. The implementation of the MUSE algorithm automated the process of optimizing the parameters in an objective and consistent manner and facilitated the extraction of predictive relationships for unified equilibrium constants. The extracted unified parameters can be implemented in reactive transport modeling in the field by either adopting the appropriate values for each surface coverage or by estimating error bounds for different conditions. The evaluation of a forward model with unified parameters successfully predicted chromate adsorption for a range of capacitance values.

The MUlti-start optimization algorithm for Surface complexation Equilibrium (MUSE) algorithm has been developed to optimize the fitting of thermodynamic constants for surface complexation modeling (SCM). Although there is a plethora of software to perform data fitting and determine intrinsic equilibrium constants, the algorithms used are highly dependent on initial values and choice of parameters. This limits their transferability to model other systems, for example, reactive transport processes. With this in mind, a hybridized optimization approach, based on a multistart algorithm combined with a local optimizer, has been developed to allow the simultaneous optimization of SCM parameters and to assess the sensitivity of these parameters to changes in the model assumptions. In this study, the CD–MUSIC formalism with a Basic Stern electrostatic model is utilized to model chromate adsorption on ferrihydrite, although the MUSE algorithm can be applied to any adsorption data set and be implemented in any mode...

Cations in natural solutions significantly impact interfacial processes, particularly dissolution and surface charge measurements for quartz and silica, which are amongst the most naturally abundant and technologically important solids. Thermodynamic parameters for cation-specific interfacial reactions have heretofore been mostly derived instead of directly measured experimentally. This work investigates the energetics of adsorption and exchange reactions of alkali metal (M+) and alkaline earth (M2+) cations with the quartz surface by flow adsorption microcalorimetry, in tandem with in-situ pH measurements. The magnitudes of the heats of adsorption and exchange were found to increase along the Hofmeister series i.e., Li+

Ferrihydrite (FH) plays an important role in controlling the fate and transport of many compounds in nature due to its large surface area and high reactivity. This study is the first attempt to build a surface complexation model using the recently proposed surface structure that incorporates tetrahedrally coordinated Fe atoms (Hiemstra, 2013). The ability of the model to describe the surface charge curves of FH with different preparation methods and Points of Zero Net Proton Charge (PZNPC) is tested. In general, FH particles that have been subject to aging are larger and have lower specific surface area (SSA) and higher PZNPCs. The structural model includes 2 types of singly coordinated (SC) oxygens that are present only on the (1-11) and (1-10) faces and 5 types of triply coordinated (TC) oxygens that are also present on the basal planes (001) and (00-1), for a total of 11 sites. The 11 - site model was able to simulate fresh FH datasets with PZNPC lower than 8.5, but could only simulate higher PZNPCs when the contribution of the more acidic basal planes was minimized. The available microscopic observations do not support this condition, which suggests TC groups on the basal planes likely have log K values higher than the macroscopic PZNPC. We attempted to test this hypothesis through three versions of simplified 3-site models, using SC and one TC on (1-10) and (1-11), with log K 8.0 (equal to fresh FH PZNPC) and one TC group on the basal planes with log K 9.5. This enables fitting of the PZNPC of aged FH datasets by adjusting the face contributions. An unresolved issue is whether this model accurately describes the relative contribution of SC and TC sites to the overall charge, which has implications for accurate description of specific ion adsorption.

Flow microcalorimetry was used to investigate the energetics associated with Rb+, K+, Na+, Cl–, and NO3– exchange at the rutile–water interface. Heats of exchange reflected differences in bulk hydration/dehydration enthalpies (Na+ > K+ > Rb+, and Cl– > NO3–) such that exchanging Na+ or Cl– from the surface was exothermic, reflecting their greater bulk hydration enthalpies. Exchange heats were measured at pH 2, 3.25, 5.8, and 11 and exhibited considerable differences as well as pH dependence. These trends were rationalized with the aid of a molecularly constrained surface complexation model (SCM) that incorporated the inner-sphere binding observed for the cations on the rutile (110) surface. Explicitly accounting for the inner-sphere binding configuration differences between Rb+, K+, and Na+, as well as accompanying differences in negative surface charge development, resulted in much better agreement with measured exchange ratios than by considering bulk hydration enthalpies alone. The observation that cal...

Here, we characterize oxalate adsorption by rutile in NaCl media (0.03 and 0.30 m) and between pH 3 and 10 over a wide temperature range which includes the near hydrothermal regime (10-150 °C). Oxalate adsorption increases with decreasing pH (as is typical for anion binding by metal oxides), but systematic trends with respect to ionic strength or temperature are absent. Surface complexation modeling (SCM) following the CD-MUSIC formalism, and as constrained by molecular modeling simulations and IR spectroscopic results from the literature, is used to interpret the adsorption data. The molecular modeling simulations, which include molecular dynamics simulations supported by free-energy and ab initio calculations, reveal that oxalate binding is outer-sphere, albeit via strong hydrogen bonds. Conversely, previous IR spectroscopic results conclude that various types of inner-sphere complexes often predominate. SCMs constrained by both the molecular modeling results and the IR spectroscopic data were developed, and both fit the adsorption data equally well. We conjecture that the discrepancy between the molecular simulation and IR spectroscopic results is due to the nature of the rutile surfaces investigated, that is, the perfect (110) crystal faces for the molecular simulations and various rutile powders for the IR spectroscopy studies. Although the (110) surface plane is most often dominant for rutile powders, a variety of steps, kinks, and other types of surface defects are also invariably present. Hence, we speculate that surface defect sites may be primarily responsible for inner-sphere oxalate adsorption, although further study is necessary to prove or disprove this hypothesis.

The zeta potential (ZP) is an oft-reported measure of the macroscopic charge state of solid surfaces and colloidal particles in contact with solvents. However, the origin of this readily measurable parameter has remained divorced from the molecular-level processes governing the underlying electrokinetic phenomena, which limits its usefulness. Here, we connect the macroscopic measure to the microscopic realm through nonequilibrium molecular dynamics simulations of electroosmotic flow between parallel slabs of the hydroxylated (110) rutile (TiO2) surface. These simulations provided streaming mobilities, which were converted to ZP via the commonly used Helmholtz-Smoluchowski equation. A range of rutile surface charge densities (0.1 to −0.4 C/m2), corresponding to pH values between about 2.8 and 9.4, in RbCl, NaCl, and SrCl2 aqueous solutions, were modeled and compared to experimental ZPs for TiO2 particle suspensions. Simulated ZPs qualitatively agree with experiment and show that “anomalous” ZP values and i...

A comprehensive set of molecular-level results, primarily from classical molecular dynamics (CMD) simulations, are used to constrain CD-MUSIC surface complexation model (SCM) parameters describing rutile powder titrations conducted in RbCl, NaCl, and NaTr (Tr = triflate, CF3SO3–) electrolyte media from 25 to 250 °C. Rb+ primarily occupies the innermost tetradentate binding site on the rutile (110) surface at all temperatures (25, 150, 250 °C) and negative charge conditions (−0.1 and −0.2 C/m2) probed via CMD simulations, reflecting the small hydration energy of this large, monovalent cation. Consequently, variable SCM parameters (Stern-layer capacitance values and intrinsic Rb+ binding constants) were adjusted relatively easily to satisfactorily match the CMD and titration data. The larger hydration energy of Na+ results in a more complex inner-sphere distribution, which shifts from bidentate to tetradentate binding with increasing negative charge and temperature, and this distribution was not matched wel...

The effect of particle size on the adsorption of Sr(2+) onto monodisperse nanometer diameter (4, 20, and 40 nm) anatase samples has been evaluated quantitatively with macroscopic experimental studies. The adsorption of Sr(2+) onto the anatase particles was evaluated by potentiometric titrations in NaCl media, at two ionic strengths (0.03 and 0.3 m), and over a wide range of pH (3-11) and surface loadings, at a temperature of 25 °C. Adsorption of Sr(2+) to the surface of the 20 and 40 nm diameter samples was similar, whereas the Sr(2+) adsorption titration curves were shallower for the 4 nm diameter samples. At high pH, the smallest particles adsorbed slightly less Sr(2+) than was adsorbed by the larger particles. At the molecular scale, density functional theory (DFT) calculations were used to evaluate the most stable Sr(2+) surface species on the (101) anatase surface (the predominant crystal face). An inner-sphere Sr-tridentate surface species was found to be the most stable. The experimental data were described with a charge distribution (CD) and multisite complexation (MUSIC) model, with a Basic Stern layer description of the electric double layer. The resulting surface complexation model explicitly incorporated the molecular-scale information from the DFT simulation results. For 20 and 40 nm diameter anatase, the CD value for the Sr-tridentate species was calculated using a bond valence interpretation of the DFT-optimized geometry. The CD value for the 4 nm sample was smaller than that for the 20 and 40 nm samples, reflecting the shallower Sr(2+) adsorption titration curves. The adsorption differences between the smallest and larger anatase particles can be rationalized by water being more highly structured near the 4 nm anatase sample and/or the Sr-tridentate surface species may require more well-developed surface terraces than are present on the 4 nm particles.

Flow microcalorimetry was used to investigate the energetics associated with Rb