Search NASA⌕ Search

SEARCH · Search NASA

Results for “Relativistic drifting particles”

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

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↗
NASA NTRS1990

Shock-drift particle acceleration in superluminal shocks - A model for hot spots in extragalactic radio sources

Shock-drift acceleration at relativistic shock fronts is investigated using a fully relativistic treatment of both the microphysics of the shock-drift acceleration and the macrophysics of the shock front. By explicitly tracing particle trajectories across shocks, it is shown how the adiabatic invariance of a particle's magnetic moment breaks down as the upstream shock speed becomes relativistic, and is recovered at subrelativistic velocities. These calculations enable the mean increase in energy of a particle which encounters the shock with a given pitch angle to be calculated. The results are used to construct the downstream electron distribution function in terms of the incident distribution function and the bulk properties of the shock. The synchrotron emissivity of the transmitted distribution is calculated, and it is demonstrated that amplification factors are easily obtained which are more than adequate to explain the observed constrasts in surface brightness between jets and hot spots.

Begelman, Mitchell C.↗
NASA NTRS2016

Outer Radiation Belt Dropout Dynamics Following the Arrival of Two Interplanetary Coronal Mass Ejections

Magnetopause shadowing and wave-particle interactions are recognized as the two primary mechanisms for losses of electrons from the outer radiation belt. We investigate these mechanisms, sing satellite observations both in interplanetary space and within the magnetosphere and particle drift modeling. Two interplanetary shocks sheaths impinged upon the magnetopause causing a relativistic electron flux dropout. The magnetic cloud (C) and interplanetary structure sunward of the MC had primarily northward magnetic field, perhaps leading to a concomitant lack of substorm activity and a 10 day long quiescent period. The arrival of two shocks caused an unusual electron flux dropout. Test-particle simulations have shown 2 to 5 MeV energy, equatorially mirroring electrons with initial values of L 5.5can be lost to the magnetosheath via magnetopause shadowing alone. For electron losses at lower L-shells, coherent chorus wave-driven pitch angle scattering and ULF wave-driven radial transport have been shownto be viable mechanisms.

Alves, L. R.↗
NASA NTRS2014

Binary Orbits as the Driver of Gamma-Ray Emission and Mass Ejection in Classical Novae

Classical novae are the most common astrophysical thermonuclear explosions, occurring on the surfaces of white dwarf stars accreting gas from companions in binary star systems. Novae typically expel about 10 (sup -4) solar masses of material at velocities exceeding 1,000 kilometers per second.However, the mechanism of mass ejection in novae is poorly understood, and could be dominated by the impulsive flash of thermonuclear energy, prolonged optically thick winds or binary interaction with the nova envelope. Classical novae are now routinely detected at giga-electronvolt gamma-ray wavelengths, suggesting that relativistic particles are accelerated by strong shocks in the ejecta. Here we report high-resolution radio imaging of the gamma-ray-emitting nova V959 Mon. We find that its ejecta were shaped by the motion of the binary system: some gas was expelled rapidly along the poles as a wind from the white dwarf, while denser material drifted out along the equatorial plane, propelled by orbital motion..At the interface between the equatorial and polar regions, we observe synchrotron emission indicative of shocks and relativistic particle acceleration, thereby pinpointing the location of gamma-ray production. Binary shaping of the nova ejecta and associated internal shocks are expected to be widespread among novae, explaining why many novae are gamma-ray emitters.

binary orbits as the driver of gamma-ray↗
DOE OSTI2024

Electron and Proton Energization in 3D Reconnecting Current Sheets in Semirelativistic Plasma with Guide Magnetic Field

