Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Free electron model”

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 811 records · Page 45

Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches

We address the calculation of microcanonical reaction rates for processes involving significant nuclear quantum effects using ring-polymer molecular dynamics (RPMD), both with and without electronically non-adiabatic transitions. After illustrating the shortcomings of the naive free-particle direct-shooting method, in which the temperature of the internal ring-polymer modes is set to the translational energy scale, we investigate alternative strategies based on the expression for the microcanonical rate in terms of the inverse Laplace transform of the thermal reaction rate. It is shown that simple application of the stationary-phase approximation (SPA) dramatically improves the performance of the microcanonical rates using RPMD, particularly in the low-energy region where tunneling dominates. Using the SPA as a Bayesian prior, numerically exact RPMD microcanonical rates are then obtained using maximum entropy inversion of the thermal reaction rates for both electronically adiabatic and non-adiabatic model systems. Finally, the direct-shooting method is revisited using the SPA-determined temperature for the internal ring-polymer modes, leading to a simple, direct-simulation method with improved accuracy in the tunneling regime. This work suggests a general strategy for the extraction of microcanonical dynamical quantities from RPMD (or other approximate thermal) simulations.

Tao, Xuecheng↗

Cryo-EM and MD infer water-mediated proton transport and autoinhibition mechanisms of V o complex

Rotary vacuolar adenosine triphosphatases (V-ATPases) drive transmembrane proton transport through a Vo proton channel subcomplex. Despite recent high-resolution structures of several rotary ATPases, the dynamic mechanism of proton pumping remains elusive. Here, we determined a 2.7-Å cryo–electron microscopy (cryo-EM) structure of yeast V o proton channel in nanodisc that reveals the location of ordered water molecules along the proton path, details of specific protein-lipid interactions, and the architecture of the membrane scaffold protein. Moreover, we uncover a state of V o that shows the c-ring rotated by ~14°. Molecular dynamics simulations demonstrate that the two rotary states are in thermal equilibrium and depict how the protonation state of essential glutamic acid residues couples water-mediated proton transfer with c-ring rotation. Our cryo-EM models and simulations also rationalize a mechanism for inhibition of passive proton transport as observed for free V o that is generated as a result of V-ATPase regulation by reversible disassembly in vivo.

59 BASIC BIOLOGICAL SCIENCES↗

Review of Solar Energetic Particle Models

Solar Energetic Particle (SEP) events are interesting from a scientific perspective as they are the product of a broad set of physical processes from the corona out through the extent of the heliosphere, and provide insight into processes of particle acceleration and transport that are widely applicable in astrophysics. From the operations perspective, SEP events pose a radiation hazard for aviation, electronics in space, and human space exploration, in particular for missions outside of the Earth’s protective magnetosphere including to the Moon and Mars. Thus, it is critical to improve the scientific understanding of SEP events and use this understanding to develop and improve SEP forecasting capabilities to support operations. Many SEP models exist or are in development using a wide variety of approaches and with differing goals. These include computationally intensive physics-based models, fast and light empirical models, machine learning-based models, and mixed-model approaches. The aim of this paper is to summarize all of the SEP models currently developed in the scientific community, including a description of model approach, inputs and outputs, free parameters, and any published validations or comparisons with data.

Kathryn Whitman↗

Vlasov Simulation of Electrostatic Solitary Structures in Multi-Component Plasmas

Electrostatic solitary structures have been observed in the Earth's magnetosheath by the Cluster spacecraft. Recent theoretical work has suggested that these solitary structures are modeled by electron acoustic solitary waves existing in a four-component plasma system consisting of core electrons, two counter-streaming electron beams, and one species of background ions. In this paper, the excitation of electron acoustic waves and the formation of solitary structures are studied by means of a one-dimensional electrostatic Vlasov simulation. The present result first shows that either electron acoustic solitary waves with negative potential or electron phase-space holes with positive potential are excited in four-component plasma systems. However, these electrostatic solitary structures have longer duration times and higher wave amplitudes than the solitary structures observed in the magnetosheath. The result indicates that a high-speed and small free energy source may be needed as a fifth component. An additional simulation of a five-component plasma consisting of a stable four-component plasma and a weak electron beam shows the generation of small and fast electron phase-space holes by the bump-on-tail instability. The physical properties of the small and fast electron phase-space holes are very similar to those obtained by the previous theoretical analysis. The amplitude and duration time of solitary structures in the simulation are also in agreement with the Cluster observation.

