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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
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↗
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 OSTI2022

Quantifying Energetic and Entropic Pathways in Molecular Systems

When examining dynamics occurring at nonzero temperatures, both energy and entropy must be taken into account to describe activated barrier crossing events. Furthermore, good reaction coordinates need to be constructed to describe different metastable states and the transition mechanisms between them. Here we use a physics-based machine learning method called state predictive information bottleneck (SPIB) to find nonlinear reaction coordinates for three systems of varying complexity. SPIB is able to correctly predict an entropic bottleneck for an analytical flat-energy double-well system and identify the entropy- and energy-dominated pathways for an analytical four-well system. Lastly, for a simulation of benzoic acid permeation through a lipid bilayer, SPIB is able to discover the entropic and energetic barriers to the permeation process. Given these results, we thus establish that SPIB is a reasonable and robust method for finding the important entropy, energy, and enthalpy barriers in physical systems, which can then be used to enhance the understanding and sampling of different activated mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2020

Energetics and Structure of Nickel Atoms and Nanoparticles on MgO(100)

The growth morphology and interfacial energetics of vapor deposited Ni on the MgO(100) surface at 300 and 100 K have been studied using single crystal adsorption calorimetry (SCAC), He + low-energy ion scattering spectroscopy (LEIS), X-ray photoelectron spectroscopy (XPS), and low-energy electron diffraction (LEED). At 300 K, the Ni atoms grow as three-dimensional nanoparticles with a saturation number density of 5 × 10 16 particles/m 2 . The differential heat of adsorption at 300 K increases rapidly with coverage, from 276 (initially) to 311 kJ/mol by 0.4 ML. Thereafter, it slowly increases asymptotically to the sublimation enthalpy of bulk Ni (430 kJ/ml) by 9 ML. The Ni 2p 3/2 XPS peak binding energy at 300 K is initially (i.e., at 0.16 ML) 1.4 eV higher than that for bulk Ni(solid), but it decreases to that value at high coverage. The Ni atoms form a metastable phase at 300 K when in nanoparticles with diameter <2.5 nm, and the adhesion energy of such Ni nanoparticles to MgO(100) was found to be 3.05 J/m 2 . At 100 K, the Ni atoms form single adatoms and then 0.17 nm thick 2D islands at low coverage with fewer Ni-Ni bonds compared to the Ni nanoparticles formed at 300 K. Thus, the initial heat (i.e., for the first ~0.03 ML) is 148 kJ/mol at 100 K, 128 kJ/mol lower than at 300 K, and remains lower for the 2D islands. With increasing coverage at 100 K, the tiny 2D Ni islands grow in size to cover nearly the entire surface before thickening. Lastly, the XPS Ni 2p 3/2 peak binding energy for 0.21 ML Ni on MgO(100) at 100 K is 2.2 eV higher than that for bulk Ni(solid), suggesting charge transfer from Ni to MgO(100) and formation of Ni 2+ at very low coverage.

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

Effect of Background Electrolyte Composition on the Interfacial Formation of Th(IV) Nanoparticles on the Muscovite (001) Basal Plane

Understanding the impact of actinide nanoparticle (NP) formation is important to assess radionuclide mobility in the environment. We combined surface X-ray diffraction (SXRD) and in situ AFM to investigate the previously reported unusual electrolyte effects on Th uptake on mica. At low [Th] (0.1 mM), interfacial structures show a broad Th electron density (~50 Å). A linear decrease of Th uptake with decreasing hydration enthalpy of the electrolyte cation (Li+, K+, NH4+, and Cs+) indicates a competitive effect between Th and the electrolyte cation. Na+ is a clear outlier from this trend. In situ AFM imaging confirms the results. Particles show a vertical size of ~1–2 nm and larger lateral dimensions of ~10–20 nm, which is typical for particles formed at interfaces (heterogeneous nucleation). At high [Th] = 1 and 3 mM, all investigated electrolytes (ACl, A = Li+, Na+, K+) show similar Th uptake, indicating a much smaller impact of electrolyte composition. The interfacial structures are dominated by a high Th loading at a distinct distance (~6.5 Å) from the surface. Therefore, the main retention mechanism at high [Th] is suggested to be the sorption of Th NPs aggregated from Th oligomers present in solution (homogeneous nucleation).

Neumann, J↗
DOE OSTI2021

Prediction of Transition-State Scaling Relationships and Universal Transition-State Vibrational and Entropic Correlations for Dehydrogenations

Linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) or transition-state scaling (TSS) relations aid with the prediction of electronic energies. However, temperature effects and pre-exponentials are often taken as constants across metal surfaces or a homologous series. Vibrational scaling relationships (VSRs) offer a way to determine such parameters. Transition-state VSRs (TSVSRs) between local minima and transition states of AH X (A = C, N, O) surface diffusions correlate with BEP relations and broaden to thermochemical property scaling. Using density functional theory, we extend TSVSRs to AHX dehydrogenation reactions on transition-metal surfaces, relating vibrational modes of local minima to transition states. We first predict the slopes of the TSS relations by incorporating bond angles using the Slater–Koster structure factors and hybridization through crystal orbital overlap population analysis and energy overlap integrals between adsorbates and metal surfaces. Additionally, we uncover universal thermochemical property scaling, enabling the estimation of entropies and temperature corrections to enthalpies across a homologous series. Here, we demonstrate both significant vibrational corrections in reactions with low intrinsic electronic barriers and considerable variation in the pre-exponential of a simple dehydrogenation reaction across metals and AH X adsorbates.

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

