Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electrons”

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 595 records · Page 33

Initial Investigation of Coupled Electron-Photon Transport with SCEPTRE and RAPTURE (Rev.1.1)

The SCEPTRE1 code models coupled photon and electron transport in material media using discrete ordinates methods. The code can use different solvers in 1D and 3D in parallel, promising users the ability to run problems that are significant both in physics and computational scale. One capability that is perhaps unique at this time in the NNSA complex is SCEPTRE’s coupled photon– electron transport physics in a discrete ordinates approach. For example, incident x–rays can generate a significant electron population, whose transport can be simulated within the same code. In this context, fully–coupled means that photons can produce electrons and electrons can produce photons; uncoupled or no–coupling means photons do not produce electrons. Partial–coupling, an intermediate option in which photons can produce electrons but electrons do not produce photons, is not considered in the present work. The MCNP code can treat this physics with a Monte Carlo approach. MCNP should be viewed as complementary to SCEPTRE in both computational method (Monte Carlo versus discrete ordinates) as well as data libraries. SCEPTRE and MCNP should provide important checks on each other, especially from complementary modeling frameworks of a S n (discrete ordinates) mesh based method versus Monte Carlo particle based method. In order to learn more about SCEPTRE, and how it might inform engineering work at the Laboratory, we used it to carry out a set of initial studies. This was a series of 1D simulations of photon beams hitting Aluminum slabs. Our studies mainly used SCEPTRE directly, and made little use of the closely associated RAPTURE code (more on this below). The goals of these studies were first, to familiarize ourselves with operating the code, and second, to begin to assess the importance of coupled electron–photon transport. To support these studies, we began to develop a set of scripts to initialize inputs, run tools, process outputs, and plot results. An initial version of these scripts has been committed to the Lightning workflow repository. This memo relates the results of those studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Energetic Nonthermal Electrons within the Above-the-looptop Regions in Solar Flares: Acceleration, Feedback, and Quasiperiodic Pulsations

Solar flares are among the most dramatic events in the solar system, releasing substantial magnetic energy and accelerating a large number of electrons to high energies. Notably, in certain events, the above-the-looptop region may contain a significant population of nonthermal electrons, both in number and energy. For the first time, we adopt a novel numerical method that combines magnetohydrodynamics with energetic particles incorporating feedback from nonthermal electrons to investigate electron acceleration and transport in solar flares. We find that a large fraction of energetic electrons are accelerated via the current sheet and termination shock regions. Most energetic electrons are concentrated in the above-the-looptop region, carrying a sizable amount of the released energy. We observe that greater feedback of nonthermal electrons leads to steeper energy spectra. The energy density of the nonthermal electrons oscillates due to the periodic impact of magnetic islands into the above-the-looptop region, which may help explain the observed quasiperiodic pulsations. Our simulations provide new insights into the origin of nonthermal electrons and associated emissions in the above-the-looptop region.

79 ASTRONOMY AND ASTROPHYSICS↗

The Electron Thermal Conductivity of Pu and Zr Substituted Gamma-Uranium

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (k) and fissile density; however, the effects of alloying elements on k remain unclear. Here, the electron thermal conductivity (k_e) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the k_e of ?-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate k_e; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the method that applies the electron relaxation time of ?-U is best, while if the elements are dissimilar, a mixed method that mixes several parameters associated with k_e from each element in the alloy is best. The introduction of all alloying elements decreases k_e; however, in binary compounds, Pu and Zr have different effects. Pu generally flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases the k_e more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel with without sacrificing k. We also note that these k_e calculation methods can be applied to non-fuel alloys that require k_e predictions, such as cladding, heat exchanger, and structural materials.

36 MATERIALS SCIENCE↗

Simulations, Modeling and Data Analysis of Parity Violating Electron Scattering Experiments

