SEARCH · Search NASA

Results for “PICS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. I. A Systematic Study

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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Synchrotron Firehose Instability

We demonstrate using linear theory and particle-in-cell (PIC) simulations that a synchrotron-cooling collisionless plasma acquires pressure anisotropy and, if the plasma beta is sufficiently high, becomes unstable to the firehose instability, in a process that we dub the synchrotron firehose instability (SFHI). The SFHI channels free energy from the pressure anisotropy of the radiating, relativistic electrons (and/or positrons) into small-amplitude, kinetic-scale, magnetic-field fluctuations, which pitch-angle scatter the particles and bring the plasma to a near-thermal state of marginal instability. The PIC simulations reveal a nonlinear cyclic evolution of firehose bursts interspersed by periods of stable cooling. We compare the SFHI for electron–positron and electron–ion plasmas. As a byproduct of the growing electron-firehose magnetic-field fluctuations, magnetized ions gain a pressure anisotropy opposite to that of the electrons. If these ions are relativistically hot, we find that they also experience cooling due to collisionless thermal coupling with the electrons, which we argue is mediated by a secondary ion-cyclotron instability. We suggest that the SFHI may be activated in a number of astrophysical scenarios, such as within ejecta from black hole accretion flows and relativistic jets, where the redistribution of energetic electrons from low to high pitch angles may cause transient bursts of radiation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

First-principles-integrated Study of Blazar Synchrotron Radiation and Polarization Signatures from Magnetic Turbulence

Blazar emission is dominated by nonthermal radiation processes that are highly variable across the entire electromagnetic spectrum. Turbulence, which can be a major source of nonthermal particle acceleration, can widely exist in the blazar emission region. The Turbulent Extreme Multi-Zone (TEMZ) model has been used to describe turbulent radiation signatures. Recent particle-in-cell (PIC) simulations have also revealed the stochastic nature of the turbulent emission region and particle acceleration therein. However, radiation signatures have not been systematically studied via first-principles-integrated simulations. In this paper, we perform combined PIC and polarized radiative transfer simulations to study synchrotron emission from magnetic turbulence in the blazar emission region. We find that the multiwavelength flux and polarization are generally characterized by stochastic patterns. Specifically, the variability timescale and average polarization degree (PD) are governed by the correlation length of the turbulence. Interestingly, magnetic turbulence can result in polarization angle swings with arbitrary amplitudes and duration, in either direction, that are not associated with changes in flux or PD. Surprisingly, these swings, which are stochastic in nature, can appear either bumpy or smooth, although large-amplitude swings (>180°) are very rare, as expected. Our radiation and polarization signatures from first-principles-integrated simulations are consistent with the TEMZ model, except that in the latter, there is a weak correlation, with zero lag, between flux and degree of polarization.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Effective Resistivity in Relativistic Collisionless Reconnection

Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing nonthermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations, we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without a guide field. In X-points, where diverging flows result in a nondiagonal thermal pressure tensor, a finite residence time for particles gives rise to a localized collisionless effective resistivity. Here, for the first time for relativistic reconnection in a fully developed plasmoid chain, we identify the mechanisms driving the nonideal electric field using a full Ohm law by means of a statistical analysis based on our PIC simulations. We show that the nonideal electric field is predominantly driven by gradients of nongyrotropic thermal pressures. We propose a kinetic physics motivated nonuniform effective resistivity model that is negligible on global scales and becomes significant only locally in X-points. It captures the properties of collisionless reconnection with the aim of mimicking its essentials in nonideal magnetohydrodynamic descriptions. This effective resistivity model provides a viable opportunity to design physically grounded global models for reconnection-powered high-energy emission.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Particle-in-cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms

Abstract Relativistic magnetic reconnection is a nonideal plasma process that is a source of nonthermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and numerical techniques for efficiently modeling such systems has been limited. With the GPU-accelerated PIC code WarpX, we explore the accuracy and potential performance benefits of two advanced Maxwell solver algorithms: a nonstandard finite-difference scheme (CKC) and an ultrahigh-order pseudo-spectral method (PSATD). We find that, for the relativistic reconnection problem, CKC and PSATD qualitatively and quantitatively match the standard Yee-grid finite-difference method. CKC and PSATD both admit a time step that is 40% longer than that of Yee, resulting in a ∼40% faster time to solution for CKC, but no performance benefit for PSATD when using a current deposition scheme that satisfies Gauss’s law. Relaxing this constraint maintains accuracy and yields a 30% speedup. Unlike Yee and CKC, PSATD is numerically stable at any time step, allowing for a larger time step than with the finite-difference methods. We found that increasing the time step 2.4–3 times over the standard Yee step still yields accurate results, but it only translates to modest performance improvements over CKC, due to the current deposition scheme used with PSATD. Further optimization of this scheme will likely improve the effective performance of PSATD.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

