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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Space charge fields in azimuthally symmetric beams: Integrated Green's function approach

Electromagnetic fields induced by the space charge in relativistic beams play an important role in Accelerator Physics. They lead to emittance growth, slice energy change, and the microbunching instability. Typically, these effects are modeled numerically since simple description exists only inthe limits of large- or small-scale current variations. In this paper we consider an axially symmetric charged beam inside a round pipe and find the solution of the space charge problem that is valid in the full range of current variations. We express the solution for the field components in terms of Green’s functions, which are fully determined by just a single function. We then find that this function is an on-axis potential from a charged disk in a round pipe, with transverse charge density ρ⊥(r), and it has a compact analythical expression. We finally provide an integrated Green’s function based approach for efficient numerical evaluation in the case when the transverse charge density stays the same along the beam

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Ultrafast quantum dynamics driven by the strong space-charge field of a relativistic electron beam

Here, we illustrate how the Coulomb field of a highly relativistic electron beam can be shaped into a broadband pulse suitable for driving ultrafast and strong-field physics. In contrast to a solid-state laser, the Coulomb field creates a pulse which can be intrinsically synchronized with an x-ray free electron laser (XFEL), can have a cutoff frequency which is broadly tunable from THz to EUV, and which acts on target systems as a “half-cycle” impulse. Explicit examples are presented to emphasize how the unique features of this excitation can be a tool for novel science at XFEL facilities like the LCLS.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Identification of Recombination Losses in CdSe/CdTe Solar Cells from Spectroscopic and Microscopic Time-Resolved Photoluminescence

Due to the lowest-cost and best reliability, CdTe solar cells are the leading thin-film photovoltaic technology. Increasing open-circuit voltage by reducing recombination represents the most promising path toward further improvements. Analysis is needed to identify limitations that cause efficiency losses. To achieve this goal for Cu-doped CdSe/CdTe solar cells, time-resolved spectroscopy and microscopy are developed and applied. Recombination lifetimes and radiative efficiency identify that defect-mediated recombination is the dominant voltage loss mechanism. When carrier lifetimes are averaged over many crystalline grains, they increase from 180 to 430 ns when Al2O3 is applied to the back contact. The quasi-Fermi-level splitting correspondingly increases from 880–905 to 906–931 mV, indicating a pathway to overcome the long-standing 900 mV voltage limitation. However, the dominant recombination losses are attributed to the absorber bulk. From microscopic carrier lifetime measurements, it is identified that space charge fields due to charged grain boundaries (GBs) lead to recombination in the CdTe absorber bulk. At high injection, GB space charge fields are screened, but that occurs above 1 Sun excitation conditions. Alloying with selenium in the near-interface CdSeTe absorber region reduces GB losses and is identified as one of the factors leading to high radiative and power conversion efficiency.

14 SOLAR ENERGY↗
DOE OSTI2024

Gas Stopper Developments for Improved Purity and Intensity of Low-Energy, Rare Isotope Ion Beams (Final Technical Report)

