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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Accelerating Finite-Temperature Kohn-Sham Density Functional Theory with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Accelerating Finite-temperature Kohn-Sham Density Functional Theory\ with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2021

Active Phase on SrCo 1–x Fe x O -δ (0 ≤ x ≤ 0.5) Perovskite for Water Oxidation: Reconstructed Surface versus Remaining Bulk

Perovskite oxides based on earth-abundant transition metals have been extensively explored as promising oxygen evolution reaction (OER) catalysts in alkaline media. The (electro)chemically induced transformation of their initially crystalline surface into an amorphous state has been reported for a few highly active perovskite catalysts. However, little knowledge is available to distinguish the contribution of the amorphized surface from that of the remaining bulk toward the OER. In this work, we utilize the promoting effects of two types of Fe modification, i.e., bulk Fe dopant and Fe ions absorbed from the electrolyte, on the OER activity of SrCoO 3-δ model perovskite to identify the active phase. Transmission electron microscopy and X-ray photoelectron spectroscopy confirmed the surface amorphization of SrCoO 3-δ as well as SrCo 0.8 Fe 0.2 O 3-δ after potential cycling in Fe-free KOH solution. By further cycling in Fe-spiked electrolyte, Fe was incorporated into the amorphized surface of SrCoO 3-δ (SrCoO 3-δ + Fe 3+ ), yielding approximately sixfold increase in activity. Despite the difference in remaining perovskites, SrCoO 3-δ + Fe 3+ and SrCo 0.8 Fe 0.2 O 3-δ exhibited remarkably similar activity. These results reflect that the in situ developed surface species are directly responsible for the measured OER activity, whereas the remaining bulk phases have little impact.

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

Electronic glasses from a broken gauge symmetry in disorder-free systems

Glass phases can be stabilized by quenched disorders, as in most spin-glass materials, or self-generated through kinetic freezing in disorder-free systems. A canonical example of the latter is structural glasses, which have been extensively studied for many decades. Yet, how the rugged energy landscape of a glass phase is spontaneously generated in disorder-free systems remains one of the key questions in glass physics. Here, in this work, we present a general electronic mechanism for the emergence of glassy phase using the example of itinerant electrons coupled to XY spins on a lattice. This model can also be viewed as the mean-field theory of a superconducting system with attractive density-density interactions. Intriguingly, the electron gauge symmetry in the strong pairing limit gives rise to a macroscopic degeneracy of XY spins. In the presence of electron hopping that breaks the gauge symmetry, the lifting of the extensive degeneracy leads to a glass phase with disordered pairings. Our findings highlight a scenario in which a glassy state originates from the breaking of quantum gauge symmetry without quenched disorders.

XY model↗
DOE OSTI2023

Simulated sulfur K-edge X-ray absorption spectroscopy database of lithium thiophosphate solid electrolytes

X-ray absorption spectroscopy (XAS) is a premier technique for materials characterization, providing key information about the local chemical environment of the absorber atom. In this work, we develop a database of sulfur K-edge XAS spectra of crystalline and amorphous lithium thiophosphate materials based on the atomic structures reported in Chem. Mater., 34, 6702 (2022). The XAS database is based on simulations using the excited electron and core-hole pseudopotential approach implemented in the Vienna Ab initio Simulation Package. Our database contains 2681 S K-edge XAS spectra for 66 crystalline and glassy structure models, making it the largest collection of first-principles computational XAS spectra for glass/ceramic lithium thiophosphates to date. This database can be used to correlate S spectral features with distinct S species based on their local coordination and short-range ordering in sulfide-based solid electrolytes. The data is openly distributed via the Materials Cloud, allowing researchers to access it for free and use it for further analysis, such as spectral fingerprinting, matching with experiments, and developing machine learning models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Phonon screening and dissociation of excitons at finite temperatures from first principles

The properties of excitons, or correlated electron–hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron–hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio GW -Bethe–Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and SrTiO 3 . We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron–hole pairs.

Science & Technology - Other Topics↗
DOE OSTI2022

Optimization of echo-enabled harmonic generation toward coherent EUV and soft X-ray free-electron laser at NSLS-II

