Search

Metal carbonates, which are ubiquitous in the near-surface mineral record, are a major product of biomineralizing organisms and serve as important targets for capturing anthropogenic CO 2 emissions. However, pathways of carbonate mineralization typically diverge from classical predictions due to the involvement of disordered precursors, such as the dense liquid phase (DLP), yet little is known about DLP formation or solidification processes. Using in situ methods we report that a highly hydrated bicarbonate DLP forms via liquid-liquid phase separation and transforms into hollow hydrated amorphous CaCO 3 particles. Acidic proteins and polymers extend DLP lifetimes while leaving the pathway and chemistry unchanged. Molecular simulations suggest that the DLP forms via direct condensation of solvated Ca² + ⋅(HCO 3 - ) 2 complexes that react due to proximity effects in the confined DLP droplets. Our findings provide insight into CaCO 3 nucleation that is mediated by liquid-liquid phase separation, advancing the ability to direct carbonate mineralization and elucidating an often-proposed complex pathway of biomineralization., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

The phenomenon of underscreening in concentrated electrolyte solutions leads to a larger decay length of the charge-charge correlation than the prediction of Debye-Hückel (DH) theory and has found a resurgence of both theoretical and experimental interest in the chemical physics community. To systematically understand and investigate this phenomenon in electrolytes requires a theory of concentrated electrolytes to describe charge-charge correlations beyond the DH theory. We review the theories of electrolytes that can transition from the DH limit to concentrations where charge correlations dominate, giving rise to underscreening and the associated Kirkwood Transitions (KTs). In this perspective, we provide a conceptual approach to a theoretical formulation of electrolyte solutions that exploits the competition between molecular-informed short-range (SR) and long-range interactions. We demonstrate that all deviations from the DH limit for real electrolyte solutions can be expressed through a single function ΣQ that can be determined both theoretically and numerically. Importantly, ΣQ can be directly related to the details of SR interactions and, therefore, can be used as a tool to understand how differences in representations of interaction can influence collective effects. The precise function form of ΣQ can be inferred through a Gaussian field theory of both the number and charge densities. The resulting formulation is validated by experiment and can accurately describe the collective phenomenon of screening in concentrated bulk electrolytes. Importantly, the Gaussian field theory predictions of the screening lengths appear to be less than ∼1 nm at concentrations above KTs., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

Understanding defect healing is necessary to control the electronic and optoelectronic performance of devices based on nanoparticle (NP) superlattices. However, a key challenge remains to understand how NP interactions and the resulting dynamics are coupled to defect self-elimination during assembly processes. Additional degrees of freedom that account for the anisotropic nature of NPs associated with rotational dynamics and torques further complicate the challenge. Here, we investigate nanocube (NC) superlattices by employing liquid-phase transmission electron microscopy, continuum theory, and molecular dynamics simulations. Our detailed analyses reveal that interparticle forces and torques due to ligand interactions dominate those from Brownian motions and van der Waals interactions. More importantly, NC translations and rotations induced by unbalanced forces and torques are transmitted to neighboring NCs, prompting "chain interactions" in a two-dimensional (2D) network and expediting self-elimination. The mechanistic understanding will further enable the design and fabrication of defect-free superlattices as well as those with tailored defects via assembly of anisotropic particles.

We establish the connection between the measured small angle x-ray scattering signal and the charge-charge correlations underlying Kirkwood transitions (KTs) in 1:1, 2:1, and 3:1 aqueous electrolytes. These measurements allow us to obtain underscreening lengths for bulk electrolytes independently verified by theory and simulations. Furthermore, we generalize the concept of KTs beyond those theoretically predicted for 1:1 electrolytes, which involves the inverse screening length, a0, and the inverse periodicity length, Q0. Above the KTs, we find a universal scaling of a0∝c-ζ/3 and Q0 ∝ c1/3 for the studied electrolyte solutions, where ζ is the ionic strength factor., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

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

Aluminum hydroxide polymorphs are of widespread importance yet their kinetics of nucleation and growth remain beyond the reach of current models. Here we attempt to unveil the reaction processes underlying the polymorphs formation at high chemical potential. We examine their formation in-situ from supersaturated alkaline sodium aluminate solutions using deuteration and time-resolved neutron pair distribution function analyses, which indicate the formation of individual Al(OD) 3 layers as an intermediate particle phase. These layers ultimately stack to form gibbsite- or bayerite-like layered heterostructures. Ex-situ characterization of the recovered precipitates using 27 Al magic angle spinning nuclear magnetic resonance spectroscopy, Raman, X-ray diffraction, and scanning electron microscopy, suggests the presence of additional intermediate states, an amorphous compound bearing both tetrahededrally- and penta-coordinated Al 3+ . These observations reveal the complex pathways to form Al(OD) 3 monolayers via either transient pentacoordinate species or amorphous-to-ordered transitions. The subsequent crystallization of admixed gibbsite/bayerite is followed by an Al(OD) 3 monolayer attachment process., (© 2024. The Author(s).)

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1 kT at separations of 10 nm and 3 kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely, dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data are relevant to another type of non-DLVO interaction, namely, the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. Our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

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

