Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Enthalpy”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Optimal Design and Techno-Economic Analysis of 3D-Printed, Intensified Packings for Absorbers and Strippers in Solvent-Based CO 2 Capture

A potential technology for the CO 2 absorption process is utilizing intensified structured packing with embedded cooling/heating channels for continuous heat exchange, which can overcome limitations of discrete methods, such as discrete intercooling and centralized reboilers, to aid in reducing energy consumption and decreasing costs. This work investigates the modeling of intensified packing (IP) for the stripper tower, extending on previous work for the absorber, which distributes heat internally within the column, improving the thermodynamics for the solvent regeneration process. The model includes submodels for steam turbine extraction to produce steam at various qualities as well as a surrogate model for calculating steam enthalpy. A cost model for a plant-scale absorption capture process was developed, allowing for the design of the plant to be optimized, subject to minimizing capture cost using two different power plant flue gas sources. In this optimization, the placement of IP in both towers is optimized to balance the trade-off between enhanced heat transfer and reduced mass transfer volume. For natural gas combined cycle flue gas, the standard process configuration had a minimum cost of $\$$65.40/tonne CO 2 , and considering IP, the minimum capture cost is reduced to $\$$62.73/tonne, with utilization in the stripper column, which reduces yearly costs by up to $\$$2.67 MM/yr. Cooling the absorber through IP, or intercoolers, was only found to be beneficial at higher capture rates, with IP in both towers having a cost of capture of $\$$68.08/tonne at 99.9% capture, a reduction of $\$$12.64/tonne when using only intercoolers at the same capture rate. When capturing from pulverized-coal power plants, the minimum cost of capture when using IP in both towers is $\$$44.18/tonne (at 97% capture), while the standard configuration with and without intercoolers was $\$$45.69 and $\$$47.22 per tonne, respectively. This results in a reduction in yearly costs of $\$$16.98 MM/yr from the base-case configuration. At this higher CO 2 concentration, cooling in the absorber from the IP becomes extremely beneficial, reducing energy consumption by up to 6%.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Comparison of Heterogeneous Langmuirian Models for Mixed-Gas Adsorption Equilibria

Mixed-gas adsorption remains a very active area of research, where newly proposed adsorption equilibrium models need to be compared to standards in the field. The generalized Langmuir (gL) isotherm, proposed by Hamid et al. to address adsorption of different size molecules on heterogeneous adsorbents without requiring evaluation of the Gibbs isotherm, has thus far been tested mainly against the adsorbed solution theory (AST). The gL isotherm was shown to be thermodynamically consistent and gave either better or similar representation than AST models. Here, the gL isotherm is contrasted with the Dual-Process Langmuir isotherm popularized by Ritter and his co-workers and the Loading Ratio Correlation popularized for simple mixtures, both of which are considered the “standard” models for process simulation. Furthermore, the comparison is based on their theoretical roots, interpretations, and representation of the equilibrium loading and the isosteric enthalpy of adsorption of unary, binary, and ternary mixtures of various gaseous adsorbates on different adsorbents over a range of experimental pressures and temperatures.

42 ENGINEERING↗
DOE OSTI2021

The Role of Water and Hydroxyl Groups in the Structures of Stetindite and Coffinite, MSiO 4 (M = Ce, U)

