Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space charge field”

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

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Separability of the Dirac equation in the STU black hole spacetime: Pairwise-equal charge case study

We study the separability of the Dirac equation in the background of four dimensional charged rotating asymptotically flat black hole solution of supergravity and the low energy string theory known as STU black hole. In particular, we analyze in detail the separability conditions in the pairwise equal charge STU black hole space-time. While in the latter case the minimally coupled Dirac equation is not separable, the introduction of a specific torsion term ensures the separability. The source of the torsion is the Kalb-Ramond field, which is an integral part of string theory, but further aspects of its properties and coupling to fermionic fields remain to be studied. To derive the torsion, two different approaches are used in conformally related frames, showing that the torsion is not unique. The correspondingly modified Dirac equations in the Einstein and string frames are shown to be separable. Furthermore, the massless radial and angular wave equations are examined; they show close similarity with corresponding equations for the standard Kerr background. A generalization of the Teukolsky equation for the pairwise equal case is conjectured. We also briefly analyze a technically sophisticated radial equation in the massive case.

Dirac equation↗
DOE OSTI2022

Improved Electrostatic Design of the Jefferson Lab 300 kV DC Photogun and the Minimization of Beam Deflection

An electron beam with high bunch charge and high repetition rate is required for electron cooling of the ion beam to achieve the high luminosity required for the proposed electron-ion colliders. An improved design of the 300 kV DC high voltage photogun at Jefferson Lab was incorporated toward overcoming the beam loss and space charge current limitation experienced in the original design. To reach the bunch charge goal of ~ few nC within 75 ps bunches, the existing DC high voltage photogun electrodes and anode-cathode gap were modified to increase the longitudinal electric field (Ez) at the photocathode. The anode-cathode gap was reduced to increase the Ez at the photocathode, and the anode aperture was spatially shifted with respect to the beamline longitudinal axis to minimize the beam deflection introduced by the geometric asymmetry of the inverted insulator photogun. The electrostatic design and beam dynamics simulations were performed to determine the required modification. Beam-based measurement from the modified gun confirmed the reduction of the beam deflection, which is presented in this contribution.

Mamun, M. A.↗
DOE OSTI2024

The science case for an intermediate energy advanced and novel accelerator linear collider facility

It is widely accepted that the next lepton collider beyond a Higgs factory would require center-of-mass energy of the order of up to 15 TeV. Since, given reasonable space and cost restrictions, conventional accelerator technology reaches its limits near this energy, high-gradient advanced acceleration concepts are attractive. Advanced and novel accelerators (ANAs) are leading candidates due to their ability to produce acceleration gradients on the order of 1–100 GV/m, leading to compact acceleration facilities. However, intermediate energy facilities (IEF) are required to test the critical technology elements on the way towards multi-TeV-class colliders. Here a science case for a 20–100 GeV center-of-mass energy ANA-based lepton collider that can be a candidate for an intermediate energy facility is presented. The IEF can provide numerous opportunities for high energy physics studies including precision Quantum Chromodynamics and Beyond the Standard Model physics measurements, investigation of charged particle interactions with extreme electromagnetic fields, and exploring muon and proton beam acceleration. Possible applications of this collider include the studies of γγ and electron beam-fixed target/beamdump collider designs. Thus, the goal of the proposed IEF is to both carry out particle physics measurements in the 20-100 GeV ranges as well as to serve as an ANA demonstrator facility.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Dynamic mode decomposition of nonequilibrium electron-phonon dynamics: accelerating the first-principles real-time Boltzmann equation

Abstract Nonequilibrium dynamics governed by electron–phonon ( e -ph) interactions plays a key role in electronic devices and spectroscopies and is central to understanding electronic excitations in materials. The real-time Boltzmann transport equation (rt-BTE) with collision processes computed from first principles can describe the coupled dynamics of electrons and atomic vibrations (phonons). Yet, a bottleneck of these simulations is the calculation of e –ph scattering integrals on dense momentum grids at each time step. Here we show a data-driven approach based on dynamic mode decomposition (DMD) that can accelerate the time propagation of the rt-BTE and identify dominant electronic processes. We apply this approach to two case studies, high-field charge transport and ultrafast excited electron relaxation. In both cases, simulating only a short time window of ~10% of the dynamics suffices to predict the dynamics from initial excitation to steady state using DMD extrapolation. Analysis of the momentum-space modes extracted from DMD sheds light on the microscopic mechanisms governing electron relaxation to a steady state or equilibrium. The combination of accuracy and efficiency makes our DMD-based method a valuable tool for investigating ultrafast dynamics in a wide range of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Testing isospin symmetry breaking in ab initio nuclear theory

