Supersaturated calcium carbonate solutions are classical

Prenucleation Cluster Pathway is Inconsistent with CaCO3 Kinetics

It has been debated whether CaCO3 nucleates classically with the attainment of a critical cluster size or nonclassically with the restructuring of a prenucleation cluster (PNC). Here, we determine from the nucleation kinetics of CaCO3 that the transition state is composed of about 10 formula units, irrespective of the supersaturation. Crucially, the size of the transition state is considerably smaller than the average PNC size estimated from experimental characterization. This size discrepancy suggests the PNCs are uninvolved in nucleation, and the kinetics indicate that if CaCO3 nucleates classically, the transition state must be an abnormally unstable (antimagic) cluster.

Organic Nucleation: Water Rearrangement Reveals the Pathway of Ibuprofen

The organic nucleation of the pharmaceutical ibuprofen is investigated, as triggered by the protonation of ibuprofen sodium salt at elevated pH. The growth and aggregation of nanoscale solution species by Analytical Ultracentrifugation and Molecular Dynamics (MD) simulations is tracked. Both approaches reveal solvated molecules, oligomers, and prenucleation clusters, their size as well as their hydration at different reaction stages. By combining surface-specific vibrational spectroscopy and MD simulations, water interacting with ibuprofen at the air-water interface during nucleation is probed. The results show the structure of water changes upon ibuprofen protonation in response to the charge neutralization. Remarkably, the water structure continues to evolve despite the saturation of protonated ibuprofen at the hydrophobic interface. This further water rearrangement is associated with the formation of larger aggregates of ibuprofen molecules at a late prenucleation stage. The nucleation of ibuprofen involves ibuprofen protonation and their hydrophobic assembly. The results highlight that these processes are accompanied by substantial water reorganization. The critical role of water is possibly relevant for organic nucleation in aqueous environments in general.

Interaction and Stability of Nanobubbles and Prenucleation Calcium Clusters during Ultrasonic Treatment of Hard Water

To investigate the stability of nanobubbles in natural hard water, a series of eight samples ranging in hardness from 0 to 332 mg/L CaCO3 were sonicated for periods of 5-45 min with an ultrasonic horn. Conductivity, temperature, ζ-potential, composition, and pH of the water were analyzed, together with the crystal structure of any calcium carbonate precipitate. Quasi-stable populations of bulk nanobubbles in Millipore and soft water are characterized by a ζ-potential of -35 to -20 mV, decaying over 60 h or more. After sonicating the hardest waters for about 10 min, they turn cloudy due to precipitation of amorphous calcium carbonate when the water temperature reaches 40 °C; the ζ-potential then jumps from -10 to +20 mV and remains positive for several days. From an analysis of the change of conductivity of the hard water before and after sonication, it is estimated that 37 ± 5% of calcium was not originally in solution but existed in nanoscale prenucleation clusters, which decorate the nanobubbles formed in the early stages of sonication. Heating and charge screening in the nanobubble colloid cause the decorated bubbles to collapse or disperse, leaving an amorphous precursor of aragonite. Sonicating the soft supernatant increases its conductivity and pH and restores the negative ζ-potential associated with bulk nanobubbles, but there is no further precipitation. Our study of the correlation between nanobubble production and calcium agglomeration spanning the hardness and composition ranges of natural waters shows that the sonication method for introducing nanobubbles is viable only for hard water if it is kept cold; the stability of the nanobubble colloid will be reduced in any case by the presence of dissolved calcium and magnesium.

Efficient Assessment of 'Instantaneous pK' Values from Molecular Dynamics Simulations

We present a molecular simulation approach to studying the role of local and momentary molecular environment for potential acid-base reactions. For this, we combine thermodynamic considerations on the pK of ionic species with rapid sampling of energy changes related to (de)protonation. Using dispersed carbonate ions in water as a reference, our approach aims at the fast assessment of the momentary protonation energy, and thus the 'instantaneous pK', of calcium-carbonate ion aggregates. The latter include transient complexes that are elusive to long sampling runs. This motivated the elaboration of approximate, yet particularly fast assessable sampling strategies. Along this line, we were able to characterize instantaneous pK values at a statistical accuracy of 0.4 pK units within sampling runs of only 10 ps duration, whereas statistical errors reduce to 0.1 pK units in 75 ps sampling runs, respectively. This readily enabled the required time resolution for the characterization of [Cax (CO3 )y ]2(x-y) aggregates with x=1,2 and y=1,2,3, respectively. In turn, the analysis of the pH-dependent nature of calcite-water interfaces and dynamically ordered liquid-like oxyanion polymers (dollop) domains is outlined at 10 ps resolution.

Characterizing the Mechanisms of Ca and Mg Carbonate Ion-Pair Formation with Multi-Level Molecular Dynamics/Quantum Mechanics Simulations

The carbonate minerals of Ca and Mg are abundant throughout the lithosphere and have recently garnered significant research interest as possible long-term carbon sinks in the sequestration of atmospheric carbon dioxide. Nonetheless, an understanding of the atomic-level processes comprising their mineralization remains limited. Here, we characterize and contrast the mechanisms of contact ion-pair formation in aqueous Ca and Mg carbonate systems, which represents the most fundamental step leading to the formation of their mineral solids. Utilizing multilevel embedded correlated wavefunction-based ab initio molecular dynamics/quantum mechanics simulations, we characterize not only the dynamics of these processes but also factors arising from the electronic structure of the involved species, revealing further details of the fundamentally different mechanisms for the interconversion between the contact ion-pairs and solvent-shared ion-pairs of Ca versus Mg carbonate.

Unlocking the mysterious polytypic features within vaterite CaCO3

Calcium carbonate (CaCO3), the most abundant biogenic mineral on earth, plays a crucial role in various fields such as hydrosphere, biosphere, and climate regulation. Of the four polymorphs, calcite, aragonite, vaterite, and amorphous CaCO3, vaterite is the most enigmatic one due to an ongoing debate regarding its structure that has persisted for nearly a century. In this work, based on systematic transmission electron microscopy characterizations, crystallographic analysis and machine learning aided molecular dynamics simulations with ab initio accuracy, we reveal that vaterite can be regarded as a polytypic structure. The basic phase has a monoclinic lattice possessing pseudohexagonal symmetry. Direct imaging and atomic-scale simulations provide evidence that a single grain of vaterite can contain three orientation variants. Additionally, we find that vaterite undergoes a second-order phase transition with a critical point of ~190 K. These atomic scale insights provide a comprehensive understanding of the structure of vaterite and offer advanced perspectives on the biomineralization process of calcium carbonate.

