Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electron density profile”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Efficient Modeling of Structural, Electronic, and Optical Properties of Silver and Gold Metal Nanoclusters and Alloys Using Optimized SCC-DFTB Parameters

Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.

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

Why Are There so Few Reports of High-Energy Electron Drift Resonances? Role of Radial Phase Space Density Gradients

Models of monochromatic Pc5 (2–7 mHz) ultralow frequency (ULF) wave interactions with high energy (greater than ~1 MeV) electrons predict drift resonant interactions that can cause rapid radial transport and acceleration. There are few reports of electron drift resonance at energies greater than ~1 MeV, in contrast to lower energies; moreover, all previous reports occur in the aftermath of interplanetary shocks. These two facts are difficult to reconcile with theory and numerical simulations predicting that greater than ~1 MeV drift resonances should occur more often and in a wider variety of driving conditions. In this study, we show that a combination of observational sampling biases and nominal radial phase space density gradients is one explanation for this discrepancy between theory and observations. In particular, we examine electron dynamics in two case studies with very similar satellite coverage, solar wind conditions, and Pc5 wave properties, yet with different radial phase space density profiles. Using global wave and particle observations, we show that the events have vastly different particle responses despite having similar wave properties. Placing these results in context with past studies, we further show here that nominal radial PSD gradients near geostationary orbit can mask the expected drift resonance particle response and explain (1) the small number of past greater than ~1 MeV drift resonance reports and (2) the restriction of these reports to interplanetary shock events. We argue that future observational studies characterizing radial transport via drift resonance should examine global particle dynamics, including observations of the radial phase space density profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Interpreting AI for fusion: An application to plasma profile analysis for tearing mode stability

Artificial intelligence models have demonstrated strong predictive capabilities for various instabilities in fusion devices such as Tokamaks, including tearing modes (TM), edge localized modes, and disruptive events, but their opaque nature raises concerns about safety and trustworthiness when applied to fusion power plants. Here, we present a physics-based interpretation framework using a TM prediction model as a demonstration that is validated through a dedicated DIII-D TM avoidance experiment. By applying Shapley analysis, we identify how profiles such as rotation, temperature, and density contribute to the model's prediction of TM stability. Our analysis shows that in our experimental scenario, core electron temperature and rotation peaking play the primary role in TM stability, while density changes have smaller effects on stability. We show that off-axis ion temperature stabilizes TMs, suggesting that off-axis neutral beam heating can further stabilize this scenario. This work presents a generalizable ML-based event prediction methodology, from training to physics-driven interpretation, bridging the gap between physics understanding and opaque ML models.

