Ion dissolution mechanism and kinetics at kink sites on NaCl surfaces

Controlled dissolution of a single ion from a salt interface

Interactions between monatomic ions and water molecules are fundamental to understanding the hydration of complex polyatomic ions and ionic process. Among the simplest and well-established ion-related reactions is dissolution of salt in water, which is an endothermic process requiring an increase in entropy. Extensive efforts have been made to date; however, most studies at single-ion level have been limited to theoretical approaches. Here, we demonstrate the salt dissolution process by manipulating a single water molecule at an under-coordinated site of a sodium chloride film. Manipulation of molecule in a controlled manner enables us to understand ion-water interaction as well as dynamics of water molecules at NaCl interfaces, which are responsible for the selective dissolution of anions. The water dipole polarizes the anion in the NaCl ionic crystal, resulting in strong anion-water interaction and weakening of the ionic bonds. Our results provide insights into a simple but important elementary step of the single-ion chemistry, which may be useful in ion-related sciences and technologies.

Ca-dimers, solvent layering, and dominant electrochemically active species in Ca(BH4)2 in THF

Divalent ions (Mg, Ca, and Zn) are being considered as competitive, safe, and earth-abundant alternatives to Li-ion electrochemistry, but present challenges for stable cycling due to undesirable interfacial phenomena. We explore the formation of electroactive species in the electrolyte Ca(BH4)2∣THF using molecular dynamics coupled with a continuum model of bulk and interfacial speciation. Free-energy analysis and unsupervised learning indicate a majority population of neutral Ca dimers and monomers with diverse molecular conformations and an order of magnitude lower concentration of the primary electroactive charged species - the monocation, CaBH[Formula: see text] - produced via disproportionation of neutral complexes. Dense layering of THF molecules within ~1 nm of the electrode surface strongly modulates local electrolyte species populations. A dramatic increase in monocation population in this interfacial zone is induced at negative bias. We see no evidence for electrochemical activity of fully-solvated Ca2+. The consequences for performance are discussed in light of this molecular-scale insight.

The elementary reactions for incorporation into crystals

The growth rates of crystals are largely dictated by the chemical reaction between solute and kinks, in which a solute molecule severs its bonds with the solvent and establishes new bonds with the kink. Details on this sequence of bond breaking and rebuilding remain poorly understood. To elucidate the reaction at the kinks we employ four solvents with distinct functionalities as reporters on the microscopic structures and their dynamics along the pathway into a kink. We combine time-resolved in situ atomic force microscopy and x-ray and optical methods with molecular dynamics simulations. We demonstrate that in all four solvents the solute, etioporphyrin I, molecules reach the steps directly from the solution; this finding identifies the measured rate constant for step growth as the rate constant of the reaction between a solute molecule and a kink. We show that the binding of a solute molecule to a kink divides into two elementary reactions. First, the incoming solute molecule sheds a fraction of its solvent shell and attaches to molecules from the kink by bonds distinct from those in its fully incorporated state. In the second step, the solute breaks these initial bonds and relocates to the kink. The strength of the preliminary bonds with the kink determines the free energy barrier for incorporation into a kink. The presence of an intermediate state, whose stability is controlled by solvents and additives, may illuminate how minor solution components guide the construction of elaborate crystal architectures in nature and the search for solution compositions that suppress undesirable or accelerate favored crystallization in industry.

Impact of Solution {Ba2+}:{SO42-} on Charge Evolution of Forming and Growing Barite (BaSO4) Crystals: A ζ-Potential Measurement Investigation

The impact of solution stoichiometry on formation of BaSO4 (barite) crystals and the development of surface charge was investigated at various predefined stoichiometries (raq = 0.01, 0.1, 1, 10, and 100, where raq = {Ba2+}:{SO42-}). Synthesis experiments and zeta potential (ζ-potential) measurements were conducted at a fixed initial degree of supersaturation (Ωbarite = 1000, where Ωbarite = {Ba2+}{SO42-}/Ksp), at circumneutral pH of ∼6, 0.02 M NaCl, and ambient temperature and pressure. Mixed-mode measurement-phase analysis light scattering (M3-PALS) showed that the particles stayed negative for raq < 1 during barite crystal formation and positive for raq > 1. At raq = 1, two populations with a positive or negative ζ-potential prevailed for ∼2.5 h before a population with a circumneutral ζ-potential (-10 to +10 mV) remained. We relate the observations of particle charge evolution to particle size and morphology evolution under the experimental conditions. Furthermore, we showed that the ζ-potential became more negative when the pH was increased for every raq. In addition, our results demonstrated that the type of monovalent background electrolyte did not influence the ζ-potential of barite crystals significantly, although NaCl showed slightly different behavior compared to KCl and NaNO3. Our results show the important role of surface charge (evolution) during ionic crystal formation under nonstoichiometric conditions. Moreover, our combined scanning electron microscopy and ζ-potential results imply that the surface charge during particle formation can be influenced by solution stoichiometry, besides the pH and ionic strength, and may aid in predicting the fate of barite in environmental settings and in understanding and improving industrial barite (surface chemistry) processes.