In the Standard Model (SM) of nuclear and particle physics, parity violation is incorporated through the representation of the weak interaction as a chiral gauge interaction. Only the left-handed components of particles and right-handed components of antiparticles participate in weak interactions in the Standard Model. This implies that parity is asymmetric for the weak interaction. Parity violating electron scattering (PVES) experiments are designed to probe the physics parameters related to the SM, with the possibility to discover physics beyond the SM (BSM) by measuring the parity violating asymmetry ¿¿¿ of longitudinally polarized electrons scattered off unpolarized targets with high precision. This dissertation will be focused on two PVES experiments, the next 208Pb Lead Radius Experiment (PREX-II), and the Measurement of a Lepton-Lepton Electroweak Reaction (MOLLER) experiment, as well as in some small sections, the Calcium Radius Experiment (CREX) and P2 experiment which are also PVES experiments). PREX-II and CREX experiments, performed in Hall A at the Thomas Jefferson National Accelerator Facility (Jefferson Lab), measured ¿¿¿ in the elastic scattering of longitudinally polarized electrons from 208Pb and 48Ca targets to provide a precise model independent determination of the neutron skin thickness of 208Pb and 48Ca nuclei, respectively. The MOLLER experiment, proposed to start in 2027 and also to be performed in Hall A at Jefferson Lab, is to measure ¿¿¿ of longitudinally polarized electrons scattered off unpolarized electrons (Møller scattering) to determine the weak charge of electrons ¿¿¿ and the weak mixing angle ¿¿ with high precision. As for the P2 experiment, which will be performed at the upcoming MESA accelerator in Mainz Germany, it is to measure the weak charge of proton ¿¿¿ using ¿¿¿ in the elastic electron-proton scattering of polarized electrons off unpolarized protons. The final results from the PREX-II experiment are presented as ¿¿¿=550±16 (¿¿¿¿)±8 (¿¿¿¿) parts-per-billion (ppb). Combining the PREX-I and PREX-II results, the neutron skin thickness from PREX experiments is determined as ¿¿-¿¿=0.283±0.071 ¿¿ in 208Pb. This thesis lists the software and computational contribution of the author to these PVES experiments, including writing scripts and software to help with the PREX-II/CREX experiments, analyzing data to provide useful information and systematic uncertainty for the PREX-II experiment, modeling and simulations for the MOLLER, and providing an alternative design of an electronic equipment for the P2 experiment.

Chen, Yufan↗

Simulations, Modeling and Data Analysis of Parity Violating Electron Scattering Experiments

In the Standard Model (SM) of nuclear and particle physics, parity violation is incorporated through the representation of the weak interaction as a chiral gauge interaction. Only the left-handed components of particles and right-handed components of antiparticles participate in weak interactions in the Standard Model. This implies that parity is asymmetric for the weak interaction. Parity violating electron scattering (PVES) experiments are designed to probe the physics parameters related to the SM, with the possibility to discover physics beyond the SM (BSM) by measuring the parity violating asymmetry ¿¿¿ of longitudinally polarized electrons scattered off unpolarized targets with high precision. This dissertation will be focused on two PVES experiments, the next 208Pb Lead Radius Experiment (PREX-II), and the Measurement of a Lepton-Lepton Electroweak Reaction (MOLLER) experiment, as well as in some small sections, the Calcium Radius Experiment (CREX) and P2 experiment which are also PVES experiments). PREX-II and CREX experiments, performed in Hall A at the Thomas Jefferson National Accelerator Facility (Jefferson Lab), measured ¿¿¿ in the elastic scattering of longitudinally polarized electrons from 208Pb and 48Ca targets to provide a precise model independent determination of the neutron skin thickness of 208Pb and 48Ca nuclei, respectively. The MOLLER experiment, proposed to start in 2027 and also to be performed in Hall A at Jefferson Lab, is to measure ¿¿¿ of longitudinally polarized electrons scattered off unpolarized electrons (Møller scattering) to determine the weak charge of electrons ¿¿¿ and the weak mixing angle ¿¿ with high precision. As for the P2 experiment, which will be performed at the upcoming MESA accelerator in Mainz Germany, it is to measure the weak charge of proton ¿¿¿ using ¿¿¿ in the elastic electron-proton scattering of polarized electrons off unpolarized protons. The final results from the PREX-II experiment are presented as ¿¿¿=550±16 (¿¿¿¿)±8 (¿¿¿¿) parts-per-billion (ppb). Combining the PREX-I and PREX-II results, the neutron skin thickness from PREX experiments is determined as ¿¿-¿¿=0.283±0.071 ¿¿ in 208Pb. This thesis lists the software and computational contribution of the author to these PVES experiments, including writing scripts and software to help with the PREX-II/CREX experiments, analyzing data to provide useful information and systematic uncertainty for the PREX-II experiment, modeling and simulations for the MOLLER, and providing an alternative design of an electronic equipment for the P2 experiment.

