Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibrium constant”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

dynamics of organic-mineral interactions at the metal oxide-solution interface as studied via binding energetics (Final report)

This project focused on addressing longstanding fundamental and experimental uncertainties on how dissolved organic substances (DOS) interact with metal oxide surface under environmentally relevant conditions. By leveraging a custom-built real-time, in-tandem flow adsorption microcalorimetry-UV-Vis/fluorescence spectroscopy platform, we characterized the binding energetics, kinetics and mechanistic pathways driving DOS-metal oxide interactions at temporal resolution on the order of 1-5 seconds. We studied a diverse suite of model organic compounds/substances – including monocarboxylates (e.g. acetate and benzoate), di-carboxylates (oxalate and succinate), amino acids, amino-based nanparticles and natural organic matter – interacting at the mineral-water interface of structurally- and/or chemically distinct metal oxides (including SiO2, boehmite, ferrihydrite, and γ-Al2O3). Our results indicated that DOS-metal oxide interactions are governed by multi-step reaction pathways, often switching between distinct, resolvable enthalpy- and entropy-driven non-electrostatic or electrostatic configurations. To quantify these interactions, we developed and implemented an analytical workflow that integrates peak deconvolution and Monte-Carlo based error propagation to determine site-specific thermodynamic and kinetic parameters for individual binding/debinding events. In addition to resolving apparent first-order rate constants of each event, we were able to quantify associated apparent equilibrium constants as well as free energy, enthalpy and entropy contribution to the activation and subsequent progression of the binding/debinding process across compounds, compound class and metal oxide surfaces. The kinetic-thermodynamic data produced in this study captured how the interplay between oxide surface reactivity and DOS molecular structure jointly drives binding-debinding dynamics. Notably, that at pH below PZC of the oxide surface, neutral species were heavily involved in monocarboxylate binding, while anionic species drove dicarboxylate binding. Also, that among amino acids 1) positional isomers show distinctive binding characteristics to each other while enantiomers show no significant differences in binding characteristics, 2) molecules that bind via outer-sphere complexation show a larger entropic shift between binding and debinding with no impact on oxide surface while 3) inner-sphere interactions increased anion exchange capacity of the oxide surface. The new insights and data from this work has great potential for improving predictive modeling of carbon dynamics and specifically organic-mineral interactions in environmental and industrial systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Microwave-driven heterogeneous catalysis for activation of dinitrogen to ammonia under atmospheric pressure

This paper presents an innovative approach to producing energy-dense, carbon–neutral liquid ammonia as a means for carrying energy. This approach synergistically integrates microwave reaction chemistry with novel heterogeneous catalysis that decouples dinitrogen activation from high-temperature and high-pressure reactions, altering reaction pathways and increasing ammonia formation rate. Results presented here demonstrate that ammonia synthesis can be conducted at 280 ? and ambient pressure to achieve ~1 mmol ammonia/g cat/h over supported ruthenium catalyst systems utilizing microwave irradiation. It is further shown that adding promoter ions such as potassium, cerium, and barium significantly improves the ammonia production rate over undoped ruthenium-based catalysts. This effect could be attributed to enhanced dielectric loss processes that lead to stronger microwave absorption by the catalyst. Measurement of the equilibrium constant under microwave conditions showed a higher ammonia yield than under thermal equilibrium conditions for both the iron- and ruthenium-based catalysts. Finally, this study also illustrates the advantages of using a variable-frequency microwave reactor for ambient-pressure ammonia synthesis. Mechanistically, investigators believed that the oscillating electric fields of the radiation can couple with adsorbed nitrogen on the surface and accelerate its dissociation. Since dinitrogen dissociation on the surface is rate limiting, this effectively accelerates the reaction.

Hu, Jianli (John)↗
DOE OSTI2022

Autoprocessing and oxyanion loop reorganization upon GC373 and nirmatrelvir binding of monomeric SARS-CoV-2 main protease catalytic domain