An exploration of machine learning models for the determination of reaction coordinates associated with conformational transitions

Determining collective variables (CVs) for conformational transitions is crucial to understanding their dynamics and targeting them in enhanced sampling simulations. Often, CVs are proposed based on intuition or prior knowledge of a system. However, the problem of systematically determining a proper reaction coordinate (RC) for a specific process in terms of a set of putative CVs can be achieved using committor analysis (CA). Identifying essential degrees of freedom that govern such transitions using CA remains elusive because of the high dimensionality of the conformational space. Various schemes exist to leverage the power of machine learning (ML) to extract an RC from CA. Here, we extend these studies and compare the ability of 17 different ML schemes to identify accurate RCs associated with conformational transitions. We tested these methods on an alanine dipeptide in vacuum and on a sarcosine dipeptoid in an implicit solvent. Our comparison revealed that the light gradient boosting machine method outperforms other methods. In order to extract key features from the models, we employed Shapley Additive exPlanations analysis and compared its interpretation with the "feature importance" approach. For the alanine dipeptide, our methodology identifies ϕ and θ dihedrals as essential degrees of freedom in the C7ax to C7eq transition. For the sarcosine dipeptoid system, the dihedrals ψ and ω are the most important for the cisαD to transαD transition. We further argue that analysis of the full dynamical pathway, and not just endpoint states, is essential for identifying key degrees of freedom governing transitions.

ATESA: An Automated Aimless Shooting Workflow

Transition path sampling methods are powerful tools for studying the dynamics of rare events in molecular simulations. However, these methods are generally restricted to experts with the knowledge and resources to properly set up and analyze the often hundreds of thousands of simulations that constitute a complete study. Aimless Transition Ensemble Sampling and Analysis (ATESA) is a new open-source software program written in Python that automates a full transition path sampling workflow based on the aimless shooting algorithm, streamlining the process and reducing the barrier to use for researchers new to this approach. This introduction to ATESA includes a demonstration of a complete transition path sampling process flow for an example reaction, including finding an initial transition state, sampling with aimless shooting, building a reaction coordinate with inertial likelihood maximization, verifying that coordinate with committor analysis, and measuring the reaction energy profile with umbrella sampling. We also describe our implementation of a termination criterion for aimless shooting based on the Godambe information calculated during model building with likelihood maximization as well as a novel approach to constraining simulations to the desired rare event pathway during umbrella sampling.

Calcite Kinks Grow via a Multistep Mechanism

The classical model of crystal growth assumes that kinks grow via a sequence of independent adsorption events where each solute transitions from the solution directly to the crystal lattice site. Here, we challenge this view by showing that some calcite kinks grow via a multistep mechanism where the solute adsorbs to an intermediate site and only transitions to the lattice site upon the adsorption of a second solute. We compute the free energy curves for Ca and CO₃ ions adsorbing to a large selection of kink types, and we identify kinks terminated both by Ca ions and by CO₃ ions that grow in this multistep way.

Obtaining Consistent Free Energies for Ion Binding at Surfaces from Solution: Pathways versus Alchemy for Determining Kink Site Stability

