Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Aqueous”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Predicting Multicomponent Aqueous Phase Adsorption Equilibria of Organic Acids Using the Generalized Brunauer–Emmett–Teller Isotherm Model

Here, to support process development of adsorptive separation of organic acids, this work presents a thermodynamic modeling methodology to predict multicomponent aqueous phase organic acid adsorption equilibria from single-component adsorption isotherms using the generalized Brunauer–Emmett–Teller isotherm model (gBET). With the organic acid fugacities rigorously accounting for the aqueous phase organic acid solution chemistry and solution nonideality, gBET precisely represents the single-component adsorption isotherms and accurately predicts the binary and ternary mixed-acid adsorption equilibria for the ranges of initial pH (∼3–7), acid concentration (100–400 mmol/L), and temperature (298.15–328.15 K) with less than 10% average absolute relative deviation. In addition, gBET with pH-independent parameters provides insights into the underlying adsorption phenomena, including the adsorbate loadings and compositions in the monolayer and subsequent layers under varying initial pH, temperature, concentration, and composition. The gBET model predictions outperform the predictions from the overloading model and the Ideal Dilute Solution Theory.

09 BIOMASS FUELS↗
DOE OSTI2017

Molecular Modeling of the Physical Properties for Aqueous Amine Solution Containing a CO 2 Hydration Catalyst

The effects of an amphiphilic CO 2 hydration catalyst (C3P) on the physical properties of aqueous monoethanolamine (MEA) solutions were studied using molecular simulations and verified experimentally. Adding 2.7–27.7 g/L of C3P in 30 wt % MEA aqueous solution did not significantly affect the solution viscosity, surface tension, or CO 2 diffusivity. These results confirm that the previously reported increase in CO 2 mass transfer by C3P is due to CO 2 hydration catalysis and not due to changes in the physical properties of the MEA solution. Additional simulations indicate that the catalyst molecules tend to aggregate in MEA solution and are preferentially adsorbed at the gas–liquid interface region. Finally, for the catalyst molecules remaining in the bulk solution, the local concentrations of CO 2 and MEA in the area immediately around the catalyst are increased while the local water concentration is decreased, relative to their concentrations in the rest of the bulk MEA solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Bridging the Transuranics with Uranium(IV) Sulfate Aqueous Species and Solid Phases

Isolating isomorphic compounds of tetravalent actinides (i.e., Th IV , U IV , Np IV , and Pu IV ) improve our understanding of the bonding behavior across the series, in addition to their relationship with tetravalent transition metals (Zr and Hf) and lanthanides (Ce). Similarities between these tetravalent metals are particularly illuminated in their hydrolysis and condensation behavior in aqueous systems, leading to polynuclear clusters typified by the hexamer [M IV 6 O 4 (OH) 4 ] 12+ building block. Prior studies have shown the predominance and coexistence of smaller species for Th IV (monomers, dimers, and hexamers) and larger species for U IV , Np IV , and Pu IV (including 38-mers and 70-mers). We show here that aqueous uranium(IV) sulfate also displays behavior similar to that of Th IV (and Zr IV ) in its isolated solidphase and solution speciation. Two single-crystal X-ray structures are described: a dihydroxide-bridged dimer (U 2 ) formulated as U 2 (OH) 2 (SO 4 ) 3 (H 2 O) 4 and a monomer-linked hexamer framework (U-U 6 ) as (U(H 2 O) 3.5 ) 2 U 6 O 4 (OH) 4 (SO 4 ) 10 (H 2 O) 9 . These structures are similar to those previously described for ThIV. Moreover, cocrystallization of monomer and dimer and of dimer and monomer–hexamer phases for both Th IV (prior) and U IV (current) indicates the coexistence of these species in solution. Because it was not possible to effectively study the sulfate-rich solutions via X-ray scattering from which U 2 and U-U 6 crystallized, we provide a parallel solution speciation study in low sulfate conditions, as a function of the pH. Raman spectroscopy, UV–vis spectroscopy, and small-angle X-ray scattering of these show decreasing sulfate binding, increased hydrolysis, increased species size, and increased complexity, with increasing pH. Furthermore, this study describes a bridge across the first half the actinide series, highlighting U IV similarities to Th IV , in addition to the previously known similarities to the transuranic elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Complexation of Lanthanides and Heavy Actinides with Aqueous Sulfur-Donating Ligands