Abstract Using 3D particle-in-cell simulation, we characterize energy conversion, as a function of guide magnetic field, in a thin current sheet in semirelativistic plasma, with relativistic electrons and subrelativistic protons. There, magnetic reconnection, the drift-kink instability (DKI), and the flux-rope kink instability all compete and interact in their nonlinear stages to convert magnetic energy to plasma energy. We compare fully 3D simulations with 2D in two different planes to isolate reconnection and DKI effects. In zero guide field, these processes yield distinct energy conversion signatures: ions gain more energy than electrons in 2Dxy(reconnection), while the opposite is true in 2Dyz(DKI), and the 3D result falls in between. The flux-rope instability, which occurs only in 3D, allows more magnetic energy to be released than in 2D, but the rate of energy conversion in 3D tends to be lower. Increasing the guide magnetic field strongly suppresses DKI, and in all cases slows and reduces the overall amount of energy conversion; it also favors electron energization through a process by which energy is first stored in the motional electric field of flux ropes before energizing particles. Understanding the evolution of the energy partition thus provides insight into the role of various plasma processes, and is important for modeling radiation from astrophysical sources such as accreting black holes and their jets.

Astronomy & Astrophysics↗
NASA NTRS1995

Simulation of proton radiation belt formation during the March 24, 1991 SSC

The rapid formation of a new proton radiation belt at L approximately = 2.5 following the March 24, 1991 Storm Sudden Commencement (SSC) observed at the Combined Release and Radiation Effects Satellite (CRRES) satellite is modeled using a relativistic guiding center test particle code. The SSC is modeled by a bipolar electric field and associated compression and relaxation in the magnetic field, superimposed on a dipole magnetic field. The source population consists of both solar and trapped inner zone protons. The simulations show that while both populations contribute to drift echoes in the 20-80 MeV range, primary conditions is from the solar protons. Proton acceleration by the SSC differs from relativistic electron acceleration in that different source populations contribute and nonrelativistic conservation of the first adiabatic invariation leads to greater energization of protons for a given decrease in L. Model drift echoes and flux distribution in L at the time of injection compare well with CRRES observations.

Hudson, M. K.↗
NASA NTRS2021

Reconnection and particle acceleration in 3D current sheet evolution in moderately-magnetized astrophysical pair plasma

Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially 2D) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in 3D, we conduct a comprehensive particle-in-cell simulation study comparing 2D and 3D evolution of long, thin current sheets in moderately-magnetized, collisionless, relativistically-hot electron- positron plasma, and find dramatic differences. We first systematically characterize this process in 2D, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths, and system sizes. We then show that 3D simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3D current-sheet evolution is not always fundamentally classical reconnection with perturbing 3D effects, but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details, and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in 2D. Intriguingly, nonthermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in 2D and 3D. Though the variety of underlying current-sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.

Gregory R. Werner↗
NASA NTRS1972

Anomalous type III bursts explained in terms of intermittent solar wind streams observed by Pioneer spacecraft.

Occasional solar radio bursts drift downward in frequency (like type III bursts), but the drift slows or ceases or it may even reverse. These anomalous forms are called L-bursts, U-bursts, etc., and are observed as low as 1 MHz by spacecraft. The conventional explanation involves guidance of relativistic exciting particles by solar magnetic fields that do not have the spiral shape expected in undisturbed conditions. A reasonable alternative explanation is possible and, indeed, it follows directly from the inferences which have been drawn from the Stanford dual-frequency experiment flown on four Pioneer spacecraft. Analysis of temporal variations in these data leads to the conclusion that narrow streams of enhanced concentration have lifetimes which are typically about one day and the streams seem to be intermittently ejected by long-enduring sources at the sun. Each enhancement is thought to lie along a spiral and to have a beginning and an end. It may happen that relativistic exciters encounter the tail of such a discrete streamer; the concentration encountered by such exciters would initially decline but, at the streamer tail, it would level off and perhaps increase. The plasma frequency versus time (the burst frequency) would then cease to drop in the typical type III profile and would level off or increase (as in a U-burst).

Croft, T. A.↗
NASA NTRS1980

Ionospheric density enhancement during relativistic electron precipitation

The temporal evolution of the ionospheric density enhancement produced by a widespread relativistic electron precipitation (REP) has been observed with the Chatanika Radar. The REP was associated with a substorm particle energization event, and both the ionospheric absorption and density perturbation exhibited an approximately 90 min periodicity associated with the particles' longitudinal drift. A 80-keV characteristic energy for the precipitating electrons is deduced from ground-based and satellite data. At the maximum of the event, electrons deposited approximately 50 ergs/sq cm per sec in the ionosphere, producing a peak density of 500,000/cu cm at 89 km altitude. At that time the radar observed densities greater than 100,000/cu cm between 70 km and 110 km altitude and riometer absorption at 30 MHz was approximately 12 db.

