Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reacting”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Atomistic mechanisms for catalytic transformations of NO to NH3, N2O, and N2 by Pd

The industrial pollutant NO is a potential threat to the environment and to human health. Thus, selective catalytic reduction of NO into harmless N2, NH3, and/or N2O gas is of great interest. Among many catalysts, metal Pd has been demonstrated to be most efficient for selectivity of reducing NO to N2. However, the reduction mechanism of NO on Pd, especially the route of N−N bond formation, remains unclear, impeding the development of new, improved catalysts. We report here the elementary reaction steps in the reaction pathway of reducing NO to NH3, N2O, and N2, based on density functional theory (DFT)-based quantum mechanics calculations. We show that the formation of N2O proceeds through an Eley-Rideal (E−R) reaction pathway that couples one adsorbed NO* with one non−adsorbed NO from the solvent or gas phase. This reaction requires high NO* surface coverage, leading first to the formation of the trans-(NO)2* intermediate with a low N−N coupling barrier (0.58 eV). Notably, trans-(NO)2* will continue to react with NO in the solvent to form N2O, that has not been reported. With the consumption of NO and the formation of N2O* in the solvent, the Langmuir-Hinshelwood (L-H) mechanism will dominate at this time, and N2O* will be reduced by hydrogenation at a low chemical barrier (0.42 eV) to form N2. In contrast, NH3 is completely formed by the L-H reaction, which has a higher chemical barrier (0.87 eV). Our predicted E-R reaction has not previously been reported, but it explains some existing experimental observations. In addition, we examine how catalyst activity might be improved by doping a single metal atom (M) at the NO* adsorption site to form M/Pd and show its influence on the barrier for forming the N−N bond to provide control over the product distribution.

Physics↗
DOE OSTI2020

Identifying crystal structures and chemical reactions at the interface of stanene on Bi 2 Te 3

Synthesizing monolayers and heterostructures is an enabling approach to extract new physical phenomena from bulk materials. Among the structures amenable to this approach is stanene, which is a monolayer of tin, similar to graphene, and has been predicted to host one-dimensional topological states at its edges. Stanene can be tuned by decorating with different adatoms, which makes it a promising platform on which to engineer topological devices. Here, we deposit Sn on Bi 2 Te 3 and characterize the growth using anomalous synchrotron x-ray scattering and x-ray photoelectron spectroscopy (XPS). X-ray diffraction data reveal the formation of epitaxial Sn-based structures, along with penetration of Sn into the Bi 2 Te 3 , with Sn intercalating between the upper 10 Bi 2 Te 3 quintuple layers. Additionally, XPS data show deposited Sn reacting to form SnTe and Bi at the Bi 2 Te 3 surface. The calculated heat of reaction for Sn and Bi2Te3 is consistent with an exothermic reaction to SnTe and Bi. Using thermodynamic calculations as a guide, we identify several candidate substrates that can stabilize the stanene phase.

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

Importance learning estimator for the site-averaged turnover frequency of a disordered solid catalyst

For disordered catalysts such as atomically dispersed “single-atom” metals on amorphous silica, the active sites inherit different properties from their quenched-disordered local environments. The observed kinetics are site-averages, typically dominated by a small fraction of highly active sites. Standard sampling methods require expensive ab initio calculations at an intractable number of sites to converge on the siteaveraged kinetics. We present a new method that efficiently estimates the site-averaged turnover frequency (TOF). The new estimator uses the same importance learning algorithm [Vandervelden et al., React. Chem. Eng. 5, 77 (2020)] that we previously used to compute the siteaveraged activation energy. We demonstrate the method by computing the site-averaged TOF for a simple disordered lattice model of an amorphous catalyst. The results show that with the importance learning algorithm, the site-averaged TOF and activation energy can now be obtained concurrently with orders of magnitude reduction in required ab initio calculations.

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

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

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

A high temperature cell for investigating interfacial structure on the molecular scale in molten salt/alloy systems

In this work, we describe the design and development of an in situ neutron reflectometry cell for high temperature investigations of structural changes occurring at the interface between inorganic salts, in their molten state up to 800 °C, and corrosion resistant alloys or other surfaces. In the cell, a molten salt is confined by an annular ring of single crystal sapphire constrained between the sample substrate and a sapphire plate using two gold O-rings, enclosing a liquid salt volume of 20 ml, along with a dynamic cell volume to accommodate expansion of the liquid with heating. Further, as a test case for the cell, we report on an in situ neutron reflectometry measurement of the interface between a eutectic salt mixture of MgCl 2 –KCl (32:68 molar ratio) and a single crystal sapphire substrate at 450 °C, resulting in the formation of a 60 Å layer having a scattering length density of 1.72 × 10 –6 Å –2 . While the origin of this layer is uncertain, it is likely to have resulted from the salt reacting with an existing impurity layer on the sapphire substrate.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

