Search NASA⌕ Search

SEARCH · Search NASA

Results for “PIC simulation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dynamic focusing of laser driven positron jets by self-generated fields

Abstract Focusing effect of laser-driven positron jets by self-generated target sheath fields has been observed for the first time experimentally and the results are supported by the computational studies. In the experiment, OMEGA EP short-pulse (0.7 ps, 500 J) irradiates mm-size gold targets with a concave back surface and reference flat-surface targets. Both targets exhibited positrons with quasi-monoenergetic energy peaks while targets with concave curvature also showed increased number of positrons at the detector. The data is consistent with hybrid-PIC simulations confirming that the time-varying electric fields driven by electrons escaping from the target significantly change the trajectories of positrons. These simulations show a small radius of curvature on the rear side increases the relative focusing effect and the positrons to electrons ratio in the escaping plasma. For the smallest radius of curvature, positron jets that are up to 10 times denser can be achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Plasma instabilities dominate radioactive transients magnetic fields: the self-confinement of leptons in Type Ia and core-collapse supernovae, and kilonovae

The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $\beta ^+$ and $\beta ^-$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material.

instabilities↗
DOE OSTI2021

Electron acceleration at oblique angles via stimulated Raman scattering at laser irradiance >10 16 W cm –2 μm 2

The generation of hot, directional electrons via laser-driven stimulated Raman scattering (SRS) is a topic of great importance in inertial confinement fusion (ICF) schemes. Little recent research has been dedicated to this process at high laser intensity, in which back, side, and forward scatter simultaneously occur in high energy density plasmas, of relevance to, for example, shock ignition ICF. In this paper, we present an experimental and particle-in-cell (PIC) investigation of hot electron production from SRS in the forward and near-forward directions from a single speckle laser of wavelength λ 0 = 1.053μm, peak laser intensities in the range I 0 = 0.2–1.0×10 17 Wcm –2 and target electron densities between ne = 0.3–1.6% n c , where n c is the plasma critical density. As the intensity and density are increased, the hot electron spectrum changes from a sharp cutoff to an extended spectrum with a slope temperature T = 34 ± 1 keV and maximum measured energy of 350 keV experimentally. Multidimensional PIC simulations indicate that the high energy electrons are primarily generated from SRS-driven electron plasma wave phase fronts with k vectors angled ~50° with respect to the laser axis. These results are consistent with analytical arguments that the spatial gain is maximized at an angle which balances the tendency for the growth rate to be larger for larger scattered light wave angles until the kinetic damping of the plasma wave becomes important. The efficiency of generated high energy electrons drops significantly with a reduction in either laser intensity or target electron density, which is a result of the rapid drop in growth rate of Raman scattering at angles in the forward direction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Transport of an intense proton beam from a cone-structured target through plastic foam with unique proton source modeling

Laser-accelerated proton beams are applicable to several research areas within high-energy density science, including warm dense matter generation, proton radiography, and inertial confinement fusion, which all involve transport of the beam through matter. Here, we report on experimental measurements of intense proton beam transport through plastic foam blocks. The intense proton beam was accelerated by the 10 ps, 700 $\textit{J}$ OMEGA EP laser irradiating a curved foil target, and focused by an attached hollow cone. The protons then entered the foam block of density 0.38 g/cm 3 and thickness 0.55 or 1.00 mm. At the rear of the foam block, a Cu layer revealed the cross section of the intense beam via proton- and hot electron-induced Cu-K α emission. Images of x-ray emission show a bright spot on the rear Cu film indicative of a forward-directed beam without major breakup. 2D fluid-PIC simulations of the transport were conducted using a unique multi-injection source model incorporating energy-dependent beam divergence. Along with postprocessed calculations of the Cu – $K_α$ emission profile, simulations showed that protons retain their ballistic transport through the foam and are able to heat the foam up to several keV in temperature. The total experimental emission profile for the 1.0 mm foam agrees qualitatively with the simulated profile, suggesting that the protons indeed retain their beamlike qualities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Stability of the modulator in a plasma-modulated plasma accelerator