Prebunching via echo-enabled harmonic generation (EEHG) is an efficient way to reduce the radiator length and improve the longitudinal coherence as well as output stability in storage-ring-based free-electron lasers. We propose a conceptual design, which uses two straight sections to seed coherent extreme-ultraviolet (EUV) and soft X-ray emission with nearly MHz repetition rate. To take the large energy spread (10 –3 ) of a storage ring into account and utilize the existing bending magnets between the two straight sections as the first chicane, we implement a special modeling tool, named EEHG optimizer. This tool has been successfully applied to maximize the prebunching with a reasonably low energy modulation, thereby generating intense coherent X-ray pulses within a short undulator length (a few meters) limited by the available space of a storage ring. Numerical simulations confirm that the optimized EEHG parameters can be directly applied to generate a 10 MW scale peak power with fully coherent ultrafast EUV to soft X-ray pulses based on the NSLS-II parameters. This method can be easily extended to other types of diffraction-limited storage rings.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

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

Ion emission from warm dense matter produced by irradiation with a soft x-ray free-electron laser

Here we report on an experiment performed at the FLASH2 free-electron laser (FEL) aimed at producing warm dense matter via soft x-ray isochoric heating. In the experiment, we focus on study of the ions emitted during the soft x-ray ablation process using time-of-flight electron multipliers and a shifted Maxwell–Boltzmann velocity distribution model. We find that most emitted ions are thermal, but that some impurities chemisorbed on the target surface, such as protons, are accelerated by the electrostatic field created in the plasma by escaped electrons. The morphology of the complex crater structure indicates the presence of several ion groups with varying temperatures. We find that the ion sound velocity is controlled by the ion temperature and show how the ion yield depends on the FEL radiation attenuation length in different materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Dihadron azimuthal asymmetry and light-quark dipole moments at the Electron-Ion Collider

Here, we propose a novel method to probe light-quark dipole moments by examining the azimuthal asymmetries between a collinear pair of hadrons in semi-inclusive deep inelastic lepton scattering off an unpolarized proton target at the Electron-Ion Collider. These asymmetries provide a means to observe transversely polarized quarks, which arise exclusively from the interference between the dipole and the Standard Model interactions, thereby depending linearly on the dipole couplings. We demonstrate that this novel approach can enhance current constraints on light-quark dipole operators by an order of magnitude, free from contamination of other new physics effects. Furthermore, it allows for a simultaneous determination of both the real and imaginary parts of the dipole couplings, offering a new avenue for investigating potential CP-violating effects at high energies.

dihadron fragmentation↗
DOE OSTI2024

Magnetotransport of Sm 2 Ir 2 O 7 across the pressure-induced quantum-critical phase boundary

Rare-earth pyrochlore iridates host two interlocking magnetic sublattices of corner-sharing tetrahedra and can harbour a unique combination of frustrated moments, exotic excitations and highly correlated electrons. They are also the first systems predicted to display both topological Weyl semimetal and axion insulator phases. We have measured the transport and magnetotransport properties of single-crystal Sm 2 Ir 2 O 7 up to and beyond the pressure-induced quantum critical point for all-in-all-out (AIAO) Ir order at p c = 63 kbar previously identified by resonant X-ray scattering and close to which Weyl semimetallic behavior has been previously predicted. Our findings overturn the accepted expectation that the suppression of AIAO order should lead to metallic conduction persisting down to zero temperature. Instead, the resistivity-minimum temperature, which tracks the decrease in the AIAO ordering temperature for pressures up to 30 kbar, begins to increase under further application of pressure, pointing to the presence of a second as-yet unidentified mechanism leading to non-metallic behavior. The magnetotransport does track the suppression of Ir magnetism, however, with a strong hysteresis observed only within the AIAO phase boundary, similar to that found for Ho 2 Ir 2 O 7 and attributed to plastic deformation of Ir domains. Around pc we find the emergence of a new type of electronic phase, characterized by a negative magnetoresistance with small hysteresis at the lowest temperatures, and hysteresis-free positive magnetoresistance above approximately 5K. The temperature dependence of our low-temperature transport data are found to be best described by a model consistent with a Weyl semimetal across the entire pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DFT (and DMFT) for DFT users

The interesting properties of f -electron systems comes from a simultaneous existence of free and discrete level electrons. The challenge we have is to describe free electrons and discrete level electrons equally well in a unified picture. The ELF index can be used to find regions in real space where discrete level physics needs to be taken into account. We have identified the HO gas as a model system that can be used to gain insight about this kind of physics. We will use the HO gas model system for creating a functional suitable for these systems via the subsystem functional scheme. With a predictive capability for the rare earths and lanthindes, we can identify new materials with new/improved properties. This is crucial for, e.g., developing material and technology substitutes for addressing the challenges associated with the use of rare earth elements and other critical materials.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Infrared excesses in early-type stars - Gamma Cassiopeiae