Orthosilicates adopt the zircon structure types ( I 4 1 /amd ), consisting of isolated SiO 4 tetrahedra joined by A-site metal cations, such as Ce and U. They are of significant interest in the fields of geochemistry, mineralogy, nuclear waste form development, and material science. Stetindite (CeSiO 4 ) and coffinite (USiO 4 ) can be formed under hydrothermal conditions despite both being thermodynamically metastable. Water has been hypothesized to play a significant role in stabilizing and forming these orthosilicate phases, though little experimental evidence exists. To understand the effects of hydration or hydroxylation on these orthosilicates, in situ high-temperature synchrotron and laboratory-based X-ray diffraction was conducted from 25 to ~850 °C. Stetindite maintains its I 4 1 /amd symmetry with increasing temperature but exhibits a discontinuous expansion along the a- axis during heating, presumably due to the removal of water confined in the [001] channels, which shrink against thermal expansion along the a -axis. Furthermore, additional in situ high-temperature Raman and Fourier transform infrared spectroscopy also confirmed the presence of the confined water. Coffinite was also found to expand nonlinearly up to 600 °C and then thermally decompose into a mixture of UO 2 and SiO 2 . A combination of dehydration and dehydroxylation is proposed for explaining the thermal behavior of coffinite synthesized hydrothermally. Additionally, we investigated high-temperature structures of two coffinite-thorite solid solutions, uranothorite (U x Th 1– x SiO 4 ), which displayed complex variations in composition during heating that was attributed to the negative enthalpy of mixing. Lastly, for the first time, the coefficients of thermal expansion of CeSiO 4 , USiO 4 , U 0.46 Th 0.54 SiO 4 , and U 0.9 Th 0.1 SiO 4 were determined to be α V = 14.49 × 10 –6 , 14.29 × 10 –6 , 17.21 × 10 –6 , and 17.23 × 10 –6 °C –1 , respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Contrasting Trivalent Lanthanide and Actinide Complexation by Polyoxometalates via Solution-State NMR

Deciphering the solution chemistry and speciation of actinides is inherently difficult due to radioactivity, rarity, and cost constraints, especially for transplutonium elements. In this context, the development of new chelating platforms for actinides and associated spectroscopic techniques is particularly important. In this study, we investigate a relatively overlooked class of chelators for actinide binding, namely, polyoxometalates (POMs). We provide the first NMR measurements on americium–POM and curium–POM complexes, using one-dimensional (1D) 31 P NMR, variable-temperature NMR, and spin-lattice relaxation time (T 1 ) experiments. The proposed POM–NMR approach allows for the study of trivalent f-elements even when only microgram amounts are available and in phosphate-containing solutions where f-elements are typically insoluble. The solution-state speciation of trivalent americium, curium, plus multiple lanthanide ions (La 3+ , Nd 3+ , Sm 3+ , Eu 3+ , Yb 3+ , and Lu 3+ ), in the presence of the model POM ligand PW 11 O 39 7– was elucidated and revealed the concurrent formation of two stable complexes, [M III (PW 11 O 39 )(H 2 O) x ] 4– and [M III (PW 11 O 39 ) 2 ] 11– . Interconversion reaction constants, reaction enthalpies, and reaction entropies were derived from the NMR data. The NMR results also provide experimental evidence of the weakly paramagnetic nature of the Am 3+ and Cm 3+ ions in solution. Furthermore, the study reveals a previously unnoticed periodicity break along the f-element series with the reversal of T 1 relaxation times of the 1:1 and 1:2 complexes and the preferential formation of the long T 1 species for the early lanthanides versus the short T 1 species for the late lanthanides, americium, and curium. Furthermore, given the broad variety of POM ligands that exist, with many of them containing NMR-active nuclei, the combined POM–NMR approach reported here opens a new avenue to investigate difficult-to-study elements such as heavy actinides and other radionuclides.

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

Metal Atom (Dis)Order and Superconductivity in YCaH n ( n = 8–20) High-Pressure Superhydrides

High-pressure superhydrides have attracted much attention due to their high superconducting critical temperatures (T c s). Our density functional theory (DFT) calculations, focusing on YCaH n (n = 8–20) compositions, found a number of nearly isoenthalpic YCaH 8 phases, differing only in the arrangement of the metal atoms, suggesting the potential stability of metal alloy superhydrides. The computed T c s of the considered YCaH 8 phases were higher than those of the isostructural I 4/ mmm MH 4 parent compounds. DFT enthalpies suggested that YCaH 12 could also be disordered; however, the T c s of the ordered variants spanned a wide range from 105 to 253 K at 200 GPa, showing that alloying could either mildly enhance or drastically reduce T c from that of the $Im\bar{3}m$ MH 6 parents. Finally, for YCaH 18 and YCaH 20 , only a single dynamically stable ordered superhydride was found, which we attribute to the differences in the structures of the most stable MH 9 and MH 10 binary hydride parents.