The role of H–H interactions and impurities on the structure and energetics of H/Pd(111)

Understanding hydrogen incorporation into palladium requires detailed knowledge of surface and subsurface structure and atomic interactions as surface hydrogen is being embedded. Using density functional theory (DFT), we examine the energies of hydrogen layers of varying coverage adsorbed on Pd(111). Here we find that H–H and H–Pd interactions promote the formation of the well-known ($\sqrt{3}$ x $\sqrt{3}$) phases but also favor an unreported (3 × 3) phase at high H coverages for which we present experimental evidence. We relate the stability of isolated H vacancies of the (3 × 3) phase to the need of H 2 molecules to access bare Pd before they can dissociate. Following higher hydrogen dosage, we observe initial steps of hydride formation, starting with small clusters of subsurface hydrogen. The interaction between H and Pd is complicated by the persistent presence of carbon at the surface. X-ray photoelectron spectroscopy experiments show that trace amounts of carbon, emerging from the Pd bulk despite many surface cleaning cycles, become mobile enough to repopulate the C-depleted surface at temperatures above 200 K. When exposed to hydrogen, these surface carbon atoms react to form benzene, as evidenced by scanning tunneling microscopy observations interpreted with DFT.

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

Synthesis and structure of carbon-doped H 3 S compounds at high pressure

Understanding of recently reported putative close-to-room-temperature superconductivity in C-S-H compounds at 267 GPa demands a reproducible synthesis protocol as well as knowledge of the compounds' structure and composition. Here, we synthesized C-S-H compounds with various carbon compositions at high pressures from elemental carbon C and methane CH 4 , sulfur S, and molecular hydrogen H 2 . Here, we focus on compounds synthesized using methane as these allow a straightforward determination of their structure and composition by combining single-crystal x-ray diffraction and Raman spectroscopy. We applied a two-stage synthesis of [(CH 4 ) x (H 2 S) (1-x) ] 2 H 2 compounds with various compositions by first reacting sulfur and mixed methane-hydrogen fluids and forming CH 4 -doped H 2 S crystals at 0.5-3 GPa and then by growing single crystals of the desired hydrogen-rich compound. Raman spectroscopy applied to this material shows the presence of CH 4 molecules incorporated into the lattice and allows the determination of the CH 4 content, while single-crystal x-ray diffraction results suggest that the methane molecules substitute H 2 S molecules. The structural behavior of these compounds is very similar to the previously investigated methane-free crystals demonstrating a transition from Al 2 Cu type I4/mcm molecular crystal to a modulated molecular structure at 20-30 GPa and back to the same basic I4/mcm structure in an extended modification with greatly modified Raman spectra. This latter phase demonstrates a distortion into a Pnma structure at 132-159 GPa and then transforms into a common Im$\bar{3}$m H 3 S phase at higher pressures; however, no structural anomaly is detected near 220 GPa, where a sharp upturn in T c has been reported.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Numerical study of pure rotational fs/ps CARS coherence beating at high pressure and for multi-species rotation-vibration non-equilibrium thermometry

Coherent anti-Stokes Raman scattering (CARS) has long been the gold standard for non-intrusively measuring gas temperature in reacting flows such as flames and plasmas. Recently, the development of ultrafast hybrid fs/ps CARS has enabled the exploitation of coherence beating between neighboring spectral lines to simultaneously measure rotational and vibrational temperatures from a single pure rotational spectrum. However, the influence of non-Boltzmann vibrational state distributions and limitations due to collisional dephasing at pressures greater than 1 atm remains unclear. In this work, we use spectral simulations to investigate the effects of non-Boltzmann vibrational state distributions and the applicability of coherence beating at pressures up to 10 atm. We show that short probe pulses can be leveraged to quantify non-Boltzmann vibrational state distributions of N2. Furthermore, we demonstrate that fs/ps CARS coherence beating can simultaneously provide sensitive measurements of rotational and vibrational temperatures of both O2 and N2 in air. A sensitivity analysis was conducted to qualitatively explain the accuracy and precision comparisons between probe delays.

Chen, Timothy Y. (ORCID:0000000330269767)↗
DOE OSTI2022

Cesium intercalation of graphene: A 2D protective layer on alkali antimonide photocathode