Chen, Yufan↗

Soluble Iron Enhances Extracellular Electron Uptake by Shewanella oneidensis MR‐1

Abstract Extracellular electron transfer (EET) is a process that microorganisms use to reduce or oxidize external insoluble electron acceptors or donors. Much of our mechanistic understanding of this process is derived from studies of transmembrane cytochrome complexes and extracellular redox shuttles that mediate outward EET to anodes and external electron acceptors. In contrast, there are knowledge gaps concerning the reverse process of inward EET from external electron donors to cells. Here, we describe a role for soluble iron (exogenous FeCl 2 ) in enhancing EET from cathodes to the model EET bacterium Shewanella oneidensis MR‐1, with fumarate serving as the intracellular electron acceptor. This iron concentration‐dependent electron uptake was eradicated upon addition of an iron chelator and occurred only in the presence of fumarate reductase, confirming an electron pathway from cathodes to this periplasmic enzyme. Moreover, S. oneidensis mutants lacking specific outer membrane and periplasmic cytochromes exhibited significantly decreased current levels relative to wild‐type. These results indicate that soluble iron can function as an electron carrier to the EET machinery of S. oneidensis .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗

Emergence of competing electronic states from non-integer nuclear charges

Understanding many-electron phenomena with competing near-degenerate electronic states is of fundamental importance to chemistry and condensed matter physics. One of the most significant challenges for exploring such many-electron phenomena is the necessity for large system sizes in order to realize competing states, far beyond those practical for first-principles methods. Here, we show how allowing non-integer nuclear charges expands the space of computationally tractable electron systems that host competing electronic states. The emergence of competing electronic states from non-integer nuclear charges is exemplified in the simple 2-electron H 2 molecule and used to examine the microscopic structure of doped quasi-1D cuprate chains, showing how non-integer nuclear charges can open a window for first-principles calculations of difficult many-electron phenomena.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Design and performance of a magnetic bottle electron spectrometer for high-energy photoelectron spectroscopy

We describe the design and performance of a magnetic bottle electron spectrometer (MBES) for high-energy electron spectroscopy. Our design features a 2 m long electron drift tube and electrostatic retardation lens, achieving sub-electronvolt (eV) electron kinetic energy resolution for high energy (several hundred eV) electrons with a close to 4π collection solid angle. A segmented anode electron detector enables the simultaneous collection of photoelectron spectra in high resolution and high collection efficiency modes. This versatile instrument is installed at the time-resolved molecular and optical sciences instrument at the Linac Coherent Light Source x-ray free-electron laser (XFEL). In this paper, we demonstrate its high resolution, collection efficiency, and spatial selectivity in measurements where it is coupled to an XFEL source. These combined characteristics are designed to enable high-resolution time-resolved measurements using x-ray photoelectron, absorption, and Auger–Meitner spectroscopy. We also describe the pervasive artifact in MBES time-of-flight spectra that arises from a periodic modulation in electron collection efficiency and present a robust analysis procedure for its removal.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Ion temperature clamping in Wendelstein 7-X electron cyclotron heated plasmas