Microfluidic investigation of pore-size dependency of barite nucleation

The understanding and prediction of mineral precipitation processes in porous media are relevant for various energy-related subsurface applications. While it is well known that thermodynamic effects can inhibit crystallization in pores with sizes <0.1 µm, the retarded observation of mineral precipitation as function of pore size is less explored. Using barite as an example and based on a series of microfluidic experiments with well-defined pore sizes and shapes, we show that retardation of observation of barite crystallite can already start in pores of 1 µm size, with the probability of nucleation scaling with the pore volume. In general, it can be expected that mineralization occurs preferentially in larger pores in rock matrices, but other parameters such as the exchange of the fluids with respect to reaction time, as well as shape, roughness, and surface functional properties of the pores may affect the crystallization process which can reverse this trend.

Scaling Solute-Solvent Distances to Improve Solubility and Ion Paring Predictions in Rigid Ion Models

Force fields based on the rigid ion model (RIM) have been developed to accurately predict the various physical and chemical properties of salts and water. However, the combined use of these models often fails to accurately predict the solubility of salts in water. To address this issue, several approaches, such as charge scaling or reparameterization, have been proposed. Nevertheless, these methods require laborious reparameterization of nonbonded force field parameters. In this article, we propose a scaling solute-solvent distance (SSSD) method to improve force fields in predicting salt solubility without changing the solute-solute and solvent-solvent interactions in the original force fields. This method can also tune the ion pairing of salt in water. One main advantage of the SSSD method is that reparameterization of the crystal and water models is not needed. We use two RIMs for the NaCl-water system (JC-SPC/E and SD-SPC/E) and the CHARMM force field for the KCl-water system to demonstrate the improved accuracy in predicting solubility by the SSSD method. Furthermore, we use the RDG-SPC/Fw force field to show that the SSSD method can also be used to tune the ion pairing of CaCO3 in water. Limitations of this method are also discussed.

Influence of the Precursor Structure on the Formation of Tungsten Oxide Polymorphs

Understanding material nucleation processes is crucial for the development of synthesis pathways for tailormade materials. However, we currently have little knowledge of the influence of the precursor solution structure on the formation pathway of materials. We here use in situ total scattering to show how the precursor solution structure influences which crystal structure is formed during the hydrothermal synthesis of tungsten oxides. We investigate the synthesis of tungsten oxide from the two polyoxometalate salts, ammonium metatungstate, and ammonium paratungstate. In both cases, a hexagonal ammonium tungsten bronze (NH4)0.25WO3 is formed as the final product. If the precursor solution contains metatungstate clusters, this phase forms directly in the hydrothermal synthesis. However, if the paratungstate B cluster is present at the time of crystallization, a metastable intermediate phase in the form of a pyrochlore-type tungsten oxide, WO3·0.5H2O, initially forms. The pyrochlore structure then undergoes a phase transformation into the tungsten bronze phase. Our studies thus experimentally show that the precursor cluster structure present at the moment of crystallization directly influences the formed crystalline phase and suggests that the precursor structure just prior to crystallization can be used as a tool for targeting specific crystalline phases of interest.

The evolution of solvation symmetry and composition in Zn halide aqueous solutions from dilute to extreme concentrations

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.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

A Tale of Two Foulants: The Coupling of Organic Fouling and Mineral Scaling in Membrane Desalination

Membrane desalination that enables the harvesting of purified water from unconventional sources such as seawater, brackish groundwater, and wastewater has become indispensable to ensure sustainable freshwater supply in the context of a changing climate. However, the efficiency of membrane desalination is greatly constrained by organic fouling and mineral scaling. Although extensive studies have focused on understanding membrane fouling or scaling separately, organic foulants commonly coexist with inorganic scalants in the feedwaters of membrane desalination. Compared to individual fouling or scaling, combined fouling and scaling often exhibits different behaviors and is governed by foulant-scalant interactions, resembling more complex but practical scenarios than using feedwaters containing only organic foulants or inorganic scalants. In this critical review, we first summarize the performance of membrane desalination under combined fouling and scaling, involving mineral scales formed via both crystallization and polymerization. We then provide the state-of-the-art knowledge and characterization techniques pertaining to the molecular interactions between organic foulants and inorganic scalants, which alter the kinetics and thermodynamics of mineral nucleation as well as the deposition of mineral scales onto membrane surfaces. We further review the current efforts of mitigating combined fouling and scaling via membrane materials development and pretreatment. Finally, we provide prospects for future research needs that guide the design of more effective control strategies for combined fouling and scaling to improve the efficiency and resilience of membrane desalination for the treatment of feedwaters with complex compositions.

A Review of Laser-Induced Crystallization from Solution

Crystallization abounds in nature and industrial practice. A plethora of indispensable products ranging from agrochemicals and pharmaceuticals to battery materials are produced in crystalline form in industrial practice. Yet, our control over the crystallization process across scales, from molecular to macroscopic, is far from complete. This bottleneck not only hinders our ability to engineer the properties of crystalline products essential for maintaining our quality of life but also hampers progress toward a sustainable circular economy in resource recovery. In recent years, approaches leveraging light fields have emerged as promising alternatives to manipulate crystallization. In this review article, we classify laser-induced crystallization approaches where light-material interactions are utilized to influence crystallization phenomena according to proposed underlying mechanisms and experimental setups. We discuss nonphotochemical laser-induced nucleation, high-intensity laser-induced nucleation, laser trapping-induced crystallization, and indirect methods in detail. Throughout the review, we highlight connections among these separately evolving subfields to encourage the interdisciplinary exchange of ideas.

Adsorption of CO32-/HCO3- on a quartz surface: cluster formation, pH effects, and mechanistic aspects