The van der Waals interaction between colloids and nanoparticles is one of the key components to understanding particle aggregation, attachment, and assembly. While the ubiquity of anisotropic particle shapes and surface roughness is well-recognized in nanocrystalline materials, the effects of both on van der Waals forces and torques have not been adequately investigated. In this study, we develop a numerical scheme to determine the van der Waals forces and torques between cubic particles with multiple configurations and relative orientations. Our results show that the van der Waals torque due to anisotropic particle shapes is appreciable at nearly all configurations and mutual angles, outcompeting Brownian torque for various materials systems and conditions. Surface roughness enhances this particle shape effect, resulting in stronger van der Waals interactions ascribed to protrusions on the surfaces. Moreover, a scaling analysis indicates that the surface roughness alters the separation dependence of the van der Waals force and, more importantly, significantly influences the dynamics of two approaching particles. Our results clearly demonstrate that surface roughness and anisotropic shape play a crucial role in the energetics and kinetics of various particle-scale and emergent phenomena, such as crystal growth by oriented attachment, nanomaterials synthesis and assembly, mud flow rheology, as well as the deposition of natural nanocrystals within the subsurface., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Battelle Memorial Institute, The Author(s). Published by Elsevier Inc. All rights reserved.)

Predicting the behavior of oxyanions in radioactive waste stored at the Department of Energy legacy nuclear sites requires the development of novel analytical methods. This work demonstrates 15 N pulsed field gradient nuclear magnetic resonance spectroscopy to quantify the diffusivity of nitrite. Experimental results, supported by molecular dynamics simulations, indicate that the diffusivity of free hydrated nitrite exceeds that of free hydrated sodium despite the greater hydrodynamic radius of nitrite. Investigations are underway to understand how the compositional and dynamical heterogeneities of the ion networks at high concentrations affect rheological and transport properties.

Understanding multiple lengthscale correlations in the pair distribution functions (PDFs) of aq. electrolytes is a persistent challenge. Here, the coordination chemistry of polyoxoanions supports an ion-network of cation-coordination polyhedra in NaNO 3(aq) and NaNO 2(aq) that induce long-range solution structure. Oxygen correlations associated with Na + -coordination polyhedra have two characteristics lengthscales; 3.5-5.5 Å and 5.5-7.5 Å, the latter solely associated oligomers. The PDF contraction between 5.5-7.5 Å observed in many electrolytes is attributed to the distinct O⋯O correlation found in dimers and dimer subunits within oligomers.

Predicting nanoparticle aggregation and attachment phenomena requires a rigorous understanding of the interplay among crystal structure, particle morphology, surface chemistry, solution conditions, and interparticle forces, yet no comprehensive picture exists. We used an integrated suite of experimental, theoretical, and simulation methods to resolve the effect of solution pH on the aggregation of boehmite nanoplatelets, a case study with important implications for the environmental management of legacy nuclear waste. Real-time observations showed that the particles attach preferentially along the (010) planes at pH 8.5 and the (101) planes at pH 11. To rationalize these results, we established the connection between key physicochemical phenomena across the relevant length scales. Starting from molecular-scale simulations of surface hydroxyl reactivity, we developed an interfacial-scale model of the corresponding electrostatic potentials, with subsequent particle-scale calculations of the resulting driving forces allowing successful prediction of the attachment modes. Finally, we scaled these phenomena to understand the collective structure at the aggregate-scale . Our results indicate that facet-specific differences in surface chemistry produce heterogeneous surface charge distributions that are coupled to particle anisotropy and shape-dependent hydrodynamic forces, to play a key role in controlling aggregation behavior.

Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.g., ReaxFF) by omitting explicit dependence on the atom coordination and employing a small parameter set based on Lennard-Jones, Gaussian, and Stillinger-Weber potentials. The model was parametrized from AIMD simulation data and is used to model aluminate reactivity in sodium hydroxide solutions with extensive validation against experimental radial distribution functions, computed free energy profiles for oligomerization, and formation energies. The model enables simulation of early stage Al(OH) 3 nucleation which has significant relevance to industrial processing of aluminum and has a computational cost that is reduced by 1 order of magnitude relative to ReaxFF.

Hypothesis: The precipitation and dissolution of aluminum-bearing mineral phases in aqueous systems often proceed via changes in both aluminum coordination number and connectivity, complicating molecular-scale interpretation of the transformation mechanism. Here, the thermally induced transformation of crystalline sodium aluminum salt hydrate, a phase comprised of monomeric octahedrally coordinated aluminate which is of relevance to industrial aluminum processing, has been studied. Because intermediate aluminum coordination states during melting have not previously been detected, it is hypothesized that the transition to lower coordinated aluminum ions occurs within ahighly disordered quasi-two-dimensional phase at the solid-solution interface., Experiments and Simulations: In situ X-ray diffraction (XRD), Raman and 27 Al nuclear magnetic resonance (NMR) spectroscopy were used to monitor the melting transition of nonasodium aluminate hydrate (NSA, Na 9 [Al(OH) 6 ] 2 ·3(OH)·6H 2 O). A mechanistic interpretation was developed based on complementary classical molecular dynamics (CMD) simulations including enhanced sampling. A reactive forcefield was developed to bridge speciation in the solution and in the solid phase., Findings: In contrast to classical dissolution, aluminum coordination change proceeds through a dynamically stabilized ensemble of intermediate states in a disordered layer at the solid-solution interface. In both melting and dissolution of NSA, octahedral, monomeric aluminum transition through an intermediate of pentahedral coordination. The intermediate dehydroxylates to form tetrahedral aluminate (Al(OH) 4 - ) in the liquid phase. This coordination change is concomitant with a breaking of the ionic aluminate-sodium ionlinkages. The solution phase Al(OH) 4 - ions subsequently polymerize into polynuclear aluminate ions. However, there are some differences between bulk melting and interfacial dissolution, with the onset of the surface-controlled process occurring at a lower temperature (∼30 °C) and the coordination change taking place more gradually as a function of temperature. This work to determine the local structure and dynamics of aluminum in the disordered layer provides a new basis to understand mechanisms controlling aluminum phase transformations in highly alkaline solutions., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