Foster, J. C.↗
DOE OSTI2021

Reconnection and particle acceleration in three-dimensional current sheet evolution in moderately magnetized astrophysical pair plasma

Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially two-dimensional) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in three dimensions, we conduct a comprehensive particle-in-cell simulation study comparing two- and three-dimensional evolution of long, thin current sheets in moderately magnetized, collisionless, relativistically hot electron–positron plasma, and find dramatic differences. We first systematically characterize this process in two dimensions, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths and system sizes. We then show that three-dimensional (3-D) simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3-D current sheet evolution is not always fundamentally classical reconnection with perturbing 3-D effects but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in two dimensions. Intriguingly, non-thermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in two and three dimensions. Although the variety of underlying current sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.

Physics↗
NASA NTRS1993

Simulation of the prompt energization and transport of radiation belt particles during the March 24, 1991 SSC

We model the rapid (about 1 min) formation of a new electron radiation belt at L about or = 2.5 that resulted from the Storm Sudden Commencement (SSC) of March 24, 1991 as observed by the Combined Release and Radiation Effects Satellite (CRRES) satellite. Guided by the observed electric and magnetic fields, we represent the time-dependent magnetospheric electric field during the SSC by an asymmetric bipolar pulse that is associated with the compression and relaxation of the Earth's magnetic field. We follow the electrons using a relativistic guiding center code. The test-particle simulations show that electrons with energies of a few MeV at L greater than 6 were energized up to 40 MeV and transported to L about or = 2.5 during a fraction of their drift period. The energization process conserves the first adiabatic invariant and is enhanced due to resonance of the electron drift motion with the time-varying electric field. Our simulation results, with an initial W(exp -8) energy flux spectra, reproduce the observed electron drift echoes and show that the interplanetary shock impacted the magnetosphere between 1500 and 1800 MLT.

Li, Xinlin↗
DOE OSTI2025

Kinetic Simulations of the Kruskal–Schwarzschild Instability in Accelerating Striped Outflows: Dynamics and Energy Dissipation

Astrophysical relativistic outflows are launched as Poynting-flux dominated, yet the mechanism governing efficient magnetic dissipation, which powers the observed emission, is still poorly understood. We study magnetic energy dissipation in relativistic "striped" jets, which host current sheets separating magnetically dominated regions with opposite field polarity. The effective gravity force g in the rest frame of accelerating jets drives the Kruskal–Schwarzschild instability (KSI), a magnetic analog of the Rayleigh–Taylor instability. By means of 2D and 3D particle-in-cell simulations, we study the linear and nonlinear evolution of the KSI. The linear stage is well described by linear stability analysis. The nonlinear stages of the KSI generate thin (skin-depth-thick) current layers, with length comparable to the dominant KSI wavelength. There, the relativistic drift-kink mode and the tearing mode drive efficient magnetic dissipation. The dissipation rate can be cast as an increase in the effective width Δ eff of the dissipative region, which follows dΔ eff /dt ≃ 0.05 $\sqrt{Δ_{eff} g}$. Our results have important implications for the location of the dissipation region in gamma-ray burst and active galactic nuclei jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Toroidal modeling of Alfvén eigenmodes excited by runaway electrons in DIII-D and ITER

The non-perturbative MHD-kinetic hybrid code MARS-K (Liu et al 2014 Phys. Plasmas 21 056105) is updated to include relativistic effects for kinetic fast particles, enabling the code to model excitation of Alfvén eigenmodes (AEs) by runaway electrons (REs) in post-disruption tokamak plasmas. Applying the updated code to RE beams in both DIII-D and ITER, a zoo of AE modes triggered by trapped REs due to precessional drift-kinetic resonances is computed while scanning the RE energy. At fixed RE energy, multiple unstable roots are also excited. These AE modes possess radially different eigenmode structures, ranging from global modes to core-localized ones. The computed mode frequency is in the Alfvén frequency range, increasing with the assumed RE energy in a staircase fashion and quantitatively matching the experimental measurement (in DIII-D). At the (more relevant) high-frequency range (above 1 MHz), the modeled eigenmodes are identified as compressional AEs (CAEs) in DIII-D and a mixture of CAE and shear Alfvén waves in ITER.

