RESUMEN

The connection between solution structure, particle forces, and emergent phenomena at solid-liquid interfaces remains ambiguous. In this case study on boehmite aggregation, we established a connection between interfacial solution structure, emerging hydration forces between two approaching particles, and the resulting structure and kinetics of particle aggregation. In contrast to expectations from continuum-based theories, we observed a nonmonotonic dependence of the aggregation rate on the concentration of sodium chloride, nitrate, or nitrite, decreasing by 15-fold in 4 molal compared to 1 molal solutions. These results are accompanied by an increase in repulsive hydration forces and interfacial oscillatory features from 0.27-0.31 nm in 0.01 molal to 0.38-0.52 nm in 2 molal. Moreover, molecular dynamics (MD) simulations indicated that these changes correspond to enhanced ion correlations near the interface and produced loosely bound aggregates that retain electrolyte between the particles. We anticipate that these results will enable the prediction of particle aggregation, attachment, and assembly, with broad relevance to interfacial phenomena.

RESUMEN

Electromagnetic fields are used in water treatment and desalination to regulate scale formation and extend the lifetime of membranes. External electric and magnetic fields can promote or suppress mineral nucleation and growth. However, the molecular-scale mechanisms of such processes remain unknown. Computing the free energies needed to form ion pairs under external fields provides important insights into understanding the elemental steps during the initial formation of mineral scales. In this paper, we used molecular dynamics combined with metadynamics simulations to investigate the free energies of forming the [Ca-CO3]0 ion pair, a fundamental building block of carbonate scales, under a range of magnetic (up to 10 T) and electric (up to 10 V m-1) fields in water. The presence of constant magnetic or electric fields favored the ion pairing reaction and lowered the free energies by up to 3% to 6%. The internal energy and entropic components of the free energy showed significant changes and exhibited non-linear behavior with increasing field strength. The [Ca-CO3]0 ion pairing is an entropy-driven process in the absence of an external field, but the mechanism shifts to an internal energy-driven process under selected external fields, suggesting possible changes in the nucleation pathways.

RESUMEN

It is established that the rates of solvent exchange at interfaces correlate with the rates of a number of mineral reactions, including growth, dissolution and ion sorption. To test if solvent exchange is limiting these rates, quasi-elastic neutron scattering (QENS) is used here to benchmark classical molecular dynamics (CMD) simulations of water bound to nanoparticulate calcite. Four distributions of solvent exchanges are found with residence times of 8.9 ps for water bound to calcium sites, 14 ps for that bound to carbonate sites and 16.7 and 85.1 ps for two bound waters in a shared calcium-carbonate conformation. By comparing rates and activation energies, it is found that solvent exchange limits reaction rates neither for growth nor dissolution, likely due to the necessity to form intermediate states during ion sorption. However, solvent exchange forms the ceiling for reaction rates and yields insight into more complex reaction pathways.

RESUMEN

Classical theories of particle aggregation, such as Derjaguin-Landau-Verwey-Overbeek (DLVO), do not explain recent observations of ion-specific effects or the complex concentration dependence for aggregation. Thus, here, we probe the molecular mechanisms by which selected alkali nitrate ions (Na+, K+, and NO3-) influence aggregation of the mineral boehmite (γ-AlOOH) nanoparticles. Nanoparticle aggregation was analyzed using classical molecular dynamics (CMD) simulations coupled with the metadynamics rare event approach for stoichiometric surface terminations of two boehmite crystal faces. Calculated free energy landscapes reveal how electrolyte ions alter aggregation on different crystal faces relative to pure water. Consistent with experimental observations, we find that adding an electrolyte significantly reduces the energy barrier for particle aggregation (∼3-4×). However, in this work, we show this is due to the ions disrupting interstitial water networks, and that aggregation between stoichiometric (010) basal-basal surfaces is more favorable than between (001) edge-edge surfaces (∼5-6×) due to the higher interfacial water densities on edge surfaces. The interfacial distances in the interlayer between aggregated particles with electrolytes (∼5-10 Å) are larger than those in pure water (a few Ångströms). Together, aggregation/disaggregation in salt solutions is predicted to be more reversible due to these lower energy barriers, but there is uncertainty on the magnitudes of the energies that lead to aggregation at the molecular scale. By analyzing the peak water densities of the first monolayer of interstitial water as a proxy for solvent ordering, we find that the extent of solvent ordering likely determines the structures of aggregated states as well as the energy barriers to move between them. The results suggest a path for developing a molecular-level basis to predict the synergies between ions and crystal faces that facilitate aggregation under given solution conditions. Such fundamental understanding could be applied extensively to the aggregation and precipitation utilization in the biological, pharmaceutical, materials design, environmental remediation, and geological regimes.

RESUMEN