The neoclassical transport optimization of the Wendelstein 7-X stellarator has not resulted in the predicted high energy confinement of gas fueled electron-cyclotron-resonance-heated (ECRH) plasmas as modelled in (Turkin et al 2011 Phys. Plasmas 18 022505) due to high levels of turbulent heat transport observed in the experiments. The electron-turbulent-heat transport appears non-stiff and is of the electron temperature gradient (ETG)/ion temperature gradient (ITG) type (Weir et al 2021 Nucl. Fusion 61 056001). As a result, the electron temperature T e can be varied freely from 1 keV–10 keV within the range of P ECRH = 1–7 MW, with electron density n e values from 0.1–1.5 × 10 20 m –3 . By contrast, in combination with the broad electron-to-ion energy-exchange heating profile in ECRH plasmas, ion-turbulent-heat transport leads to clamping of the central ion temperature at T i ~ 1.5 keV ± 0.2 keV. In a dedicated ECRH power scan at a constant density of $\langle n_{e} \rangle$ = 7 × 10 19 m –3 , an apparent 'negative ion temperature profile stiffness' was found in the central plasma for (r/a < 0.5), in which the normalized gradient ∇T i /T i decreases with increasing ion heat flux. The experiment was conducted in helium, which has a higher radiative density limit compared to hydrogen, allowing a broader power scan. This 'negative stiffness' is due to a strong exacerbation of turbulent transport with an increasing ratio of T e /T i in this electron-heated plasma. This finding is consistent with electrostatic microinstabilities, such as ITG-driven turbulence. Theoretical calculations made by both linear and nonlinear gyro-kinetic simulations performed by the GENE code in the W7-X three-dimensional geometry show a strong enhancement of turbulence with an increasing ratio of T e /T i . The exacerbation of turbulence with increasing T e /T i is also found in tokamaks and inherently enhances ion heat transport in electron-heated plasmas. This finding strongly affects the prospects of future high-performance gas-fueled ECRH scenarios in W7-X and imposes a requirement for turbulence-suppression techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extracting correlation effects from momentum-resolved electron energy loss spectroscopy: Synergistic origin of the dispersion kink in Bi 2.1 Sr 1.9 CaCu 2 O 8 + x

We employ momentum-resolved electron energy loss spectroscopy (M-EELS) on Bi 2.1 Sr 1.9 CaCu 2 O 8 + x to resolve the issue of the kink feature in the electron dispersion widely observed in the cuprates. To this end, we utilize the GW approximation to relate the density response function measured in M-EELS to the self-energy, isolating contributions from phonons, electrons, and the momentum dependence of the effective interaction to the decay rates. The phononic contributions, present in the M-EELS spectra due to electron-phonon coupling, lead to kink features in the corresponding single-particle spectra at energies between 40 and 80 meV, independent of the doping level. We find that a repulsive interaction constant in momentum space is able to yield the kink attributed to phonons in ARPES. Hence, our analysis of the M-EELS spectra points to local repulsive interactions as a factor that enhances the spectroscopic signatures of electron-phonon coupling in cuprates. We conclude that the strength of the kink feature in cuprates is determined by the combined action of electron-phonon coupling and electron-electron interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Enhanced Monte Carlo Simulations for Electron Energy Loss Mitigation in Real-Space Nanoimaging of Thick Biological Samples and Microchips

High-resolution imaging using Transmission Electron Microscopy (TEM) is essential for applications such as grain boundary analysis, microchip defect characterization, and biological imaging. However, TEM images are often compromised by electron energy spread and other factors. In TEM mode, where the objective and projector lenses are positioned downstream of the sample, electron–sample interactions cause energy loss, which adversely impacts image quality and resolution. This study introduces a simulation tool to estimate the electron energy loss spectrum (EELS) as a function of sample thickness, covering electron beam energies from 300 keV to 3 MeV. Leveraging recent advances in MeV-TEM/STEM technology, which includes a state-of-the-art electron source with 2-picometer emittance, an energy spread of 3 × 10 -5 , and optimized beam characteristics, we aim to minimize energy spread. By integrating EELS capabilities into the BNL Monte Carlo (MC) simulation code for thicker samples, we evaluate electron beam parameters to mitigate energy spread resulting from electron–sample interactions. Based on our simulations, we propose an experimental procedure for quantitively distinguishing between elastic and inelastic scattering. The findings will guide the selection of optimal beam settings, thereby enhancing resolution for nanoimaging of thick biological samples and microchips.

36 MATERIALS SCIENCE↗

Mechanistic implications of excited high-spin states, spin–spin coupling, and differential [2Fe–2S] + cluster temperature relaxations in the electron-bifurcating NfnABC from Thermococcus sibiricus