The separation of trivalent lanthanides and actinides is challenging because of their similar sizes and charge densities. S-donating extractants have shown significant selectivity for trivalent actinides over lanthanides, with single-stage americium/ lanthanide separation efficiencies for some thiol-based extractants reported at >99.999%. While such separations could transform the nuclear waste management landscape, these systems are often limited by the hydrolytic and radiolytic stability of the extractant. Progress away from thiol-based systems is limited by the poorly understood and complex interactions of these extractants in organic phases, where molecular aggregation and micelle formation obfuscates assessment of the metal–extractant coordination environment. Because S-donating thioethers are generally more resistant to hydrolysis and oxidation and the aqueous phase coordination chemistry is anticipated to lack complications brought on by micelle formation, we have considered three thioethers, 2,2'-thiodiacetic acid (TDA), (2R,5S)- tetrahydrothiophene-2,5-dicarboxylic acid, and 2,5-thiophenedicarboxylic acid (TPA), as possible trivalent actinide selective reagents. Formation constants, extended X-ray absorption fine structure spectroscopy, and computational studies were completed for thioether complexes with a variety of trivalent lanthanides and actinides including Nd, Eu, Tb, Am, Cm, Bk, and Cf. TPA was found to have moderately higher selectivity for the actinides because of its ability to bind actinides in a different manner than lanthanides, but the utility of TPA is limited by poor water solubility and high rigidity. While significant competition with water for the metal center limits the efficacy of aqueous-based thioethers for separations, the characterization of these solution-phase, S-containing lanthanide and actinide complexes is the most comprehensively available in the literature to date. Here, this is due to the breadth of lanthanides and actinides considered as well as the techniques deployed and serves as a platform for the further development of Scontaining reagents for actinide separations. Additionally, this paper reports on the first bond lengths for Cf and Bk with a neutral S donor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Transient Radiation-Induced Berkelium(III) and Californium(III) Redox Chemistry in Aqueous Solution

Despite the significant impact of radiation-induced redox reactions on the accessibility and lifetimes of actinide oxidation states, fundamental knowledge of aqueous actinide metal ion radiation chemistry is limited, especially for the late actinides. A quantitative understanding of these intrinsic radiation-induced processes is essential for investigating the fundamental properties of these actinides. We present here a picosecond electron pulse reaction kinetics study into the radiation-induced redox chemistry of trivalent berkelium (Bk(III)) and californium (Cf(III)) ions in acidic aqueous solutions at ambient temperature. New, and first of a kind, second-order rate coefficients are reported for the transient radical-induced reduction of Bk(III) and Cf(III) by the hydrated electron (eaq−) and hydrogen atom (H•), demonstrating significant reactivity (up to 1011M-1s-1) indicative of a preference of these metals to adopt divalent states. Additionally, we reportthe first-ever second-order rate coefficients for the transient radical-induced oxidation of these elements by reaction with hydroxyl (•OH) and nitrate (NO3•) radicals, which also exhibited fast reactivity (ca. 108 M-1s-1). Transient Cf(II), Cf(IV), and Bk(IV) absorption spectra are also reported. Overall, the presented data highlight the existence of rich, complex, intrinsic late actinide radiation-induced redox chemistry that has the potential to influence the findings of other areas of actinide science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Complexation of NpO 2 + with Amine-Functionalized Diacetamide Ligands in Aqueous Solution: Thermodynamic, Structural, and Computational Studies

Complexation of Np(V) with three structurally related amine-functionalized diacetamide ligands, including 2,2'-azanediylbis( N, N'-dimethylacetamide) (ABDMA), 2,2'-(methylazanediyl)bis( N, N'-dimethylacetamide) (MABDMA), and 2,2'-(benzylazanediyl)bis( N, N'-dimethylacetamide) (BnABDMA), in aqueous solutions was investigated. The stability constants of two successive complexes, namely, NpO 2 L + and NpO 2 L 2 + , where L stands for the ligands, were determined by absorption spectrophotometry. The results suggest that the stability constants of corresponding Np(V) complexes follow the trend: MABDMA > ABDMA ≈ BnABDMA. The data are discussed in terms of the basicity of the ligands and compared with those for the complexation of Np(V) with an ether oxygen-linked diacetamide ligand. Extended X-ray absorption fine structure data indicate that, similar to the complexation with Nd 3+ and UO 2 2+ , the ligands coordinate to NpO 2 + in a tridentate mode through the amine nitrogen and two oxygen atoms of the amide groups. Furthermore, computational results, in conjunction with spectrophotometric data, verify that the 1:2 complexes (NpO 2 (L) 2 + ) in aqueous solutions are highly symmetric with Np at the inversion center, so that the f-f transition of Np(V) is forbidden and NpO 2 (L) 2 + does not display significant absorption in the near-IR region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Complexation of Light Trivalent Lanthanides with N-(2-Hydroxyethyl)ethylenediamine-N,N',N'-triacetic Acid in Aqueous Solutions: Thermodynamic Analysis and Coordination Modes