We explore the regime of operation of the modulator stage of a recently proposed laser-plasma accelerator scheme [Phys. Rev. Lett. 127, 184801 (2021)], dubbed the plasma-modulated plasma accelerator (P-MoPA). The P-MoPA scheme offers a potential route to high-repetition-rate, GeV-scale plasma accelerators driven by picosecond-duration laser pulses from, for example, kilohertz thin-disk lasers. The first stage of the P-MoPA scheme is a plasma modulator in which a long, high-energy “drive” pulse is spectrally modulated by copropagating in a plasma channel with the low-amplitude plasma wave driven by a short, low-energy “seed” pulse. The spectrally modulated drive pulse is converted to a train of short pulses, by introducing dispersion, which can resonantly drive a large wakefield in a subsequent accelerator stage with the same on-axis plasma density as the modulator. In this paper we derive the 3D analytic theory for the evolution of the drive pulse in the plasma modulator and show that the spectral modulation is independent of transverse coordinate, which is ideal for compression into a pulse train. We then identify a transverse mode instability (TMI), similar to the TMI observed in optical fiber lasers, which sets limits on the energy of the drive pulse for a given set of laser-plasma parameters. We compare this analytic theory with particle-in-cell (PIC) simulations and find that even higher energy drive pulses can be modulated than those demonstrated in the original proposal.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). Here, the formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Investigation on laser absorption and x-ray radiation in microstructured titanium targets heated by short-pulse relativistic laser pulses

The enhancement effect of a microstructured surface on laser absorption and characteristic K α emission has been investigated by measuring K-shell emission from titanium (Ti) targets irradiated with high-intensity ( ∼ 10 20 W cm − 2 ), subpicosecond (500 fs) laser pulses. The experimental results indicate a modest enhancement ( 1.6 × ) of K α emission from microstructured targets compared to flat foils, but with a similar intensity and profile of He α and Li-like satellites. Particle-in-cell (PIC) simulations are implemented to further understand the underlying physical processes in the laser interaction with both targets, interpreting the mechanisms responsible for the K α enhancement. The reasons for the lower-than-expected enhancement of K α emission are discussed. The rapid heating of the bulk plasma might result in the premature shutdown of K α emission before the thermalization of hot electrons or even the end of laser pulses, suggesting that the use of K α emission as a diagnostic of the hot-electron yield or relaxation could lead to a misinterpretation. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Expansion-Driven Self-Magnetization of High-Energy-Density Plasmas

Understanding plasma self-magnetization is one of the fundamental challenges in both laboratory and astrophysical plasmas. Self-magnetization can modify plasma transport properties, altering the dynamical evolution of plasmas. Multiple high-energy-density (HED) experiments have observed the formation of ion-scale magnetic filaments of megagauss strength, though their origin remains debated. Here, in this study, we conduct 2D collisional particle-in-cell (PIC) simulations with a laser ray-tracing module for a fully self-consistent simulation of the plasma ablation, expansion, and magnetization. The simulations use a planar geometry, effectively suppressing the Biermann magnetic fields, to focus on anisotropy-driven instabilities. The laser intensity is varied between 10 13 and 10 14 W/cm 2 , which is relevant to HED and inertial fusion experiments where collisions must be considered. We find that, above a critical intensity, the plasma rapidly self-magnetizes via an expansion-driven Weibel process, producing a plasma beta of 100 (𝛽 = 8⁢𝜋⁢𝑘 𝐵 ⁢𝑛 𝑒 ⁢𝑇 𝑒 /𝐵 2 ) and Hall parameter 𝜔 ce ⁢𝜏 𝑒 >1 within the first few hundred picoseconds. The magnetic field is sufficiently strong to modify plasma heat transport, and simulations with an artificially suppressed magnetic field show noticeably different temperature profiles.

Lezhnin, K. V. [Princeton Plasma Physics Laborator↗
DOE OSTI2024

Theoretical Assessment of the Transition Between Electron Emission Mechanisms for Nonplanar Diodes

Theoretically and computationally describing the operation of nanodiodes requires characterizing the transitions between multiple electron emission mechanisms for nanodiodes with complicated geometries. This motivates our development of techniques to determine when simplified theories for individual mechanisms suffice compared to more complete, but more computationally expensive, models. Leveraging recent theories that define a canonical gap distance to translate planar theory to nonplanar diodes, we derive the conditions for the transitions among thermal emission, field emission, and space-charge-limited current density (SCLCD) in vacuum and with collisions for non-Cartesian coordinate systems, including spherical, cylindrical, and prolate spheroidal coordinate systems. Particle-in-cell (PIC) simulations of the current density as a function of applied voltage for a tip-to-plate geometry in vacuum agreed qualitatively with the asymptotes for thermal emission at low voltage and SCLCD at higher voltage using the canonical gap distance. As a result, this demonstrates the utility of this approach for guiding system design and suggests future extensions to save simulation time for more realistic geometries that are more computationally expensive.

Conformal mapping↗
DOE OSTI2021

First Experiments and Radiographs on the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) Dense Plasma Focus