Electron bifurcation (EB) is a mechanism of biological energy transduction in which multiple oxidation–reduction (redox) reactions are thermodynamically coupled within a single enzyme, enabling the enzyme to harness the excess free energy from an exergonic process to drive an endergonic process. Because of this unprecedented chemistry, there is interest to translate EB principles to artificial and bioengineered systems, but a hurdle is that knowledge pertaining to the fundamental design principles of EB enzymes remains scarce. Here, we investigated the fundamental physical and electronic properties of electron transfer sites in a spectroscopically uncharacterized member of the BfuABC family of EB enzymes, the NADH-dependent reduced-ferredoxin:NADP + oxidoreductase from Thermococcus sibiricus (Tsi NfnABC). Cryo-EM structures of Tsi NfnABC previously demonstrated that it contains twelve redox cofactors: two flavins (one FAD and one FMN), eight [4Fe–4S] clusters, and two [2Fe–2S] clusters. The FMN, one [4Fe–4S] cluster, and one [2Fe–2S] cluster comprise the bifurcating active site termed the electron-bifurcating flavobicluster (BF-FBC), which is found in all BfuABC family members. By using electron paramagnetic resonance spectroscopy, we identified spectral signatures originating from interactions between the FMN radical and [4Fe–4S] + cluster in the BF-FBC and observed temperature dependent behavior of the BF-FBC's [2Fe–2S] + cluster indicative of moderately slow spin–lattice relaxation. Additionally, we uncovered numerous spectral features corresponding to half-integer, S > ½ spin states of [4Fe–4S] + clusters, including one attributable to the consequences of lysine-ligation of a [4Fe–4S] cluster unique to NfnABC. We contextualize these findings to electron transfer theory and NfnABC's structure. Our insights further the understanding of how enzymes are designed to exert control over electron transfer to conduct thermodynamically challenging reactions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Near-Monochromatic Tuneable Cryogenic Niobium Electron Field Emitter

We report creating, manipulating, and detecting coherent electrons is at the heart of future quantum microscopy and spectroscopy technologies. Leveraging and specifically altering the quantum features of an electron beam source at low temperatures can enhance its emission properties. Here, we describe electron field emission from a monocrystalline, superconducting niobium nanotip at a temperature of 5.9 K. The emitted electron energy spectrum reveals an ultranarrow distribution down to 16 meV due to tunable resonant tunneling field emission via localized band states at a nanoprotrusion’s apex and a cutoff at the sharp low-temperature Fermi edge. This is an order of magnitude lower than for conventional field emission electron sources. The self-focusing geometry of the tip leads to emission in an angle of 3.7°, a reduced brightness of 3.8 × 10 8 A/(m 2 srV), and a stability of hours at 4.1 nA beam current and 69 meV energy width. This source will decrease the impact of lens aberration and enable new modes in low-energy electron microscopy, electron energy loss spectroscopy, and high-resolution vibrational spectroscopy.

36 MATERIALS SCIENCE↗

A Case for Electron-Astrophysics

Abstract The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the research theme of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. In addition, plasma electrons often play an important role for the spatial transfer of thermal energy due to the high heat flux associated with their velocity distribution. The regulation of this electron heat flux is likewise not understood. By focussing on these and other fundamental electron processes, the research theme of electron-astrophysics links outstanding science questions of great importance to the fields of space physics, astrophysics, and laboratory plasma physics. In this White Paper, submitted to ESA in response to the Voyage 2050 call, we review a selection of these outstanding questions, discuss their importance, and present a roadmap for answering them through novel space-mission concepts.

79 ASTRONOMY AND ASTROPHYSICS↗

Molecular Geometry and Electronic Structure of Copper Corroles