N-(2-hydroxyethyl)ethylenediamine-N,N',N'-triacetic acid (HEDTA, denoted as H 3 L) is a strong chelating ligand that is widely used in separation of f-elements as relevant to the nuclear fuel cycle. There is much to be known about the structure and composition of the coordination sphere of the complexes of HEDTA with lanthanides. The complexation of HEDTA with light lanthanides (La 3+ , Nd 3+ and Eu 3+ ) was investigated thermodynamically and structurally in aqueous solutions. Potentiometry and microcalorimetry were performed to acquire the complexation constants (25-70°C) and enthalpies (25°C), respectively, at I = 1.0mol·L -1 NaClO 4 . Coordination modes of the complexes were analyzed by luminescence spectroscopy and nuclear magnetic resonance spectroscopy. The results indicate that there are two successive Ln 3+ /HEDTA complexes, LnL(aq) and Ln 2 (H -1 L) 2 2- (Ln 3+ referred to La 3+ , Nd 3+ and Eu 3+ ; H -1 L 4- referred to the deprotonation of the hydroxyl group) during the titration. The hydroxyl group of HEDTA is coordinated in Ln 3+ /HEDTA complex. The dinuclear Ln 2 (H -1 L) 2 2- complex is present as a carboxyl-bridged centrosymmetric dimer, and two carboxyl groups in bridging positions are coordinated to two adjacent Ln 3+ cations. Complexation of the NdL(aq) is exothermic, while formation of the hydrolytic complex Nd 2 (H -1 L) 2 2- is endothermic. Both NdL aq and Nd 2 (H -1 L) 2 2- complexes are driven by entropic force. These data will help to predict the behavior of lanthanides in the separation process where HEDTA is used as the aqueous complexant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Importance of Nuclear Quantum Effects on Aqueous Electrolyte Transport under Confinement in Ti 3 C 2 MXenes

Protons display a high chemical activity and strongly affect the charge storage capability in confined interlayer spaces of two-dimensional (2D) materials. As such, an accurate representation of proton dynamics under confinement is important for understanding and predicting charge storage dynamics in these materials. While often ignored in atomistic-scale simulations, nuclear quantum effects (NQEs), e.g., tunneling, can be significant under confinement even at room temperature. Here, using the thermostatted ring polymer molecular dynamics implementation of path integral molecular dynamics (PIMD) in conjunction with the ReaxFF force field, density functional tight binding (DFTB), and NequIP neural network potential simulations, we investigate the role of NQEs on proton and water transport in bulk water and aqueous electrolytes under confinement in Ti 3 C 2 MXenes. Although overall NQEs are relatively small, especially in bulk, we find that they can alter both quantitative values and qualitative trends on both proton transport and water self-diffusion under confinement relative to classical MD predictions. Therefore, our results suggest the need for NQEs to be considered to simulate aqueous systems under confinement for both qualitative and quantitative accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Ten-Fold Solvent Kinetic Isotope Effect for the Nonradiative Relaxation of the Aqueous Ferrate(VI) Ion

Hypervalent iron intermediates have been invoked in the catalytic cycles of many metalloproteins, and thus, it is crucial to understand how the coupling between such species and their environment can impact their chemical and physical properties in such contexts. In this work, we take advantage of the solvent kinetic isotope effect (SKIE) to gain insight into the nonradiative deactivation of electronic excited states of the aqueous ferrate-(VI) ion. Here we observe an exceptionally large SKIE of 9.7 for the nanosecond-scale relaxation of the lowest energy triplet ligand field state to the ground state. Proton inventory studies demonstrate that a single solvent O-H bond is coupled to the ion during deactivation, likely due to the sparse vibrational structure of ferrate-(VI). Such a mechanism is consistent with that reported for the deactivation of f-f excited states of aqueous trivalent lanthanides, which exhibit comparably large SKIE values. This phenomenon is ascribed entirely to dissipation of energy into a higher overtone of a solvent acceptor mode, as any impact on the apparent relaxation rate due to a change in solvent viscosity is negligible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Aqueous Photochemistry of Secondary Organic Aerosol of α-Pinene and α-Humulene Oxidized with Ozone, Hydroxyl Radical, and Nitrate Radical