The monomeric catalytic domain (residues 1–199) of SARS-CoV-2 main protease (MPro 1-199 ) fused to 25 amino acids of its flanking nsp4 region mediates its autoprocessing at the nsp4-MPro 1-199 junction. We report the catalytic activity and the dissociation constants of MPro 1-199 and its analogs with the covalent inhibitors GC373 and nirmatrelvir (NMV), and the estimated monomer-dimer equilibrium constants of these complexes. Mass spectrometry indicates the presence of the accumulated adduct of NMV bound to MPro WT and MPro 1-199 and not of GC373. A room temperature crystal structure reveals a native-like fold of the catalytic domain with an unwound oxyanion loop (E state). In contrast, the structure of a covalent complex of the catalytic domain-GC373 or NMV shows an oxyanion loop conformation (E* state) resembling the full-length mature dimer. These results suggest that the E-E* equilibrium modulates autoprocessing of the main protease when converting from a monomeric polyprotein precursor to the mature dimer.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Tumor‐pH‐value responsive non‐peripheral substituted phthalocyanines: Synthesis, investigation of photophysical and photochemical properties

Cancer is one of the diseases with the highest mortality rate worldwide. Although PDT has recently produced encouraging outcomes, there are still many areas that need to be improved. The first of these is the negative consequences faced by patients treated with the PDT method when exposed to sunlight. For this reason, a new PDT method has been developed in recent years, and it is aimed at using photosensitizer molecules that can be active in acidic conditions. Since the pH values of tumor tissues are more acidic than normal tissues, preparing molecules that act effectively in acidic conditions will allow for more effective results in treating cancer with PDT. In this context, within the scope of this study, 3‐(4‐propionylphenoxy)phthalonitrile ( 1 ) and its non‐peripheral tetra‐substituted phthalocyanine derivatives [( 2 ), ( 3 ), and ( 4 )] were prepared. With these phthalocyanine derivatives, the novel compounds ( 5 ), ( 6 ), and ( 7 ) were synthesized for the first time. The aggregation tendencies of newly synthesized phthalocyanines ( 5–7 ) were investigated in solvent media. The effects of pH changes upon UV–Vis and fluorescence spectra were performed. The electronic and emission spectra of synthesized phthalocyanine derivatives are highly sensitive to pH changes. Formation constant (Log K) values of mono‐ and di‐protonated phthalocyanine forms were calculated by the Henderson–Hasselback equation. The mono‐ and di‐protonated species' equilibrium constants (logK 1 and logK 2 ) were calculated as ~5.0. This value may be promising for pH‐sensitizing photosensitizers. Also, the photophysical and photochemical properties of synthesized metallophthalocyanine derivatives ( 2 ) and ( 5 ) were studied at different pH values. The singlet oxygen quantum yield of ( 2 ) and ( 5 ) was calculated to be 0.78 and 0.81 in DMSO, respectively. When pH = 6.4, that is, tumor‐pH‐values, this value for ( 5 ) has increased to 0.92. The newly synthesized phthalocyanines are suitable photosensitizers for PDT applications, especially with high singlet oxygen quantum yield at pH 6.4.

Genc Bilgicli, Hayriye↗
DOE OSTI2021

Electrochemical Modeling of Iodide Oxidation in Metal-Halide Molten Salts

The oxidation of iodide in NaI-AlBr 3 , NaI-AlCl 3 , and NaI-GaCl 3 molten salts was analyzed using simulation software to extract relevant kinetic parameters. The experimental oxidation potentials were ordered AlCl 3 < AlBr 3 < GaCl 3 , with higher oxidation potentials correlating with softer Lewis acidity of the metal halide. An iodide oxidation and metal halide speciation model was developed and simulated to fit the electrochemical response, enabling determination of electrochemical charge transfer parameters and chemical equilibrium constants. NaI-AlBr 3 displayed the fastest electron transfer rates yet showed the lowest current densities. All salts revealed smaller than expected current densities, explained by equilibrium between various species, where some are not electrochemically active at the studied potentials. These equilibrium reactions are due to the various metal halide species, controlling the reactant concentration of iodide and the resultant current. We hypothesize the electrochemically active iodide species, present as a metal halide monomer (MX 3 I − ), is decreased dramatically from the expected concentration, sequestered as a more stable metal halide dimer species (M 2 X 6 I − ) with a higher oxidation potential. Traditional Tafel analysis of the experimental data supports the validity of the simulations. These results increase understanding of iodide oxidation in low-temperature Lewis acidic molten salts and inform task-specific molten salt design.