anions↗
DOE OSTI2020

Dynamics of Cation-Induced Conformational Changes in Nanometer-Sized Uranyl Peroxide Clusters

Conformational changes of the pyrophosphate (Pp)-functionalized uranyl peroxide nanocluster [(UO 2 ) 24 (O 2 ) 24 (P 2 O 7 ) 12 ] 48– ({U 24 Pp 12 }), dissolved as a Li/Na salt, can be induced by the titration of alkali cations into solution. The most symmetric conformer of the molecule has idealized octahedral ( O h ) molecular symmetry. One-dimensional 31 P NMR experiments provide direct evidence that both K + and Rb + ions trigger an O h -to- D 4 h conformational change within {U 24 Pp 12 }. Variable-temperature 31 P NMR experiments conducted on partially titrated {U 24 Pp 12 } systems show an effect on the rates; increased activation enthalpy and entropy for the D 4 h -to- O h transition is observed in the presence of Rb + compared to K + . Two-dimensional, exchange spectroscopy 31 P NMR revealed that magnetization transfer links chemically unique Pp bridges that are present in the D 4 h conformation and that this magnetization transfer occurs via a conformational rearrangement mechanism as the bridges interconvert between two symmetries. The interconversion is triggered by the departure and reentry of K (or Rb) cations out of and into the cavity of the cluster. This rearrangement allows Pp bridges to interconvert without the need to break bonds. Cs ions exhibit unique interactions with {U 24 Pp 12 } clusters and cause only minor changes in the solution 31 P NMR signatures, suggesting that O h symmetry is conserved. Single-crystal X-ray diffraction measurements reveal that the mixed Li/Na/Cs salt adopts D 2 h molecular symmetry, implying that while solvated, this cluster is in equilibrium with a more symmetric form. Thus, these results highlight the unusually flexible nature of the actinide-based {U 24 Pp 12 } and its sensitivity to countercations in solution.

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

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to stop unphysical rearrangements during geometry relaxation. A simple fixed-atom or “coordinate-lock” approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

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

QC-GN 2 oMS 2 : a Graph Neural Net for High Resolution Mass Spectra Prediction

Predicting the mass spectrum of a molecular ion is often accomplished via three generalized approaches: rules-based methods for bond breaking, deep learning, or quantum chemical (QC) modeling. Rules-based approaches are often limited by the conditions for different chemical subspaces and perform poorly under chemical regimes with few defined rules. QC modeling is theoretically robust but requires significant amounts of computational time to produce a spectrum for a given target. Among deep learning techniques, graph neural networks (GNNs) have performed better than previous work with fingerprint-based neural networks in mass spectra prediction. To explore this technique further, we investigate the effects of including quantum chemically derived information as edge features in the GNN to increase predictive accuracy. The models we investigated include categorical bond order, bond force constants derived from extended tight-binding (xTB) quantum chemistry, and acyclic bond dissociation energies. Throughout this work, we evaluated these models against a control GNN with no edge features in the input graphs. Bond dissociation enthalpies yielded the best improvement with a cosine similarity score of 0.462 relative to the baseline model (0.437). In this work we also apply dynamic graph attention which improves performance on benchmark problems and supports the inclusion of edge features. Between implementations, we investigate the nature of the molecular embedding for spectra prediction and discuss the recognition of fragment topographies in distinct chemistries for further development in tandem mass spectrometry prediction.

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

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

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

Anion–Cation–Anion Ion Triplet Characterization by Computation and Photoelectron Spectroscopy