magnetohydrodynamics (MHD)↗

Magnetic structure of Fe 16 N 2 determined by polarized neutron diffraction on thin-film samples

Strain-free Fe 16 N 2 thin films were obtained on MgO (001) single-crystal substrates with Cr seed layers. Additionally, x-ray diffraction and transmission electron microscopy confirmed that the film is textured in all three spatial directions with crystallographic orientation determined by epitaxial constraints from the substrate. Average grain size was estimated to be in the range of tens of nanometers. Ten pieces of 40-nm-thick Fe 16 N 2 thin-film samples were stacked together for half-polarized neutron diffraction study at room temperature. A total of six reflections were usable, from which the magnetic structure of Fe 16 N 2 was determined and compared to existing models. Contrary to most electronic structure calculations, the present results suggest that the magnetic moment of Fe atoms at the 8h site is large and larger than that of the 4d site. Our results support the giant magnetization of Fe 16 N 2 .

36 MATERIALS SCIENCE↗

Measurement of deeply virtual Compton scattering off 4 He with the CEBAF Large Acceptance Spectrometer at Jefferson Lab

We report on the measurement of the beam spin asymmetry in the deeply virtual Compton scattering off $^4$He using the CEBAF Large Acceptance Spectrometer (CLAS) at Jefferson Lab using a 6 GeV longitudinally polarized electron beam incident on a pressurized $^4$He gaseous target. We detail the method used to ensure the exclusivity of the measured reactions, in particular the upgrade of CLAS with a radial time projection chamber to detect the low-energy recoiling $^4$He nuclei and an inner calorimeter to extend the photon detection acceptance at forward angles. Our results confirm the theoretically predicted enhancement of the coherent ($e^4$He$~\to~e'$$^4$He$'\gamma'$) beam spin asymmetries compared to those observed on the free proton, while the incoherent ($e^4$He$~\to~e'$p$'\gamma'$X$'$) asymmetries exhibit a 30$\%$ suppression. From the coherent data, we were able to extract, in a model-independent way, the real and imaginary parts of the only $^4$He Compton form factor, $\cal H_A$, leading the way toward 3D imaging of the partonic structure of nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

Microwave emission above steady and moving sunspots

Two-dimensional maps of radio brightness temperature and polarization, computed assuming thermal emission with free-free and gyroresonance absorption, are compared with observations of active region 2502, performed at Westerbork at lambda = 6.16 cm during a period of 3 days in June 1980. The computation is done assuming a homogeneous model in the whole field of view and a force-free extrapolation of the photospheric magnetic field observed at MSFC with a resolution of 2.34 arcsec. The mean results are the following: (1) a very good agreement is found above the large leading sunspot of the group, assuming a potential extrapolation of the magnetic field and a constant conductive flux in the transition region ranging from .2 x 10 to the 6th to 10 to the 7th erg/sq cm 5; (2) a strong radio source, associated with a new-born moving sunspot, cannot be ascribed to thermal emission. It is suggested that this source may be due to synchrotron radiation by mildly relativistic electrons accelerated by resistive instabilities occurring in the evolving magnetic configuration. An order-of-magnitude computation of the expected number of accelerated particles seems to confirm this hypothesis.

Drago, F. Chiuderi↗

Solar Wind Strahl Observations and Their Implication to the Core-Halo Formation due to Scattering

A study of the kinetic properties of the strahl electron velocity distribution functions (VDF?s) in the solar wind is presented. This study focuses on the mechanisms that control and regulate the electron VDF?s and the stability of the strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. These parameters are used to investigate the stability of the strahl population. The analysis check for whether the strahl electrons are constrained by some instability (e.g., the whistler or KAW instabilities), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling of the 3D-VDF?s from which the moments properties of the various populations are obtained. The results of this study have profound implication on the current hypothesis about the probable formation of the solar wind halo electrons produced from the scattering of the strahl population. This hypothesis is strengthened by direct observations of the strahl electrons being scattered into the core-halo in an isolated event. The observation implies that the scattering of the strahl is not a continuous process but occurs in bursts in regions where conditions for wave growth providing the scattering are optimum. Sometimes, observations indicate that the strahl component is anisotropic (Tper/Tpal approx. 2). This provides a possible free energy source for the excitation of whistler waves as a possible scattering mechanism, however this condition is not always observed. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer.

Vinas, Adolfo F.↗

X-ray photoelectron spectroscopy and rotating disk electrode measurements of smooth sputtered Fe-N-C films

Electrocatalysts for the oxygen reduction reaction (ORR) based on complexes of iron and nitrogen in a carbon matrix (Fe-N-C) are a promising alternative to platinum group metal (PGM) based catalysts in polymer electrolyte membrane (PEM) fuel cells. Further improvements of Fe-N-C catalysts would benefit from model thin film studies of activity and stability of catalytic sites, but synthesis of Fe-N-C model thin films is challenging. Here we report on synthesis and characterization of Fe-N-C thin films produced by co-sputtering iron and carbon in a reactive nitrogen atmosphere onto removable glassy carbon rotating disk electrode (RDE) tips. Scanning electron microscopy (SEM) measurements indicate that the Fe-N-C films deposited at high temperature are smoother than the films annealed at high temperature. Electrocatalytic activity measured on the thin Fe-N-C films is greater for both high-temperature samples than for the room-temperature sample. From the analysis of X-ray photoelectron spectroscopy (XPS) data, exposure of the films to high temperatures results in increased graphitization of the carbon with the Fe-N-C films, and increased relative amount of graphitic and hydrogenated nitrogen species. Overall the results of this study demonstrate the feasibility of a thin film model system approach for studying active sites in PGM-free catalysts.

14 SOLAR ENERGY↗

Fantastic n = 4: Ce5Co4+xGe13−ySny of the A n +1M n X3 n +1 homologous series

Ce-based intermetallics are of interest due to the potential to study the interplay of localized magnetic moments and conduction electrons. Our work on Ce-based germanides led to the identification of a new homologous series An+1MnX3n+1 (A = rare earth, M = transition metal, X = tetrels, and n = 1–6). This work presents the single-crystal growth, structure determination, and anisotropic magnetic properties of the n = 4 member of the Cen+1ConGe3n+1 homologous series. Ce5Co4+xGe13−ySny consists of three Ce sites, three Co sites, seven Ge sites, and two Sn sites, and the crystal structure is best modeled in the orthorhombic space group Cmmm where a = 4.3031(8) Å, b = 45.608(13) Å, and c = 4.3264(8) Å, which is in close agreement with the previously reported Sn-free analog where a = 4.265(1) Å, b = 45.175(9) Å, and c = 4.293(3) Å. Anisotropic magnetic measurements show Kondo-like behavior and three magnetic transitions at 6, 4.9, and 2.4 K for Ce5Co4+xGe13−ySny.

Weiland, Ashley (ORCID:0000000171983559)↗

Grid-Free 2D Plasma Simulations of the Complex Interaction Between the Solar Wind and Small, Near-Earth Asteroids

We present results from a new grid-free 2D plasma simulation code applied to a small, unmagnetized body immersed in the streaming solar wind plasma. The body was purposely modeled as an irregular shape in order to examine photoemission and solar wind plasma flow in high detail on the dayside, night-side, terminator and surface-depressed 'pocket' regions. Our objective is to examine the overall morphology of the various plasma interaction regions that form around a small body like a small near-Earth asteroid (NEA). We find that the object obstructs the solar wind flow and creates a trailing wake region downstream, which involves the interplay between surface charging and ambipolar plasma expansion. Photoemission is modeled as a steady outflow of electrons from illuminated portions of the surface, and under direct illumination the surface forms a non-monotonic or ''double-sheath'' electric potential upstream of the body, which is important for understanding trajectories and equilibria of lofted dust grains in the presence of a complex asteroid geometry. The largest electric fields are found at the terminators, where ambipolar plasma expansion in the body-sized night-side wake merges seamlessly with the thin photoelectric sheath on the dayside. The pocket regions are found to be especially complex, with nearby sunlit regions of positive potential electrically connected to unlit negative potentials and forming adjacent natural electric dipoles. For objects near the surface, we find electrical dissipation times (through collection of local environmental solar wind currents) that vary over at least 5 orders of magnitude: from 39 Micro(s) inside the near-surface photoelectron cloud under direct sunlight to less than 1 s inside the particle-depleted night-side wake and shadowed pocket regions

Asteroids↗

Far Outer Galaxy H II Regions

We have made a multifrequency (6, 3.6, and 2 cm), high-resolution (3"-6"), radio continuum survey of IRAS selected sources from the catalogue of Wouterloot & Brand (1989) to search for and study H II regions in the far outer Galaxy. We identified 31 sources in this catalog with well determined galactocentric distances, and with R approx.. greater than 15 kpc and L(sub FIR) approx.greater than 10(exp 4) solar luminosity, indicating the presence of high-mass star-formation. We have observed 11 of these sources with the Very Large Array (VLA). We observed the sources at 6 and 2 cm using "scaled arrays", making possible a direct and reliable comparison of the data at these two wavelengths for the determination of spectral indices. We detected a total of 12 radio sources, of which 10 have spectral indices consistent with optically-thin free-free emission from H II regions. Combined with previous VLA observations by other investigators, we have data on a total of 15 H II regions at galactocentric distances of 15 to 18.2kpc, among the most remote H II regions found in our Galaxy. The sizes of the H II regions range from approx. less than 0.10 to 2.3 pc. Using the measured fluxes and sizes, we determine the electron densities, emission measures, and excitation parameters of the H II regions, as well as the fluxes of Lyman continuum photons needed to keep the nebulae ionized. The sizes and electron densities are consistent with most of the sources detected in this survey being compact or ultracompact H II regions. Seven of the fifteen H II regions have sizes approx. less than 0.20 pc. Assuming simple pressure-driven expansion of the H II regions, these sizes indicate ages approx. less than 5 x 10(exp 4) yr, or only 1% of the lifetime of an O star, which implies an unlikely overabundance of O stars in the outer Galaxy. Thus, the large number of compact H II regions suggests that the time these regions spend in a compact phase must be much longer than their dynamical expansion times. Five of the fifteen H II regions have cometary shapes; the remainder are spherical or unresolved. Comparison of the radio continuum data with molecular line maps suggests that the cometary shape of the two H II regions in S 127 may be due to pressure confinement of the expanding ionized gas, as in the "blister" or "champagne flow" models of H II regions. Comparison of the radio continuum data with the IRAS far-infrared data indicates that the five most luminous H II regions are consistent with a single 0 or B star exciting a dust-free H II region. Subject headings: stars: formation - ISM: H II regions - ISM: individual objects: S 127 radio continuum: interstellar

Rudolph, A. L.↗

Criegee Intermediates Compete Well with OH as a Cleaning Agent for Atmospheric Amides

Elucidating the chemical kinetics of stabilized Criegee intermediates (sCIs) in the atmosphere is critically important for climate modeling. Here, we report a class of very rapid bimolecular reactions of two sCIs, namely, CH 2 OO and syn-CH 3 CHOO, with amides. We used electronic structure calculations and kinetics calculations to elucidate a universal mechanism by which the oxygen atom of the carbonyl group in an amide is added to the carbon atom of the COO group in the sCI with simultaneous transfer of the amide hydrogen to the terminal oxygen atom of the carbonyl oxide of the sCI. The barriers for the mechanism are submerged below reactants by ∼9 kcal/mol, which means that the tight transition state is not the rate-determining step, and the rate constants are determined by loose free energy bottlenecks between the reactants and the precursor complexes. We calculate the rate constants due to these loose and barrierless transition states by variable-reaction-coordinate variational transition-state theory. We find that bimolecular reactions of sCIs with amides are very rapid, with rate constants 1 to 5 × 10 –10 cm 3 molecule –1 s –1 , which is a factor of 2 faster than the rate constants of carboxylic acids with sCIs and 10 2 faster than the rate constants of OH reactions with amides. This has the consequence that, under a wide range of conditions, sCIs are the major sink for atmospheric amides, in contrast to the usual assumption that OH radicals are the dominant sinks.

Amides↗

Understanding the Role of Borohydride Doping in Electrochemical Stability of Argyrodite Li 6 PS 5 Cl Solid‐State Electrolyte

This work elucidates the mechanism by which lithium borohydride (LiBH 4 ) doping into argyrodite-type Li 6 PS 5 Cl (LBH-LPSCl) solid-state electrolyte (SSE) enhances electrochemical stability. State-of-the-art electrochemical performance is achieved with 5 wt% borohydride. Symmetric cells achieve critical current density (CCD) of 7.3 mA cm −2 , versus 2.6 mA cm −2 for baseline-LPSCl. All solid-state batteries (ASSBs) employing lithium metal and NMC811 cathode are stable over 400 cycles at 0.5C, with capacity retention of 83%. An anode-free ASSB (AF-ASSB) is stable over 600 cycles, with capacity loss of 0.04% per cycle. 5LBH-LPSCl allows for enhanced low temperature operation, down to −14 °C. Yet the difference in electrolytes’ bulk microstructures and hardnesses are minimal, while ionic conductivity is incrementally improved (≈50%). Theoretical modeling indicates limited effect of substitution on thermodynamic stability of PS 4 3- units, which decompose when contacting Li. Instead, enhanced electrochemical stability is site-specific kinetic effect: In situ electrodeposition experiments using X-ray photoelectron spectroscopy (XPS) and time-of-flight secondary ion mass spectrometry (TOF-SIMS) reveal tri-layer SEI based predominately on Li 3 P/LiBH 4 /Li 2 S that blocks electrons while facilitating ion transport. This SEI manifests reduced interface resistance and accelerated nucleation and growth of metallic Li. With baseline-LPSCl the SEI based on Li 3 P/Li 2 S is substantially thicker, generating localized stresses that promote interfacial cracking while cycling.

argyrodite↗

A reduced-model (nSOLT) simulation of neutral recycling effects on plasma turbulence in the divertor region of MAST-U

The 2D scrape-off-layer turbulence code (nSOLT), which includes 1D kinetic neutral–plasma interactions, is applied to study effects of neutral recycling on plasma turbulence for parameters illustrative of the MAST-U divertor region. Neutral recycling is modeled by injecting a fraction of the parallel plasma flux to the divertor back into the simulation domain as a source of Franck–Condon-distributed neutrals. Stationary sources, concentrated at the magnetic separatrix, model plasma streaming into the divertor region from the upstream scrape-off-layer and sustain plasma turbulence absent neutral recycling. Starting from one such no-neutrals equilibrium, we initiate recycling in a numerical experiment designed to diagnose and identify the effects of various neutral–plasma interactions on the divertor plasma, divertor turbulence, and plasma exhaust. The onset of recycling triggers an initial burst of enhanced cross field plasma transport that is quelled by ionization cooling and charge–exchange (CX) friction, with growing neutral pressure, leading to a quiescent, turbulence-free state. Diagnosis of this transient burst reveals that (1) the sudden increase in plasma density due to ionization dominates the onset of the burst, (2) electron cooling due to ionization increases collisionality and disconnects blob filaments from the sheath, and (3) CX friction drives tripole polarization of a blob that can dominate the curvature-driven dipole polarization, leading to the stagnation of blob propagation and reduced radial turbulent transport. It is shown that CX friction is negligible compared to sheath physics in determining equilibrium mean flow shearing rates, for parameters considered herein (specifically a short connection length to the divertor target), while it can significantly reduce interchange-instability growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Solving the “Coloring Problem” in InPd 3– x Ag x ( x = 0–0.7) by Phase Diagrams Modeling and Diffraction Experiments

Here, a series of InPd 3–x Ag x (x = 0–1) compositions were synthesized by conventional high-temperature synthesis, and as-synthesized samples were characterized by powder X-ray diffraction experiments. Up to x = 0.7, InPd 3–x Ag x adopts the ternary substitutional variant of the InPd 3 structure (TiAl 3 -type), when x > 0.7, elemental Ag starts to segregate along with the main phase. Accurate structural characterization in InPd 3–x Ag x faces a critical challenge due to the narrow X-ray scattering contrast among constituents In, Pd, and Ag and nearly identical neutron scattering lengths of Pd and Ag. To overcome this “coloring problem”, a combination of calculation of phase diagrams modeling (CALPHAD) and diffraction techniques (X-ray and neutron) was employed. In the compositional range 0 ≤ x ≤ 0.7, InPd 3–x Ag x presents a ternary variant of the TiAl 3 -type structure, where Ag atoms selectively substitute one (the 2b Wyckoff site) of the two Pd sites in InPd 3 . Notably, in contrast to the isologous InPd 3–x Cu x (x = 0–1) system, Ag substitution does not form an ordered VRh 2 Sn-type structure at the limiting composition. The distinct site preference in InPd 3–x Ag x is elucidated by charge population analysis, electronic structure calculations, and orbital-resolved chemical bonding investigations, and the extent of substitution is supported by formation free energy calculations.

36 MATERIALS SCIENCE↗