We report a dense plasma focus (DPF) is a relatively compact coaxial plasma gun, which completes its discharge as a Z-pinch. These devices are designed to operate at a variety of scales to produce short (<100 ns) pulses of ions, X-rays, and/or neutrons. LLNL recently constructed and brought into operation a new device, the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) DPF, which is designed for radiography and high-yield operations. This device has been commissioned and has achieved neutron yields of up to 3.8E11 neutrons/pulse at 2.5-MA peak current while operating at up to 1 MJ of stored energy in its original pulsed power configuration. MJOLNIR is equipped with a wide range of diagnostics, including nuclear activation detectors, neutron time-of-flight (nToF) detectors, a fast-framing camera, optical light gates, and a time-gated neutron and X-ray imager. LLNL also runs unique particle-in-cell (PIC) simulations of DPF discharges in the Chicago code and has gained significant insight into the various physical factors that influence neutron yield. MJOLNIR is one of the first DPFs whose design and continual upgrades are heavily influenced by model predictions. In this article, we describe insights from modeling, device operation, and recent results. Comparisons between modeling predictions and measurements, as well as X-ray and neutron images are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Direct Measurement of Backstreaming Ions from the Axis 2 Septum Dump and RF Generation

The DARHT-II accelerator produces a 2-kA, 17-MeV beam in a 1600-ns pulse. After exiting the accelerator, the pulse is sliced into four short pulses by a kicker and quadrupole septum and then transported for several meters to a high Z material target for conversion to x-rays for radiography. The un-kicked beam is diverted to a graphite dump. The interaction of this beam with the dump produces ions that propagate back to the target beamline and affect the beam properties. Coupling of the electron beam and ions to the beam dump vacuum chamber can also produce significant amounts of radiofrequency fields (RF) that disrupt the beam transport to the target. This requires a change in the nominal tune in order to minimize the heating of the carbon beam dump and suppress desorption of ions. Initially, ions are generated through stimulated desorption and ionization of the monolayers of water vapor stuck to the beam dump surface. These ions propagate upstream at high speeds and act as an electrostatic focusing lens by reducing the beam space charge. As the beam size on the dump decreases, the dump is sufficiently heated so that a surface plasma is formed and much greater quantities of ions are generated through temperature limited thermionic emission. The interaction of these ions with the beam creates RF which causes significant disruption of the beam. The results of a 3D PIC simulation show that ions can propagate back to the target beam line in about 500 ns. A beam position monitor (BPM) located just upstream of the septum was modified to include four capacitive pickups in addition to four of the usual inductive pickups. The inductive and capacitive pickups provide a measure of the beam current. The capacitive pickups are very sensitive to low velocity ions. As ions from the beam dump pass through the BPM the inductive and capacitive pickup will record two different measures of the beam current. Initial results from this BPM will be discussed along with other measurements which indicate the presence of backstreaming ions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Plasma Photonics for the Control of Kilojoule Lasers

The main goal of this project was to develop plasma photonic devices (plasma lenses or ellipsoidal plasma mirrors (EPMs)) for controlling the focus of high-energy (kJ-class) laser systems. In order to defocus a kJclass laser beam, a plasma optic must produce a very specific plasma density profile. For a plasma lens, this profile must be peaked on axis and fall off parabolicly with radius. We developed four plasma-lens target designs that were predicted to yield the ideal plasma profile required for defocusing kJ-class laser beams. Particle-in-cell (PIC) simulations using the code OSIRIS were used to explore the design. The resulting plasma lens must also be underdense, which sets an upper limit on the solid density of the plasma-lens targets to ~ 1 mg/cm 3 or less. This low of a density is highly challenging to fabricate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Collaborative Research: Explosive reconnection in relativistic magnetically-dominated plasmas (Final Report)

Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Investigation on laser absorption and x-ray radiation in microstructured titanium targets heated by short-pulse relativistic laser pulses