Lee, Rose Y. (ORCID:0000000226587220)↗
DOE OSTI2023

On the Thermodynamic Condition for Adsorption Azeotropes

Adsorption azeotropy is a common phenomenon in mixed-gas adsorption equilibria. Since the first literature report of binary adsorption azeotropes in 1933, numerous studies investigated the thermodynamic conditions for adsorption azeotropes but failed to reach definitive conclusions. Based on the generalized Langmuir isotherm model for multicomponent adsorption equilibria which takes into account the vacant site as part of the adsorbed phase, this study presents the thermodynamic condition for adsorption azeotropes as derived from the generalized Langmuir isotherm to be the equality of the ratios of adsorbed phase activity coefficient γi and adsorption equilibrium constant $K$$^{o}_{i}$ for the two adsorbates in the binary 1–2 adsorption system, i.e., γ 1 /$K$$^{o}_{1}$ = γ 2 /$K$$^{o}_{2}$. This adsorption azeotropic condition is analogous to the vapor–liquid equilibrium azeotropic condition, i.e., the equality of the products of liquid phase activity coefficient and saturation vapor pressure Pisat for the two components, i.e., γ 1 $P$$^{sat}_{1}$ = γ 2 $P$$^{sat}_{2}$. We validated the thermodynamic condition for adsorption azeotropes with 14 azeotrope-forming adsorption systems in this study. As a result, we investigated the effects of the pressure, temperature, and adsorbed phase nonideality on azeotrope formation.

42 ENGINEERING↗
DOE OSTI2022

Re-evaluation of rate constants for the reaction N 2 H 4 (+ M) $^→_←$ NH 2 + NH 2 (+ M)

Rate constants for the dissociation/recombination reaction N 2 H 4 (+ M) $^→_←$ NH 2 + NH 2 (+ M) are determined by a combination of quantum-chemical calculations and statistical unimolecular rate theory. Between 1100 and 2500 K, limiting low-pressure rate constants for hydrazine dissociation in the bath gas Ar of k 0 = [Ar] 6.1 × 10 20 (T/1000 K) –7.3 exp (-34490 K/T) cm 3 mol –1 s –1 , limiting high-pressure rate constants of k ∞ = 7.6 × 10 16 (T/1000 K) –1.0 exp(-33600 K/T) s –1 , and center broadening factors of the falloff curves (between 1100 and 1600 K in Ar) of F cent = 0.71 exp(- T/1460 K) + 0.29 exp(- T/21 K) + exp(-13400 K/T) were calculated. Using equilibrium constants K C = 1.7 × 10 3 (T/1000 K) –1.5 exp(-33,460 K/T) mol cm –3 , between 300 and 600 K limiting low-pressure rate constants for the reverse recombination of NH 2 radicals in the bath gas N 2 of k rec,0 = [N 2 ] 4.4 × 10 20 (T/400 K) –6.9 exp (-1630 K/T) cm 6 mol –2 s –1 , limiting high-pressure rate constants of k rec,∞ = 4.4 × 10 13 (T/1000 K) 0.44 exp(-140 K/T) cm 3 mol -1 s –1 , and center broadening factors (between 300 and 600 K in N 2 ) of F cent = exp(- T/1130 K) + exp (-10340 K/T) were obtained. A comparison with experimental results for hydrazine dissociation from the literature suggests incomplete falloff extrapolations toward k ∞ and experimental problems in the determination of k 0 at temperatures above about 1600 K. Finally, implications of the present re-evaluated rate constants for the modeling of high temperature ammonia oxidation kinetics are discussed, showing an only small influence of their precise values on the overall properties of the process.

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