Ion triplets of the chloride salts of two commonly used weakly coordinating cations are reported (i.e., Cl – ·NMe 4 + Cl – (1 – ) and Cl – ·PPh 4 + Cl – (2 – )). Negative ion photoelectron spectra at 20 K afford vertical and adiabatic detachment energies of 5.18 and 5.00 eV (1 – ) and 5.03 and 4.70 eV (2 – ), respectively. These results are well reproduced by coupled cluster calculations with single, double, and perturbative triple excitations (CCSD(T)) whereas M06-2X is systematically too small by ~0.3 eV (i.e., 7 kcal mol –1 ). The structures of both 1 – and 2 – have five or six C–H···Cl – interactions that stabilize these cluster anions by 32 (1 – ) and 27 (2 – ) kcal mol –1 as given by their chloride dissociation enthalpies. In conclusion, these values drop to 7.4 and 3.8 kcal mol –1 in dichloromethane based up conductor-like polarizable continuum model calculations and suggest that X – ·M + X – ion triplets with a weakly coordinating cation maybe the reactive form of salts under some conditions.

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

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗
DOE OSTI2022

IMPROVING THE ACCURACY OF COMPOSITE METHODS: A G4MP2 METHOD WITH G4-LIKE ACCURACY AND IMPLICATIONS FOR MACHINE LEARNING

G4MP2 theory has proven to be a reliable and accurate quantum chemical composite method for the calculation of molecular energies using an approximation based on second-order perturbation theory to lower computational costs compared to G4 theory. However, it has been found to have significantly increased errors when applied to larger organic molecules with 10 or more nonhydrogen atoms. We report here on an investigation of the cause of the failure of G4MP2 theory for such larger molecules. One source of error is found to be the "higher-level correction (HLC)", which is meant to correct for deficiencies in correlation contributions to the calculated energies. This is because the HLC assumes that the contribution is independent of the element and the type of bonding involved, both of which become more important with larger molecules. We address this problem by adding an atom-specific correction, dependent on atom type but not bond type, to the higher-level correction. We find that a G4MP2 method that incorporates this modification of the higher-level correction, referred to as G4MP2A, becomes as accurate as G4 theory (for computing enthalpies of formation) for a test set of molecules with less than 10 nonhydrogen atoms as well as a set with 10-14 such atoms, the set of molecules considered here, with a much lower computational cost. The G4MP2A method is also found to significantly improve ionization potentials and electron affinities. Finally, we implemented the G4MP2A energies in a machine learning method to predict molecular energies.

Dandu, Naveen↗
DOE OSTI2024

Capturing CO 2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO 2 ) n , n = 1,2

Pyrimidine has two in-plane CH(δ+)/N̈(δ–)/CH(δ+) binding sites that are complementary to the (δ–/2δ+/δ–) quadrupole moment of CO 2 . For this study, we recorded broadband microwave spectra over the 7.5–17.5 GHz range for pyrimidine-(CO 2 ) n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14 N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO 2 complex is planar, with CO 2 occupying one of the quadrupolar binding sites, forming a structure in which the CO 2 is stabilized in the plane by interactions with the C–H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO 2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191–195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C 2v symmetry in which the second CO 2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO 2 is –13.7 kJ mol –1 , including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO 2 ) 2 is almost exactly double that size, indicating little interaction between the two CO 2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.

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

Hydrolysis Reaction Pathways of Thorium Oxide Nanoclusters

Density functional theory benchmarked by correlated molecular orbital theory is used to develop a fundamental and predictive understanding of the interaction of thorium oxide nanoclusters with gas phase water to provide insight into nuclear-waste storage, production of thorium nuclear reactor fuels, and reprocessing of spent fuel. The structures of Th n O 2n (n = 3 – 6) clusters and their interactions with water have been studied at the B3LYP, MP2, and CCSD(T) levels. Hydrolysis is initiated by the formation of Lewis acid-base adducts, with relative H 2 O binding energies (physisorption) ranging from −15 kcal/mol to −22 kcal/mol. The initial H 2 O physisorption energy is ca. −21 kcal/mol regardless of the cluster size and is consistent with the experimentally obtained initial adsorption energy on a thorium dioxide surface. The physisorption enthalpies for additional water molecules can be affected by the presence of terminal groups OH groups generated by proton transfer to a Th-O near the site of adsorption. The hydrolysis products (chemisorption) form either bridging or terminal hydroxides. More exothermic hydrolysis steps were predicted for the formation of terminal hydroxides as compared to the formation of bridging hydroxides. Here, the calculated transition state barriers for transfer of protons from bound water complexes to form the chemisorption products are very low. Overall, water readily reacts with thorium oxide clusters preferring hydroxide products over hydrated complexes. First and second order fits were predicted for the combined physisorption and chemisorption energies for the hydrolysis of thorium oxide clusters. Finally, ionization energies and electron affinities were calculated as were HOMO-LUMO gaps to provide additional insights into the properties of the thorium oxide and hydroxide clusters.