Soluble inorganic carbon is an important component of a soil carbon pool, and its fate in soils, sediments, and underground water environments has great effects on many physiochemical and geological processes. However, the dynamical processes, behaviors and mechanism of their adsorption by soil active components, such as quartz, are still unclear. The aim of this work is to systematically address the anchoring mechanism of CO₃^2- and HCO₃^- on a quartz surface at different pH values. Three pH values (pH 7.5, pH 9.5 and pH 11) and three carbonate salt concentrations (0.07, 0.14 and 0.28 M) are considered, and molecular dynamics methods are used. The results indicate that the pH value regulates the adsorption behavior of CO₃^2- and HCO₃^- on the quartz surface by affecting the CO₃^2-/HCO₃^- ratio and the surface charge of quartz. In general, both HCO₃^- and CO₃^2- ions were able to adsorb on the quartz surface and the adsorption capacity of CO₃^2- is higher than that of HCO₃^-. HCO₃^- ions tended to uniformly distribute in an aqueous solution and contact the quartz surface in the form of single molecules instead of clusters. In contrast, CO₃^2- ions were mainly adsorbed as clusters which became larger as the concentration increased. Na⁺ ions were essential for the adsorption of HCO₃^- and CO₃^2-, because some of the Na⁺ and CO₃^2- ions spontaneously associated together to form clusters, promoting the clusters to be adsorbed on the quartz surface through cationic bridges. The local structures and dynamics trajectory of CO₃^2- and HCO₃^- showed that the anchoring mechanism of carbonate solvates on quartz involved H-bonds and cationic bridges, which changed in relation to the concentration and pH values. However, the HCO₃^- ions mainly adsorbed on the quartz surface via H-bonds while the CO₃^2- ions tended to be adsorbed through cationic bridges. These results may help in understanding the geochemical behavior of soil inorganic carbon and further the processes of the Earth's carbon chemical cycle.

Chitosan as a Canvas for Studies of Macromolecular Controls on CaCO3 Biological Crystallization

A mechanistic understanding of how macromolecules, typically as an organic matrix, nucleate and grow crystals to produce functional biomineral structures remains elusive. Advances in structural biology indicate that polysaccharides (e.g., chitin) and negatively charged proteoglycans (due to carboxyl, sulfate, and phosphate groups) are ubiquitous in biocrystallization settings and play greater roles than currently recognized. This review highlights studies of CaCO₃ crystallization onto chitinous materials and demonstrates that a broader understanding of macromolecular controls on mineralization has not emerged. With recent advances in biopolymer chemistry, it is now possible to prepare chitosan-based hydrogels with tailored functional group compositions. By deploying these characterized compounds in hypothesis-based studies of nucleation rate, quantitative relationships between energy barrier to crystallization, macromolecule composition, and solvent structuring can be determined. This foundational knowledge will help researchers understand composition-structure-function controls on mineralization in living systems and tune the designs of new materials for advanced applications.

Local Molecular Field Theory for Coulomb Interactions in Aqueous Solutions

Coulomb interactions play a crucial role in a wide array of processes in aqueous solutions but present conceptual and computational challenges to both theory and simulations. We review recent developments in an approach addressing these challenges─local molecular field (LMF) theory. LMF theory exploits an exact and physically suggestive separation of intermolecular Coulomb interactions into strong short-range and uniformly slowly varying long-range components. This allows us to accurately determine the averaged effects of the long-range components on the short-range structure using effective single particle fields and analytical corrections, greatly reducing the need for complex lattice summation techniques used in most standard approaches. The simplest use of these ideas in aqueous solutions leads to the short solvent (SS) model, where both solvent-solvent and solute-solvent Coulomb interactions have only short-range components. Here we use the SS model to give a simple description of pairing of nucleobases and biologically relevant ions in water.

Molecular cluster analysis using local order parameters selected by machine learning

Accurately extracting local molecular structures is essential for understanding the mechanisms of phase and structural transitions. A promising method to characterize the local molecular structure is defining the value of the local order parameter (LOP) for each particle. This work develops the Molecular Assembly structure Learning package for Identifying Order parameters (MALIO), a machine learning package that can propose an optimal (set of) LOP(s) quickly and automatically for a huge number of LOP species and various methods of selecting neighboring particles for the calculation. We applied this package to distinguish between the nematic and smectic phases of uniaxial liquid crystal molecules, and selected candidate LOPs that could be used to precisely observe the nematic-smectic phase transition. The LOP candidates were used to observe the nucleation and subsequent percolation transition, and the effect of the choice of LOP species and neighboring particles on the statistics of local molecular structures (clusters) was examined. The procedure revealed the time evolution of the number of clusters and the dependence of the percolation curve on the number of neighboring particles for each LOP species. The LOP species with the lowest dependence on the number of neighboring particles was the best-performing LOP species in the MALIO screening strategy. These results not only show that machine learning can powerfully screen a huge number of LOP species and suggest only a few promising candidates, but also indicate that MALIO can select the best LOP species.

New insights into the early stage nucleation of calcium carbonate gels by reactive molecular dynamics simulations

The precipitation of calcium carbonate (CaCO₃) is a key mechanism in carbon capture applications relying on mineralization. In that regard, Ca-rich cementitious binders offer a unique opportunity to act as a large-scale carbon sink by immobilizing CO₂ as calcium carbonate by mineralization. However, the atomistic mechanism of calcium carbonate formation is still not fully understood. Here, we study the atomic scale nucleation mechanism of an early stage amorphous CaCO₃ gel based on reactive molecular dynamics (MD) simulations. We observe that reactive MD offers a notably improved description of this reaction as compared to classical MD, which allows us to reveal new insights into the structure of amorphous calcium carbonate gels and formation kinetics thereof.

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na2SO4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Herein, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

Small-Molecular-Weight Additives Modulate Calcification by Interacting with Prenucleation Clusters on the Molecular Level

Small-molecular-weight (MW) additives can strongly impact amorphous calcium carbonate (ACC), playing an elusive role in biogenic, geologic, and industrial calcification. Here, we present molecular mechanisms by which these additives regulate stability and composition of both CaCO₃ solutions and solid ACC. Potent antiscalants inhibit ACC precipitation by interacting with prenucleation clusters (PNCs); they specifically trigger and integrate into PNCs or feed PNC growth actively. Only PNC-interacting additives are traceable in ACC, considerably stabilizing it against crystallization. The selective incorporation of potent additives in PNCs is a reliable chemical label that provides conclusive chemical evidence that ACC is a molecular PNC-derived precipitate. Our results reveal additive-cluster interactions beyond established mechanistic conceptions. They reassess the role of small-MW molecules in crystallization and biomineralization while breaking grounds for new sustainable antiscalants.

Formation mechanisms and supervisory prediction of scaling in water supply pipelines: A review