Copper corroles are known for their unique multiconfigurational electronic structures in the ground state, which arise from the transfer of electrons from the π orbitals of the corrole to the d-orbital of copper. While density functional theory (DFT) provides reasonably good molecular geometries, the determination of the ground spin state and the associated energetics is heavily influenced by functional choice, particularly the percentage of the Hartree–Fock exchange. Using extended multireference perturbation theory methods (XMS-CASPT2), the functional choice can be assessed. The molecular geometries and electronic structures of both the unsubstituted and the meso-triphenyl copper corroles were investigated. A minimal active space was employed for structural characterization, while larger active spaces are required to examine the electronic structure. The XMS-CASPT2 investigations conclusively identify the ground electronic state as a multiconfigurational singlet (S o ) with three dominant electronic configurations in its lowest energy and characteristic saddled structure. In contrast, the planar geometry corresponds to the triplet state (T o ), which is approximately 5 kcal/mol higher in energy compared to the S o state for both the bare and substituted copper corroles. Notably, the planarity of the T o geometry is reduced in the substituted corrole compared with that in the unsubstituted one. By analyzing the potential energy surface (PES) between the S o and T o geometries using XMS-CASPT2, the multiconfigurational electronic structure is shown to transition toward a single electron configuration as the saddling angle decreases (i.e., as one approaches the planar geometry). Despite the ability of the functionals to reproduce the minimum energy structures, only the TPSSh-D 3 PES is reasonably close to the XMS-CASPT2 surface. Significant deviations along the PES are observed with other functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating the Chemical Reactivity of DFT-Simulated Liquid Water with Hydrated Electrons via the Dual Descriptor

Modeling the various properties of liquid water, particularly its reactivity, has been a longstanding challenge for simulation methods. Recently, ab initio simulations based on density functional theory (DFT) have come to the fore as tenable methods for calculating the properties and reactivity of water, with varying degrees of success for different exchange-correlation functionals. In particular, hybrid-GGA and meta-GGA functionals have been shown to reproduce many of the structural, dynamical, and energetic properties of water to a high degree of accuracy relative to their computational cost. Here, we show that the dual descriptor (DD) measure of nucleophilicity and electrophilicity, which is sometimes used to elucidate organic chemistry reaction mechanisms, can also be used to characterize the reactivity of DFTsimulated liquid water. The DD is especially apt for understanding the reactivity of excess electrons with water as its calculation explicitly involves adding and removing an excess electron from a reference system. We use the DD to explore the reactivity of water simulated using three different DFT functionals: the LDA functional (LDA), a hybrid-GGA functional (PBE0), and a hybrid meta- GGA functional (SCAN0). Using the DD, we show that the SCAN0 functional with the standard 25% Hartree–Fock exchange produces simulated liquid water with many regions that are far more reactive than either PBE0 or LDA. To understand the implications of these highly reactive regions, we then add a strong nucleophile in the form of an excess electron and find that although PBE0 and LDA predict stable hydrated electrons, the excess electron reacts nearly instantaneously with SCAN0 water via proton abstraction to form a hydrogen atom and hydroxide ion. We show that the DD provides the ability to not only predict whether or not liquid water will react with a hydrated electron but also which particular waters will be involved solely from analyzing pure water configurations generated with each functional. We rationalize this result in terms of the known trap-seeking behavior of injected hydrated electrons, which are able to find the most electronegative region in bulk water. These results highlight the utility of the dual descriptor as a fast and interpretable method for investigating condensed-phase reactivity with excess electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-Time Time-Dependent Nuclear–Electronic Orbital Approach: Dynamics beyond the Born–Oppenheimer Approximation

The quantum mechanical treatment of both electrons and nuclei is crucial in nona- diabatic dynamical processes such as proton-coupled electron transfer. The nuclear electronic orbital (NEO) method provides an elegant framework for including nuclear quantum effects beyond the Born-Oppenheimer approximation. To enable the study of nonequilibrium properties, we derive and implement a real-time NEO (RT-NEO) approach based on time-dependent Hatree-Fock or density functional theory, in which both the electronic and nuclear degrees of freedom are propagated in a time-dependent variational framework. Nuclear and electronic spectral features can be resolved from the time-dependent dipole moment computed using the RT-NEO method. The test cases demonstrate the dynamical interplay between the quantum nuclei and the electrons through vibronic coupling. Moreover, vibrational excitation in the RT-NEO approach is demonstrated by applying a resonant driving eld, and electronic excitation is demonstrated by simulating excited state intramolecular proton transfer. This work shows that the RT-NEO approach is a promising tool to study nonadiabatic quantum dynamical processes within a time-dependent variational description for the coupled electronic and nuclear degrees of freedom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