Formation of secondary organic aerosols (SOA) from biogenic volatile organic compounds 13 (BVOC) occurs via O 3 - and OH-initiated reactions during the day and reactions with NO 3 during the 14 night. We explored the effect of these three oxidation conditions on the molecular composition and 15 aqueous photochemistry of model SOA prepared from two common BVOC. A common monoterpene, α- 16 pinene, and sesquiterpene, α-humulene, were used to form SOA in a smog chamber via BVOC + O 3 , 17 BVOC + NO 3 , and BVOC + OH + NOx oxidation. Samples of SOA were collected, extracted in water, 18 and photolyzed in an aqueous solution in order to simulate the photochemical cloud processing of SOA. 19 The extent of change in the molecular level composition of SOA over 4 hours of photolysis (roughly 20 equivalent to 64 hours of photolysis under ambient conditions) was assessed with high-resolution 21 electrospray ionization mass spectrometry. The analysis revealed significant differences in the molecular 22 composition between monoterpene and sesquiterpene SOA formed by the different oxidation pathways. 23 The composition further evolved during photolysis with the most notable change corresponding to the 24 nearly-complete removal of nitrogen-containing organic compounds. Hydrolysis of SOA compounds also 25 occurred in parallel with photolysis. The preferential loss of larger SOA compounds during photolysis 26 and hydrolysis made the SOA compounds more volatile on average. This study suggests that cloud- and 27 fog-processing may under certain conditions lead to a reduction in the SOA loading as opposed to an 28 increase in SOA loading commonly assumed in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Aerosol Fragmentation Driven by Coupling of Acid–Base and Free-Radical Chemistry in the Heterogeneous Oxidation of Aqueous Citric Acid by OH Radicals

A key uncertainty in the heterogeneous oxidation of carboxylic acids by hydroxyl radicals (OH) in aqueous-phase aerosol is how the free-radical reaction pathways might be altered by acid-base chemistry. In particular, if acid-base reactions occur concurrently with acyloxy radical formation and unimolecular decomposition of alkoxy radicals, there is a possibility that differences in reaction pathways impact the partitioning of organic carbon between the gas and aqueous phases. To examine these questions, a kinetic model is developed for the OH-initiated oxidation of citric acid aerosol at high relative humidity. The reaction scheme, containing both free-radical and acid-base elementary reaction steps with physically validated rate coefficients, accurately predicts the experimentally observed molecular composition, particle size, and average elemental composition of the aerosol upon oxidation. The difference between the two reaction channels centers on the reactivity of carboxylic acid groups. Free-radical reactions mainly add functional groups to the carbon skeleton of neutral citric acid, because carboxylic acid moieties deactivate the unimolecular fragmentation of alkoxy radicals. In contrast, the conjugate carboxylate groups originating from acid-base equilibria activate both acyloxy radical formation and carbon-carbon bond scission of alkoxy radicals, leading to the formation of low molecular weight, highly oxidized products such as oxalic and mesoxalic acid. Subsequent hydration of carbonyl groups in the oxidized products increases the aerosol hygroscopicity and accelerates the substantial water uptake and volume growth that accompany oxidation. These results frame the oxidative lifecycle of atmospheric aerosol: it is governed by feedbacks between reactions that first increase the particle oxidation state, then eventually promote water uptake and acid-base chemistry. When coupled to free-radical reactions, acid-base channels lead to formation of low molecular weight gas-phase reaction products and decreasing particle size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Effect of Ions on the Optical Absorption Spectra of Aqueously Solvated Chromophores

In the condensed phase, ions often create heterogeneous local environments around a solute, which may impart chemical reactivity or perturbations to physico-chemical properties. Although the former has been the subject of some study, the latter - particularly as is pertains to optical absorption spectroscopy - is much less understood. As such, in this work, the computed UV-Vis absorption spectrum is examined for the aqueously solvated chromophore anion of green fluorescent protein for different local ion configurations. The strong ability of water to screen the ions from the chromophore results in little change in excitation energy compared to a purely aqueous environment. However, upon forming a contact ion pair with a sodium ion at either of the two electronegative oxygen sites of the chromophore, there is a spectral shift to either higher or lower energies. Surprisingly, our analysis suggests that the cause of the spectral shift is dominated not by the electrostatic presence of the ion, but instead by ion disruption of the hydrogen bond network at the oxygen contact ion pair site.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Resolving Heterogeneous Dynamics of Excess Protons in Aqueous Solution with Rate Theory

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 life-time 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 life-time for the Zundel species is 0.8 ps and 1.3 ps for 2M and 8M HCl, respectively. In this paper, 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 counter-ion in acid solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Concentration-Dependent Solvation Structure and Dynamics of Aqueous Sulfuric Acid Using Multinuclear NMR and DFT