The scaling problem in the water supply pipeline will increase the resistance coefficient of the pipeline and the pressure of the water supply pipeline, which will not only affect the operation safety of the water supply pipeline, but also cause energy waste. The scale in the pipeline will also enrich heavy metal ions and pathogenic microorganisms, affecting the safety of water supply water quality and causing secondary pollution of water quality. At present, a lot of research has been done on the composition structure and crystallization process of the scale. The study found that calcite is the main component of the scale; the scale process is a heterogeneous nucleation process induced by heavy metal particles and their corrosion products in the pipeline. The introduction of electrochemical detection technology, density functional theory and molecular dynamics simulation has greatly improved the accuracy and timeliness of water scaling conditions detection and realized the visualization of scaling mechanism. In this paper, the measurement methods of the scale in the water supply pipeline and the corresponding material composition and crystal structure characteristics are reviewed, and the mechanism of the scale and the water quality conditions are summarized. At the end of this paper, based on summarizing the existing water quality scaling tendency evaluation methods, it is proposed to establish a water quality potential scaling risk assessment framework based on Puckorius scaling index (PSI) and electrochemical impedance spectroscopy (EIS) in the future.

Physical insight into the entropy-driven ion association

The ion association is widely believed to be dominated by the favorable entropy change arising from the release of water molecules from ion hydration shells. However, no direct thermodynamic evidence exists to validate the reliability and suitability of this view. Herein, we employ complicated free energy calculations to rigorously split the free energy including its entropic and enthalpic components into the water-induced contributions and ion-ion interaction terms for several ion pairs from monatomic to polyatomic ions, spanning the size range from small kosmotropes to large chaotropes (Na⁺ , Cs⁺ , Ca²⁺ , F^- , I^- , CO₃ ^2- , and HPO₄ ^2- ). Our results successfully reveal that though ion associations are indeed determined by a delicate balance between the favorable entropy variation and the repulsive enthalpy change, the entropy gain dominated by the solvent occurs only for the monatomic ion pairing. The water-induced entropic contribution significantly goes against the ion pairing between polyatomic anion and cation, which is, alternatively, dominated by the favorable entropy from the ion-ion interaction term, due to the configurational arrangement of polyatomic anions involved in ion association. The structural and dynamic analysis demonstrates that the entropy penalty from the water phase is primarily ascribed to the enhanced stability of water molecules around the cation imposed by the incoming anion. Our study successfully provides a fundamental understanding of water-mediated ion associations and highlights disparate lengthscale dependencies of the dehydration thermodynamics on the specific types of ions.

Crystal Nucleation and Growth of Inorganic Ionic Materials from Aqueous Solution: Selected Recent Developments, and Implications

In this review article, selected, latest theoretical, and experimental developments in the field of nucleation and crystal growth of inorganic materials from aqueous solution are highlighted, with a focus on literature after 2015 and on non-classical pathways. A key point is to emphasize the so far underappreciated role of water and solvent entropy in crystallization at all stages from solution speciation through to the final crystal. While drawing on examples from current inorganic materials where non-classical behavior has been proposed, the potential of these approaches to be adapted to a wide-range of systems is also discussed, while considering the broader implications of the current re-assessment of pathways for crystallization. Various techniques that are suitable for the exploration of crystallization pathways in aqueous solution, from nucleation to crystal growth are summarized, and a flow chart for the assignment of specific theories based on experimental observations is proposed.

Light Controlled 3D Crystal Morphology for Droplet Evaporative Crystallization on Photosensitive Hydrophobic Substrate

Controlling crystal morphology is crucial in analytical chemistry and smart materials synthesis, etc. However, flexible manipulation of 3D crystal morphology still remains challenging. Herein, we present a novel and facile light strategy for droplet evaporative crystallization to manipulate macroscopic crystal morphology on photosensitive hydrophobic substrate possessing photothermal conversion property. We demonstrate that the spherical coronal shell and alms bowl-like crystal skeletons can be achieved on smooth photosensitive hydrophobic substrate, depending on the salt concentration. Rough photosensitive hydrophobic substrate further creates a bubble-assisted light strategy, by which a cylindrical shell-like crystal skeleton with a directionally controllable cavity is achieved. Amazingly, the proper additive endows droplet evaporative crystallization to form a closed crystal skeleton with the solution wrapped inside. The present study provides new ideas for designing a novel optical droplet microfluidic platform for controlling crystal morphology.

Computational approaches for XANES, VtC-XES, and RIXS using linear-response time-dependent density functional theory based methods

The emergence of state-of-the-art X-ray light sources has paved the way for novel spectroscopies that take advantage of their atomic specificity to shed light on fundamental physical, chemical, and biological processes both in the static and time domains. The success of these experiments hinges on the ability to interpret and predict core-level spectra, which has opened avenues for theory to play a key role. Over the last two decades, linear-response time-dependent density functional theory (LR-TDDFT), despite various theoretical challenges, has become a computationally attractive and versatile framework to study excited-state spectra including X-ray spectroscopies. In this context, we focus our discussion on LR-TDDFT approaches for the computation of X-ray Near-Edge Structure (XANES), Valence-to-Core X-ray Emission (VtC-XES), and Resonant Inelastic X-ray Scattering (RIXS) spectroscopies in molecular systems with an emphasis on Gaussian basis set implementations. We illustrate these approaches with applications and provide a brief outlook of possible new directions.

Process-Specific Effects of Sulfate on CaCO3 Formation in Environmentally Relevant Systems

Additives, such as ions, small molecules, and macromolecules, have been found to regulate the formation of CaCO₃ and control its morphologies and properties. However, a single additive usually affects dominantly one process in CaCO₃'s formation and is seldom found to significantly affect multiple CaCO₃ formation processes. Here, we used in situ grazing incidence X-ray techniques to observe the heterogeneous formation of CaCO₃ and found that a series of formation processes (i.e., nucleation, growth, and Ostwald ripening) were modulated by sulfate. In the nucleation process, increased interfacial free energy and bulk free energy cooperatively increased the nucleation barrier and decreased nucleation rates. In the growth process, sulfate reduced the electrostatic repulsion between CaCO₃ precursors and nuclei, promoting CaCO₃ growth. This influence on the growth counteracted the inhibition effect in the nucleation process, causing a nearly 100% increase in the volume of heterogeneously formed CaCO₃. Meanwhile, adsorbed sulfate on CaCO₃ nuclei may poison the surface of smaller CaCO₃ nuclei, inhibiting Ostwald ripening. These revealed sulfate's active roles in controlling CaCO₃ formation advance our understanding of sulfate-incorporated biomineralization and scaling phenomena in natural and engineered aquatic environments.