The enhancement effect of a microstructured surface on laser absorption and characteristic Kα emission has been investigated by measuring K-shell emission from titanium (Ti) targets irradiated with high-intensity (~10 20 W cm -2 ), subpicosecond (500 fs) laser pulses. The experimental results indicate a modest enhancement (1.6×) of Kα emission from microstructured targets compared to flat foils, but with a similar intensity and profile of Heα and Li-like satellites. Particle-in-cell (PIC) simulations are implemented to further understand the underlying physical processes in the laser interaction with both targets, interpreting the mechanisms responsible for the Kα enhancement. The reasons for the lower-than-expected enhancement of Kα emission are discussed. The rapid heating of the bulk plasma might result in the premature shutdown of Kα emission before the thermalization of hot electrons or even the end of laser pulses, suggesting that the use of Kα emission as a diagnostic of the hot-electron yield or relaxation could lead to a misinterpretation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

New model for the ion collection by cylindrical probes over a wide range of collisionality

Langmuir probes remain one of the most important diagnostic tools for plasma processing applications. Modern probe analysis usually relies on the electron current part of the Langmuir probe characteristic using the Druyvesteyn method. However, for electronegative plasmas or for discharges containing dust the analysis of the ion current attracted by the probe can be desirable to determine the ion density. But, even at low pressures of a few Pa, the ion current is affected by collisions due to the large cross section for charge exchange. Available theories for collisional or collision-enhanced ion currents onto probes are complex and not well validated. Thus, in this contribution, we compare available collisional probe theories for the ion current to results of particle-in-cell (PIC) simulations. To this end, the probe surrounded by a semi-infinite plasma is simulated using a modified version of the open-source code EDIPIC. A dataset of currents for different neutral gas pressures is obtained and compared to the different theories from the literature. Based on these results, we propose a simpler and more intuitive model for the ion current collected by the probe, based on the model of Gatti and Kortshagen (Phys. Rev. E 78, 046402, 2008), developed for the charging of dust particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴ speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously previously impossible with conventional tools while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗
DOE OSTI2021

Ion Acceleration and the Development of a Power-law Energy Spectrum in Magnetic Reconnection

How charged particles are accelerated efficiently and form a power-law energy spectrum in magnetic reconnection is a problem that is not well understood. In a previous paper, it was shown that the electron Kelvin–Helmholtz instability (EKHI) in force-free magnetic reconnection generates fast-expanding vortices that can accelerate electrons in a few tens of ion gyroperiods (less than 1 ms in the solar corona) to form a power-law energy distribution. In this paper, we present a particle-in-cell (PIC) simulation study of ion acceleration in force-free magnetic reconnection in the presence of the EKHI-induced turbulence. We find that ions are not significantly accelerated by the EKHI-induced stochastic electric field until the magnetic vortices expand to sizes comparable to the ion gyroradius. The Alfvén waves generated by the EKHI couple with the magnetic vortices, leading to resonance between the ions inside the magnetic vortices and Alfvén waves and enhanced ion heating. The induced Alfvén wave resonance results in a broken power-law energy spectrum with a breakpoint at $\sim {m}_{i}{v}_{A}^{2}$, where v A is the Alfvén velocity. Here, we show that the process that forms the nonthermal tail is a second-order Fermi mechanism and the mean spectral index is α = (1 + 4 a 2 D/R) /2, where D is the spatial scale of the inductive electric field, R is that of vortices, and a = B g /B 0 , with ratio of guide field B g and asymptotic B 0 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

On the Short-scale Spatial Variability of Electron Inflows in Electron-only Magnetic Reconnection in the Turbulent Magnetosheath Observed by MMS

We investigate the detailed properties of electron inflow in an electron-only reconnection event observed by the four Magnetospheric Multiscale (MMS) spacecraft in the Earth's turbulent magnetosheath downstream of the quasi-parallel bow shock. The lack of ion coupling was attributed to the small-scale sizes of the current sheets, and the observed bidirectional super-Alfvénic electron jets indicate that the MMS spacecraft crossed the reconnecting current sheet on both sides of an active X-line. Remarkably, the MMS spacecraft observed the presence of large asymmetries in the two electron inflows, with the inflows (normal to the current sheet) on the two sides of the reconnecting current layer differing by as much as a factor of four. Furthermore, even though the four MMS spacecraft were separated by less than seven electron skin depths, the degree of inflow asymmetry was significantly different at the different spacecraft. The asymmetry in the inflow speeds was larger with increasing distances downstream from the reconnection site, and the asymmetry was opposite on the two sides of the X-line. We compare the MMS observations with a 2D kinetic particle-in-cell (PIC) simulation and find that the asymmetry in the inflow speeds stems from in-plane currents generated via the combination of reconnection-mediated inflows and parallel flows along the magnetic separatrices in the presence of a large guide field.

79 ASTRONOMY AND ASTROPHYSICS↗