Ion incorporation or removal from a solid at the interface with solution is a fundamental part of crystal growth. Despite this, there have been few quantitative determinations of the thermodynamics for such processes from atomistic molecular dynamics due to the associated technical challenges. In this study, we compute the free energies for ion removal from kink sites at the interface between NaCl and water as an illustrative example. To examine the influence of the free energy technique used, we compare methods that follow an explicit pathway for dissolution with those that focus on the thermodynamics of the initial and final states using metadynamics and free energy perturbation, respectively. While the initial results of the two approaches are found to be completely different, it is demonstrated that the thermodynamics can be reconciled with appropriate corrections for the standard states, thus illustrating the need for caution in interpreting raw free energy curves for ion binding as widely found in the literature. In addition, a new efficient approach is introduced to correct for the system size dependence of kink site energies both due to the periodic interaction of charges in an inhomogeneous dielectric system and due to the dipolar interactions between pairs of kinks along a row. Ultimately, it is shown that with suitable care, both classes of free energy techniques are capable of producing kink site stabilities that are consistent with the solubility of the underlying bulk solid. However, the precise values for individual kink sites exhibit a small systematic offset, which can be ascribed to the contribution of the interfacial potential to the pathway-based results. For the case of NaCl, the free energies of the kink sites relative to a 1 M aqueous solution for Na⁺ and Cl^- are found to be surprisingly different and of opposite sign, despite the ions having very similar hydration free energies.

Solubility of Organic Salts in Solvent-Antisolvent Mixtures: A Combined Experimental and Molecular Dynamics Simulations Approach

We combine molecular dynamics simulations with experiments to estimate solubilities of an organic salt in complex growth environments. We predict the solubility by simulations of the growth and dissolution of ions at the crystal surface kink sites at different solution concentrations. Thereby, the solubility is identified as the solution's salt concentration, where the energy of the ion pair dissolved in solution equals the energy of the ion pair crystallized at the kink sites. The simulation methodology is demonstrated for the case of anhydrous sodium acetate crystallized from various solvent-antisolvent mixtures. To validate the predicted solubilities, we have measured the solubilities of sodium acetate in-house, using an experimental setup and measurement protocol that guarantees moisture-free conditions, which is key for a hygroscopic compound like sodium acetate. We observe excellent agreement between the experimental and the computationally evaluated solubilities for sodium acetate in different solvent-antisolvent mixtures. Given the agreement and the rich data the simulations produce, we can use them to complement experimental tasks, which in turn will reduce time and capital in the design of complicated industrial crystallization processes of organic salts.

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.

Anomalous Phase Behaviors of Monolayer NaCl Aqueous Solutions Induced by Effective Coulombic Interactions within Angstrom-Scale Slits

Interests in subnanofluidic devices have called for molecular dynamics (MD) simulation studies of the thermodynamic behavior of monolayer salt solution within angstrom-scale slits. However, it still remains a grand challenge to accurately describe the Coulombic interactions by incorporating the effects of charge transfer and electronic dielectric screening. Herein, by using the electronic continuum model, where the effective ion charges are fine-tuned with a scaling factor of λ, we present simulation evidence that the effective Coulombic interactions among Na⁺/Cl^- ions can strongly affect the behavior of monolayer ionic aqueous solution. Our microsecond-scale MD simulations show that only the counterions with moderate effective charges (0.3 ≤ λ ≤ 0.8) can dissolve in monolayer water, whereas the high effective charges (λ ≥ 0.85) induce ions to assemble into monolayer nanocrystals, and ions with the low effective charges (λ ≤ 0.2) exhibit gas-like nanobubble. These findings could provide deeper insights into the physical chemistry behind subnanofluidic iontronic devices.

Explaining reaction coordinates of alanine dipeptide isomerization obtained from deep neural networks using Explainable Artificial Intelligence (XAI)

A method for obtaining appropriate reaction coordinates is required to identify transition states distinguishing product and reactant in complex molecular systems. Recently, abundant research has been devoted to obtaining reaction coordinates using artificial neural networks from deep learning literature, where many collective variables are typically utilized in the input layer. However, it is difficult to explain the details of which collective variables contribute to the predicted reaction coordinates owing to the complexity of the nonlinear functions in deep neural networks. To overcome this limitation, we used Explainable Artificial Intelligence (XAI) methods of the Local Interpretable Model-agnostic Explanation (LIME) and the game theory-based framework known as Shapley Additive exPlanations (SHAP). We demonstrated that XAI enables us to obtain the degree of contribution of each collective variable to reaction coordinates that is determined by nonlinear regressions with deep learning for the committor of the alanine dipeptide isomerization in vacuum. In particular, both LIME and SHAP provide important features to the predicted reaction coordinates, which are characterized by appropriate dihedral angles consistent with those previously reported from the committor test analysis. The present study offers an AI-aided framework to explain the appropriate reaction coordinates, which acquires considerable significance when the number of degrees of freedom increases.

Selective desolvation in two-step nucleation mechanism steers crystal structure formation