The Effect of Inverse Compton Losses on Particle Acceleration in Three-dimensional Relativistic Reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived “free phase” where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a “trapped phase,” where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection, including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor γrad exceeds the magnetization σ. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is dN free /dγ ∝ γ −1 , as in the uncooled case. The spectrum of trapped electrons above the cooling break γcool (in the range γ cool < γ < γ rad ) is dN/dγ ∝ γ −3 , steeper than the scaling dN/dγ ∝ γ −2 of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by Zhang et al. for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. I. Curvilinear Special Relativistic Module

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.

Hakobyan, Hayk [Flatiron Institute, New York, NY (↗

DOE OSTI2022

Origin of Intense Electron Heating in Relativistic Blast Waves

The modeling of gamma-ray burst afterglow emission bears witness to strong electron heating in the precursor of Weibel-mediated, relativistic collisionless shock waves propagating in unmagnetized electron–ion plasmas. In this Letter, we propose a theoretical model, which describes electron heating via a Joule-like process caused by pitch-angle scattering in the decelerating, self-induced microturbulence and the coherent charge-separation field induced by the difference in inertia between electrons and ions. The emergence of this electric field across the precursor of electron–ion shocks is confirmed by large-scale particle-in-cell (PIC) simulations. Integrating the model using a Monte Carlo-Poisson method, we compare the main observables to the PIC simulations to conclude that the above mechanism can indeed account for the bulk of electron heating.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

The effect of inverse Compton losses on particle acceleration in three-dimensional relativistic reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived "free phase" where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a "trapped phase", where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor $γ_{\rm rad}$ exceeds the magnetization $σ$. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is $dN_{\rm free}/dγ\propto γ^{-1}$, as in the uncooled case. The spectrum of trapped electrons above the cooling break $γ_{\rm cool}$ (in the range $γ_{\rm cool}<γ<γ_{\rm rad}$) is $dN/dγ\propto γ^{-3}$, steeper than the scaling $dN/dγ\propto γ^{-2}$ of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by arXiv:2302.12269 for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure, that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

FOS: Physical sciences↗

NASA NTRS1978

Preliminary design study of an alternate heat source assembly for a Brayton isotope power system

Results are presented for a study of the preliminary design of an alternate heat source assembly (HSA) intended for use in the Brayton isotope power system (BIPS). The BIPS converts thermal energy emitted by a radioactive heat source into electrical energy by means of a closed Brayton cycle. A heat source heat exchanger configuration was selected and optimized. The design consists of a 10 turn helically wound Hastelloy X tube. Thermal analyses were performed for various operating conditions to ensure that post impact containment shell (PICS) temperatures remain within specified limits. These limits are essentially satisfied for all modes of operation except for the emergency cooling system for which the PICS temperatures are too high. Neon was found to be the best choice for a fill gas for auxiliary cooling system operation. Low cycle fatigue life, natural frequency, and dynamic loading requirements can be met with minor modifications to the existing HSA.

Strumpf, H. J.↗

NASA NTRS1985

Self-consistent simulation of plasma interactions with secondary-emitting insulators

A cylindrical particle-in-cell (PIC) plasma simulation code applicable to plasma densities encountered in low Earth orbit (LEO) is described. The simulated geometries include that of a plain disk and a disk surrounded by a dielectric. Both configurations are mounted upon a ground plate in contact with a plasma environment. Techniques allowing simulation of dielectric charging using PIC time scales are discussed. Current versus voltage characteristic curves are calculated and the results are compared to experimental data.

Brandon, S. T.↗

NASA NTRS1985

The gaseous component of the disk around Beta Pictoris

Optical spectra of alpha Lyr, alpha PsA, and beta Pic have been obtained at a velocity resolution of 3 km/s. No circumstellar absorption lines of Ca II or Na I are detected toward alpha Lyr or alpha PsA at sensitive limits. In the favorable case of beta Pic, where the circumstellar disk imaged by Smith and Terrile (1984) is seen nearly edge-on, a strong, narrow, circumstellar Ca II K absorption line previously reported by Slettebak (1982) and weaker, still narrower circumstellar Na I D lines are detected. Negative results of high sensitivity also are obtained for the Ca I 4226 A and CH(+) 4232 A lines, along with upper limits on the Zn II 2026, 2062 A doublet from archival IUE spectra. Under assumptions which agree with other well-established observations of the gaseous abundances of calcium and zinc, the total gaseous column density of hydrogen along a radius of the circumstellar disk is between 10 to the 18th and 4 x 10 to the 20th/sq cm. Within the boundaries of the dust disk detected by Smith and Terrile (1984) the total gaseous mass then is less than about 2, or less than 1 percent of the mass of the planetary system. A simplified model of the density distribution in the gaseous disk yields a characteristic total density n(H) of about 100,000/cu cm, which exceeds that of all interplanetary gas at earth's position by a factor of about 10,000.

Hobbs, L. M.↗

NASA NTRS1986

Observations of gaseous circumstellar disks. III

Echelle spectra obtained at the Ca II K line and the Na I D lines are reported for seven stars which show excess fluxes in at least one of the IRAS bands at 25, 60, and 100 microns. Detectable circumstellar absorption lines arising from the gaseous components of the presumed circumstellar dust disks are found in none of these cases, in contrast to the previously reported case of Beta Pic. Conversely, six A or F stars which were known before the launch of IRAS to show gaseous circumstellar disks do not show excess fluxes in any of the IRAS bands. Similar echelle spectra of these stars show Ca II circumstellar absorption lines comparable in four cases to that of Beta Pic.

Hobbs, L. M.↗

NASA NTRS1989

Optical polarimetry and thermal imaging of the disc around Beta Pictoris

The dust disk around the main-sequence A5 star Beta Pic is characterized on the basis of recently obtained data. Optical polarimetry indicates that the disk has polarization of about 15 percent and polarization vectors normal to the scattering plane, indicating that the features observed are indeed a dust disk made visible by scattering. Thermal images of the central 13-arcsec-diameter region, obtained at 10 and 20 microns by Telesco et al. (1988) using the NASA Marshall 20-pixel bolometer array, suggest that the dust grains in the disk have a mean diameter of 1 micron or smaller. Also discussed are (1) unsuccessful searches for other nearby stars with disks and (2) the potential value of Hubble Space Telescope images for understanding the structure of Beta Pic.

Wolstencroft, R. D.↗

NASA NTRS1991

A search for circumstellar dust around HR 10, a proposed Beta Pictoris star

Archival IRAS data and recent IUE observations of the A shell star HR 10 (HD 256) are presented. HR 10 is detected in the IRAS Faint Source Survey and ADDSCAN data at 12 and 25 microns and marginally detected at 60 microns. The IUE data indicate a spectral type of A2 with 0.1 mag of UV selective extinction following a galactic average extinction law. HR 10 shows an infrared excess, compared to normal A stars of the same spectral type, beginning at 12 microns and extending to longer wavelengths. The IR excess is consistent with free-free emission and a suggestion of CS dust. The IR colors for this star differ significantly from those of Beta Pic and are closer to those of Vega. The absence of circumstellar selective extinction shortward of 3200 A suggests that the circumstellar grain size distribution does not resemble that of the interstellar medium, but may be more like the grains around Vega and Beta Pic.

Cheng, Kwang-Ping↗

NASA NTRS1991

Monte Carlo simulation of ion-neutral charge exchange collisions and grid erosion in an ion thruster

A combined particle-in-cell (PIC)/Monte Carlo simulation model has been developed in which the PIC method is used to simulate the charge exchange collisions. It is noted that a number of features were reproduced correctly by this code, but that its assumption of two-dimensional axisymmetry for a single set of grid apertures precluded the reproduction of the most characteristic feature of actual test data; namely, the concentrated grid erosion at the geometric center of the hexagonal aperture array. The first results of a three-dimensional code, which takes into account the hexagonal symmetry of the grid, are presented. It is shown that, with this code, the experimentally observed erosion patterns are reproduced correctly, demonstrating explicitly the concentration of sputtering between apertures.

Peng, Xiaohang↗

NASA NTRS1991

The circumstellar disks of Beta Pictoris analogs

A survey using data from the International Ultraviolet Explorer (IUE) of previously known B and A shell stars with IRAS detections has resulted in the identification of three stars, HD 93563, Sigma Her, and 51 Oph, which have spectral signatures of infalling circumstellar plasma similar to Beta Pic. Two of these systems have infrared flux distributions indicating the presence of circumstellar dust disks, while the other, HD 93563, has an infrared excess consistent with free-free emission from the plasma envelope. With the identification of three such systems, it is clear that infalling circumstellar plasma is more common than previously anticipated among late-type B shell stars. The absence of dust in one system, HD 93563, suggests that infalling plasma in these stars, and possibly also in Beta Pic itself, may not be due to either erosion of a dust disk or to high cometary bombardment rates, but may instead be linked to stellar activity.

Grady, Carol A.↗