This final technical report summarizes the work of the Michigan State University (MSU) team supported by grant # DE-SC0021423 awarded by the Office of Nuclear Physics, Department of Energy. Objectives: The successful fulfillment of the FRIB science mission hinges on ensuring the availability of fast, stopped, and reaccelerated beams consisting of rare isotopes. This project's research and development focus was dedicated to supporting the advancement and creation of a cutting-edge linear gas stopper. The primary aim is to efficiently convert the high-intensity fast beams of rare isotopes provided by FRIB into high-quality, low-energy beams. These beams are essential for conducting stopped beam experiments or for subsequent reacceleration. The overarching goal is to advance technology, aiming to increase the beam rate capability of the linear gas stopper for medium-to-heavy-mass rare isotopes by more than tenfold compared to the currently most effective gas stopper in operation, and to improve the purity of the delivered beams. Project Description: The existing technology employed in gas stopping devices designed for low-energy, rare-isotope beams presents limitations in both the purity of extracted beams and the intensities of injected beams. These limitations are incompatible with the requirements of the recently commissioned rare isotope beam facility, FRIB. Our research and development efforts, aligned with the previously outlined objectives, focused on addressing the most critical aspects for enhancing beam-rate capability and purity. Specifically, advanced particle-in-cell simulations were developed and integrated into a simulation pipeline to explore the efficacy of multi-layer RF carpets on increasing ion transport efficiency with high incoming beam rates that generate space charge fields which can limit it. We also explored the possibility of using a collision-induced-dissociation (CID) gas cell to break up molecular contaminant ions that are generated during the stopping process. A prototype CID gas cell was constructed and tested with beams from an offline ion source, validating the concept with the successful demonstration of breaking of molecular ions. The outcome of this research enabled the formulation a conceptual design for a next-generation linear gas stopping device specifically tailored for FRIB. This device is envisioned to deliver rare-isotope-ion beams at a rate of 10 8 particles per second or higher, accompanied by advancements in purity. Methods employed: This project leverages advancements in technologies initially designed for the Advanced Cryogenic Gas Stopper (ACGS), the current state-of-the-art linear gas stopper, through the use of new simulations and beam purification via collision-induced-dissociation. The methods include: 1. Development of a prototype low-energy, low-pressure CID gas stopper. This prototype features a thin, approximately 20 nm, Si 3 N 4 entrance window designed for dissociating stable and rare-isotope molecular ions. The goal is to enhance beam purification and overall efficiency. 2. Creation of Particle-in-cell (PIC) simulations to assess the advantages of multi-layer RF carpets and multi-point extraction for ion transport efficiency. These simulations rely on the 3DCylPIC package, specifically designed for studying devices of this nature. The goal is to quantify and mitigate ion transport losses due to space charge generated in the stopping process of large numbers of ions. 3. Perform ion transport simulations across an RF carpet using an 8-phase travelling wave and evaluate its performance. Compared to the 4-phase RF carpets used in ACGS, the 8-phase carpets will double the wavelength of the generated traveling wave allowing for larger maximum RF amplitudes that could result in improved ion transport efficiency for high-intensity incoming beams when large space charge fields are present. Impact: Tackling the primary challenges associated with transforming high-energy projectile fragment beams into low-energy beams—specifically, addressing efficiency and purity—holds significant promise for advancing FRIB science. This advancement will particularly impact precision mass measurements, laser spectroscopy of short-lived nuclei, and studies in astrophysics and nuclear reactions using reaccelerated beams. These domains play a crucial role in addressing key questions outlined in the 2023 NSAC long-range plan, spanning nuclear structure, nuclear astrophysics, and fundamental symmetries. Additionally, they contribute to addressing 10 out of the 17 benchmarks identified by the NRC RISAC. The development of a next-generation gas stopping device capable of delivering low-energy, rare-isotope beams at a rate of 10 8 particles per second, or more, with high purity holds the potential to unlock experiments that would otherwise be unfeasible. Furthermore, it is expected to reduce the time required for experiments at FRIB, thereby maximizing scientific output. The research and development activities performed as part of this project bolstered essential competencies at FRIB in beam physics and ion source technologies, provided valuable training opportunities for junior scientists.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Poisson Equation for a (General) Homogeneous d-Dimensional Ellipsoid with Applications to Beam Envelope Tracking

This note describes the solution of the free-space Poisson equation in the interior of a $d$-dimensional homogeneous ellipsoid, and the associated space charge fields. An explicit formula (\ref{Sformula}) is provided that relates the $d\times d$ matrix describing the space charge (quadratic) potential to the $d\times d$ covariance matrix of the ellipsoid. For the cases $d=2$ and $d=3$, this result is used to determine the linear map corresponding to a space charge kick, that may be used to push the beam $6\times 6$ covariance matrix during envelope tracking. The treatment of upright ellipsoids for $d=2$ and $d=3$ is well-represented in the literature. However, the approach taken here emphasizes a general ellipsoid with arbitrary correlations in any dimension. The Appendix provides a general solution of the free-space Poisson equation in dimension $d$ for a source distribution with ellipsoidal symmetry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Space-Charge Aberrations in Single-Shot Time-Resolved Transmission Electron Microscopy

In this paper, we discuss effects of space-charge fields on imaging performance in a single-shot time-resolved transmission electron microscope. Using a Green’s function perturbation method, we derive analytical estimates for the effects of space-charge nonlinearity on the image formation process and the associated aberration coefficients. The dependence of these coefficients on the initial beam phase space distribution is elucidated. The results are validated by particle tracking simulations and provide fundamental scaling laws for the trade-off between temporal and spatial resolutions in single-shot time-resolved transmission electron microscopy.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Scanning ultrafast electron microscopy reveals photovoltage dynamics at a deeply buried $p^- \mathrm{Si/SiO_2}$ interface