A primary means to generate hydrated electrons in laboratory experiments is excitation to the charge-transfer-to-solvent (CTTS) state of a solute such as I - (aq), but this initial step in the genesis of e - (aq) has never been simulated directly using ab initio molecular dynamics. We report the first such simulations, combining ground- and excited-state simulations of I - (aq) with a detailed analysis of fluctuations in the Coulomb potential experienced by the nascent solvated electron. What emerges is a two-step picture of the evolution of e - (aq) starting from the CTTS state: I - (aq) + hν → I -* (aq) → I • (aq) + e - (aq). Notably, the equilibrated ground state of e - (aq) evolves from I -* (aq) without any nonadiabatic transitions, simply as a result of solvent reorganization. The methodology used here should be applicable to other photochemical electron transfer processes in solution, an important class of problems directly relevant to photocatalysis and energy transfer.

Ion trajectory simulation in mass spectrometry systems from injection to detection is technically challenging but very important for better understanding the ion dynamics in instrument development. Here, we present SimELIT ( Sim ulator of E ulerian and L agrangian I on T rajectories), a novel ion trajectory simulation platform. SimELIT is built upon a suite of multiphysics solvers compiled into OpenFOAM (an open-source numerical solver library particularly used for computational mechanics), with a simple web-based graphical user interface (GUI) allowing users to define the details of OpenFOAM cases and run simulations. SimELIT is a modular program and can provide extensions of physics (e.g., gas flows, electrodynamic fields) and thus enable ion trajectory simulations from the ion source to detector. The current version (SimELIT) provides two numerical solvers for ion trajectory simulations─(1) a Lagrangian particle tracker in vacuum and (2) a Eulerian ion density solver in background gas in the presence of electric fields. Here, we describe the architecture of SimELIT, including its use of Docker and the React Framework, and demonstrate the computation of ion trajectories of multiple m / z values in a static/linear voltage drop in vacuum (across a 1 m long flight tube). Further, the drift motion of ions under 1 Torr pressure conditions in a static background (N 2 ) gas through a 20 V/cm static electric field is shown. The results produced from SimELIT were compared with SIMION and theoretical estimates. In addition, we report the computation of ion trajectories in electrodynamic fields within a planar FAIMS device operating at atmospheric pressure.

Although nanometer-sized aluminum hydroxide clusters (i.e., ϵ-Al 13 , [Al 13 O 4 (OH) 24 (H 2 O) 12 ] 7+ ) command a central role in aluminum ion speciation and transformations between minerals, measurement of their translational diffusion is often limited to indirect methods. Here, 27 Al pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has been applied to the AlO 4 core of the ϵ-Al 13 cluster with complementary theoretical simulations of the diffusion coefficient and corresponding hydrodynamic radii from a boundary element-based calculation. The tetrahedral AlO 4 center of the ϵ-Al 13 cluster is symmetric and exhibits only weak quadrupolar coupling, which results in favorable T 1 and T 2 27 Al NMR relaxation coefficients for 27 Al PFGSTE NMR studies. Stokes-Einstein relationship was used to relate the 27 Al diffusion coefficient of the ϵ-Al 13 cluster to the hydrodynamic radius for comparison with theoretical simulations, dynamic light scattering from literature, and previously published 1 H PFGSTE NMR studies of chelated Keggin clusters. This first-of-its-kind observation proves that 27 Al PFGSTE NMR diffusometry can probe symmetric Al environments in polynuclear clusters of greater molecular weight than previously considered., (© 2021 John Wiley & Sons, Ltd.)

Aluminate salts precipitated from caustic alkaline solutions exhibit a correlation between the anionic speciation and the identity of the alkali cation in the precipitate, with the aluminate ions occurring either in monomeric (Al(OH) 4 - ) or dimeric (Al 2 O(OH) 6 2- ) forms. The origin of this correlation is poorly understood as are the roles that oligomeric aluminate species play in determining the solution structure, prenucleation clusters, and precipitation pathways. Characterization of aluminate solution speciation with vibrational spectroscopy results in spectra that are difficult to interpret because the ions access a diverse and dynamic configurational space. To investigate the Al(OH) 4 - and Al 2 O(OH) 6 2- anions within a well-defined crystal lattice, inelastic neutron scattering (INS) and Raman spectroscopic data were collected and simulated by density functional theory for K 2 [Al 2 O(OH) 6 ], Rb 2 [Al 2 O(OH) 6 ], and Cs[Al(OH) 4 ]·2H 2 O. These structures capture archetypal solution aluminate species: the first two salts contain dimeric Al 2 O(OH) 6 2- anions, while the third contains the monomeric Al(OH) 4 - anion. Comparisons were made to the INS and Raman spectra of sodium aluminate solutions frozen in a glassy state. In contrast to solution systems, the crystal lattice of the salts results in well-defined vibrations and associated resolved bands in the INS spectra. The use of a theory-guided analysis of the INS of this solid alkaline aluminate series revealed that differences were related to the nature of the hydrogen-bonding network and showed that INS is a sensitive probe of the degree of completeness and strength of the bond network in hydrogen-bonded materials. Results suggest that the ionic size may explain cation-specific differences in crystallization pathways in alkaline aluminate salts.