HYPOTHESIS: Understanding the mechanisms of proton transfer on quartz surfaces in water is critical for a range of processes in geochemical, environmental, and materials sciences. The wide range of surface acidities (>9 pKa units) found on the ubiquitous mineral quartz is caused by the structural variations of surface silanol groups. Molecular scale simulations provide essential tools for elucidating the origin of site-specific surface acidities. SIMULATIONS: We used density-functional tight-binding-based molecular dynamics combined with rare-event metadynamics simulations to probe the mechanisms of deprotonation reactions from ten representative surface silanol groups found on both pristine and defect-rich quartz (101) surfaces with Si vacancies. FINDINGS: The results show that deprotonation is a highly dynamic process where both the surface hydroxyls and bridging oxygen atoms serve as the proton acceptors, in addition to water. Deprotonation of embedded silanols through intrasurface proton transfer exhibited lower pKa values with less H-bond participation and higher energy barriers, suggesting a new mechanism to explain the bimodal acidity observed on quartz surface. Defect sites, recently shown to comprise a significant portion of the quartz (101) surface, diversify the coordination and local H-bonding environments of the surface silanols, changing both the deprotonation pathways and energetics, leading to a wider range of pKa values (2.4 to 11.5) than that observed on pristine quartz surface (10.4 and 12.1).

RESUMEN

The emergence of cation-anion species, or contact ion pairs, is fundamental to understanding the physical properties of aqueous solutions when moving from the ideal, low-concentration limit to the manifestly non-ideal limits of very high solute concentration or constituent ion activity. We focus here on Zn halide solutions both as a model system and also as an exemplar of the applications spanning from (i) electrical energy storage via the paradigm of water in salt electrolyte (WiSE) to (ii) the physical chemistry of brines in geochemistry to (iii) the long-standing problem of nucleation. Using a combination of experimental and theoretical approaches we quantify the halide coordination number and changing coordination geometry without embedded use of theoretical equilibrium constants. These results and the associated methods, notably including the use of valence-to-core X-ray emission spectroscopy, provide new insights into the Zn halide system and new research directions in the physical chemistry of concentrated electrolytes.

RESUMEN

Nucleation and growth of minerals has broad implications in the geological, environmental and materials sciences. Recent developments in fast X-ray nanotomography have enabled imaging of crystal growth in solutions in situ with a resolution of tens of nanometres, far surpassing optical microscopy. Here, a low-cost, custom-designed aqueous flow cell dedicated to the study of heterogeneous nucleation and growth of minerals in aqueous environments is shown. To gauge the effects of radiation damage from the imaging process on growth reactions, radiation-induced morphological changes of barite crystals (hundreds of nanometres to ∼1 µm) that were pre-deposited on the wall of the flow cell were investigated. Under flowing solution, minor to major crystal dissolution was observed when the tomography scan frequency was increased from every 30 min to every 5 min (with a 1 min scan duration). The production of reactive radicals from X-ray induced water radiolysis and decrease of pH close to the surface of barite are likely responsible for the observed dissolution. The flow cell shown here can possibly be adopted to study a wide range of other chemical reactions in solutions beyond crystal nucleation and growth where the combination of fast flow and fast scan can be used to mitigate the radiation effects.

RESUMEN

The interactions of heavy metals with minerals influence the mobility and bioavailability of toxic elements in natural aqueous environments. The sorption of heavy metals on covalently bonded minerals is generally well described by surface complexation models (SCMs). However, understanding sorption on sparingly soluble minerals is challenging because of the dynamically evolving chemistry of sorbent surfaces. The interpretation can be even more complicated when multiple metal ions compete for sorption. In the present study, we observed synergistically enhanced uptake of lead and selenate on the barite (001) surface through two sorption mechanisms: lattice incorporation that dominates at lower coverages and two-dimensional monolayer growth that dominates at higher coverages. We also observed a systematic increase in the sorption affinity with increasing co-sorbed ion coverages, different from the assumption of invariant binding constants for individual adsorption processes in classical SCMs. Computational simulations showed thermodynamically favorable co-incorporation of lead and selenate by simultaneously substituting for barium and sulfate in neighboring sites, resulting in the formation of molecular clusters that locally match the net dimension of the substrate lattice. These results emphasize the importance of ion-ion interactions at mineral-water interfaces that control the fate and transport of contaminants in the environment.

Asunto(s)

RESUMEN

We characterize a concentrated 7.3 m CaCl2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca2+-Cl- interactions were primarily "hidden" in the Ca-NDIS signal due to overlapping with Ca2+-Ow and Ca2+-Hw contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models "on average," the question "which model is the best" has not been answered but replaced by the question "which model is better for a given research." An overall "good" model can still be inappropriate or, in some instances, "bad" and, unfortunately, produce erroneous results. The accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

Asunto(s)