skip to main content

Title: Reactive force fields for aqueous and interfacial magnesium carbonate formation
We develop Mg/C/O/H ReaxFF parameter sets for two environments: an aqueous force field for magnesium ions in solution and an interfacial force field for minerals and mineral–water interfaces. Since magnesium is highly ionic, we choose to fix the magnesium charge and model its interaction with C/O/H through Coulomb, Lennard-Jones, and Buckingham potentials. We parameterize the forcefields against several crystal structures, including brucite, magnesite, magnesia, magnesium hydride, and magnesium carbide, as well as Mg 2+ water binding energies for the aqueous forcefield. Then, we test the forcefield for other magnesium-containing crystals, solvent separated and contact ion-pairs and single-molecule/multilayer water adsorption energies on mineral surfaces. We also apply the forcefield to the forsterite–water and brucite–water interface that contains a bicarbonate ion. We observe that a long-range proton transfer mechanism deprotonates the bicarbonate ion to carbonate at the interface. Free energy calculations show that carbonate can attach to the magnesium surface with an energy barrier of about 0.22 eV, consistent with the free energy required for aqueous Mg–CO 3 ion pairing. Also, the diffusion constant of the hydroxide ions in the water layers formed on the forsterite surface are shown to be anisotropic and heterogeneous. These findings can help explain the experimentally observed more » fast nucleation and growth of magnesite at low temperature at the mineral–water–CO 2 interface in water-poor conditions. « less
Award ID(s):
Publication Date:
Journal Name:
Physical Chemistry Chemical Physics
Page Range or eLocation-ID:
23106 to 23123
Sponsoring Org:
National Science Foundation
More Like this
  1. We exploit gas-phase cluster ion techniques to provide insight into the local interactions underlying divalent metal ion-driven changes in the spectra of carboxylic acids at the air–water interface. This information clarifies the experimental findings that the CO stretching bands of long-chain acids appear at very similar energies when the head group is deprotonated by high subphase pH or exposed to relatively high concentrations of Ca 2+ metal ions. To this end, we report the evolution of the vibrational spectra of size-selected [Ca 2+ ·RCO 2 − ] + ·(H 2 O) n =0 to 12 and RCO 2 − ·(H 2 O) n =0 to 14 cluster ions toward the features observed at the air–water interface. Surprisingly, not only does stepwise hydration of the RCO 2 − anion and the [Ca 2+ ·RCO 2 − ] + contact ion pair yield solvatochromic responses in opposite directions, but in both cases, the responses of the 2 (symmetric and asymmetric stretching) CO bands to hydration are opposite to each other. The result is that both CO bands evolve toward their interfacial asymptotes from opposite directions. Simulations of the [Ca 2+ ·RCO 2 − ] + ·(H 2 O) n clusters indicate thatmore »the metal ion remains directly bound to the head group in a contact ion pair motif as the asymmetric CO stretch converges at the interfacial value by n = 12. This establishes that direct metal complexation or deprotonation can account for the interfacial behavior. We discuss these effects in the context of a model that invokes the water network-dependent local electric field along the C–C bond that connects the head group to the hydrocarbon tail as the key microscopic parameter that is correlated with the observed trends.« less
  2. Abstract
    Site description. This data package consists of data obtained from sampling surface soil (the 0-7.6 cm depth profile) in black mangrove (Avicennia germinans) dominated forest and black needlerush (Juncus roemerianus) saltmarsh along the Gulf of Mexico coastline in peninsular west-central Florida, USA. This location has a subtropical climate with mean daily temperatures ranging from 15.4 °C in January to 27.8 °C in August, and annual precipitation of 1336 mm. Precipitation falls as rain primarily between June and September. Tides are semi-diurnal, with 0.57 m median amplitudes during the year preceding sampling (U.S. NOAA National Ocean Service, Clearwater Beach, Florida, station 8726724). Sea-level rise is 4.0 ± 0.6 mm per year (1973-2020 trend, mean ± 95 % confidence interval, NOAA NOS Clearwater Beach station). The A. germinans mangrove zone is either adjacent to water or fringed on the seaward side by a narrow band of red mangrove (Rhizophora mangle). A near-monoculture of J. roemerianus is often adjacent to and immediately landward of the A. germinans zone. The transition from the mangrove to the J. roemerianus zone is variable in our study area. An abrupt edge between closed-canopy mangrove and J. roemerianus monoculture may extend for up to several hundred metersMore>>
  3. Decavanadate (V 10 O 28 6− or V10) is a paradigmatic member of the polyoxidometalate (POM) family, which has been attracting much attention within both materials/inorganic and biomedical communities due to its unique structural and electrochemical properties. In this work we explored the utility of high-resolution electrospray ionization (ESI) mass spectrometry (MS) and ion exclusion chromatography LC/MS for structural analysis of V10 species in aqueous solutions. While ESI generates abundant molecular ions representing the intact V10 species, their isotopic distributions show significant deviations from the theoretical ones. A combination of high-resolution MS measurements and hydrogen/deuterium exchange allows these deviations to be investigated and interpreted as a result of partial reduction of V10. While the redox processes are known to occur in the ESI interface and influence the oxidation state of redox-active analytes, the LC/MS measurements using ion exclusion chromatography provide unequivocal evidence that the mixed-valence V10 species exist in solution, as extracted ion chromatograms representing V10 molecular ions at different oxidation states exhibit distinct elution profiles. The spontaneous reduction of V10 in solution is seen even in the presence of hydrogen peroxide and has not been previously observed. The susceptibility to reduction of V10 is likely to be shared bymore »other redox active POMs. In addition to the molecular V10 ions, a high-abundance ionic signal for a V 10 O 26 2− anion was displayed in the negative-ion ESI mass spectra. None of the V 10 O 26 cations were detected in ESI MS, and only a low-abundance signal was observed for V 10 O 26 anions with a single negative charge, indicating that the presence of abundant V 10 O 26 2− anions in ESI MS reflects gas-phase instability of V 10 O 28 anions carrying two charges. The gas-phase origin of the V 10 O 26 2− anion was confirmed in tandem MS measurements, where mild collisional activation was applied to V10 molecular ions with an even number of hydrogen atoms (H 4 V 10 O 28 2− ), resulting in a facile loss of H 2 O molecules and giving rise to V 10 O 26 2− as the lowest-mass fragment ion. Water loss was also observed for V 10 O 28 anions carrying an odd number of hydrogen atoms ( e.g. , H 5 V 10 O 28 − ), followed by a less efficient and incomplete removal of an OH˙ radical, giving rise to both HV 10 O 26 − and V 10 O 25 − fragment ions. Importantly, at least one hydrogen atom was required for ion fragmentation in the gas phase, as no further dissociation was observed for any hydrogen-free V10 ionic species. The presented workflow allows a distinction to be readily made between the spectral features revealing the presence of non-canonical POM species in the bulk solution from those that arise due to physical and chemical processes occurring in the ESI interface and/or the gas phase.« less
  4. Information resulting from a comprehensive investigation into the intrinsic strengths of hydrated divalent magnesium clusters is useful for elucidating the role of aqueous solvents on the Mg2+ ion, which can be related to those in bulk aqueous solution. However, the intrinsic Mg–O and intermolecular hydrogen bond interactions of hydrated magnesium ion clusters have yet to be quantitatively measured. In this work, we investigated a set of 17 hydrated divalent magnesium clusters by means of local vibrational mode force constants calculated at the ωB97X-D/6-311++G(d,p) level of theory, where the nature of the ion–solvent and solvent–solvent interactions were interpreted from topological electron density analysis and natural population analysis. We found the intrinsic strength of inner shell Mg–O interactions for [Mg(H2O)n]2+ (n = 1–6) clusters to relate to the electron density at the bond critical point in Mg–O bonds. From the application of a secondary hydration shell to [Mg(H2O)n]2+ (n = 5–6) clusters, stronger Mg–O interactions were observed to correspond to larger instances of charge transfer between the lp(O) orbitals of the inner hydration shell and the unfilled valence shell of Mg. As the charge transfer between water molecules of the first and second solvent shell increased, so did the strength of theirmore »intermolecular hydrogen bonds (HBs). Cumulative local vibrational mode force constants of explicitly solvated Mg2+, having an outer hydration shell, reveal a CN of 5, rather than a CN of 6, to yield slightly more stable configurations in some instances. However, the cumulative local mode stretching force constants of implicitly solvated Mg2+ show the six-coordinated cluster to be the most stable. These results show that such intrinsic bond strength measures for Mg–O and HBs offer an effective way for determining the coordination number of hydrated magnesium ion clusters.« less
  5. All-solid-state batteries (ASSBs) have garnered increasing attention due to the enhanced safety, featuring nonflammable solid electrolytes as well as the potential to achieve high energy density. 1 The advancement of the ASSBs is expected to provide, arguably, the most straightforward path towards practical, high-energy, and rechargeable batteries based on metallic anodes. 1 However, the sluggish ion transmission at the cathode-electrolyte (solid/solid) interface would result in the high resistant at the contact and limit the practical implementation of these all solid-state materials in real world batteries. 2 Several methods were suggested to enhance the kinetic condition of the ion migration between the cathode and the solid electrolyte (SE). 3 A composite strategy that mixes active materials and SEs for the cathode is a general way to decrease the ion transmission barrier at the cathode-electrolyte interface. 3 The active material concentration in the cathode is reduced as much as the SE portion increases by which the energy density of the ASSB is restricted. In addition, the mixing approach generally accompanies lattice mismatches between the cathode active materials and the SE, thus providing only limited improvements, which is imputed by random contacts between the cathode active materials and the SE during the mixingmore »process. Implementing high-pressure for the electrode and electrolyte of ASSB in the assembling process has been verified is a but effective way to boost the ion transmission ability between the cathode active materials and the SE by decreasing the grain boundary impedance. Whereas the short-circuit of the battery would be induced by the mechanical deformation of the electrolyte under high pressure. 4 Herein, we demonstrate a novel way to address the ion transmission problem at the cathode-electrolyte interface in ASSBs. Starting from the cathode configuration, the finite element method (FEM) was employed to evaluate the current concentration and the distribution of the space charge layer at the cathode-electrolyte interface. Hierarchical three-dimensional (HTD) structures are found to have a higher Li + transfer number (t Li+ ), fewer free anions, and the weaker space-charge layer at the cathode-electrolyte interface in the resulting FEM simulation. To take advantage of the HTD structure, stereolithography is adopted as a manufacturing technique and single-crystalline Ni-rich (SCN) materials are selected as the active materials. Next, the manufactured HTD cathode is sintered at 600 °C in an N 2 atmosphere for the carbonization of the resin, which induces sufficient electronic conductivity for the cathode. Then, the gel-like Li 1.4 Al 0.4 Ti 1.6 (PO 4 ) 3 (LATP) precursor is synthesized and filled into the voids of the HTD structure cathode sufficiently. And the filled HTD structure cathodes are sintered at 900 °C to achieve the crystallization of the LATP gel. Scanning transmission electron microscopy (STEM) is used to unveil the morphology of the cathode-electrolyte interface between the sintered HTD cathode and the in-situ generated electrolyte (LATP). A transient phase has been found generated at the interface and matched with both lattices of the SCN and the SE, accelerating the transmission of the Li-ions, which is further verified by density functional theory calculations. In addition, Electron Energy Loss Spectroscopy demonstrates the preserved interface between HTD cathode and SEs. Atomic force microscopy is employed to measure the potential image of the cross-sectional interface by the peak force tapping mode. The average potential of modified samples is lower than the sample that mix SCN and SEs simply in the 2D planar structure, which confirms a weakened space charge layer by the enhanced contact capability as well as the ion transmission ability. To see if the demonstrated method is universally applicable, LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NCM811) is selected as the cathode active material and manufactured in the same way as the SCN. The HTD cathode based on NCM811 exhibits higher electrochemical performance compared with the reference sample based on the 2D planar mixing-type cathode. We believe such a demonstrated universal strategy provides a new guideline to engineer the cathode/electrolyte interface by revolutionizing electrode structures that can be applicable to all-solid-state batteries. Figure 1. Schematic of comparing of traditional 2D planar cathode and HTD cathode in ASSB Tikekar, M. D. , et al. , Nature Energy (2016) 1 (9), 16114 Banerjee, A. , et al. , Chem Rev (2020) 120 (14), 6878 Chen, R. , et al. , Chem Rev (2020) 120 (14), 6820 Cheng, X. , et al. , Advanced Energy Materials (2018) 8 (7) Figure 1« less