Alkali antimonide photocathodes have wide applications in free-electron lasers and electron cooling. The short lifetime of alkali antimonide photocathodes necessitates frequent replacement of the photocathodes during a beam operation. Furthermore, exposure to mediocre vacuum causes loss of photocathode quantum efficiency due to the chemical reaction with residual gas molecules. Theoretical analyses have shown that covering an alkali antimonide photocathode with a monolayer graphene or hexagonal boron nitride protects it in a coarse vacuum environment due to the inhibition of chemical reactions with residual gas molecules. Alkali antimonide photocathodes require an ultra-high vacuum environment, and depositing a monolayer 2D material on it poses a serious challenge. In the present work, we have incorporated a novel method known as intercalation, in which alkali atoms pass through the defects of a graphene thin film to create a photocathode material underneath. Initially, Sb was deposited on a Si substrate, and a monolayer graphene was transferred on top of the Sb film. Heat cleaning around 550–600 °C effectively removed the Sb oxides, leaving metallic Sb underneath the graphene layer. Depositing Cs on top of a monolayer graphene enabled the intercalation process. Atomic force microscopy, Raman spectroscopy, x-ray photoelectron spectroscopy, low energy electron microscopy, and x-ray diffraction measurements were performed to evaluate photocathode formation underneath the monolayer graphene. Our analysis shows that Cs penetrated the graphene and reacted with Sb and formed Cs 3 Sb.

2D materials↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

Chapman–Jouguet deflagration criteria and compressibility dynamics of turbulent fast flames for turbulence-induced deflagration-to-detonation transition

This work characterizes the compressibility dynamics in turbulent fast flames for a range of turbulent flame speeds. These turbulent fast flames experience increased effects of compressibility through the formation of strong shocks and may develop a runaway acceleration combined with a pressure buildup that leads to turbulence induced deflagration-to-detonation transition (tDDT). Simultaneous high-speed particle image velocimetry, OH* chemiluminescence, schlieren, and pressure measurements are used to examine the reacting flow field and flame dynamics. We examine flames with turbulent flame speeds ranging from 100 to 600 m/s. At lower turbulent flame speeds, the flame is not able to produce favorable background conditions for deflagration-to-detonation transition (DDT) onset, and thus flame compressibility and turbulence amplification are less dominant, resulting in a weaker acoustic coupling between the flame and compressed region. As the turbulent burning velocities exceed the Chapman–Jouguet deflagration speed, favorable background conditions are produced, as we observe flame-generated shocks and flame-generated turbulence with higher turbulent velocities and larger turbulent scales. At this regime, the flame is categorized to be at the runaway transition regime that leads to tDDT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Room temperature colossal superparamagnetic order in aminoferrocene–graphene molecular magnets

Intensive studies are published for graphene-based molecular magnets due to their remarkable electric, thermal, and mechanical properties. However, to date, most of all produced molecular magnets are ligand based and subject to challenges regarding the stability of the ligand(s). The lack of long-range coupling limits high operating temperature and leads to a short-range magnetic order. Herein, we introduce an aminoferrocene-based graphene system with room temperature superparamagnetic behavior in the long-range magnetic order that exhibits colossal magnetocrystalline anisotropy of 8 × 10 5 and 3 × 10 7 J/m 3 in aminoferrocene and graphene-based aminoferrocene, respectively. These values are comparable to and even two orders of magnitude larger than pure iron metal. Aminoferrocene [C 10 H 11 FeN] + is synthesized by an electrophilic substitution reaction. It was then reacted with graphene oxide that was prepared by the modified Hammers method. The phase structure and functionalization of surface groups were characterized and confirmed by XRD, FT-IR, and Raman spectroscopy. To model the behavior of the aminoferrocene between two sheets of hydroxylated graphene, we have used density functional theory by placing the aminoferrocene molecule between two highly ordered hydroxylated sheets and allowing the structure to relax. Here, the strong bowing of the isolated graphene sheets suggests that the charge transfer and resulting magnetization could be strongly influenced by pressure effects. In contrast to strategies based on ligands surface attachment, our present work that uses interlayer intercalated aminoferrocene opens routes for future molecular magnets as well as the design of qubit arrays and quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The equations of state of statically compressed palladium and rhodium

Here, the pressure–volume equations of state of palladium and rhodium statically compressed in neon are presented. Vinet fits give parameters for palladium: V0=58.678(73) Å 3 , B 0 =189.3(30) GPa, B' 0 =5.473(63), and rhodium: V 0 =55.062(63) Å3, B 0 =241.3(65) GPa, B' 0 =5.34(24). Both metals are observed to react with hydrocarbons under pressure to form hydrides. Existing equations of state are discussed with regard to potential inadvertent hydrogen contamination as a source for discrepancies and anomalous fitted parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Introduction to the special topic on inertial confinement fusion diagnostics