The two-step nucleation (TSN) theory and crystal structure prediction (CSP) techniques are two disjointed yet popular methods to predict nucleation rate and crystal structure, respectively. The TSN theory is a well-established mechanism to describe the nucleation of a wide range of crystalline materials in different solvents. However, it has never been expanded to predict the crystal structure or polymorphism. On the contrary, the existing CSP techniques only empirically account for the solvent effects. As a result, the TSN theory and CSP techniques continue to evolve as separate methods to predict two essential attributes of nucleation - rate and structure. Here we bridge this gap and show for the first time how a crystal structure is formed within the framework of TSN theory. A sequential desolvation mechanism is proposed in TSN, where the first step involves partial desolvation to form dense clusters followed by selective desolvation of functional groups directing the formation of crystal structure. We investigate the effect of the specific interaction on the degree of solvation around different functional groups of glutamic acid molecules using molecular simulations. The simulated energy landscape and activation barriers at increasing supersaturations suggest sequential and selective desolvation. We validate computationally and experimentally that the crystal structure formation and polymorph selection are due to a previously unrecognized consequence of supersaturation-driven asymmetric desolvation of molecules.

Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics

Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the "spectral gap optimization of order parameters (SGOOP)" approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. We believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.

Two-dimensional monolayer salt nanostructures can spontaneously aggregate rather than dissolve in dilute aqueous solutions

It is well known that NaCl salt crystals can easily dissolve in dilute aqueous solutions at room temperature. Herein, we reported the first computational evidence of a novel salt nucleation behavior at room temperature, i.e., the spontaneous formation of two-dimensional (2D) alkali chloride crystalline/non-crystalline nanostructures in dilute aqueous solution under nanoscale confinement. Microsecond-scale classical molecular dynamics (MD) simulations showed that NaCl or LiCl, initially fully dissolved in confined water, can spontaneously nucleate into 2D monolayer nanostructures with either ordered or disordered morphologies. Notably, the NaCl nanostructures exhibited a 2D crystalline square-unit pattern, whereas the LiCl nanostructures adopted non-crystalline 2D hexagonal ring and/or zigzag chain patterns. These structural patterns appeared to be quite generic, regardless of the water and ion models used in the MD simulations. The generic patterns formed by 2D monolayer NaCl and LiCl nanostructures were also confirmed by ab initio MD simulations. The formation of 2D salt structures in dilute aqueous solution at room temperature is counterintuitive. Free energy calculations indicated that the unexpected spontaneous salt nucleation behavior can be attributed to the nanoscale confinement and strongly compressed hydration shells of ions.

Aqueous solutions under nanoscale confinement exhibit interesting physicochemical properties. This work reports evidence on the spontaneous formation of two-dimensional alkali chloride crystalline/non-crystalline nanostructures in dilute aqueous solution under nanoscale confinement by computer simulations.

A survey on fractionation: the optimal control of distilling in batch and semibatch configurations

Since the middle of the last century, discussion about the operation of discontinuous fractionation to meet multifarious goals, such as product purity and recovery rate, by monitoring process variables including reflux or/and heat duty, is been on. The engineering practice showed intolerable events to occur; hereof the operation must be supervised, which makes it difficult to be in agreement with the batch distillation objectives. Hence, to uphold the effectuation of new operating policies into the industrial “know-how” techniques, different optimal control strategies can be conceived. The objective of this work is to offer a literature survey on the investigations of optimal control functioning for selected simple distillation column configurations employed in batch/semibatch distillation of homogeneous/reactive mixtures, as well as the approaches used in this regard. Available optimal control schemes have been reviewed in detail, emphasizing its major assets.

Stability of Two-Dimensional Ionic Clusters at Solid-Liquid Interfaces

The stability of two-dimensional clusters (2DCs) at the interface between ionic crystals and their solutions was investigated by molecular dynamics simulations. We found that 2DCs show a remarkable feature of odd-even alternation in stability. In NaCl and NaBr systems, the clusters containing an odd number of ions are more stable than those with an even number of ions, while in KCl systems, it is the other way round. Accordingly, the stability of water molecules in the first hydration shell of 2DCs also shows an odd-even alternation, which is consistent with the associated 2DCs. The odd-even alternation is discussed based on a competition mechanism between two factors: the Coulomb repulsion in charged 2DCs and the interaction between charges and water dipoles. Our discussion indicates that this odd-even alternation should be a universal feature in similar systems and would be important for understanding the nucleation and crystallization of solutions on ionic crystal surfaces.