In this work we present the first steps towards benchmarking isospin symmetry breaking in ab initio nuclear theory for calculations of superallowed Fermi β decay. Using the valence-space in-medium similarity renormalization group, we calculate b and c coefficients of the isobaric multiplet mass equation, starting from two different Hamiltonians constructed from chiral effective field theory. We compare results to experimental measurements for all T = 1 isobaric analog triplets of relevance to superallowed β decay for masses A = 10 to A = 74 and find an overall agreement within approximately 250 keV of experimental data for both b and c coefficients. A greater level of accuracy, however, is obtained by a phenomenological Skyrme interaction or a classical charged-sphere estimate. Lastly, we show that evolution of the valence-space operator does not meaningfully improve the quality of the coefficients with respect to experimental data, which indicates that higher-order many-body effects are likely not responsible for the observed discrepancies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Direct observation of elemental fluctuation and oxygen octahedral distortion-dependent charge distribution in high entropy oxides

The enhanced compositional flexibility to incorporate multiple-principal cations in high entropy oxides (HEOs) offers the opportunity to expand boundaries for accessible compositions and unconventional properties in oxides. Attractive functionalities have been reported in some bulk HEOs, which are attributed to the long-range compositional homogeneity, lattice distortion, and local chemical bonding characteristics in materials. However, the intricate details of local composition fluctuation, metal-oxygen bond distortion and covalency are difficult to visualize experimentally, especially on the atomic scale. Here, we study the atomic structure-chemical bonding-property correlations in a series of perovskite-HEOs utilizing the recently developed four-dimensional scanning transmission electron microscopy techniques which enables to determine the structure, chemical bonding, electric field, and charge density on the atomic scale. The existence of compositional fluctuations along with significant composition-dependent distortion of metal-oxygen bonds is observed. Consequently, distinct variations of metal-oxygen bonding covalency are shown by the real-space charge-density distribution maps with sub-ångström resolution. The observed atomic features not only provide a realistic picture of the local physico-chemistry of chemically complex HEOs but can also be directly correlated to their distinctive magneto-electronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Redesign of the Jefferson Lab -300 kV DC Photo-Gun for High Bunch Charge Operations

Production of high bunch charge beams for the Electron-Ion Collider (EIC) is a challenging task. High bunch charge (a few nC) electron beam studies at Jefferson Lab using an inverted insulator DC high voltage photo-gun showed evidence of space charge limitations starting at 0.3 nC, limiting the maximum delivered bunch charge to 0.7 nC for beam at -225 kV, 75 ps (FWHM) pulse width, and 1.64 mm (rms) laser spot size. The low extracted charge is due to the modest longitudinal electric field (Ez) at the photocathode leading to beam loss at the anode and downstream beam pipe. To reach the few nC high bunch charge goal, and to correct the beam deflection exerted by the non-symmetric nature of the inverted insulator photo-gun the existing photo-gun was modified. This contribution discusses the electrostatic design of the modified photo-gun obtained using CST Studio Suite’s electromagnetic field solver. Beam dynamics simulations performed using General Particle Tracer (GPT) with the resulting electrostatic field map obtained from the modified electrodes confirmed the validity of the new design.

Wijethunga, S.A.K.↗
DOE OSTI2025

Bridging the gap between molecules and materials in quantum chemistry with localized active spaces

The number of materials that “bridge the gap” between single molecules and extended solids, such as metal-organic frameworks and organic semiconductors, has been increasing. Consequently, there is a growing need for modeling approaches that effectively integrate the real-space molecular perspective employed by computational chemists and the reciprocal-space dispersive perspective employed by computational physicists. Here, we propose the localized active space (LAS) approach as a promising method to successfully bridge this gap. The LAS approach extends the active space concept from multiconfigurational methods such as complete active space self-consistent field theory to multiple molecular fragments via a product-form wave function ansatz. Here, we apply this method to solid state phenomena by treating each unit cell as a fragment with different sets of local quantum numbers (e.g., charge and excitation number). State interaction between these LAS states (LASSI) thus provides a comprehensive basis for the study of charge and energy transfer, meeting and surpassing the capabilities of single-reference fragmentation approaches such as constrained density functional theory (cDFT). Most centrally, we show how combining this LASSI approach with multiconfigurational pair-density functional theory (MC-PDFT) provides an elegant and efficient method to compute band structures that capture multiconfigurational character. We apply the LASSI band structure approach to the computation of band gaps in stretched hydrogen chain, polyacetylene, and bulk nickel oxide (NiO), finding good or excellent quantitative agreement with reference values in all cases. Additionally, we use the LAS basis in one-dimensional model systems to demonstrate its ability to treat difficult solid-state phenomena such as exciton transfer and excitation at p-n junctions.

method development↗
DOE OSTI2022