The geometric organization and force networks of 3D dense suspensions that exhibit both shear thinning and thickening have been examined as a function of varying strength of interparticle attractive interactions using lubrication flow discrete element simulations. Significant rearrangement of the geometric topology does not occur at either the local or global scale as these systems transition across the shear thinning and shear thickening regimes. In contrast, massive rearrangements in the balance of attractive, lubrication, and contact forces are observed with interesting behavior of network growth and competition. In agreement with prior work, in shear thinning regions the attractive force is dominant, however as the shear thickening region is approached there is growth of lubrication forces. Lubrication forces oppose the attraction forces, but as viscosity continues to increase under increasing shear stress, the lubrication forces are dominated by contact forces that also resist attraction. Contact forces are the dominant interactions during shear thickening and are an order of magnitude higher than their values in the shear-thinning regime. At high attractive interaction strength, contact networks can form even under shear thinning conditions, however high shear stress is still required before contact networks become the driving mechanism of shear thickening. Analysis of the contact force network during shear thickening generally indicates a uniformly spreading network that rapidly forms across empty domains; however the growth patterns exhibit structure that is significantly dependent upon the strength of interparticle interactions, indicating subtle variations in the mechanism of shear thickening.

Amongst the challenges for a variety of research fields are the visualization of solid-liquid interfaces and understanding how they are affected by the solution conditions such as ion concentrations, pH, ligands, and trace additives, as well as the underlying crystallography and chemistry. In this context, three-dimensional fast force mapping (3D FFM) has emerged as a promising tool for investigating solution structure at interfaces. This capability is based on atomic force microscopy (AFM) and allows the direct visualization of interfacial regions in three spatial dimensions with sub-nanometer resolution. Here we provide a detailed description of the experimental protocol for acquiring 3D FFM data. The main considerations for optimizing the operating parameters depending on the sample and application are discussed. Moreover, the basic methods for data processing and analysis are discussed, including the transformation of the measured instrument observables into tip-sample force maps that can be linked to the local solution structure. Finally, we shed light on some of the outstanding questions related to 3D FFM data interpretation and how this technique can become a central tool in the repertoire of surface science.

Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined.In this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion-ion, ion-water, and water-water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions.DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air-water interface, ions in bulk, and at the air-water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

The molecular speciation of aluminum (Al 3+ ) in alkaline solutions is fundamental to its precipitation chemistry within a number of industrial applications that include ore refinement and industrial processing of Al wastes. Under these conditions, Al 3+ is predominantly Al(OH) 4 - , while at high [Al 3+ ] dimeric species are also known to form. To date, the mechanism of dimer formation remains unclear and is likely influenced by complex ion···ion interactions. In the present work, we investigate a suite of potential dimerization pathways and the role of ion pairing on energetics using static DFT calculations and DFT and density functional tight binding molecular dynamics. Specific cation effects imparted by the background electrolyte cations Na + , Li + , and K + have been examined. Our simulations predict that, when the Al species are ion-paired with either cation, the formation of the oxo-bridged Al 2 O(OH) 6 2- is favored with respect to the dihydroxo-bridged Al 2 (OH) 8 2- , in agreement with previous spectroscopic work. The formation of both dimers first proceeds by bridging of two monomeric units via one hydroxo ligand, leading to a labile Al 2 (OH) 8 2- isomer. The effect of contact ion pairing of Li + and K + on the dimerization energetics is distinctly more favorable than that of Na + , which may have an effect on further oligomerization.

We report the heterolysis of molecular hydrogen under ambient conditions by the crystalline frustrated Lewis pair (FLP) 1-{2-[bis(pentafluorophenyl)boryl]phenyl}-2,2,6,6-tetramethylpiperidine (KCAT). The gas-solid reaction provides an approach to prepare the solvent-free, polycrystalline ion pair KCATH2 through a single crystal to single crystal transformation. The crystal lattice of KCATH2 increases in size relative to the parent KCAT by approximately 2%. Microscopy was used to follow the transformation of the highly colored red/orange KCAT to the colorless KCATH2 over a period of 2 h at 300 K under a flow of H 2 gas. There is no evidence of crystal decrepitation during hydrogen uptake. Inelastic neutron scattering employed over a temperature range from 4-200 K did not provide evidence for the formation of polarized H 2 in a precursor complex within the crystal at low temperatures and high pressures. However, at 300 K, the INS spectrum of KCAT transformed to the INS spectrum of KCATH2. Calculations suggest that the driving force is more favorable in the solid state compared to the solution or gas phase, but the addition of H 2 into the KCAT crystal is unfavorable. Ab Initio methods were used to calculate the INS spectra of KCAT, KCATH2, and a possible precursor complex of H 2 in the pocket between the B and N of crystalline KCAT. Ex-situ NMR showed that the transformation from KCAT to KCATH2 is quantitative and our results suggest that the hydrogen heterolysis process occurs via H 2 diffusion into the FLP crystal with a rate-limiting movement of H 2 from inactive positions to reactive sites.

Hypothesis: Understanding the stability and rheological behavior of suspensions composed of anisotropic particles is challenging due to the complex interplay of hydrodynamic and colloidal forces. We propose that orientationally-dependent interactions resulting from the anisotropic nature of non-spherical sub-units strongly influences shear-induced particle aggregation/fragmentation and suspension rheological behavior., Experiments: Wide-, small-, and ultra-small-angle X-ray scattering experiments were used to simultaneously monitor changes in size and fractal dimensions of boehmite aggregates from 6 to 10,000 Å as the sample was recirculated through an in-situ capillary rheometer. The latter also provided simultaneous suspension viscosity data. Computational fluid dynamics modeling of the apparatus provided a more rigorous analysis of the fluid flow., Findings: Shear-induced aggregation/fragmentation was correlated with a complicated balance between hydrodynamic and colloidal forces. Multi-scale fractal aggregates formed in solution but the largest could be fragmented by shear. Orientationally-dependent interactions lead to a relatively large experimental suspension viscosity when the hydrodynamic force was small compared to colloidal forces. This manifests even at low boehmite mass fractions., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020. Published by Elsevier Inc.)