Anion-Assisted Delivery of Multivalent Cations to Inert Electrodes

To understand and control key electrochemical processes-metal plating, corrosion, intercalation, etc.-requires molecular-scale details of the active species at electrochemical interfaces and their mechanisms for desolvation from the electrolyte. Using free energy sampling techniques we reveal the interfacial speciation of divalent cations in ether-based electrolytes and mechanisms for their delivery to an inert graphene electrode interface. Surprisingly, we find that anion solvophobicity drives a high population of anion-containing species to the interface that facilitate the delivery of divalent cations, even to negatively charged electrodes. Our simulations indicate that cation desolvation is greatly facilitated by cation-anion coupling. We propose anion solvophobicity as a molecular-level descriptor for rational design of electrolytes with increased efficiency for electrochemical processes limited by multivalent cation desolvation.

Absolute chemical potentials for complex molecules in fluid phases: A centroid reference for predicting phase equilibria

Solid-fluid phase equilibria are difficult to predict in simulations because bound degrees of freedom in the crystal phase must be converted to free translations and rotations in the fluid phase. Here, we avoid the solid-to-fluid transformation step by starting with chemical potentials for two reference systems, one for the fluid phase and one for the solid phase. For the solid, we start from the Einstein crystal and transform to the fully interacting molecular crystal. For the fluid phase, we introduce a new reference system, the "centroid," and then transform to gas phase molecules. We illustrate the new calculations by predicting the sublimation vapor pressure of succinic acid in the temperature range of 300 K-350 K.

CrystalGrower: a generic computer program for Monte Carlo modelling of crystal growth

A Monte Carlo crystal growth simulation tool, CrystalGrower, is described which is able to simultaneously model both the crystal habit and nanoscopic surface topography of any crystal structure under conditions of variable supersaturation or at equilibrium. This tool has been developed in order to permit the rapid simulation of crystal surface maps generated by scanning probe microscopies in combination with overall crystal habit. As the simulation is based upon a coarse graining at the nanoscopic level features such as crystal rounding at low supersaturation or undersaturation conditions are also faithfully reproduced. CrystalGrower permits the incorporation of screw dislocations with arbitrary Burgers vectors and also the investigation of internal point defects in crystals. The effect of growth modifiers can be addressed by selective poisoning of specific growth sites. The tool is designed for those interested in understanding and controlling the outcome of crystal growth through a deeper comprehension of the key controlling experimental parameters.

Generic in silico methodology – CrystalGrower – for simulating crystal habit and nanoscopic surface topology to determine crystallisation free energies.

Olanzapine crystal symmetry originates in preformed centrosymmetric solute dimers

The symmetries of a crystal are notoriously uncorrelated to those of its constituent molecules. This symmetry breaking is typically thought to occur during crystallization. Here we demonstrate that one of the two symmetry elements of olanzapine crystals, an inversion centre, emerges in solute dimers extant in solution prior to crystallization. We combine time-resolved in situ scanning probe microscopy to monitor the crystal growth processes with all-atom molecular dynamics simulations. We show that crystals grow non-classically, predominantly by incorporation of centrosymmetric dimers. The growth rate of crystal layers exhibits a quadratic dependence on the solute concentration, characteristic of the second-order kinetics of the incorporation of dimers, which exist in equilibrium with a majority of monomers. We show that growth by dimers is preferred due to overwhelming accumulation of adsorbed dimers on the crystal surface, where it is complemented by dimerization and expedites dimer incorporation into growth sites.

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.

Ion Distribution and Hydration Structure at Solid-Liquid Interface between NaCl Crystal and Its Solution

The interface structure between NaCl crystal and its solution has been investigated at the saturated concentration of 298 K by molecular dynamics simulations. We have found that there are many fine structures at this complex interface. Near the surface of crystal, most of Na⁺ only coordinate with water molecules, while almost all Cl^- coordinate with Na⁺ in addition to water molecules. An ion coordinating with more water molecules is farther away from the epitaxial position of lattice. As approaching to the interface, the first hydration shell of ions has the tendency of being ordered, while the orientation of dipole of water molecules in the first hydration shell becomes more disordered than that in the solution. Generally, the first hydration shell of Na⁺ is less affected by nearest Cl^-, whereas the first hydration shell of Cl^- is significantly affected by nearest Na⁺.