The feasibility of electron beam slip stacking in EIC Rapid Cycling Synchrotron

Announced in 2020, the Electron Ion Collider will be constructed at Brookhaven National Lab (BNL) by modifying the existing Relativistic Heavy Ion Collider (RHIC) facility. Among the alterations include the installation of an electron storage ring (ESR) and a Rapid Cycling Synchrotron (RCS) in the existing RHIC tunnel. The RCS will accelerate polarized electrons from 400MeV to 5GeV,10GeV and 18 GeV and transfer it to the ESR. The LINAC source produces 7nC per shot meanwhile the RCS requires 28nC per bunch of charge. Currently, it is planned to inject 4 bunches into adjacent RF buckets and then merge them into a single RF bucket using fields from several RF cavities. However, this requires a kicker with an exceptionally fast rise and fall time to inject into adjacent RF bucket without disturbing the existing circulating bunch and this requires harmonic RF kickers. Our study focuses on determining the feasibility of an alternate particle injection method that will utilize slip stacking, a method that accommodates the circulation of multiple beams with different momentum in the same azimuthal space of the beampipe, thus permitting the accumulation of charge from several shots into a single RF bucket. We find the momentum aperture required for slip stacking using current RCS beam parameters exceeds the RCS’s momentum deviation limit,dp/p = 1.5%. We determine the appropriate dimensions the beam must have to enable slip stacking and conduct a literature review of slip stacking schemes in other facilities. The literature review tentatively indicates that slip stacking has not been attempted on electron beams at a major accelerator facility.. The slip stacking approach has the potential to save the cost of the harmonic RF kicker and bunch merging cavities.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Beam physics research with the IOTA electron lens

The electron lens in the Fermilab Integrable Optics Test Accelerator (IOTA) will enable new research in nonlinear integrable optics, space-charge compensation, electron cooling, and the stability of intense beams. This research addresses scientific questions on high-brightness beams and operational challenges of high-power accelerators for nuclear and particle physics. We review the roles that electron lenses play in this field and the physical principles behind their applications. The design criteria and specifications for the IOTA storage ring and electron lens are then discussed. We conclude with a description of the components of the apparatus.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Mapping hidden space-charge distributions across crystalline metal oxide/group IV semiconductor interfaces

The electronic structures of semiconducting heterojunctions are critically dependent on composition including the presence and concentrations of dopants, both intended and unintended. Dopant profiles in the interfacial region can have major effects on band energies which in turn drive transport properties. Here, in this study, we use core-level photoelectron line shapes excited with hard x rays to extract information about electric fields resulting from internal charge transfer in epitaxial La 0.03 Sr 0.97 Zr x Ti 1–x O 3 /Ge(001) (0.1≤x≤0.7) heterostructures. Experiments were carried out for heterojunctions involving both n- and p-type Ge substrates. These heterojunctions were not amenable to electronic characterization of all regions by transport measurements because the doped substrates act as electrical shunts, precluding probing the more resistive films and masking interface conductivity. However, the core-level line shapes were found to be a rich source of information on built-in potentials that exist throughout the heterostructure, and yielded valuable insight into the impact of band bending on band alignment at the buried interfaces. The electronic effects expected for Ge with uniform n- and p-type doping are eclipsed by those of unintended oxygen dopants in the Ge near the interface. This study illustrates the power of hard x-ray photoemission spectroscopy and related modeling to determine electronic structure in material systems for which insight from traditional transport measurements is limited.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Degradation-related defect level in weathered silicon heterojunction modules characterized by deep level transient spectroscopy

Commercial silicon heterojunction photovoltaic modules, known as amorphous-silicon-based heterojunction with intrinsic thin-film layer (HIT) modules, show average degradation after 10 years in the field. HIT modules weathered outdoors in Colorado and Florida display mostly uniform decreases in intensity when mapped with photoluminescence (PL) imaging compared to a control module. Flash-table-based current-voltage curves show that degradation is dominated by voltage loss. Samples are cored from each of the modules, and deep level transient spectroscopy (DLTS) detects three electron-trap defect states in all modules with activation energies of electron emission from the defects of 0.07, 0.16, and 0.50 eV. DLTS measurements on the weathered modules show an additional deep-level, electron-trap defect state with an activation energy of 0.51 eV and a trap density of approximately 10 12 cm –3 . The capture rate is measured using varying short filling pulse times, and the resulting capture cross section is estimated to be 1.1x10 –16 cm 2 . The development of the weathering-related defect level correlates to decreases in carrier lifetime, PL intensity, and module voltage. Various depths of the space charge region are probed with increments in applied reverse bias and filling-pulse bias. Furthermore, this DLTS depth profiling shows a trend of trap density increasing with less applied reverse bias, suggesting that the weathering-related defect increases carrier recombination toward the interface between the bulk silicon wafer and the junction-forming amorphous-silicon passivation layers.