Elucidating the Influence of Metal Surface Composition on Organic Adsorbate Binding Using Active Particle Dynamics

The adsorption strengths of organic compounds on metal surfaces are sensitive to the metal composition, and they play a central role in many catalytic reactions, helping to control the coverage of the reactant and altering the overall reaction rate. While adsorption energies are straightforward to measure and calculate in vacuum and gas-phase environments, adsorption energetics can be dramatically altered by the presence of solvent in liquid-phase reactions. However, the effects of metal composition on binding strengths in a liquid environment are less well understood, primarily due to the difficulty of accurate in situ measurements of organic binding on metal surfaces in the liquid phase. Here, we utilize the motion of active particles in water to probe the adsorption energies of an organic adsorbate (furfural) on a range of metal surfaces (pure Pd, pure Pt, and four PdAu alloy compositions) to elucidate the effect of metal composition. Janus particles with catalytic caps of particular metal compositions all exhibited active motion resulting from consumption of H 2 O 2 ; adsorbate binding was inferred through the decrease in the velocity of active motion and was modeled by a Langmuir adsorption isotherm. The measured adsorption affinities were used to extract the adsorption enthalpy of furfural on the different metals. The Pd surface was found to bind furfural more strongly than the Pt surface by some 10 kJ/mol. Furthermore, the adsorption of furfural on the alloys was found to increase monotonically in magnitude with Pd content. As a result, the data reported herein aid the development of accurate understanding of organic adsorption in the presence of solvent and the role of the metal surface in tuning adsorption strengths to optimize catalytic processes in the liquid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Polyethylene in Dead-End Silica Nanopores: Forces and Mobility from Non-Equilibrium Statistical Mechanics and Exchange Spectroscopy Nuclear Magnetic Resonance

Billions of tons of plastic have been produced, and only a small fraction of this has been recycled. Tennakoon et al. [Nature Catalysis 3, 893 (2020)] developed a catalyst that repeatedly cleaves C 10 –C 30 hydrocarbons from the end of a polyethylene chain. The reaction occurs at a Pt nanoparticle at the base of a cylindrical silica mesopore with a diameter of 2 nm and a length of 110 nm. Portions of the polymer situated inside the pore can be differentiated from those outside using 13 C nuclear magnetic resonance (NMR), allowing the dynamics and extent of polymer threading to be monitored using two-dimensional (2D) exchange spectroscopy NMR. Here we construct a Fokker–Planck equation for the polymer dynamics by assuming a reptation diffusivity and a graduated adsorption free energy that depends linearly on the depth of polymer penetration in the pore. The solutions allow us to predict the intensities of the 2D NMR resonances as a function of time. We use the solutions to extract a polymer diffusivity at each temperature and estimate the per-segment desorption free energy, enthalpy, and entropy. Random and systematic errors are examined to test key assumptions in the theory and interpretation of the experiments.

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

Ethylene Hydrogenation Molecular Mechanism on MoC y Nanoparticles

Ethylene hydrogenation catalyzed by MoC y nanoparticles has been studied by means of density functional theory methods and several models. These include MetCar (Mo 8 C 12 ), Nanocube (Mo 14 C 13 ), and Mo 12 C 12 nanoparticles as representatives of experimental MoC y nanostructures. The effect of hydrogen coverage has been studied in detail by considering low-, intermediate-, and high-hydrogen regimes. The calculated enthalpy and energy barriers show that ethylene hydrogenation is feasible on the MetCar, Mo 12 C 12 , and Nanocube but at low, medium, and high hydrogen coverages, respectively. An additional step, related to the H* migration from a Mo to a C site in the nanoparticle, has been found to be the key to establishing the best hydrogenation system. In most cases, the reactions are exothermic, featuring low hydrogenation energy barriers, especially for the Nanocube at high hydrogen coverage. In addition, the calculated adsorption Gibbs free energy shows that, for this system, the C 2 H 4 adsorption is feasible in the 300–400 K temperature range and pressures from 10 –10 to 2 atm. For the hydrogenation steps, calculated transition state theory rates show that the overall process is limited by the first hydrogenation step (C 2 H 4 → C 2 H 5 ) at temperatures of 330–400 K. However, at the lower temperatures of 300–320 K, the reaction rates are comparable for the two steps. The present results indicate that the Mo 14 C 13 Nanocube models of MoC y nanoparticles exhibit appropriate thermodynamic and kinetic features to catalyze ethylene hydrogenation at a high-hydrogen-coverage regime. The present findings provide a basis for understanding the chemistry of active MoC y catalysts, suggest appropriate working conditions for the reaction to proceed, and provide a basis for future experimental studies.

03 NATURAL GAS↗