Solvent fluctuations in the solvation shell determine the activation barrier for crystal growth rates

Solution crystallization is a common technique to grow advanced, functional crystalline materials. Supersaturation, temperature, and solvent composition are known to influence the growth rates and thereby properties of crystalline materials; however, a satisfactory explanation of how these factors affect the activation barrier for growth rates has not been developed. We report here that these effects can be attributed to a previously unrecognized consequence of solvent fluctuations in the solvation shell of solute molecules attaching to the crystal surface. With increasing supersaturation, the average hydration number of the glutamic acid molecule decreases and can reach an asymptotic limit corresponding to the number of adsorption sites on the molecule. The hydration number of the glutamic acid molecule also fluctuates due to the rapid exchange of solvent in the solvation shell and local variation in the supersaturation. These rapid fluctuations allow quasi-equilibrium between fully solvated and partially desolvated states of molecules, which can be used to construct a double-well potential and thereby to identify the transition state and the required activation barrier. The partially desolvated molecules are not stable and can attach spontaneously to the crystal surface. The activation barrier versus hydration number follows the Evans-Polanyi relation. The predicted absolute growth rates of the α-glutamic acid crystal at lower supersaturations are in reasonable agreement with the experimental observations.

Reconsidering Calcium Dehydration as the Rate-Determining Step in Calcium Mineral Growth

The dehydration of cations is generally accepted as the rate-limiting step in many processes. Molecular dynamics (MD) can be used to investigate the dynamics of water molecules around cations, and two different methods exist to obtain trajectory-based water dehydration frequencies. Here, these two different post-processing methods (direct method versus survival function) have been implemented to obtain calcium dehydration frequencies from a series of trajectories obtained using a range of accepted force fields. None of the method combinations reproduced the commonly accepted experimental water exchange frequency of 10-8.2 s-1. Instead, our results suggest much faster water dynamics, comparable with more accurate ab initio MD simulations and with experimental values obtained using neutron scattering techniques. We obtained the best agreement using the survival function method to characterize the water dynamics, and we show that different method combinations significantly affect the outcome. Our work strongly suggests that the fast water exchange kinetics around the calcium ions is not rate-limiting for reactions involving dissolved/solvated calcium. Our results further suggest that, for alkali and most of the earth alkali metals, mechanistic rate laws for growth, dissolution, and adsorption, which are based on the principle of rate-limiting cation dehydration, need careful reconsideration.

Inhibition of Spiral Growth and Dissolution at the Brushite (010) Interface by Chondroitin 4-Sulfate

Modulation of mineralization and demineralization of calcium phosphates (Ca-Ps) with organic macromolecules is a critical process which prevents human kidney stone disease. As a long unbranched polysaccharide of urinary glycosaminoglycans, chondroitin 4-sulfate (Ch4S) has been shown to play an essential role in inhibiting the formation of kidney stones. However, the mechanism of the role of Ch4S remains poorly understood. Here, we used in situ atomic force microscopy to observe the growth and dissolution of spirals on brushite (CaHPO₄·2H₂O) (010) surfaces. The results show that Ch4S preferentially inhibits the [101]_Cc step growth/dissolution by step pinning. This step-specific effect appears to be related to specific binding of Ch4S to Ca sites, as the observed inhibition is not seen in other crystallographic directions where there are fewer Ca terminations. Moreover, Ch4S promotes an increase in the terrace width of [101̅]_Cc by the modification of the interfacial energies of the step edge. These in vitro direct observations of Ch4S modulating brushite mineralization and demineralization reveal a dual control of both step kinetics and interfacial energy.

Free-Energy Landscape of the Dissolution of Gibbsite at High pH

The individual elementary reactions involved in the dissolution of a solid into solution remain mostly speculative due to a lack of direct experimental probes. In this regard, we have applied atomistic simulations to map the free-energy landscape of the dissolution of gibbsite from a step edge as a model of metal hydroxide dissolution. The overall reaction combines kink formation and kink propagation. Two individual reactions were found to be rate-limiting for kink formation, that is, the displacement of Al from a step site to a ledge adatom site and its detachment from ledge/terrace adatom sites into the solution. As a result, a pool of mobile and labile adsorbed species, or adatoms, exists before the release of Al into solution. Because of the quasi-hexagonal symmetry of gibbsite, kink site propagation can occur in multiple directions. Overall, our results will enable the development of microscopic mechanistic models of metal oxide dissolution.