14 SOLAR ENERGY↗
DOE OSTI2020

Two-fluid and kinetic transport physics of Kelvin–Helmholtz instabilities in nonuniform low-beta plasmas

Hall-magnetohydrodynamic (Hall-MHD) theory, two-fluid simulations, and kinetic simulations are used in this work to investigate the cross-field transport properties of Kelvin–Helmholtz instabilities in nonuniform low-beta collisionless plasmas. Hall-MHD analysis shows how the linear properties of the instability are modified by density gradients and magnetization. High-order accurate two-fluid and kinetic simulations, with complete dynamics of finite-mass electrons and ions, are applied to a suite of parameter cases to systematically assess the effects of diamagnetic drift, magnetization, charge separation, and finite Larmor motion. Initialization of exact two-species kinetic equilibria facilitates the study of isolated physical effects and enables detailed cross-comparisons between two-fluid and kinetic simulations, including for cases where ion gyroradii are comparable to gradient scale lengths. For nonuniform plasmas with significant space charge, the results of two-fluid and kinetic simulations are found to disagree with Hall-MHD predictions. Kelvin–Helmholtz instability growth rates, per unit shear, are shown to be smaller when ion diamagnetic drift and $E \times B$ drift are parallel and larger when the two drifts are antiparallel. The effect is attributed to polarization drift in the shear layer, which leads to redistribution of charge, alters the electric field that drives plasma advection, and consequently modifies growth rates. Instability-induced mass transport for different parameters is characterized in terms of the flux across the shear layer and a simplified diffusion model. Distribution functions from kinetic simulations are shown to deviate substantially from Maxwellian reconstructions, indicating the importance of kinetic physics during the nonlinear phase of the instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Direct deflection of Millicharged radiation

Millicharged particles are generic in theories of dark sectors. A cosmic or local abundance of them may be produced by the early universe, stellar environments, or the decay or annihilation of dark matter/dark energy. Furthermore, if such particles are light, these production channels result in a background of millicharged radiation. We show that light-shining-through-wall experiments employing superconducting RF cavities can also be used as “direct deflection” experiments to search for this relativistic background. The millicharged plasma is first subjected to an oscillating electromagnetic field of a driven cavity, which causes charge separation in the form of charge and current perturbations. In turn, these perturbations can propagate outwards and resonantly excite electromagnetic fields in a well-shielded cavity placed nearby, enabling detection. We estimate that future versions of the existing Dark SRF experiment can probe orders of magnitude of currently unexplored parameter space, including millicharges produced from the Sun, the cosmic neutrino background, or other mechanisms that generate a thermal abundance with energy density as small as ~ 10$^{−4}$ that of the cosmic microwave background.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗
DOE OSTI2021

Comparison of Feschenko BSM and Fast Faraday Cup with Low Energy Ion Beams

A comparison between Fast Faraday Cup and Feschenko longitudinal bunch shape detectors was recently performed at HELIAC Advanced Demonstrator beamline at GSI. Feschenko bunch shape monitor (BSM) uses the time to space conversion by means of secondary electrons emitted from a wire correlated to a rf deflector, while the fast Faraday cup (FFC) measures the deposited charge in a cup geometry matched to 50 Ohm. The FFC design aims to minimize the bunch shape dilution due to field polarization and secondary electrons produced on irradiation. An He¹⁺ with 100 µA average current and 1.4 MeV/u kinetic energy is utilized for this comparison. A buncher upstream of the detectors was operated to focus the beam longitudinally. The results are discussed in this contribution.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Repulsive black holes and higher-derivatives

In two-derivative theories of gravity coupled to matter, charged black holes are self-attractive at large distances, with the force vanishing at zero temperature. However, in the presence of massless scalar fields and four-derivative corrections, zero-temperature black holes no longer need to obey the no-force condition. In this paper, we show how to calculate the long-range force between such black holes. We develop an efficient method for computing the higher-derivative corrections to the scalar charges when the theory has a shift symmetry, and compute the resulting force in a variety of examples. We find that higher-derivative corrected black holes may be self-attractive or self-repulsive, depending on the value of the Wilson coefficients and the VEVs of scalar moduli. Indeed, we find black hole solutions which are both superextremal and self-attractive. Furthermore, we present examples where no choice of higher-derivative coefficients allows for self-repulsive black hole states in all directions in charge space. This suggests that, unlike the Weak Gravity Conjecture, which may be satisfied by the black hole spectrum alone, the Repulsive Force Conjecture requires additional constraints on the spectrum of charged particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