The understanding and control of charge carrier interactions with defects at buried insulator/semiconductor interfaces is essential for achieving optimum performance in modern electronics. In this work, we report on the use of scanning ultrafast electron microscopy (SUEM) to remotely probe the dynamics of excited carriers at a Si surface buried below a thick thermal oxide. Our measurements illustrate a previously unidentified SUEM contrast mechanism, whereby optical modulation of the space-charge field in the semiconductor modulates the electric field in the thick oxide, thus affecting its secondary electron yield. By analyzing the SUEM contrast as a function of time and laser fluence we demonstrate the diffusion mediated capture of excited carriers by interfacial traps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Photorefraction-induced Bragg scattering in cryogenic lithium niobate ring resonators

We report intracavity Bragg scattering induced by the photorefractive (PR) effect in high-Q lithium niobate ring resonators at cryogenic temperatures. We show that when a cavity mode is strongly excited, the PR effect imprints a long-lived periodic space-charge field. This residual field in turn creates a refractive index modulation pattern that dramatically enhances the back scattering of an incoming probe light, and results in selective and reconfigurable mode splittings. This PR-induced Bragg scattering effect, despite being undesired for many applications, could be utilized to enable optically programmable photonic components.

42 ENGINEERING↗
DOE OSTI2024

EIC hadron ring longitudinal stability simulations at injection

At the EIC injection energy (23.8 GeV for protons), the hadron beam experiences a relatively strong space charge field. The corresponding wakefield, in conjunction with the wakes from the RF cavities (i.e. the geometrical wakes), may become the cause of single-bunch instabilities, both longitudinal and transverse. In this report, we focus on the longitudinal stability aspect of this picture, using a simplified model for the ring.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Broadband hard X-ray attosecond pulses from extremely chirped electron beams

Attosecond pulses from free-electron lasers have opened the doors to atomic site-specific studies of bound electronic dynamics on their natural, sub-femtosecond timescales. Key to their success has been electron beam shaping techniques enabling the generation of sub-femtosecond current spikes with peak currents on the order of 10 kA. We demonstrate in an RF linac the generation of current spikes with extreme chirps on the order of 350 MeV/micron, directly competitive with the chirps expected from beam-driven plasma wakefield accelerators. Leveraging chirp-taper compensation, we use these highly chirped beams to generate hard X-ray attosecond pulses with bandwidths exceeding 30 eV, a factor of two beyond previous demonstrations. We simultaneously present the first explicit experimental evidence of chirp-taper compensation in an attosecond XFEL, finding that optimal tapering improves the bandwidth and pulse energy by factors of two and five, respectively, for our conditions. In addition to the immediate utility of such broadband hard X-ray pulses, electron beams with such extreme chirps can be utilized for unique new experimental modalities by performing further compression after the undulators. Such post-lasing compression can enable subsequent superradiant light emission at longer wavelengths, or direct excitation of quantum systems with the beam's intense space-charge field for unique attosecond pump-probe possibilities.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Coherent transverse instability of bunched beams in circular accelerators

Coherent transverse oscillations of a bunch in a ring accelerator are considered. Three different models of the bunch are used: a hollow bunch in a square potential well, a square bunch in a parabolic potential well, and a parabolic bunch in a parabolic well. The bunch space-charge field, wakefield, and chromaticity are taken into account in all the cases. The bunch eigenmodes and various instabilities are studied. The results are compared and analyzed to establish the range of applicability of the models. Role of the Landau damping due to dependence of the SC field on the longitudinal coordinate is discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Plasma-induced collective effects in muon ionization cooling

Muon ionization cooling is a critical process that determines the achievable beam brightness and overall performance of a muon collider. In this study, we investigate collective effects induced by high-intensity muon beams inside the absorber, effects that have not been included in previous ionization cooling simulations. We show that beam-induced plasma can neutralize the beam space-charge field, while the self-generated azimuthal magnetic field can introduce additional beam focusing through a z-pinch-like mechanism. We also demonstrate that sufficiently high plasma densities may modify the stopping-power rate. The implications of these plasma-induced phenomena for ionization cooling performance are discussed.

