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

Stabilizing High–Nickel Cathodes with High–Voltage Electrolytes

Electrolytes connect the two electrodes in a lithium battery by providing Li + transport channels between them. Advanced electrolytes are being explored with high-nickel cathodes and the lithium-metal anode to meet the high energy density and cycle life goals, but the origin of the performance differences with different electrolytes is not fully understood. Here, the mechanisms involved in protecting the high-capacity, cobalt-free cathode LiNiO 2 with a model high-voltage electrolyte (HVE) are delineated. The kinetic barrier posed by a thick surface degradation layer with poor Li + -ion transport is found to be the major contributor to the fast capacity fade of LiNiO 2 with the conventional carbonate electrolyte. In contrast, HVE reduces the side reactions between the electrolyte and the electrodes, leading to a thinner nano-interphase layer comprised of more beneficial species. Crucially, the HVE leads to a different surface reorganization pathway involving the formation of a thinner nanoscale LiNi 2 O 4 spinel phase on the LiNiO 2 surface. Finally, with a high 3D Li + -ion and electronic conductivity, the spinel LiNi 2 O 4 reorganization nanolayer preserves fast Li + transport across the cathode–electrolyte interface, reduces reaction heterogeneity in the electrode and alleviates intergranular cracking within secondary particles, resulting in superior long-term cycle life.

25 ENERGY STORAGE↗
DOE OSTI2023

Probing the Mechanical Performance of Micro-architected Porous Structures Through In Situ Characterization and Analysis

Micropores play critical roles in both natural and man-made materials. Such pores take on a variety of shapes and sizes ranging from spherical to irregular sphere-like voids with diameters spanning from the nanometer to millimeter scales. When porous structures are mechanically loaded, the pores direct the stress around their free surfaces, altering the material’s mechanical response relative to fully dense materials. In this study, for the first time, we create micro-architected porous samples using nanolithography and investigate the role of pore morphology by conducting a series of in situ micropillar compression tests in scanning electron microscopy. The findings demonstrate that porosity is the primary factor influencing the mechanical response of these micro-architected materials, as often seen on a macroscopic level. Additionally, we observed that pore geometry had a significant impact on Young’s modulus, yield stress, and strain energy density as a secondary parameter. Then, the extracted Young’s modulus was compared to macroscopic empirical models and determined the analytical models sufficiently described the impact of porosity in the microscopic scale but failed to capture the impact of second-order parameters. Furthermore, these results suggest how porous materials can be tailored to achieve desired mechanical properties based on the engineering applications of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pairing instability on a Luttinger surface: A non-Fermi liquid to superconductor transition and its Sachdev-Ye-Kitaev dual

Superconductivity results from an instability of the Fermi surface—contour of poles of the single-particle propagator—to an infinitesimally small attraction between electrons. Here, we instead discuss the analogous problem on a model Luttinger surface, or contour of zeros of the Green function. At zero temperature ($β → \infty$) and a critical interaction strength ($u_{c∞}$) characterized by the residue of self-energy pole, we find that the pair susceptibility diverges leading to a superconducting instability. We evaluate the pair fluctuation partition function and find that the spectral density in the normal state has an interaction-driven, power-law $\frac{1}{\sqrt{ω}}$ type, van-Hove singularity (vHS), indicating non-Fermi liquid (NFL) physics. Crucially, in the strong coupling limit ($βu \gg 1$), the leading-order fluctuation free-energy terms in the normal state of this NFL-SC transition resembles the equivalent [$\textit{O}(1)$] terms of the Sachdev-Ye-Kitaev (SYK) model. This free-energy contribution takes a simple form $–βF = βu_{c∞} – γ \text{ln}(βu_{c∞}$), where $\gamma$ is a constant equal to $\frac{1}{2}$. Weak impurity scattering ($τ \gg β^{–1}$) leaves the low-energy spectral density unaffected but leads to an interaction-driven enhancement of superconductivity. Our results in this work shed light on the role played by order-parameter fluctuations in providing the key missing link between Mott physics and strongly coupled toy models exhibiting gravity duals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Light-front holography model of the EMC effect

A new two-component model of the EMC effect based on Light-Front Holographic QCD (LFHQCD) is presented. The model suggests the EMC effect is the result of the nuclear potential breaking SU(6) symmetry. The model separates the $F_{2}^{A}$ nuclear structure function into two parts: a free contribution, involving the addition of proton and neutron structure functions weighted by the number of protons and neutrons respectively, and a nuclear/medium modified contribution that involves a universal function for all nuceli. Further, the model displays a connection with the correlation between the size of the EMC effect and the SRC pair density, a 2 - extracted from kinematic plateaus at around x > 1 in inclusive quasi-elastic (QE) scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

First-principles ionized-impurity scattering and charge transport in doped materials

Scattering of carriers with ionized impurities governs charge transport in doped semiconductors. However, electron interactions with ionized impurities cannot be fully described with quantitative first-principles calculations, so their understanding relies primarily on simplified models. Here we show an ab initio approach to compute the interactions between electrons and ionized impurities or other charged defects. It includes the short- and long-range electron-defect (e-d) interactions on equal footing and allows for efficient interpolation of the e-d matrix elements. Here we combine the e-d and electron-phonon interactions in the Boltzmann transport equation to compute the carrier mobilities in doped silicon over a wide range of temperature and doping concentrations, seamlessly spanning the defect- and phonon-limited transport regimes. The individual contributions of the defect- and phonon-scattering mechanisms to the carrier relaxation times and mean-free paths are analyzed. Our method provides a powerful tool to study electronic interactions in doped materials. It broadens the scope of first-principles transport calculations, enabling studies of a wide range of doped semiconductors and oxides with application to electronics, energy, and quantum technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

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↗
DOE OSTI2020

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↗
DOE OSTI2020

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↗
DOE OSTI2021

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↗
DOE OSTI2026

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↗
DOE OSTI2020

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↗
DOE OSTI2021

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)↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
DOE OSTI2023

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↗
DOE OSTI2025

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↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