Machine Learning in Environmental Chemistry: Application to Surface Complexation Modeling

Environmental chemistry – or biogeochemistry – is the scientific discipline typically invoked when examining and quantifying groundwater or surface water contamination, and nutrient cycling in the environment. Over the last three decades, there have been significant advances in mechanistic model development to describe and predict these complex biogeochemical processes. In particular, surface complexation models (SCMs) have been developed to describe the rock/soil surface reactions of metals and radionuclides, and their partitioning between various mobile species in the aqueous phase or immobile species sorbed on solid surfaces. Often represented by a simplified linear isotherm constant – Kd – in reactive transport models, these reactions play a critical role in many environmental science applications; particularly in contamination risk assessments and nuclear waste disposal performance assessments. In the past several decades, efforts by various institutions across the world have focused on developing SCMs based on datasets from laboratory measurements, including the identification of key parameters such as equilibrium constants.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Kinetics and Energetics of Electron Transfer to Dimer Radical Cations

In this work, spectra of the dimer cations naphthalene (Nap 2 •+ ) and ethene (Ethene 2 •+ ) were measured in liquid dichloromethane (DCM). The spectra peak at very different energies, 1.2 and 3.3 eV. In DCM dimerization stabilizes Nap 2 •+ by Δ$G_\text{d}°$(Nap 2 •+ ) = -218 meV relative to the monomer Nap •+ as determined from the dimerization equilibrium constant. Both dimers can transfer a positive charge to hole acceptor molecules, but for both the rate constants rise more gradually with reaction energetics than do many charge transfer reactions previously studied. A striking observation finds that the rate constant for hole transfer from the Nap 2 •+ dimer to phenanthrene is smaller by two decades than that from biphenyl•+ monomer to Nap, although both reactions have the same –Δ$\textit{G}°$=0.05 eV. A plausible interpretation for these observations is the presence of an energy of reorganization, λ(M 2 ), for the dimer that involves movement apart of the two partners in the dimer. While the dimerization equilibrium cannot be measured for Ethene 2 •+ the charge transfer data imply that both Δ$G_\text{d}°$(Ethene 2 •+ ) and λ(Ethene 2 •+ ), are considerably larger, perhaps by factors of 2-4 than for Nap 2 •+ .

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

Reaction path model of the formation of abiotic immiscible hydrocarbon fluids in subducted carbonated serpentinites, Lanzo Massif (Western Italian Alps)

Fluids generated from subducted slabs participate in the cycling of deep carbon in the crust and upper mantle. In these fluids, aqueous carbon species vary in oxidation state between +IV and -IV depending on whether the fluids are oxidizing or reducing, respectively. Most studies of subduction-zone fluids have focused on oxidized carbon species. However, recent studies of natural samples have demonstrated the occurrence of deep, reducing fluids, generated in both subducted oceanic upper mantle and crustal rocks. CH 4 -H 2 -rich fluid inclusions in subducted carbonated serpentinites have demonstrated the existence of abiotic, immiscible, hydrocarbon fluids at upper mantle conditions. To investigate the formation of such immiscible hydrocarbon fluids during the evolution of subducted carbonated serpentinites, we used equilibrium constants from the Deep Earth Water model to carry out predictive chemical mass transfer modeling to simulate the alteration reactions. A novel feature of the models was the inclusion of an immiscible hydrocarbon fluid containing six components (CH 4,f , C 2 H 6,f , C 3 H 8,f , isoC 4 H 10,f , CO 2,f , H 2,f ). This feature enabled prediction of the formation of a separate immiscible fluid in equilibrium with aqueous species and minerals. We developed a predictive reaction path model of invasive H 2,f reacting with carbonated serpentinites and interstitial aqueous fluids for comparison with the natural samples from the Lanzo Massif, western Italian Alps. Over a range of temperatures and pressures, immiscible hydrocarbon fluids formed in association with altered mineral assemblages. Reaction progress caused the transformation of carbonated serpentinites and the formation of clinopyroxene, brucite, graphite, and hydrocarbon fluids, along with changes of pH, logfO 2 , and aqueous species. CH 4,f was the most abundant hydrocarbon species in all the models. The overall results at 2.0 GPa and 400 to 450 °C were consistent with the natural samples from the Lanzo Massif. Interestingly, large amounts of H 2 O formed due to oxidation of H 2 . More hydrocarbons and H 2 O formed in models with lower fluid/rock mass ratios or with more reactant H 2 . Models at different pressure and temperature conditions showed similar results with some variation in the relative stabilities of aragonite, graphite and olivine solid solution, and associated differences in mineral sequences, hydrocarbon fluids, values of aqueous species, and the final log fO 2 and pH. As a result, our models strongly support the laboratory and field evidence that reduction of carbonated serpentinites by infiltrating H 2 fluids can cause the formation of immiscible, abiotic hydrocarbon fluids in subduction zones.