In the first decade of the 21st century, the United States began operating three leading Inertial Confinement Fusion (ICF) and High-Energy Density Physics (HEDP) facilities. The National Ignition Facility (NIF) began its quest to first surpass the Lawson Criterion and then reach thermonuclear ignition and thermonuclear gain. The Omega EP (Extended Performance) laser provides unprecedented energy, reliability, and consistency for petawatt (PW)-class laser experiments. The Omega laser continues to support over 1000 ICF and HEDP experiments each year. The Z pulsed power machine is a Z-pinch plasma confinement device that was refurbished to increase current drive and improve reliability. The extreme states of matter created in these new facilities required advancement of many technologies. The most important among these was the significant improvement in instrumentation capabilities needed to diagnose plasma conditions and resulting nuclear products. The instruments needed a better temporal, spatial, and energy resolution. They also needed to be hardened against EMP and radiation. Typical events in an ICF implosion include compressing the deuterium–tritium (DT) fuel to greater than 60 g/cm 3 and heating the DT to a temperature greater than 10 keV (110 000 000 K). The DT atoms react to produce 14 MeV neutrons from a volume with a radius of ~30 μm. When the fuel ignites, the DT burn prop agates into cold fuel and the neutron-emitting region can reach sizes of over 200 μm. Furthermore, the burn typically lasts less than 150 ps, but the burn duration will decrease as the DT yield increases. Burn durations as short as 90 ps have been recorded and are predicted to decrease still further to 20 ps. X-ray imaging diagnostics are needed to measure the acceleration and velocity of the capsule enclosing the DT, where the implosion velocities can reach over 400 km/s.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Surface chlorination of IrO2(110) by HCl

The ability to controllably chlorinate metal-oxide surfaces can provide opportunities for designing selective oxidation catalysts. In the present study, we investigated the surface chlorination of IrO2(110) by HCl using temperature programmed reaction spectroscopy (TPRS), x-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. We find that exposing IrO2(110) to HCl, followed by heating to 650 K in ultrahigh vacuum, produces nearly equal quantities of on-top and bridging Cl atoms on the surface, Clt and Clbr, where the Clbr atoms replace O-atoms that are removed from the surface by H2O formation. After HCl adsorption at 85 K, only H2O desorbs at low Cl coverages during TPRS, but HCl begins to desorb in increasing yields as the Cl coverage is increased above about 0.5 monolayer (ML). The desorption of Cl2 was not observed under any conditions, in good agreement with the high barrier for this reaction predicted by DFT. A maximum Cl coverage of 1 ML, with nearly equal coverages of Clt and Clbr atoms, could be generated by reacting HCl with IrO2(110) in UHV. Our results suggest that a kinetic competition between recombinative HCl and H2O desorption under the conditions studied limits the saturation Cl coverage to a value less than the 2 ML maximum predicted by thermodynamics. XPS further shows that the partitioning of Cl between the Clt and Clbr states can be altered by subjecting partially chlorinated IrO2(110) to reductive or oxidative treatments, demonstrating that the Cl site population can change dynamically in response to the gas environment. Our results provide insights for understanding the chlorination of IrO2(110) by HCl and can enable future experimental studies to determine how Cl-modification alters the surface chemical reactivity of IrO2(110) and potentially enhances selectivity toward partial oxidation chemistry.

Chemistry↗
DOE OSTI2024

Low-temperature formation of Ti2AlN during post-deposition annealing of reactive multilayer systems

M n+1 AXn-phase Ti 2 AlN thin-films were synthesized using reactive sputtering-based methods involving the deposition of single-layer TiAlN, and Ti/AlN and TiN/TiAl multilayers of various modulation periods at ambient temperature and subsequent annealing at elevated temperatures. Ex situ and in situ x-ray diffraction measurements were used to characterize the Ti 2 AlN formation temperature and phase fraction. During annealing, Ti/AlN multilayers yielded Ti 2 AlN at a significantly lower in situ temperature of 650 °C compared to TiN/TiAl multilayers or single-layer TiAlN (750 °C). The results suggest a reactive multilayer mechanism whereby distinct Ti and AlN layers react readily to release exothermic energy resulting in lower phase transition temperatures compared to TiN and TiAl layers or mixed TiAlN. With a modulation period of 5 nm, however, Ti/AlN multilayers yielded Ti 2 AlN at a higher temperature of 750 °C, indicating a disruption of the reactive multilayer mechanism due to a higher fraction of low-enthalpy interfacial TiAlN within the film.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2025

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