Adsorption↗
DOE OSTI2026

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

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

Solvation Energy of Ions in a Stockmayer Fluid

We calculate the solvation energy of monovalent and divalent ions in various liquids with coarse-grained molecular dynamics simulations. Our theory treats the solvent as a Stockmayer fluid, which accounts for the intrinsic dipole moment of molecules and the rotational dynamics of the dipoles. Despite the simplicity of the model, we obtain qualitative agreement between the simulations and experimental data for the free energy and enthalpy of ion solvation, which indicates that the primary contribution to the solvation energy arises mainly from the first and possibly second solvation shells near the ions. Our results suggest that a Stockmayer fluid can serve as a reference model that enables direct comparison between theory and experiment and may be invoked to scale up electrostatic interactions from the atomic to the molecular length scale.

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

Complex Formation and Dissociation Dynamics on Amorphous Silica Surfaces

Benzene complex formation and dissociation dynamics with silanols on the amorphous silica surfaces of nanoporous SiO 2 , from a benzene/carbon tetrachloride solution, were measured by the growth of off-diagonal peaks in the two dimensional infrared (2D IR) chemical exchange spectrum of the isolated Si-OD stretch. The presence of two types of isolated silanols, referred to as type I and II, was revealed, with dissociation time constants of 82 ps and 4.0 ps for type I and II complexes, respectively. The type I silanols are associated with the main IR absorption feature in the Si-OD stretching region, while the type II silanols give rise to a broader shoulder to lower frequency. Polarization selective pump-probe (PSPP) measurements provided the vibrational lifetimes and orientational relaxation rates of the two silanols in the CCl 4 (free) and benzene (complex) environments. The type II silanols constitute roughly 30% of the isolated silanol population and exhibit a substantially faster rate of vibrational relaxation, making the type I dynamics the dominant contribution to the PSPP and 2D IR signals. From the measured dissociation times, the enthalpies of formation for the two surface complexes were obtained; the formation of the type I complex being significantly more exothermic. Furthermore, as the type II site is preferentially removed from the amorphous silica surface with increasing activation temperature, the results provide a reasonable explanation for the increased exothermicity of benzene adsorption on silica with increasing activation temperature in previous calorimetry experiments.

14 SOLAR ENERGY↗
DOE OSTI2021

Role of a Multivalent Ion–Solvent Interaction on Restricted Mg 2+ Diffusion in Dimethoxyethane Electrolytes

The diffusion behavior of Mg 2+ in electrolytes is not as readily accessible as that from Li + or Na + utilizing PFG NMR, due to the low sensitivity, poor resolution, and rapid relaxation encountered when attempting 25 Mg NMR. In MgTFSI 2 /DME solutions, “bound” DME (coordinating to Mg 2+ ) and “free” DME (bulk) are distinguishable from 1 H NMR. With the exchange rates between them obtained from 2D 1 H EXSY NMR, we can extract the self-diffusivities of free DME and bound DME (which are equal to that of Mg 2+ ) before the exchange occurs using PFG diffusion NMR measurements coupled with analytical formulas describing diffusion under two-site exchange. Further, the high activation enthalpy for exhange (65–70 kJ/mol) can be explained by the structural change of bound DME as evidenced by its reduced C–H bond length. Comparison of the diffusion behaviors of Mg 2+ , TFSI – , DME, and Li + reveals a relative restriction to Mg 2+ diffusion that is caused by the long-range interaction between Mg 2+ and solvent molecules, especially those with suppressed motions at high concentrations and low temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