Probing the effects of polymers on the early stages of calcium carbonate formation by stoichiometric co-titration

Potentiometric titrations are a powerful tool to study the early stages of the precipitation of minerals such as calcium carbonate and were used among others for the discovery and characterisation of key precursors like prenucleation clusters. Here we present a modified procedure for conducting such titration experiments, in which the reactants (i.e. calcium and (bi)carbonate ions) are added simultaneously in stoichiometric amounts, while both the amount of free calcium and the optical transmission of the solution are monitored online. Complementarily, the species occurring at distinct stages of the crystallisation process were studied using cryogenic transmission electron microscopy. This novel routine was applied to investigate CaCO₃ nucleation in the absence and presence of polymeric additives with different chemical functionalities. The obtained results provide new insights into the critical steps underlying nucleation and subsequent ripening, such as the role of liquid mineral-rich phases and their transformation into solid particles. The studied polymers proved to interfere at multiple stages along the complex mineralisation pathway of calcium carbonate, with both the degree and mode of interaction depending on the chosen polymer chemistry. In this way, the methodology developed in this work allows the mechanisms of antiscalants - or crystallisation modifiers in general - to be elucidated at an advanced level of detail.

X-ray compatible microfluidics for in situ studies of chemical state, transport and reaction of light elements in an aqueous environment using synchrotron radiation

This paper presents an X-ray compatible microfluidic platform for in situ characterization of chemical reactions at synchrotron light sources. We demonstrate easy to implement techniques to probe reacting solutions as they first come into contact, and study the very first milliseconds of their reaction in real-time through X-ray absorption spectroscopy (XAS). The devices use polydimethylsiloxane (PDMS) microfluidic channels sandwiched between ultrathin, X-ray transparent silicon nitride observation windows and rigid substrates. The new approach has three key advantages: i) owing to the assembly techniques employed, the devices are suitable for both high energy and tender (1-5 keV) X-rays; ii) they can operate in a vacuum environment (a must for low energy X-rays) and iii) they are robust enough to survive a full 8 hour shift of continuous scanning with a micro-focused beam, providing higher spatial and thus greater time resolution than previous studies. The combination of these opens new opportunities for in situ studies. This has so far not been possible with Kapton or glass-based flow cells due to increased attenuation of the low energy beam passing through these materials. The devices provide a well-defined mixing region to collect spatial maps of spatially stable concentration profiles, and XAS point spectra to elucidate the chemical structure and characterize the chemical reactions. The versatility of the approach is demonstrated through in situ XAS measurements on the mixing of two reactants in a microfluidic laminar flow device, as well as a segmented droplet based system for time resolved analysis.

Complete inhibition of a polyol nucleation by a micromolar biopolymer additive

Preventing spontaneous crystallization of supersaturated solutions by additives is of critical interest to successful process design and implementation, with numerous applications in chemical, pharmaceutical, medical, pigment, and food industries, but challenges remain in laboratory and industry settings and fundamental understanding is lacking. When copresented with antifreeze proteins (AFPs), otherwise spontaneously crystallizing osmolytes are maintained at high supersaturations for months in over-wintering organisms. Thus, we here explore the inhibition phenomenon by AFPs, using persistent crystallization of a common sugar alcohol, D-mannitol, as a case study. We report experimentally that DAFP1, an insect AFP, completely inhibits D-mannitol nucleation. Computer simulations reveal a new mechanism for crystallization inhibition where the population of the crystal-forming conformers are selectively bound and randomized in solution by hydrogen bonding to the protein surface. These results highlight the advantages of using natural polymers to address crystallization inhibition challenges and suggest new strategies in controlling the nucleation processes.

Nested Formation of Calcium Carbonate Polymorphs in a Bacterial Surface Membrane with a Graded Nanoconfinement: An Evolutionary Strategy to Ensure Bacterial Survival

It is the intention of this study to elucidate the nested formation of calcium carbonate polymorphs or polyamorphs in the different nanosized compartments. With these observations, it can be concluded how the bacteria can survive in a harsh environment with high calcium carbonate supersaturation. The mechanisms of calcium carbonate precipitation at the surface membrane and at the underlying cell wall membrane of the thermophilic soil bacterium Geobacillus stearothermophilus DSM 13240 have been revealed by high-resolution transmission electron microscopy and atomic force microscopy. In this Gram-positive bacterium, nanopores in the surface layer (S-layer) and in the supporting cell wall polymers are nucleation sites for metastable calcium carbonate polymorphs and polyamorphs. In order to observe the different metastable forms, various reaction times and a low reaction temperature (4 °C) have been chosen. Calcium carbonate polymorphs nucleate in the confinement of nanosized pores (⌀ 3–5 nm) of the S-layer. The hydrous crystalline calcium carbonate (ikaite) is formed initially with [110] as the favored growth direction. It transforms into the anhydrous metastable vaterite by a solid-state transition. In a following reaction step, calcite is precipitated, caused by dissolution of vaterite in the aqueous solution. In the larger pores of the cell wall (⌀ 20–50 nm), hydrated amorphous calcium carbonate is grown, which transforms into metastable monohydrocalcite, aragonite, or calcite. Due to the sequence of reaction steps via various metastable phases, the bacteria gain time for chipping the partially mineralized S-layer, and forming a fresh S-layer (characteristic growth time about 20 min). Thus, the bacteria can survive in solutions with high calcium carbonate supersaturation under the conditions of forced biomineralization.

Liquid-Liquid Phase Separation in Nucleation Process of Biomineralization

Biomineralization is a typical interdisciplinary subject attracting biologists, chemists, and geologists to figure out its potential mechanism. A mounting number of studies have revealed that the classical nucleation theory is not suitable for all nucleation process of biominerals, and phase-separated structures such as polymer-induced liquid precursors (PILPs) play essential roles in the non-classical nucleation processes. These structures are able to play diverse roles biologically or pathologically, and could also give inspiring clues to bionic applications. However, a lot of confusion and dispute occurred due to the intricacy and interdisciplinary nature of liquid precursors. Researchers in different fields may have different opinions because the terminology and current state of understanding is not common knowledge. As a result, our team reviewed the most recent articles focusing on the nucleation processes of various biominerals to clarify the state-of-the-art understanding of some essential concepts and guide the newcomers to enter this intricate but charming field.