Accurately predicting the molecular structure of solutions is a fundamental scientific challenge. Using quantum mechanical density functional theory (DFT) to make these predictions is hindered by significant variation depending on which DFT functional is used. Here, we present a simple metric that can determine the reliability of a DFT functional for predicting solvation structure. We then show that including a simple interaction term to correct this metric leads to quantitative agreement with experimental measurements of liquid structure. We demonstrate the utility of this method by using it to accurately describe the hydration structure around the Na + and K + ions as well as the structural properties of pure water with a computationally cheap functional.

This paper describes a formalism for extracting spatially varying transport coefficients from simulations of a molecular fluid in a nanochannel. This approach is applied to self-diffusion of a Lennard-Jones fluid confined between two parallel surfaces. A numerical grid is laid over the domain confining the fluid, and fluid properties are projected onto the grid cells. The time correlation functions between properties in different grid cells are calculated and can be used as the basis for a fitting procedure for extracting spatially varying diffusion coefficients from the simulation. Results for the Lennard-Jones system show that transport behavior varies sharply near the liquid-solid boundary and that the changes depend on the details of the liquid-solid interaction. A quantitative difference between the reduced and detailed models is discussed. It is found that the difference could be associated with assumptions about the form of the transport equations at molecular scales in lieu of problems with the method itself. The study suggests that this approach to fitting molecular simulations to continuum equations may guide the development of appropriate coarse-grained equations to model transport phenomena at nanometer scales.

We study the prototypical S N 2 reaction Cl - + CH 3 Cl → CH 3 Cl + Cl - in water using quantum mechanics/molecular mechanics (QM/MM) computer simulations with transition path sampling and inertial likelihood maximization. We have identified a new solvent coordinate to complement the original atom-exchange coordinate used in the classic analysis by Chandrasekhar, Smith, and Jorgensen [J. Am. Chem. Soc. 107, 154 (1985)]. The new solvent coordinate quantifies instantaneous solvent-induced polarization relative to the equilibrium average charge density at each point along the reaction pathway. On the basis of likelihood scores and committor distributions, the new solvent coordinate improves upon the description of solvent dynamical effects relative to previously proposed solvent coordinates. However, it does not increase the transmission coefficient or the accuracy of a transition state theory rate calculation.

Rate theories have found great utility across the chemical sciences by providing a physically transparent way to analyze dynamical processes. Here we demonstrate the benefits of using transition state theory and Marcus theory to study the rate of proton transfer in HCl solutions. By using long ab initio molecular dynamics simulations, we show that good agreement is obtained between these two different formulations of rate theory and how they can be used to study the pathways and lifetime of proton transfer in aqueous solution. Since both rate theory formulations utilize identical sets of molecular data, this provides a self-consistent theoretical picture of the rates of aqueous phase proton transfer. Specifically, we isolate and quantify the rates of proton transfer, ion-pair dissociation, and solvent exchange in concentrated HCl solutions. Our analysis predicts a concentration dependence to both proton transfer and ion-pairing. Moreover, our estimate of the lifetime for the Zundel species is 0.8 and 1.3 ps for 2 M and 8 M HCl, respectively. We demonstrate that concentration effects can indeed be quantified through the combination of state-of-the-art simulation and theory and provides a picture of the important correlations between the cation (hydronium) and the counterion in acid solutions.

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

The ability to reproduce the experimental structure of water around the sodium and potassium ions is a key test of the quality of interaction potentials due to the central importance of these ions in a wide range of important phenomena. Here, we simulate the Na + and K + ions in bulk water using three density functional theory functionals: (1) the generalized gradient approximation (GGA) based dispersion corrected revised Perdew, Burke, and Ernzerhof functional (revPBE-D3) (2) the recently developed strongly constrained and appropriately normed (SCAN) functional (3) the random phase approximation (RPA) functional for potassium. We compare with experimental X-ray diffraction (XRD) and X-ray absorption fine structure (EXAFS) measurements to demonstrate that SCAN accurately reproduces key structural details of the hydration structure around the sodium and potassium cations, whereas revPBE-D3 fails to do so. However, we show that SCAN provides a worse description of pure water in comparison with revPBE-D3. RPA also shows an improvement for K + , but slow convergence prevents rigorous comparison. Finally, we analyse cluster energetics to show SCAN and RPA have smaller fluctuations of the mean error of ion-water cluster binding energies compared with revPBE-D3.

Tetrahedrally coordinated aluminate Al(OH) 4 - and dialuminate Al 2 O(OH) 6 2- anions are considered to be major species in aluminum-rich alkaline solutions. However, their relative abundance remains difficult to spectroscopically quantify due to local structure similarities and poorly understood effects arising from extent of polymerization and counter-cations. To help unravel these relationships here we report detailed characterization of three solid-phase analogues as structurally and compositionally well-defined reference materials. We successfully synthesized a cesium salt of the aluminate monomer, CsAl(OH) 4 ·2H 2 O, for comparison to potassium and rubidium salts of the aluminate dimer, K 2 Al 2 O(OH) 6 , and Rb 2 Al 2 O(OH) 6 , respectively. Single crystal and powder X-ray diffraction methods clearly reveal the structure and purity of these materials for which a combination of 27 Al MAS-NMR, Al K-edge X-ray absorption and Raman/IR spectroscopies was then used to fingerprint the two major tetrahedrally coordinated Al species. The resulting insights into the effect of Al-O-Al bridge formation between aluminate tetrahedra on spectroscopic features may also be generalized to the many materials that are based on this motif.

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O) 6 3+ , accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH) 4 - requires an understanding of how these errors cancel in order to define the limits of accurately predicting Al27 chemical shielding in Al(OH) 4 - . In this work, we estimate the absolute shielding of the Al27 nucleus in Al(OH) 4 - at the coupled cluster level (515.1 ± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O) 6 3+ and Al(OH) 4 - ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the Al27 chemical shift of Al(OH) 4 - from first principles.