Spectrophotometry of the classical Be star Gamma Cas (1-4 microns, with about 2% spectral resolution) is presented. These data, together with existing broad-band observations, are accurately described by simple isothermal LTE models for the IR excess which differ from most previously published work in three ways: (1) hydrogenic bound-free emission is included; (2) the attenuation of the star by the shell is included; and (3) no assumption is made that the shell contribution is negligible in some bandpass. It is demonstrated that the bulk of the IR excess consists of hydrogenic bound-free and free-free emission from a shell of hot ionized hydrogen gas, although a small thermal component cannot be ruled out. The bound-free emission is strong, and the Balmer, Paschen, and Brackett discontinuities are correctly represented by the shell model with physical parameters as follows: a shell temperature of approximately 18,000 K, an optical depth (at 1 micron) of about 0.5, an electron density of approximately 1 trillion per cu cm, and a size of about 2 trillion cm. Phantom shells (i.e., ones which do not alter the observed spectrum of the underlying star) are discussed.

Scargle, J. D.↗
NASA NTRS1995

Chromospheric evaporation and decimetric radio emission in solar flares

We have discovered decimetric signatures of the chromospheric evaporation process. Evidence for the radio detection of chromospheric evaporation is based on the radio-inferred values of (1) the electron density, (2) the propagation speed, and (3) the timing, which are found to be in good agreement with statistical values inferred from the blueshifted Ca XIX soft X-ray line. The physical basis of our model is that free-free absorption of plasma emission is strongly modified by the steep density gradient and the large temperature increase in the upflowing flare plasma. The steplike density increase at the chromospheric evaporation front causes a local discontinuity in the plasma frequency, manifested as almost infinite drift rate in decimetric type III bursts. The large temperature increase of the upflowing plasma considerably reduces the local free-free opacity (due to the T(exp -3/2) dependence) and thus enhances the brightness of radio bursts emitted at the local plasma frequency near the chromospheric evaporation front, while a high-frequency cutoff is expected in the high-density regions behind the front, which can be used to infer the velocity of the upflowing plasma. From model calculations we find strong evidence that decimetric bursts with a slowly drifting high-frequency cutoff are produced by fundamental plasma emission, contrary to the widespread belief that decimetric bursts are preferentially emitted at the harmonic plasma level. We analyze 21 flare episodes from 1991-1993 for which broadband (100-3000 MHz) radio dynamic spectra from Pheonix, hard X-ray data from (BATSE/CGRO) and soft X-ray data from Burst and Transient Source Experiment/Compton Gamma Ray Observatory (GOES) were available.

Aschwanden, Markus J.↗
NASA NTRS1999

Inter-Diffusion in the Presence of Free Convection

Because of their technological importance, establishment of the precise values of interdiffusion coefficients is important in multicomponent fluid systems. Such values are not available because diffusion is influenced by free convection due to compositionally induced density variations. In this project, earth based diffusion experiments are being performed in a viscous fluid system PbO-SiO2 at temperatures between 500-1000 C. This system is chosen because it shows a large variation in density with small changes in composition and is expected to show a large free convection effect. Infinite diffusion couples at different temperatures and times are being studied with different orientations with respect to gravity. Composition fields will be measured using an Electron Microprobe Analyzer and will be compared with the results of a complementary modeling study to extract the values of the true diffusion coefficient from the measured diffusion profiles.

Gupta, Prabhat K.↗
NASA NTRS1996

The Positronium Radiative Combination Spectrum: Calculation in the Limit of Thermal Positrons and Low Densities

We calculate the intensities of the positronium de-excitation lines for two processes: (1) the radiative combination of free thermal electrons and positrons for transitions with principal quantum number n less than 20, and (2) charge exchange between free positrons and hydrogen and helium atoms, restricting our evaluation to the Lyman-alpha line. We consider a low-density medium modeled by the case A assumption of Baker & Menzel and use the "nL method" of Pengelly to calculate the absolute intensities. We also evaluate the positronium fine and hyperfine intensities and show that these transitions are in all cases much weaker than positronium de-excitation lines in the same wavelength range. We also extrapolate our positronium de-excitation intensities to the submillimeter, millimeter, and centimeter wavelengths. Our results favor the search of infrared transitions of positronium lines for point sources when the visual extinction A, is greater than approx. 5.

Wallyn, P.↗