Pressure-induced crystallization and densification of amorphized calcium carbonate hexahydrate controlled by interfacial water

Amorphous calcium carbonate (ACC) is widely known as a metastable precursor in the formation of crystalline calcium carbonate biominerals. However, the exact role of water during the crystallization of ACC remains elusive. Here, a novel ACC with high specific surface area and nanopores is synthesized by solvent-induced dehydration and amorphization of crystalline calcium carbonate hexahydrate (ikaite), denoted as I-ACC. Comparing I-ACC and typical spherical ACC (S-ACC) nanoparticles, it reveals that the crystallization pathways of ACC under heating or pressure are not dictated by the total amount of water in ACC as reported, but rather the interfacial water that is released from ACC bulk and adsorbed on the surface of the particles. We show that the crystallization pathways of I-ACC to calcite single crystal with high specific surface area or vaterite can be easily controlled by tuning the release of water during heating. In addition, densely packed pure vaterite can be obtained via pressured-induced transformation of I-ACC at room temperature, which is otherwise difficult to form using S-ACC. These insights contribute to the understanding of the biological control of mineral formation via amorphous precursors and offer new opportunities to bioprocess inspired fabrication of strong bulk material at room temperature.

Structure and phase analysis of calcium carbonate powder prepared by a simple solution method

This paper focused on the analysis of the crystal structure and phase transformation of CaCO₃ synthesized by simple solution method from 0.5 M Ca(NO₃)₂ precursor and 0.5 M Na₂CO₃ precipitant at ambient temperature (300 K). The pH of the sample solution at various reaction times of 5, 10, 15, and 30 min were measured and correlated with the supersaturating condition in the presence of the Na₂CO₃ which is responsible for vaterite phase formation. The formation of the polymorph structure of obtained CaCO₃ powders was characterized using powder X-ray diffraction patterns and their crystal structure and phase transformation were evaluated using the Rietveld refinement method. Moreover, the qualitative analysis of the CaCO₃ powders phase was conducted by Fourier Transform Infrared (FTIR) spectroscopy to evaluate the effect of reaction time correlated with their crystal formation. The XRD analysis showed that the vaterite formation was 89 % at a reaction time of 15 min and confirmed also by FTIR that the amount of vaterite increased due to the effect of increasing reaction time. The crystallite size of vaterite was stable at 36 nm at the reaction time of 15 and 30 min. The morphology of the CaCO₃ powders obtained from Scanning Electron Microscope (SEM) was spherical with sizes of 2–5 μm. It was highlighted that the supersaturating condition started occurred at a reaction time of 15 min at pH 7.88 which was responsible for vaterite formation took place. It was concluded that the amount of precipitant (Na₂CO₃) and reaction times play an important role to determine the saturation of carbonate source to allow vaterite phase formation of CaCO₃ powders to occur.

•

The engineered crystal structure and phase transformation of CaCO₃ were achieved at room temperature using Ca(NO₃)₂ and Na₂O₃ precursor.

•

The vaterite phase formed at 15 min and 30 min reaction times with pH ~ 7.9 shows the spherical morphologies (2–5 µm in diameter).

•

Different reaction times and pH vs. formation and transformation of various CaCO₃ phase was demonstrated experimentally.

Barite precipitation in porous media: Impact of pore structure and surface charge on ionic diffusion

Several scientific fields such as global carbon sequestration, deep geological radioactive waste disposal, and oil recovery/fracking encounter safety assessment issues originating from pore-scale processes such as mineral precipitation and dissolution. These processes occur in situations where the pore solution contains chemical complexity (such as pH, ionic strength, redox chemistry, etc.…) and the porous matrix contains physical complexity (such as pore size distribution, surface charge, surface roughness, etc.…). Thus, to comprehend the participation of each physicochemical phenomenon on governing mineral precipitation, it is essential to investigate the precipitation behavior of a given mineral in different confined volumes. In this study, a counter-diffusion approach was used to investigate barite precipitation in two porous materials: micritic chalk and compacted kaolinite. The two materials present similar water and anionic tracer diffusivities and total accessible porosities but distinct pore size distributions with pore throats of c.a. 660 nm in chalk versus c.a. 35 nm in kaolinite. X-ray tomography results obtained on the two materials showed a distinct distribution of barite precipitates: a 500 μm-thick homogeneous layer in chalk versus spherical clusters spread in a thickness of 2 mm in kaolinite. Mass balance calculations showed that barite precipitation led to a porosity decrease in the chalk reacted zone from 45% to 12% and in the kaolinite reacted zone from 36% to 34.5%. In contrast, water tracer diffusion experiments showed that diffusivity decreased by a factor of 28 in chalk and by a factor of 1000 in kaolinite. Such a discrepancy was attributed to the difference in the pore size distribution that would lead to the distinct barite precipitation patterns, capable of altering in a very different manner the connectivity within the reacted zone of the two selected porous media. Such local alterations in connectivity linked to pore volume reduction would also magnify surface charge effects on ionic transport, as indicated by chloride diffusion experiments and electrophoric tests using zeta potential measurements. Indeed, ³⁶Cl^- was strongly more hindered than water, when diffused in reacted materials, with a diffusivity decrease by a factor of 450 in chalk and a total restriction of ³⁶Cl^- in kaolinite. These experiments clearly provide an insight of how local pore structure properties combined with mineral reactivity could help in predicting the evolution of pore scale clogging and its impact on water and ionic diffusive transport.

Multistep nucleation of anisotropic molecules

Phase transition of anisotropic materials is ubiquitously observed in physics, biology, materials science, and engineering. Nevertheless, how anisotropy of constituent molecules affects the phase transition dynamics is still poorly understood. Here we investigate numerically the phase transition of a simple model system composed of anisotropic molecules, and report on our discovery of multistep nucleation of nuclei with layered positional ordering (smectic ordering), from a fluid-like nematic phase with orientational order only (no positional order). A trinity of molecular dynamics simulation, machine learning, and molecular cluster analysis yielding free energy landscapes unambiguously demonstrates the dynamics of multistep nucleation process involving characteristic metastable clusters that precede supercritical smectic nuclei and cannot be accounted for by the classical nucleation theory. Our work suggests that molecules of simple shape can exhibit rich and complex nucleation processes, and our numerical approach will provide deeper understanding of phase transitions and resulting structures in anisotropic materials such as biological systems and functional materials.