Farre-Kaga, Hiro J. [Princeton Univ., NJ (United S↗
DOE OSTI2022

Boron nitride on SiC(0001)

Here, in the field of van der Waals heterostructures, the twist angle between stacked two-dimensional layers has been identified to be of utmost importance for the properties of the heterostructures. In this context, we previously reported the growth of a single layer of unconventionally oriented epitaxial graphene that forms in a surfactant atmosphere [F. C. Bocquet et al., Phys. Rev. Lett. 125, 106102 (2020)]. The resulting G-R0° layer is aligned with the SiC lattice, and hence represents an important milestone towards high-quality twisted bilayer graphene, a frequently investigated model system in this field. Here, we focus on the surface structures obtained in the same surfactant atmosphere, but at lower preparation temperatures at which a boron nitride template layer forms on SiC(0001). In a comprehensive study based on complementary experimental and theoretical techniques, we find—in contrast to the literature—that this template layer is a hexagonal B x N y layer, but not high-quality hBN. It is aligned with the SiC lattice and gradually replaced by low-quality graphene in the 0° orientation of the B x N y template layer upon annealing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Kinetic modeling of neutral transport for a continuum gyrokinetic code

In this work we present the first-of-its-kind coupling of a continuum full- f gyrokinetic turbulence model with a 6D continuum model for kinetic neutrals, carried out using the Gkeyll code. Our objective is to improve the first-principle understanding of the role of neutrals in plasma fueling, detachment, and their interaction with edge plasma profiles and turbulence statistics. Our model includes only atomic hydrogen and incorporates electron-impact ionization, charge exchange, and wall recycling. These features have been successfully verified with analytical predictions and benchmarked with the DEGAS2 Monte Carlo neutral code. We carry out simulations for a scrape-off layer (SOL) with simplified geometry and National Spherical Torus Experiment parameters. We compare these results to a baseline simulation without neutrals and find that neutral interactions reduce the normalized density fluctuation levels and associated skewness and kurtosis, while increasing auto-correlation times. A flatter density profile is also observed, similar to the SOL density shoulder formation in experimental scenarios with high fueling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Quantum Spin Torque Driven Transmutation of an Antiferromagnetic Mott Insulator

The basic model of spin-transfer torque (STT) in antiferromagnetic spintronics considers the exchange of angular momentum between quantum spins of flowing electrons and noncollinear-to-them localized spins treated as classical vectors. These vectors are assumed to realize Néel order in equilibrium, ↑ ↓ ∙∙∙ ↑ ↓, and their STT-driven dynamics is described by the Landau-Lifshitz-Gilbert (LLG) equation. However, many experimentally employed materials (such as archetypal NiO) are strongly electron-correlated antiferromagnetic Mott insulators (AFMIs) whose localized spins form a ground state quite different from the unentangled Néel state | ↑ ↓ ∙∙∙ ↑ ↓$\rangle$. The true ground state is entangled by quantum spin fluctuations, leading to the expectation value of all localized spins being zero, so that LLG dynamics of classical vectors of fixed length rotating due to STT cannot even be initiated. Instead, a fully quantum treatment of both conduction electrons and localized spins is necessary to capture the exchange of spin angular momentum between them, denoted as quantum STT. We use a recently developed time-dependent density matrix renormalization group approach to quantum STT to predict how injection of a spin-polarized current pulse into a normal metal layer coupled to an AFMI overlayer via exchange interaction and possibly small interlayer hopping—mimicking, e.g., topological-insulator/NiO bilayer employed experimentally—will induce a nonzero expectation value of AFMI localized spins. This new nonequilibrium phase is a spatially inhomogeneous ferromagnet with a zigzag profile of localized spins. The total spin absorbed by AFMI increases with electron-electron repulsion in AFMIs, as well as when the two layers do not exchange any charge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Characterization of an aerosolized nanoparticle beam beyond the diffraction limit through strong field ionization

Abstract The study of nanomaterials is an active area of research for technological applications as well as fundamental science. A common method for studying properties of isolated nanoparticles is by an in-vacuum particle beam produced via an aerodynamic lens. Despite being common practice, characterization of such beams has proven difficult as light scattering detection techniques fail for particles with sizes beyond the diffraction limit. Here we present a new technique for characterizing such nanoparticle beams using strong field ionization. By focusing an ultrafast, mJ-level laser into the particle beam, a nanoparticle within the laser focus is ionized and easily detected by its ejected electrons. This method grants direct access to the nanoparticle density at the location of the focus, and by scanning the focus through the transverse and longitudinal profiles of the particle beam we attain the 3-dimensional particle density distribution for a cylindrically symmetric beam. Further, we show that strong field ionization is effective in detecting spherical nanoparticles as small as 10 nm in diameter. Additionally, this technique is an effective tool in optimizing the particle beam for specific applications. As an example we show that the particle beam density and width can be manipulated by restricting the gas flow into the aerodynamic lens.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Refined interpretation of electron temperature response to neutral beam injection at DIII-D

Accurate particle and power deposition profiles of neutral beam injection (NBI) are essential to transport studies, and that information is usually acquired through Monte Carlo simulations with a given collisional model. The deposition process of the energetic beam particles leads to the informative electron temperature (T e ) evolution trajectory, which can be captured by electron cyclotron emission (ECE) system due to its good spatial and temporal resolution. Previously, some work has been done to interpret the T e responses to the pulsed NBI as a linear heating source with Fourier-based techniques, although that approach fell short when the fast ion slowing-down time becomes significant (~100 ms). It has been observed in DIII-D that the modulated NBI pulses (10–50 Hz) reduce local core T e values ~0.1 keV through cold electron dilution in high-T e (>2 keV) plasmas alongside accumulative heating. Here, a novel approach to interpret the T e response to NBI was developed by linearizing and modeling the detailed T e evolution trajectory using coherently averaged ECE data based on the different time scales of the terms in the local power and particle balance equations. The technique does not require absolute calibrations of ECE and is independent of collisional models. The resulting beam deposition profiles show good consistency and reasonable agreement with Monte Carlo calculations based on the atomic data from the Atomic Data and Analysis Structure (ADAS). Local electron density response measured by Thomson scattering (TS) also suggests the same features when the beam pulse is large enough for that diagnostic to resolve. The remaining discrepancies are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Constraining solar electron number density via neutrino flavor data at Borexino

Understanding the physics of the deep solar interior, and the more exotic environs of core-collapse supernovae (CCSN) and binary neutron-star (NS) mergers, is of keen interest in many avenues of research. To date, this physics is based largely on simulations via forward integration. While these simulations provide valuable constraints, it could be insightful to adopt the "inverse approach" as a point of comparison. Within this paradigm, parameters of the solar interior are not output based on an assumed model, but rather are inferred based on real data. Here, we take the specific case of solar electron number density, which historically is taken as output from the standard solar model. We show how one may arrive at an independent constraint on that density profile based on available neutrino flavor data from the Earth-based Borexino experiment. The inference technique's ability to offer a unique lens on physics can be extended to other datasets, and to analogous questions for CCSN and NS mergers, albeit with simulated data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A Study of a Surrogate Nuclear Thermal Propulsion Fuel Element

We report a surrogate fuel, Zirconium Carbide/Titanium Nitride has been prepared using a simple method to obtain a more uniform distribution of titanium nitride in the zirconium carbide matrix. Powders were sintered using Spark Plasma Sintering at 1750oC and 1850oC using a modified sintering profile. This resulted in a fairly uniform distribution of titanium nitride in the zirconium carbide matrix with a density of 95% of theoretical at 1850oC. Scanning electron microscopy, indentation hardness and x-ray diffraction results are given.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First Evidence of Local E × B Drift in the Divertor Influencing the Structure and Stability of Confined Plasma near the Edge of Fusion Devices

The structure of the edge plasma in a magnetic confinement system has a strong impact on the overall plasma performance. Here we uncover for the first time a magnetic-field-direction dependent density shelf, i.e., local flattening of density radial profile near the magnetic separatrix, in high confinement plasmas with low edge collisionality in the DIII-D tokamak. The density shelf is correlated with a doubly peaked density profile near the divertor target plate, which tends to occur for operation with the ion BX∇B drift direction away from the X-point, as currently employed for DIII-D advanced tokamak scenarios. This double-peaked divertor plasma profile is connected via the EXB drifts, arising from a strong radial electric field induced by the radial electron temperature gradient near the divertor target. The drifts lead to the reversal of the poloidal flow above the divertor target, resulting in the formation of the density shelf. The edge density shelf can be also enhanced at higher heating power, preventing large, periodic bursts of the plasma, i.e. Edge-localized Modes, in the edge region, consistent with ideal magnetohydrodynamics calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Synthesis and Characterization of Pyridine Dipyrrolide Uranyl Complexes

The first actinide complexes of the pyridine dipyrrolide (PDP) ligand class, ( Mes PDP Ph )UO 2 (THF) and ( Cl2Ph PDP Ph )UO 2 (THF), are reported as the U VI uranyl adducts of the bulky aryl substituted pincers ( Mes PDP Ph ) 2– and ( Cl2Ph PDP Ph ) 2– (derived from 2,6-bis(5-(2,4,6-trimethylphenyl)-3-phenyl-1H-pyrrol-2-yl)pyridine (H 2 Mes PDP Ph , Mes = 2,4,6-trimethylphenyl), and 2,6-bis(5-(2,6-dichlorophenyl)-3-phenyl-1H-pyrrol-2-yl)pyridine (H 2 Cl2Ph PDP Ph , Cl 2 Ph = 2,6-dichlorophenyl), respectively). Following the in situ deprotonation of the proligand with lithium hexamethyldisilazide to generate the corresponding dilithium salts (e.g., Li 2 Ar PDP Ph , Ar = Mes of Cl 2 Ph), salt metathesis with [UO 2 Cl 2 (THF) 2 ] 2 afforded both compounds in moderate yields. The characterization of each species has been undertaken by a combination of solid- and solution-state methods, including combustion analysis, infrared, electronic absorption, and NMR spectroscopies. In both complexes, single-crystal X-ray diffraction has revealed a distorted octahedral geometry in the solid state, enforced by the bite angle of the rigid meridional ( Ar PDP Ph ) 2– pincer ligand. The electrochemical analysis of both compounds by cyclic voltammetry in tetrahydrofuran (THF) reveals rich redox profiles, including events assigned as U VI /U V redox couples. Here, a time-dependent density functional theory study has been performed on ( Mes PDP Ph )UO 2 (THF) and provides insight into the nature of the transitions that comprise its electronic absorption spectrum.

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

Longitudinal tapering in gas jets for increased efficiency of 10-GeV class laser plasma accelerators

Modern laser plasma accelerators often require plasma waveguides tens of centimeters long to propagate a high-intensity drive laser pulse. Tapering the longitudinal gas density profile in 10 cm scale gas jets could allow for single stage laser plasma acceleration well beyond 10 GeV with current petawatt-class laser systems. Via simulation and interferometry measurements, we show density control by longitudinally adjusting the throat width and jet angle. Density profiles appropriate for tapering were calculated analytically and via particle-in-cell simulations and were matched experimentally. Further, these simulations show that tapering can increase electron beam energy using 19 J laser energy from ∼9 GeV to > 12 GeV in a 30 cm plasma and the accelerated charge by an order of magnitude.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Influence of conical wire array geometry on plasma flow and temperature profiles of radiatively cooled jets

The influence of conical wire array geometry on the formation and dynamics of pulsed-power driven plasma jets is investigated. In the experiments, the jet becomes isolated from the inflows as it passes through an aperture, allowing the study of its intrinsic evolution for different conical angles. Here, our results show that, regardless of the array opening angle, the jets are supersonic, highly collisional, and exhibit an exponential axial density decay with a characteristic scale length of 𝐿 𝑛 ≈ 3 mm, significantly shorter than the overall length of the jet. In contrast, axial velocity systematically increases with larger array opening angles. The near invariance of the density profile is attributed to a compensating mechanism between geometric divergence and axial acceleration, consistent with an asymptotically steady flow regime. Additionally, temperature measurements reveal ion-electron thermal decoupling near the base of the jet, with equilibration downstream. Analysis of radiative cooling and collisional energy exchange timescales indicates that the plasma evolves in a radiatively stable regime and that the two-temperature profile is fully described by the competition between both methods. These findings provide valuable information into the internal structure and evolution of pulsed-power plasma jets in regimes dominated by geometric shaping and radiative cooling, demonstrating experimental control over flow acceleration and collimation relevant to both laboratory and astrophysical contexts.

Physics - Plasma physics↗
DOE OSTI2021

Global Theory of Microtearing Modes in the Tokamak Pedestal

The pedestal of H-mode tokamaks displays strong magnetic fluctuations correlated with the evolution of the electron temperature. The microtearing mode (MTM), a temperature-gradient-driven instability that alters magnetic topology, is compatible with these observations. Here we extend the conventional theory of the MTM to include the global variation of the temperature and density profiles. The offset between the rational surface and the location of the pressure gradient maximum (μ) emerges as a crucial parameter for MTM stability. In conclusion, the extended theory matches observations on the Joint European Torus tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Evidence and modeling of turbulence bifurcation in L-mode confinement transitions on Alcator C-Mod

Analysis and modeling of rotation reversal hysteresis experiments show that a single turbulent bifurcation is responsible for the Linear to Saturated Ohmic Confinement (LOC/SOC) transition and concomitant intrinsic rotation reversal on Alcator C-Mod. Plasmas on either side of the reversal exhibit different toroidal rotation profiles and therefore different turbulence characteristics despite the profiles of density and temperature, which are indistinguishable within measurement uncertainty. Elements of this bifurcation are also shown to persist for auxiliary heated L-modes. The deactivation of subdominant (in the linear growth rate and contribution to heat transport) ion temperature gradient and trapped electron mode instabilities is identified as the only possible change in turbulence within a reduced quasilinear transport model across the reversal, which is consistent with the measured profiles and inferred heat and particle fluxes. Experimental constraints on a possible change from strong to weak turbulence, outside the description of the quasilinear model, are also discussed. These results indicate an explanation for the LOC/SOC transition that provides a mechanism for the hysteresis through the dynamics of subdominant modes and changes in their relative populations and does not involve a change in the most linearly unstable ion-scale drift-wave instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Terahertz Nano-Imaging of Electronic Strip Heterogeneity in a Dirac Semimetal

Emerging topological semimetals offer promise of realizing topological electronics enabled by terahertz (THz) current persistent against impurity scattering. Yet most fundamental issues remain on how to image nanoscale conductivity inhomogeneity. Here in this paper we show noninvasive and contactless conductivity mapping at THz-nm limit of electronic heterogeneity and nanostrip junctions in a Dirac material ZrTe 5 . A clear Dirac Fermion density transition, manifested as the exclusive THz conductivity contrast, is quantitatively analyzed and profiled on both sides of the junction. This also allows the determination of variable junction width of ~25–220 nm, depending on the THz conductivity contrast of adjacent strips. The unique THz-nm contrast is absent in mid-infrared nano-imaging measurements since topological semimetals with small Fermi pockets exhibit a better matching of their plasma frequency and scattering rate to the THz spectral region. The first-principles calculations provide two compelling implications: the conductivity nanocontrast can be induced by a small anisotropic strain, even less than 0.5%, due to an extreme strain sensitivity in ZrTe 5 ; A nanoscale topological phase transition is realized across some junctions induced by the strain, between strong topological insulators (TIs) and weak TIs/Dirac semimetals (DSMs).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