58 GEOSCIENCES↗
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
DOE OSTI2021

Solvated Electron in Acetonitrile: Radiation Yield, Absorption Spectrum, and Equilibrium between Cavity- and Solvent-Localized States

The equilibrium between a solvent cavity-localized electron, ecav-, and a dimeric solvent anion, (CH 3 CN)$_{2}^{•-}$, which are the two lowest energy states of the solvated electron in acetonitrile, has been investigated by pulse radiolysis at 233-353 K. The enthalpy and entropy for the e$_{cav}^{-}$ to (CH 3 CN)$_{2}^{•-}$ conversion amount to -11.2 ± 0.3 kcal/mol and -39.3 ± 1.2 cal/(mol K), corresponding to a 0.44 ± 0.35 equilibrium constant at 25 °C. The radiation yield of the solvated electron has been quantified using a Co(II) macrocycle that scavenges electrons with a 1.55 × 10 11 M -1 s -1 rate constant. The apparent yield increases without saturation over the attainable scavenger concentration range, reaching 2.8 per 100 eV; this value represents the lower limit for the acetonitrile ionization yield in pulse radiolysis. The apparent molar absorption coefficient of (20.8 ± 1.5) × 10 3 M -1 cm -1 at 1450 nm and 20 °C for the solvated electron and individual vis-near-infrared (NIR) absorption spectra of e$_{cav}^{-}$ and (CH 3 CN)$_{2}^{•-}$ are derived from the data. Variances with previous reports are thoroughly discussed. All together, these results resolve several controversies concerning the solvated electron properties in acetonitrile and furnish requisite data for quantitative pulse radiolysis investigations in this commonly used solvent.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Methanediol from cloud-processed formaldehyde is only a minor source of atmospheric formic acid

Atmospheric formic acid is severely underpredicted by models. A recent study proposed that this discrepancy can be resolved by abundant formic acid production from the reaction ( 1 ) between hydroxyl radical and methanediol derived from in-cloud formaldehyde processing and provided a chamber-experiment-derived rate constant, k 1 = 7.5 × 10 −12 cm 3 s −1 . High-level accuracy coupled cluster calculations in combination with E,J -resolved two-dimensional master equation analyses yield k 1 = (2.4 ± 0.5) × 10 −12 cm 3 s −1 for relevant atmospheric conditions ( T = 260–310 K and P = 0–1 atm). We attribute this significant discrepancy to HCOOH formation from other molecules in the chamber experiments. More importantly, we show that reversible aqueous processes result indirectly in the equilibration on a 10 min. time scale of the gas-phase reaction HCHO + H 2 O ⇌ HOCH 2 OH (2) with a HOCH 2 OH to HCHO ratio of only ca . 2%. Although HOCH 2 OH outgassing upon cloud evaporation typically increases this ratio by a factor of 1.5–5, as determined by numerical simulations, its in-cloud reprocessing is shown using a global model to strongly limit the gas-phase sink and the resulting production of formic acid. Based on the combined findings in this work, we derive a range of 1.2–8.5 Tg/y for the global HCOOH production from cloud-derived HOCH 2 OH reacting with OH. The best estimate, 3.3 Tg/y, is about 30 times less than recently reported. The theoretical equilibrium constant K eq (2) determined in this work also allows us to estimate the Henry’s law constant of methanediol (8.1 × 10 5 M atm −1 at 280 K).

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