In Situ X-ray Absorption Spectroscopy and Droplet-Based Microfluidics: An Analysis of Calcium Carbonate Precipitation

Droplet-based microfluidic systems are ideally suited for the investigation of nucleation and crystallization processes. To best leverage the features of such platforms (including exquisite time resolution and high-throughput operation), sensitive and in situ detection schemes are needed to extract real-time chemical information about all species of interest. In this regard, the extension of conventional (UV, visible, and infrared) optical detection schemes to the X-ray region of the electromagnetic spectrum is of high current interest, as techniques such as X-ray absorption spectroscopy (XAS) provide for the element-specific investigation of the local chemical environment. Accordingly, herein, we report for the first time the integration of millisecond droplet-based microfluidics with XAS. Such a platform allows for the sensitive acquisition of X-ray absorption data from picoliter-volume droplets moving at high linear velocities. Significantly, the high-temporal resolution of the droplet-based microfluidic platform enables unprecedented access to the early stages of the reaction. Using such an approach, we demonstrate in situ monitoring of calcium carbonate precipitation by extracting XAS spectra at the early time points of the reaction with a dead time as low as 10 ms. We obtain insights into the kinetics of the formation of amorphous calcium carbonate (ACC) as a first species during the crystallization process by monitoring the proportion of calcium ions converted into ACC. Within the confined and homogeneous environment of picoliter-volume droplets, the ACC content reaches 60% over the first 130 ms. More generally, the presented method offers new opportunities for the real-time monitoring of fast chemical and biological processes.

Uncovering the Role of Bicarbonate in Calcium Carbonate Formation at Near-Neutral pH

Mechanistic pathways relevant to mineralization are not well‐understood fundamentally, let alone in the context of their biological and geological environments. Through quantitative analysis of ion association at near‐neutral pH, we identify the involvement of HCO₃⁻ ions in CaCO₃ nucleation. Incorporation of HCO₃⁻ ions into the structure of amorphous intermediates is corroborated by solid‐state nuclear magnetic resonance spectroscopy, complemented by quantum mechanical calculations and molecular dynamics simulations. We identify the roles of HCO₃⁻ ions as being through (i) competition for ion association during the formation of ion pairs and ion clusters prior to nucleation and (ii) incorporation as a significant structural component of amorphous mineral particles. The roles of HCO₃⁻ ions as active soluble species and structural constituents in CaCO₃ formation are of fundamental importance and provide a basis for a better understanding of physiological and geological mineralization.

Enhancing the analysis of disorder in X-ray absorption spectra: application of deep neural networks to T-jump-X-ray probe experiments

Many chemical and biological reactions, including ligand exchange processes, require thermal energy for the reactants to overcome a transition barrier and reach the product state. Temperature-jump (T-jump) spectroscopy uses a near-infrared (NIR) pulse to rapidly heat a sample, offering an approach for triggering these processes and directly accessing thermally-activated pathways. However, thermal activation inherently increases the disorder of the system under study and, as a consequence, can make quantitative interpretations of structural changes challenging. In this Article, we optimise a deep neural network (DNN) for the instantaneous prediction of Co K-edge X-ray absorption near-edge structure (XANES) spectra. We apply our DNN to analyse T-jump pump/X-ray probe data pertaining to the ligand exchange processes and solvation dynamics of Co2+ in chlorinated aqueous solution. Our analysis is greatly facilitated by machine learning, as our DNN is able to predict quickly and cost-effectively the XANES spectra of thousands of geometric configurations sampled from ab initio molecular dynamics (MD) using nothing more than the local geometric environment around the X-ray absorption site. We identify directly the structural changes following the T-jump, which are dominated by sample heating and a commensurate increase in the Debye-Waller factor.

Controlling CaCO3 Particle Size with {Ca2+}:{CO3 2-} Ratios in Aqueous Environments

The effect of stoichiometry on the new formation and subsequent growth of CaCO₃ was investigated over a large range of solution stoichiometries (10^-4 < r _aq < 10⁴, where r _aq = {Ca²⁺}:{CO₃ ^2-}) at various, initially constant degrees of supersaturation (30 < Ω_cal < 200, where Ω_cal = {Ca²⁺}{CO₃ ^2-}/K _sp), pH of 10.5 ± 0.27, and ambient temperature and pressure. At r _aq = 1 and Ω_cal < 150, dynamic light scattering (DLS) showed that ion adsorption onto nuclei (1-10 nm) was the dominant mechanism. At higher supersaturation levels, no continuum of particle sizes is observed with time, suggesting aggregation of prenucleation clusters into larger particles as the dominant growth mechanism. At r _aq ≠ 1 (Ω_cal = 100), prenucleation particles remained smaller than 10 nm for up to 15 h. Cross-polarized light in optical light microscopy was used to measure the time needed for new particle formation and growth to at least 20 μm. This precipitation time depends strongly and asymmetrically on r _aq. Complementary molecular dynamics (MD) simulations confirm that r _aq affects CaCO₃ nanoparticle formation substantially. At r _aq = 1 and Ω_cal ≫ 1000, the largest nanoparticle in the system had a 21-68% larger gyration radius after 20 ns of simulation time than in nonstoichiometric systems. Our results imply that, besides Ω_cal, stoichiometry affects particle size, persistence, growth time, and ripening time toward micrometer-sized crystals. Our results may help us to improve the understanding, prediction, and formation of CaCO₃ in geological, industrial, and geo-engineering settings.

Role of Water in CaCO3 Biomineralization

Biomineralization occurs in aqueous environments. Despite the ubiquity and relevance of CaCO₃ biomineralization, the role of water in the biomineralization process has remained elusive. Here, we demonstrate that water reorganization accompanies CaCO₃ biomineralization for sea urchin spine generation in a model system. Using surface-specific vibrational spectroscopy, we probe the water at the interface of the spine-associated protein during CaCO₃ mineralization. Our results show that, while the protein structure remains unchanged, the structure of interfacial water is perturbed differently in the presence of both Ca²⁺ and CO₃^2– compared to the addition of only Ca²⁺. This difference is attributed to the condensation of prenucleation mineral species. Our findings are consistent with a nonclassical mineralization pathway for sea urchin spine generation and highlight the importance of protein hydration in biomineralization.