The intermediate-range (1-4 nm) correlation of cations, anions, and water in aqueous alkaline earth salt solutions is measured using synchrotron X-ray diffraction. We differentiate from the entire solution structure factor, S X ( Q ), a separate region at low Q (<1.5 Å -1 ) containing local diffraction maxima (prepeaks) that indicate nanometer-scale oscillatory behavior. These features are quantitatively reproduced with classical molecular dynamics simulations. At high concentrations, the prepeaks emerge from correlations arising from the existence of small quasi close-packed lattice-like structures comprised of cation hydration spheres. We also analyze the concentration dependence of the prepeak and discuss the overall results in light of the rich literature dealing with intermediate-range correlations underlying universal phenomena in concentrated electrolytes.

When hydrolyzable cations such as aluminum interact with solid-water interfaces, macroscopic interfacial properties ( e.g ., surface charge and potential) and interfacial phenomena ( e.g ., particle adhesion) become tightly linked with the microscopic details of ion adsorption and speciation. We use in situ atomic force microscopy to directly image individual aluminum ions at a mica-water interface and show how adsorbate populations change with pH and aluminum activity. Complementary streaming potential measurements then allow us to build a triple layer model (TLM) that links surface potentials to adsorbate populations, via equilibrium binding constants. Our model predicts that hydrolyzed species dominate the mica-water interface, even when unhydrolyzed species dominate the solution. Ab initio molecular dynamics (AIMD) simulations confirm that aluminum hydrolysis is strongly promoted at the interface. The TLM indicates that hydrolyzed adsorbates are responsible for surface-potential inversions, and we find strong correlations between hydrolyzed adsorbates and particle-adhesion forces, suggesting that these species mediate adhesion by chemical bridging.

The interplay between crystal and solvent structure, interparticle forces and ensemble particle response dynamics governs the process of crystallization by oriented attachment (OA), yet a quantitative understanding is lacking. Using ZnO as a model system, we combine in situ TEM observations of single particle and ensemble assembly dynamics with simulations of interparticle forces and responses to relate experimentally derived interparticle potentials to the underlying interactions. We show that OA is driven by forces and torques due to a combination of electrostatic ion-solvent correlations and dipolar interactions that act at separations well beyond 5 nm. Importantly, coalignment is achieved before particles reach separations at which strong attractions drive the final jump to contact. The observed barrier to attachment is negligible, while dissipative factors in the quasi-2D confinement of the TEM fluid cell lead to abnormal diffusivities with timescales for rotation much less than for translation, thus enabling OA to dominate.

Confinement of water in sub-nanometer pores strongly alters its vibrational dynamics from that of bulk water. The effect of confinement can, furthermore, be finely tuned by small changes in the size and symmetry of the confining pore. Using inelastic neutron scattering (INS), we recently studied the dynamics of water confined in the channels of beryl and cordierite in which, at low temperatures, water shows similar behavior, indicating an absence of hydrogen bonds acting on the water molecule and a shallow water potential in the direction perpendicular to the channels. In addition, we observed multiple tunneling modes (between 0.66 and 14.7 meV) in the INS spectra of beryl due to transitions between the split ground-state of the water protons. Here, we present a study of (i) the effect of pressure on the dynamics of water in beryl, (ii) the dynamics of water in beryl containing alkali metals (which results in changing the orientation of the water molecule in the crystal), and (iii) the dynamics of water in cordierite at low energies. We found a shift in the tunneling and vibrational modes of water in beryl to higher energies at 22 kbar relative to 1 bar. No tunneling modes were observed for water in cordierite and type-II water in beryl. Therefore, we conclude that very small differences in the size and structure of the pores and the orientation of the water molecule in these minerals result in changes in the potential of the water protons and drastic changes in the confined water dynamics.

A systematic analysis of the hydration structure of Cs + ions in solution is derived from simulations carried out using a series of molecular models built upon a hierarchy of approximate representations of many-body effects in ion-water interactions. It is found that a pairwise-additive model, commonly used in biomolecular simulations, provides poor agreement with experimental X-ray spectra, indicating an incorrect description of the underlying hydration structure. Although the agreement with experiment improves in simulations with a polarizable model, the predicted hydration structure is found to lack the correct sequence of water shells. Progressive inclusion of explicit many-body effects in the representation of Cs + -water interactions as well as accounting for nuclear quantum effects is shown to be necessary for quantitatively reproducing the experimental X-ray spectra. Besides emphasizing the importance of many-body effects, these results suggest that molecular models rigorously derived from many-body expansions hold promise for realistic simulations of aqueous solutions.

Highly concentrated particle suspensions (also called slurries) can undergo a sharp increase in viscosity, or shear thickening, under applied stress. Understanding the fundamental features leading to such rheological change is crucial to optimize flow conditions or to design flow modifiers for slurry processing. While local changes to the particle environment under applied shear can be related to changes in viscosity, there is a broader need to connect the shear thickening transition to the fundamental organization of particle-interaction forces which lead to long-range organization. In particular, at a high volume fraction of particles, recent evidence indicates frictional forces between contacting particles is of importance. Herein, the network of frictional contact forces is analyzed within simulated two-dimensional shear thickening suspensions. Two topological metrics are studied to characterize the response of the contact force network (CFN) under varying applied shear stress. The metrics, geodesic index and the void parameter, reflect complementary aspects of the CFN: One is the connectedness of the contact network and the second is the distribution of spatial areas devoid of particle-particle contacts. Considered in relation to the variation of the viscosity, the topological metrics show that the network grows homogeneously at large scales but with many local regions devoid of contacts, indicating clearly the role of CFN growth in causing the large change in the rheological response at the shear thickening transition.