Enhanced prediction of Cs removal by CST from Hanford tank waste with K accountability

Abstract The treatment of Hanford tank waste is one of the most challenging environmental cleanup activities to date. To expedite the processing of liquid waste stored in underground tanks in Washington State it is necessary to remove the significant dose contributor, 137 Cs. Crystalline silicotitanate ion exchanger is currently used to remove 137 Cs from the aqueous phase of Hanford tank wastes in preparation for vitrification at the Waste Treatment and Immobilization Plant (WTP). Improving the understanding of potassium impacts on ion exchange behavior of Cs will help in the operation of a critical component of one of the most complex treatment processes in the world today. Optimization of this process can result in significant cost savings and less waste production. Toward this effort, a series of batch contact tests varied in potassium concentration were conducted to look at the impact of potassium concentration on Cs distribution. Experimental distribution ratios ( K d ) were compared to the distribution ratios predicted using the ZAM model. A significant underprediction of Cs capacity in the presence of potassium was seen with the existing model. A revision of the equilibrium constants was determined and provided a statistically better fit for determining the Cs K d values in tank waste matrices.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Structural Manipulation of a Zirconocene-Based Porous Coordination Cage Using External and Host–Guest Stimuli

Porous coordination cages (PCCs) constructed from a zirconocene-based cluster have drawn great attention for their high stability and structural tunability. These features have enabled these cages to be utilized for a wide array of applications. Here this work demonstrates the structural flexibility and stimulus-responsive features of two interconvertible zirconocene-based PCCs. The subjects, PCC-20c (capsular) and PCC-20t (tetrahedral), are constructed from the same cluster and ligand but exhibit different shapes and properties. In the synthetic study of the cages, it is found that higher temperatures of synthesis yield PCC-20c in several different crystalline phases while lower temperatures of synthesis yield pure phase PCC-20t. Thermodynamic studies on the interconversion of the two cages show that the conversion from PCC-20t into PCC-20c is an entropically driven process and can be largely affected by the solvent environment. Additionally, the process of PCC-20t/c conversion can be mediated by the counter anions and light irradiation as evidenced by the change in equilibrium constants under these influences. The current studies are conducted to decipher how the PCC-20t/c phase can be controllably synthesized, interconverted, and their properties can be modified under the use of external stimuli: heat, light, and guest molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Zavitsas’ hydration model for electrolytes only stable in the presence of another: NaAl(OH) 4 in aqueous NaOH solution

Many aqueous electrolytes are only stable in the presence of another electrolyte, such as electrolytes that are only soluble in strong acid or base. An example is sodium aluminate [NaAl(OH) 4 ], which is only stable in aqueous NaOH. This complicates developing thermodynamic parameters because it is difficult to separate the contributions of individual electrolytes in multicomponent solutions to measured bulk thermodynamic properties. The present study develops a method to determine the liquid phase parameters from solubility data for Zavitsas’ Hydration model, a model that incorporates hydration parameters into the activities of dissolved species. Furthermore, this study uses gibbsite [Al(OH) 3 ] solubility data in aqueous NaOH solution to develop Zavitsas' model parameters for NaAl(OH) 4 simultaneously with the equilibrium constants. A good fit of the solubility data was found, showing that Zavitsas’ model is effective for this system.

Bayer process↗