Variation in Properties of Pre-Nucleation Calcium Carbonate Clusters Induced by Aggregation: A Molecular Dynamics Study

Numerous studies have speculated calcium carbonate (CaCO3) nucleation induced by pre-nucleation clusters (PNCs) aggregation. However, it is challenging for experiments to directly obtain the relationship between PNCs aggregation and nucleation. Herein, we employ molecular dynamics simulations to explore the variation during PNCs aggregation, which can describe the beginning stage of CaCO3 nucleation induced by PNCs aggregation in supersaturated solutions. The results reveal that the formation of CaCO3 nucleus consists of PNCs spontaneous growth, PNCs solubility equilibrium, and aggregation of PNCs inducing nucleation. The PNCs aggregation, accompanied by the variation in the configuration and stability of CaCO3 aggregate, breaks the solubility equilibrium of PNCs and creates conditions for the formation of the more stable nucleus. Besides, the CaCO3 nucleus with the higher coordination number and the lower hydration number form when decreasing the CaCO3 concentration or increasing the temperature. This work not only sheds light on the formation of the CaCO3 nucleus but also contributes to the explanation for CaCO3 polymorphism.

Multiple Pathways in NaCl Homogeneous Crystal Nucleation

NaCl crystal nucleation from metastable solutions has long been considered to occur according to a single-step mechanism where the growth in the size and crystalline order of the emerging nuclei...

Nanoscale kinetics of amorphous calcium carbonate precipitation in H2O and D2O

Calcium carbonate (CaCO₃) is one of the most well-studied and abundant natural materials on Earth. Crystallisation of CaCO₃ is often observed to proceed via an amorphous calcium carbonate (ACC) phase, as a precursor to more stable crystalline polymorphs such as vaterite and calcite. Despite its importance, the kinetics of ACC formation have proved difficult to study, in part due to rapid precipitation at moderate supersaturations, and the instability of ACC with respect to all other polymorphs. However, ACC can be stabilised under confinement conditions, such as those provided by a nanopipette. This paper demonstrates electrochemical mixing of a Ca²⁺ salt (CaCl₂) and a HCO₃^- salt (NaHCO₃) in a nanopipette to repeatedly and reversibly precipitate nanoparticles of ACC under confined conditions, as confirmed by scanning transmission electron microscopy (STEM). Measuring the current as a function of applied potential across the end of the nanopipette and time provides millisecond-resolved measurements of the induction time for ACC precipitation. We demonstrate that under conditions of electrochemical mixing, ACC precipitation is extremely fast, and highly pH sensitive with an apparent third order dependence on CO₃^2- concentration. Furthermore, the rate is very similar for the equivalent CO₃^2- concentrations in D₂O, suggesting that neither ion dehydration nor HCO₃^- deprotonation represent significant energetic barriers to the formation of ACC. Finite element method simulations of the electrochemical mixing process enable the supersaturation to be estimated for all conditions and accurately predict the location of precipitation.

Probing Sulfur Chemical and Electronic Structure with Experimental Observation and Quantitative Theoretical Prediction of Kα and Valence-to-Core Kβ X-ray Emission Spectroscopy

An extensive experimental and theoretical study of the Kα and Kβ high-resolution X-ray emission spectroscopy (XES) of sulfur-bearing systems is presented. This study encompasses a wide range of organic and inorganic compounds, including numerous experimental spectra from both prior published work and new measurements. Employing a linear-response time-dependent density functional theory (LR-TDDFT) approach, strong quantitative agreement is found in the calculation of energy shifts of the core-to-core Kα as well as the full range of spectral features in the valence-to-core Kβ spectrum. The ability to accurately calculate the sulfur Kα energy shift supports the use of sulfur Kα XES as a bulk-sensitive tool for assessing sulfur speciation. The fine structure of the sulfur Kβ spectrum, in conjunction with the theoretical results, is shown to be sensitive to the local electronic structure including effects of symmetry, ligand type and number, and, in the case of organosulfur compounds, to the nature of the bonded organic moiety. This agreement between theory and experiment, augmented by the potential for high-access XES measurements with the latest generation of laboratory-based spectrometers, demonstrates the possibility of broad analytical use of XES for sulfur and nearby third-row elements. The effective solution of the forward problem, i.e., successful prediction of detailed spectra from known molecular structure, also suggests future use of supervised machine learning approaches to experimental inference, as has seen recent interest for interpretation of X-ray absorption near-edge structure (XANES).

Hierarchical phenomena in multicomponent liquids: simulation methods, analysis, chemistry

Complex, multicomponent, solutions have often been studied solely through the lens of specific applications of interest. Yet advances to both simulation methodologies (enhanced sampling, etc.) and analysis techniques (network analysis algorithms and others), are creating a trove of data that reveal transcending characteristics across vast compositional phase space. This perspective discusses technical considerations of the reliable and accurate simulations of complex solutions, followed by the advances to analysis algorithms that elucidate coupling of different length and timescale behavior (hierarchical phenomena). The different manifestations of hierarchical phenomena are presented across an array of solution environments, emphasizing fundamental and ongoing science questions. With a more advanced molecular understanding in hand, a quintessential application (solvent extraction) is discussed, where significant opportunities exist to re-imagine the technical scope of an established technology.

Amino Acid and Oligopeptide Effects on Calcium Carbonate Solutions

Biological organisms display sophisticated control of nucleation and crystallization of minerals. In order to mimic living systems, deciphering the mechanisms by which organic molecules control the formation of mineral phases from solution is a key step. We have used computer simulations to investigate the effects of the amino acids arginine, aspartic acid, and glycine on species that form in solutions of calcium carbonate (CaCO₃) at lower and higher levels of supersaturation. This provides net positive, negative, and neutral additives. In addition, we have prepared simulations containing hexapeptides of the amino acids to consider the effect of additive size on the solution species. We find that additives have limited impact on the formation of extended, liquid-like CaCO₃ networks in supersaturated solutions. Additives control the amount of (bi)carbonate in solution, but more importantly, they are able to stabilize these networks on the time scales of the simulations. This is achieved by coordinating the networks and assembled additive clusters in solutions. The association leads to subtle changes in the coordination of CaCO₃ and reduced mobility of the cations. We find that the number of solute association sites and the size and topology of the additives are more important than their net charge. Our results help to understand why polymer additives are so effective at stabilizing dense liquid CaCO₃ phases.