Reverse current induced by injection of a relativistic electron beam into a pinched plasma
Pinched plasma reversed current induction due to relativistic electron beam injection
SEARCH · Engineering Papers
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.
Pinched plasma reversed current induction due to relativistic electron beam injection
Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special relativistic module. Entity is the first PIC code designed to solve the Vlasov–Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 ns per particle per time step. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.
Abstract It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multiwavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multiwavelength flares. Nevertheless, previous works have not directly evaluated γ -ray signatures from first-principles simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multiwavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and Fermi-LAT γ -ray bands as well as closely correlated optical and γ -ray flares. The swings in optical polarization angle are also accompanied by γ -ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self-Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast γ -ray flares or orphan γ -ray flares under the leptonic scenario, complementary to the frequently considered minijet scenario. It may also imply neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.
Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) are responsible for particle acceleration in relativistic pair jets. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic pair jet propagating through a pair plasma. Simulations show that the Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. Simulation results show that this instability generates and amplifies highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The "jitter' I radiation from deflected electrons can have different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. The growth rate of the Weibel instability and the resulting particle acceleration depend on the magnetic field strength and orientation, and on the initial particle distribution function. In this presentation we explore some of the dependencies of the Weibel instability and resulting particle acceleration on the magnetic field strength and orientation, and the particle distribution function.
Modifications of large transverse momentum single-hadron, dihadron, and γ -hadron spectra in relativistic heavy-ion collisions are direct consequences of parton-medium interactions in the quark-gluon plasma (QGP). The interaction strength and underlying dynamics can be quantified by the jet transport coefficient q ̂ . We carry out the first global constraint on q ̂ using a next-to-leading order pQCD parton model with higher-twist parton energy loss and combining world experimental data on single-hadron, dihadron, and γ -hadron suppression at both RHIC and LHC energies with a wide range of centralities. The global Bayesian analysis using the information field (IF) priors provides the most stringent constraint on q ̂ ( T ) . We demonstrate in particular the progressive constraining power of the IF Bayesian analysis on the strong temperature dependence of q ̂ using data from different centralities and colliding energies. We also discuss the advantage of using both inclusive and correlation observables with different geometric biases. As a verification, the obtained q ̂ ( T ) is shown to describe data on single-hadron anisotropy at high transverse momentum well. Predictions for future jet quenching measurements in oxygen-oxygen collisions are also provided. Published by the American Physical Society 2024
In nuclear and particle physics, reconciling sophisticated simulations with experimental data is vital for understanding complex systems like the Quark Gluon Plasma (QGP) generated in heavy ion collisions. However, computational demands pose challenges, motivating using Gaussian Process emulators for efficient parameter extraction via Bayesian calibration. We conduct a comparative analysis of Gaussian Process emulators in heavy-ion physics to identify the most adept emulator for parameter extraction with minimal uncertainty. Furthermore, our study contributes to advancing computational techniques in heavy-ion physics, enhancing our ability to interpret experimental data and understand QGP properties.
Photon pair-production processes in marginally relativistic (transrelativistic) accretion plasmas are investigated in both Planckian and non-Planckian cases. Pair production in a plasma with an equilibrium (Planck) spectrum is reviewed, and pair-concentration calculations are performed for three general non-Planckian situations most relevant to black-hole accretion scenarios: steady-state transrelativistic plasmas of relatively high density characterized by a pure bremsstrahlung spectrum, a comptonized bremsstrahlung spectrum, and an unsaturated Compton scattering spectrum. The results obtained indicate that for transrelativistic temperatures (600 million to 6 billion K) photon pair production is not generally a dominant process for a plasma with a pure bremsstrahlung spectrum, but becomes dominant for plasmas where comptonization is important. It is also shown that photon pair-creation processes in a transrelativistic bremsstrahlung-radiating plasma that is more than marginally optically thick to Compton scattering significantly alter the plasma's spectrum by forcing it to become black-body before it reaches relativistic temperatures. Pair production and instabilities in unsteady-state plasmas are briefly considered
Abstract Intense lasers can accelerate electrons to very high energy over a short distance. Such compact accelerators have several potential applications including fast ignition, high energy physics, and radiography. Among the various schemes of laser-based electron acceleration, vacuum laser acceleration has the merits of super-high acceleration gradient and great simplicity. Yet its realization has been difficult because injecting free electrons into the fast-oscillating laser field is not trivial. Here we demonstrate free-electron injection and subsequent vacuum laser acceleration of electrons up to 20 MeV using the relativistic transparency effect. When a high-contrast intense laser drives a thin solid foil, electrons from the dense opaque plasma are first accelerated to near-light speed by the standing laser wave in front of the solid foil and subsequently injected into the transmitted laser field as the opaque plasma becomes relativistically transparent. It is possible to further optimize the electron injection/acceleration by manipulating the laser polarization, incident angle, and temporal pulse shaping. Our result also sheds light on the fundamental relativistic transparency process, crucial for producing secondary particle and light sources.
In this paper, we develop a first-principles model for the relativistic magnetic reconnection rate in strongly magnetized pair plasmas. By considering the energy budget and required current density near the x-line, we analytically show that in the magnetically dominated relativistic regime, the x-line thermal pressure is significantly lower than the upstream magnetic pressure due to the extreme energy needed to sustain the current density, consistent with kinetic simulations. This causes the upstream magnetic field lines to collapse in, producing the open outflow geometry which enables fast reconnection. The result is important for understanding a wide range of extreme astrophysical environments, where fast reconnection has been evoked to explain observations such as transient flares and nonthermal particle signatures.
Lasers with intensity >1.4x10 18 W.cm -2 irradiating solid targets couple 10-50% of its energy to “hot electrons” near the critical surface, with temperature kT > mc 2 . When these hot electrons impact high-Z target ions (e.g. Au, Pt), they emit copious bremsstrahlung gamma-rays and create e + e - pairs via the Trident process. For targets thicker than ~0.1mm, the gamma-rays create secondary pairs via the Bethe-Heitler processes. The concept of using lasers to create pairs was first studied by Shearer et al (1973). Motivated by the rapid advance of short-pulse ultra-intense lasers based on chirped-pulse amplification, Liang (1994) first proposed irradiating solid Au targets with laser intensity ≥10 20 W.cm -2 to create dense e + e - pairs. Using particle-in-cell (PIC, Birdsall and Langdon, 1991) simulations, Liang et al (1998, 2002) estimated that high-energy PW lasers can in principle achieve ultra-high in-situ pair densities. Subsequent studies supported this idea and found that for gold foils thicker than ~50 μm, the pair yield is dominated by the Bethe-Heitler (BH, Heitler 1954) process. Cowan et al (1999, 2000) using the LLNL Nova PW-laser to irradiate 125μm gold foils first observed pair creation, followed by Chen et al (2009) who demonstrated copious BH pair creation using the Titan, Omega-EP and Orion lasers. However, the emerging e + /e - ratio of these experiments was ≤ few %, and the pair density was only 10 13 /cc so that the pair jet transverse size R(~mm) was < pair skin depth c/ω + (=8πn+e 2 /m) 1/2 ). This was insufficient to qualify as a bona fide ”pair plasma”.
Not Available
Blazars are relativistic magnetized plasma outflows from supermassive black holes that point very close to our line of sight. Their emission is nonthermal-dominated and highly variable across the entire electromagnetic spectrum. Relativistic magnetic reconnection has been proposed as the driver of particle acceleration during blazar flares. While recent particle-in-cell (PIC) simulations have self-consistently studied the evolution of magnetic reconnection and particle acceleration therein, the resulting radiation signatures have not been systematically explored. In particular, the polarization signatures, which directly reflect the characteristic strongly dynamical evolution of magnetic field during reconnection, have not been carefully investigated. In this paper, we present a systematic study of radiation and polarization signatures arising from magnetic reconnection in blazars, based on combined PIC and polarized radiation transfer simulations with various physical parameters. We identify a harder-when-brighter trend in the spectral evolution. Moreover, higher-frequency bands (ultraviolet to X-ray) tend to flare earlier than lower-frequency bands (infrared to optical) in the synchrotron spectral component. Most importantly, polarization signatures appear more variable with higher frequencies. We find that the variation in temporal polarization depends strongly on the guide field strength. Specifically, reconnection with a significant guide field component leads to a very high polarization degree that contradicts typical blazar observations, while large polarization angle rotations are unique signatures of magnetic reconnection between nearly antiparallel magnetic field lines. These rotations are at least 90° and can extend to >180°, and they may be in either direction. Furthermore, these results imply that blazars that have shown large polarization angle rotations intrinsically have more nearly antiparallel magnetic field morphology.
Stellar wind is of great interest, because many of the astrophysical systems possess wind like structures. For a long time only nonrelativistic winds have been studied. However, recently it was proposed by Kennel et al. that the relativistic pulsar wind with the plasma, consisting of electrons and positrons, can be responsible for the observed features of the Crab nebula. The study of Kennel et al. revealed the inconsistency of the assumption of the wind zero temperature with the observational data. It has been shown that only for high relativistic temperatures can high Mach numbers be reached, which allows the possibility of a shock formation. The analysis of Kennel et al. was extended into the 2-D case by using conventional magnetohydrodynamic equations for a relativistic plasma with an isotropic relativistic temperature. The state equation is assumed polytropic.
Not Available
Magnetic field amplification by relativistic streaming plasma instabilities is central to a wide variety of high-energy astrophysical environments as well as to laboratory scenarios associated with intense lasers and electron beams. We report on a new secondary nonlinear instability which arises for relativistic dilute electron beams after the saturation of the linear Weibel instability. This instability grows due to the transverse magnetic pressure associated with the beam current filaments, which cannot be quickly neutralized due to the inertia of background ions. We show that it can amplify the magnetic field strength and spatial scale by orders of magnitude, leading to large-scale plasma cavities with strong magnetic field and to very efficient conversion of the beam kinetic energy into magnetic energy. The instability growth rate, saturation level, and scale length are derived analytically and shown to be in good agreement with fully-kinetic simulations.
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.
We describe here efforts to create and study magnetized electron–positron pair plasmas, the existence of which in astrophysical environments is well-established. Laboratory incarnations of such systems are becoming ever more possible due to novel approaches and techniques in plasma, beam and laser physics. Traditional magnetized plasmas studied to date, both in nature and in the laboratory, exhibit a host of different wave types, many of which are generically unstable and evolve into turbulence or violent instabilities. This complexity and the instability of these waves stem to a large degree from the difference in mass between the positively and the negatively charged species: the ions and the electrons. The mass symmetry of pair plasmas, on the other hand, results in unique behaviour, a topic that has been intensively studied theoretically and numerically for decades, but experimental studies are still in the early stages of development. A levitated dipole device is now under construction to study magnetized low-energy, short-Debye-length electron–positron plasmas; this experiment, as well as a stellarator device that is in the planning stage, will be fuelled by a reactor-based positron source and make use of state-of-the-art positron cooling and storage techniques. Relativistic pair plasmas with very different parameters will be created using pair production resulting from intense laser–matter interactions and will be confined in a high-field mirror configuration. We highlight the differences between and similarities among these approaches, and discuss the unique physics insights that can be gained by these studies.
We present a systematic study of $\Lambda$ hyperon's polarization observables using event-by-event (3+1)D relativistic hydrodynamics. The effects of initial hot spot size and QGP's specific shear viscosity on the polarization observables are quantified. We examine the effects of the two formulations of the thermal shear tensor on the polarization observables using the same hydrodynamic background. With event-by-event simulations, we make predictions for the Fourier coefficients of $\Lambda$'s longitudinal polarization $P^z$ with respect to the event planes of different orders of anisotropic flow. We propose new correlations among the Fourier coefficients of $P^z$ and charged hadron anisotropic flow coefficients to further test the mapping from fluid velocity gradients to hyperon's polarization. Finally, we present a system size scan with Au+Au, Ru+Ru, and O+O collisions at $\sqrt{s_\mathrm{NN}} = 200$ GeV to study the system size dependence of polarization observables at the Relativistic Heavy-ion Collider.