Alfvén eigenmodes↗
DOE OSTI2025

A Particle-in-cell Method for Plasmas with A Generalized Momentum Formulation, Part III: A family of Gauge Conserving Methods

In this paper, we introduce a new family of spatially co-located field solvers for particle-in-cell applications which evolve the potential formulation of Maxwell’s equations under the Lorenz gauge. Our recent work [2] introduced the concept of time-consistency, which connects charge conservation to the preservation of the gauge at the semi-discrete level. It will be shown that there exists a large family of time discretizations which satisfy this property. Additionally, it will be further shown that for large classes of time marching methods, the satisfaction of the gauge condition automatically implies the satisfaction of Gauss’s law for electricity, with the potential formulation ensuring that that Gauss’s law for magnetism is satisfied by definition. We focus on popular time marching methods including centered differences, backward differences, and diagonally-implicit Runge-Kutta methods, which are coupled to a spectral discretization in space. We demonstrate the theory by testing the methods on a relativistic Weibel instability and a drifting cloud of electrons.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Electron acceleration to relativistic energies by traveling interplanetary shocks

Between 1978 and 1985, the Pioneer 10/11 measurements identified 11 traveling interplanetary (IP) shocks in the outer heliosphere, which accelerated ions to at least 11-20 MeV/nucleon and some electrons to 7 MeV. The measurements of the particle, plasma, and magnetic-field characteristics at shock fronts over the radial range 7 to 28 AU were used to determine the conditions which must exist at a traveling IP shock for both ion and relativistic electron acceleration to occur. Based on these conditions, a model was developed, which uses the shock-drift mechanism in conjunction with Fermi acceleration to explain the simultaneous acceleration of electrons and ions in the outer heliosphere.

Lopate, C.↗
NASA NTRS1981

Dielectric tensor of a weakly relativistic, nonequilibrium, and magnetized plasma

The dielectric tensor of a magnetized plasma is investigated in order to obtain analytic expressions which can be used to study wave-plasma interactions for the case of wave frequencies lying close to the cyclotron harmonics of ions or electrons. Work on this topic by Trubnikov (1959) and Shkarovsky (1966) is extended to the nonequilibrium cases, in which the particles may have drift motions, thermal anisotropy, and loss-cone-type distributions. Account is taken of these features, and a set of analytic expressions in terms of an integral is derived. The discussion presented is pertinent to the absorption and emission of radiation near the electron cyclotron harmonics.

Tsai, S. T.↗
DOE OSTI2024

Conditional guided generative diffusion for particle accelerator beam diagnostics

Abstract Advanced accelerator-based light sources such as free electron lasers (FEL) accelerate highly relativistic electron beams to generate incredibly short (10s of femtoseconds) coherent flashes of light for dynamic imaging, whose brightness exceeds that of traditional synchrotron-based light sources by orders of magnitude. FEL operation requires precise control of the shape and energy of the extremely short electron bunches whose characteristics directly translate into the properties of the produced light. Control of short intense beams is difficult due to beam characteristics drifting with time and complex collective effects such as space charge and coherent synchrotron radiation. Detailed diagnostics of beam properties are therefore essential for precise beam control. Such measurements typically rely on a destructive approach based on a combination of a transverse deflecting resonant cavity followed by a dipole magnet in order to measure a beam’s 2D time vs energy longitudinal phase-space distribution. In this paper, we develop a non-invasive virtual diagnostic of an electron beam’s longitudinal phase space at megapixel resolution (1024 × 1024) based on a generative conditional diffusion model. We demonstrate the model’s generative ability on experimental data from the European X-ray FEL.

43 PARTICLE ACCELERATORS↗