The US government currently spends significant resources managing the legacies of the Cold War, including 300 million liters of highly radioactive wastes stored in hundreds of tanks at the Hanford (WA) and Savannah River (SC) sites. The materials in these tanks consist of highly radioactive slurries and sludges at very high pH and salt concentrations. The solid particles primarily consist of aluminum hydroxides and oxyhydroxides (gibbsite and boehmite), although many other materials are present. These form complex aggregates that dramatically affect the rheology of the solutions and, therefore, efforts to recover and treat these wastes. In this paper, we have used a combination of transmission and cryo-transmission electron microscopy, dynamic light scattering, and X-ray and neutron small and ultrasmall-angle scattering to study the aggregation of synthetic nanoboehmite particles at pH 9 (approximately the point of zero charge) and 12, and sodium nitrate and calcium nitrate concentrations up to 1 m. Although the initial particles form individual rhombohedral platelets, once placed in solution they quickly form well-bonded stacks, primary aggregates, up to ∼1500 Å long. These are more prevalent at pH = 12. Addition of calcium nitrate or sodium nitrate has a similar effect as lowering pH, but approximately 100 times less calcium than sodium is needed to observe this effect. These aggregates have fractal dimension between 2.5 and 2.6 that are relatively unaffected by salt concentration for calcium nitrate at high pH. Larger aggregates (>∼4000 Å) are also formed, but their size distributions are discrete rather than continuous. The fractal dimensions of these aggregates are strongly pH-dependent, but only become dependent on solute at high concentrations.

The (meta)stability of low water activity sodium hydroxide/aluminate (Na + OH - /Al(OH) 4 - ) electrolytes dictates kinetics in the Bayer process for aluminum refining and high-level nuclear waste processing. We utilized quasi-elastic neutron scattering (QENS) and proton nuclear magnetic resonance spectroscopy ( 1 H NMR) in extremely concentrated sodium aluminate solutions to investigate the picosecond (ps) to microsecond (ms) timescale motions of H-bearing species (Al(OH) 4 - monomers/clusters, OH - and H 2 O). In the QENS data, in contrast to typical liquids, no short-time translational diffusion was observed at 293 K, but two types of localized motions were found: (i) local backbone tumbling or a formation of large hydrated ion clusters on the order of 40-60 ps; and (ii) much slower, complex, and collective dynamics of the ensemble of H-bearing species on the order of 350-750 ps. Variable temperature, pulsed field gradient, diffusion-ordered 1 H NMR was used to determine the ensemble translational motion along with relaxometry to calculate rotational correlation coefficients. The ensemble rotational correlation times were on the order of 184-300 ps from 1 H NMR, which is consistent with the timescale of the QENS components. Complementary molecular dynamics simulation of NaOH solutions exhibit extensive ion networks potentially responsible for the observed dynamical coupling of water with the motion of large hydrated ion clusters. Understanding these collective motions will aid in predicting the behavior of complex solutions during aluminum production and during nuclear waste processing.

Pulsed field gradient nuclear magnetic resonance (PFG-NMR) measurements were successfully applied to the 27 Al ( I = 5/2) nucleus in concentrated electrolytes to investigate the diffusion of aluminate ions [Al(OH) 4 - ] in simulant high-level nuclear waste (3 M NaOH) between 25 and 85 °C. The temperature-dependent diffusion coefficients obtained from 1 H, 23 Na, and 27 Al PFG-NMR were well fit by a Vogel-Fulcher-Tammann model and a power law equation. Comparison of 27 Al diffusion coefficients of 0.1 M Al(OH) 4 - in ∼3 M MOH (where M = Na + , K + , (CH 3 ) 4 N + ) at room temperature varied in agreement with the expected changes in solution viscosity via Stokes-Einstein relationship, confirming that the dominant Al species at these conditions are Al(OH) 4 - monomers. This 27 Al PFG-NMR study extends an established methodology to a previously unexplored nucleus enabling this experimental technique to be leverage for exploring ion transport, speciation, and solution structure in concentrated electrolytes.

Although oriented aggregation of particles is a widely recognized mechanism of crystal growth, the impact of many fundamental parameters, such as crystallographically distinct interfacial structures, solution composition, and nanoparticle morphology, on the governing mechanisms and assembly kinetics are largely unexplored. Thus, the collective dynamics of systems exhibiting OA has not been predicted. In this context, we investigated the structure and dynamics of boehmite aggregation as a function of solution pH and ionic strength. Cryogenic transmission electron microscopy shows that boehmite nanoplatelets assemble by oriented attachment on (010) planes. The coagulation rate constants obtained from dynamic light scattering during the early stages of aggregation span 7 orders of magnitude and cross both the reaction-limited and diffusion-limited regimes. Combining a simple scaling analysis with calculations for stability ratios and rotational/translational diffusivities of irregular particle shapes, the effects of orientation for irregular-shaped particles on the early stages of aggregation are understood via angular dependencies of van der Waals, electrostatic, and hydrodynamic interactions. Using Monte Carlo simulations, we found that a simple geometric parameter, namely, the contact area between two attaching nanoplatelets, presents a useful tool for correlating nanoparticle morphologies to the emerging larger-scale aggregates, hence explaining the unusually high fractal dimensions measured for boehmite aggregates. Our findings on nanocrystal transport and interactions provide insights toward the predictive understanding of nanoparticle growth, assembly, and aggregation, which will address critical challenges in developing synthesis strategies for nanostructured materials, understanding the evolution of geochemical reservoirs, and addressing many environmental problems.