Yonehara, Katsuya [Fermilab] (ORCID:00000002544041↗
DOE OSTI2023

Dimensional analysis on microscale gas breakdown with electric field nonuniformity and positive space charge effects

A dimensional method was employed to evaluate the microscale gas breakdown characteristics at atmospheric pressure, resulting in a universal breakdown curve applicable to different types of gases (e.g., Ar, Xe, Ne, and N2). As the gap distance decreases, the breakdown mode transitions from ion-induced secondary electron emission to the field emission regime. In the field emission regime, the positive space charge effect becomes more significant. We discovered that incorporating the positive space charge effect in the field emission regime can be achieved by modifying the local electric field enhancement factor β. Consequently, we propose an effective electric field enhancement factor, βeff, which scales linearly with β, to accurately reproduce the breakdown curve while considering the positive space charge effect. This proposed approach significantly simplifies the numerical model. Additionally, we examined the effects of gas pressure, gap distance, cathode properties (e.g., work function and secondary electron emission coefficient), and electric field nonuniformity.

Physics↗
DOE OSTI2026

Calculating the space-charge-limited current density for nonplanar geometries by simulating the charge-free electric field

Calculating the space-charge-limited-current density (SCLCD) for a complicated diode geometry often requires computationally expensive particle-in-cell (PIC) codes. Here, this paper addresses this issue by using the charge-free electric field $E_0$ calculated using COMSOL Multiphysics to determine local and global SCLCD. The SCLCD obtained by using the surface average of $|E_0|^2$ on the cathode recovers theoretical results for one-dimensional (1D) planar, cylindrical, and tip-to-tip geometries in appropriate limits. We next compared tip-to-tip calculations with the SCLCD obtained using the PIC code Empire. The SCLCD calculated using COMSOL agreed well with Empire for flatter 1D tip-to-tip geometries and diverged with increasing sharpness. Physically, Empire predicts lower SCLCD than COMSOL because the electrons spread due to concentrated space-charge at the tip, whereas theory assumes that the electrons follow the charge-free electric field lines. We further assess the behavior of the SCLCD for tips protruding from the centers of flat, circular plates of various areas. Larger plate areas with constant tip size recover the 1D planar SCLCD globally and 1D tip-to-tip SCLCD locally, while reducing the difference between Empire and COMSOL calculations since larger plates capture more of the emitted electrons, reducing SCLCD suppression due to beam spreading. These results show that charge-free electric field simulations can be used to determine the SCLCD without needing to simulate particle dynamics in PIC.

Wright, Jack K. [Purdue Univ., West Lafayette, IN ↗
DOE OSTI2026

pyEF: A Python Framework for QM and QM/MM Atom-Wise Electric Field Analysis

We introduce pyEF, a software package for computing molecular electric fields, electrostatic interaction energies, and electrostatic potentials from quantum mechanical (QM) atom-centered multipole expansions with atom-wise decomposable contributions. We demonstrate the computational efficiency and accuracy of this QM-derived electric field evaluation tool through several tests. To assess the influence of the underlying QM method and charge partitioning scheme on these electrostatic quantities, we analyze over 250 configurations of an acetone solute molecule in five solvents of variable polarity. We find that electric field calculations are highly sensitive to the choice of charge partitioning method. Even among real-space charge schemes, acetone Stark tuning rates differ by up to a factor of 2. Benchmarking computed solvent dipole moments against experimental bulk values, we conclude that the CM5, ADCH, and Hirshfeld-I charge schemes most reliably capture solvent electrostatics and therefore provide a more faithful foundation for computing electric fields. When constructed from these real-space charges, electric fields are nearly insensitive to basis set size and monotonically increase in magnitude with higher Fock exchange. We also demonstrate efficient convergence of QM electrostatics when more distant molecules are represented solely by MM point charges, reducing computational overhead. Leveraging these findings, we demonstrate the use of pyEF to deduce environmental effects on a transition metal complex from a Ga 4 L 6 12– nanocage and quantify the dominant role of organic linkers in orchestrating electrostatic preorganization.

electric fields↗
DOE OSTI2023

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