Sulfuric acid is a ubiquitous compound for industrial processes and aqueous sulfate solutions also play a critical role in electrolytes for many prominent battery chemistries. While the thermodynamic literature for it is quite well-developed, comprehensive studies of the solvation structure, and particularly molecular-scale dynamical and transport properties, are less available. This study applies a multinuclear NMR approach to the elucidation of the solvation structure and dynamics over a wide temperature (-10 – 50?C) and concentration (0 – 18 M) range, combining 17O shift, linewidth, and T1 relaxation measurements, 33S shift and linewidth measurements, and 1H PFG-NMR measurements of the proton self-diffusivity. In conjunction, these results indicate a crossover between two regimes of solvation structure and dynamics, occurring above the concentration associated with the deep eutectic point (~ 4.5 M), with the high-concentration regime dominated by strong water-sulfate correlation. This description was borne out in detail by the activation energy trends with increasing concentration derived from the relaxation of both the H2O/H3O+ and H2SO4/HSO4-/SO42- 17O resonances and the 1H self-diffusivity. However, the 17O chemical shift difference between the H2O/H3O+ and H2SO4/HSO4-/SO42- resonances across the entire temperature range is nevertheless strikingly linear. A computational approach coupling MD simulations and DFT NMR shift calculations to reproduce this trend is presented, which will be the subject of further development. This combination of multinuclear, dynamical NMR and computational methods, and the results furnished by this study, will provide a platform for future studies on battery electrolytes where aqueous sulfate chemistry plays a central role in the solution structure.

Bazak, Jonathan D.↗
DOE OSTI2023

Local Molecular Field Theory for Coulomb Interactions in Aqueous Solutions

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Radiolytic Evaluation of 3,4,3-LI(1,2-HOPO) in Aqueous Solutions

We report the octadentate hydroxypyridinone ligand 3,4,3-LI(1,2-HOPO) (abbreviated as HOPO) has been identified as a promising candidate for both chelation and f-element separation technologies, two applications that require optimal performance in radiation environments. However, the radiation robustness of HOPO is currently unknown. Here, we employ a combination of time-resolved (electron pulse) and steady-state (alpha self-radiolysis) irradiation techniques to elucidate the basic chemistry of HOPO and its f-element complexes in aqueous radiation environments. Chemical kinetics were measured for the reaction of HOPO and its Nd(III) ion complex ([Nd III (HOPO)] - ) with key aqueous radiation-induced radical transients (eaq - , H · atom, and · OH and NO 3 · radicals). The reaction of HOPO with eaq - is believed to proceed via reduction of the hydroxypyridinone moiety, while transient adduct spectra indicate that reactions with the H · atom and · OH and NO 3 · radicals proceeded by addition to HOPO's hydroxypyridinone rings, potentially allowing for the generation of an extensive suite of addition products. Complementary steady-state 241 Am(III)-HOPO complex ([ 241 Am III (HOPO)] - ) irradiations showed the gradual release of 241 Am(III) ions with increasing alpha dose up to 100 kGy, although complete ligand destruction was not observed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Furthermore, our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Molecular Dynamics Simulation Study of the Early Stages of Nucleation of Iron Oxyhydroxide Nanoparticles in Aqueous Solutions

Nucleation is a fundamental step in crystal growth. Of environmental and materials relevance are reactions that lead to nucleation of iron oxyhydroxides in aqueous solutions. These reactions are difficult to study experimentally due to their rapid kinetics. Here, we used classical molecular dynamics simulations to investigate nucleation of iron hydroxide/oxyhydroxide nanoparticles in aqueous solutions. Results show that in a solution containing ferric ions and hydroxyl groups, iron–hydroxyl molecular clusters form by merging ferric monomers, dimers, and other oligomers, driven by strong affinity of ferric ions to hydroxyls. When deprotonation reactions are not considered in the simulations, these clusters aggregate to form small iron hydroxide nanocrystals with a six-membered ring-like layered structure allomeric to gibbsite. By comparison, in a solution containing iron chloride and sodium hydroxide, the presence of chlorine drives cluster assembly along a different direction to form long molecular chains (rather than rings) composed of Fe–O octahedra linked by edge sharing. Further, in chlorine-free solutions, when deprotonation reactions are considered, the simulations predict ultimate formation of amorphous iron oxyhydroxide nanoparticles with local atomic structure similar to that of ferrihydrite nanoparticles. Overall, our simulation results reveal that nucleation of iron oxyhydroxide nanoparticles proceeds via a cluster aggregation-based nonclassical pathway.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