Aluminum hydroxide (Al(OH) 3 , gibbsite) dissolution and precipitation processes in alkaline environments play a commanding role in aluminum refining and nuclear waste processing, yet mechanistic aspects underlying sluggish kinetics during crystallization have remained obscured due to a lack of in situ probes capable of isolating incipient ion pairs. At a molecular level Al is cycling between tetrahedral ( T d ) coordination in solution to octahedral ( O h ) in the solid. We explored dissolution of Al(OH) 3 that was used to produce variably saturated aluminate (Al(OH) 4 - )-containing solutions under alkaline conditions (pH >13) with in situ 27 Al magic angle spinning (MAS)-nuclear magnetic resonance (NMR) spectroscopy, and interrogated the results with ab initio molecular dynamics (AIMD) simulations complemented with chemical shift calculations. The collective results highlight the overall stability of the solvation structure for T d Al in the Al(OH) 4 - oxyanion as a function of both temperature and Al concentration. The observed chemical shift did not change significantly even when the Al concentration in solution became supersaturated upon cooling and limited precipitation of the octahedral Al(OH) 3 phase occurred. However, subtle changes in Al(OH) 4 - speciation correlated with the dissolution/precipitation reaction were found. AIMD-informed chemical shift calculations indicate that measurable perturbations should begin when the Al(OH) 4 - ···Na + distance is less than 6 Å, increasing dramatically at shorter distances, coinciding with appreciable changes to the electrostatic interaction and reorganization of the Al(OH) 4 - solvation shell. The integrated findings thus suggest that, under conditions incipient to and concurrent with gibbsite crystallization, nominally expected contact ion pairs are insignificant and instead medium-range (4-6 Å) solvent-separated Al(OH) 4 - ···Na + pairs predominate. Moreover, the fact that these medium-range interactions bear directly on resulting gibbsite characteristics was demonstrated by detailed microscopic and X-ray diffraction analysis and by progressive changes in the fwhm of the O h resonance, as measured by in situ NMR. Sluggish gibbsite crystallization may arise from the activation energy associated with disrupting this robust medium-range ion pair interaction.

The aqueous solvation structure of the Na + ion is studied both by molecular dynamics (MD) simulations and K-edge X-ray absorption near edge spectroscopy (XANES). We present a systematic study contrasting the differences in the predicted XANES spectra, using molecular configurations generated from classical and quantum MD. Using a set of suitable order parameters, we elucidate how the spectroscopic features are influenced by the local solvation structure around the Na + ion. Our findings suggest that XANES is sensitive to fluctuations between the first and second solvation shells and can distinguish between classical and quantum interaction representations used to generate molecular ensembles.

The characterization of prenucleation species is essential to understand crystallization mechanisms across many chemical systems and often involves the use of vibrational spectroscopy. Nowhere is this more evident than in the development of "green" aluminum processing technologies, where detailed understanding of the speciation of aluminum and its polynuclear analogues in highly alkaline, low water solutions is elusive. The aluminate anion Al(OH) 4 - predominates in alkaline conditions, yet equilibrium with dimeric species, either μ-oxo Al 2 O(OH) 6 2- or di-μ-hydroxo Al 2 (OH) 8 2- , can be assumed. Using ab initio molecular dynamics with full solvation and the presence of counterions, this work reconciles previous contradictory studies that had concluded only a single species under relevant solution conditions. We reveal that the two dimers are energetically separated by 2 kcal/mol in pure water but that the stability of each can be reversed by ion pairing expected in saturated salt solutions. Simulated Raman and IR spectra for each species (accounting for anharmonicity and the fluctuating solvating environment) provide the first proof that the considered species are "spectroscopic siblings", whose multiple overlapping bands prevent definitive assertions in terms of speciation when compared to the experimental spectra. These observations are likely to hold in higher order aluminate oligomers and as such present a massive challenge toward understanding the crystallization mechanisms relevant to aluminum processing.

Studies on ion pairing at interfaces have been intensified recently because of their importance in many chemical reactive phenomena, such as ion-ion interactions that are affected by interfaces and their influence on kinetic processes. In this study, we performed simulations to examine the thermodynamics and kinetics of small polarizable sodium iodide ions in the bulk and near the water liquid-vapor interface. Using classical transition state theory, we calculated the dissociation rates and corrected them with transmission coefficients obtained from the reactive flux formalism and Grote-Hynes theory. Our results show that in addition to affecting the free energy of ions in solution, the interfacial environments significantly influence the kinetics of ion pairing. The results on the relaxation time obtained using the reactive flux formalism and Grote-Hynes theory present an unequivocal picture that the interface suppresses ion dissociation. The effects of the use of molecular models on the ion interactions as well as the ion-pair configurations at the interface are also quantified and discussed.

Understanding the nature of ionic hydration at a fundamental level has eluded scientists despite intense interest for nearly a century. In particular, the microscopic origins of the asymmetry of ion solvation thermodynamics with respect to the sign of the ionic charge remains a mystery. Here, we determine the response of accurate quantum mechanical water models to strong nanoscale solvation forces arising from excluded volumes and ionic electrostatic fields. This is compared to the predictions of two important limiting classes of classical models of water with fixed point changes, differing in their treatment of "lone pair" electrons. Using the quantum water model as our standard of accuracy, we find that a single fixed classical treatment of lone pair electrons cannot accurately describe solvation of both apolar and cationic solutes, emphasizing the need for a more flexible description of local electronic effects in solvation processes. However, we explicitly show that all water models studied respond to weak long-ranged electrostatic perturbations in a manner that follows macroscopic dielectric continuum models, as would be expected. We emphasize the importance of these findings in the context of realistic ion models, using density functional theory and empirical models, and discuss the implications of our results for quantitatively accurate reduced descriptions of solvation in dielectric media.