Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Power Law”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Transition to Petschek Reconnection in Subrelativistic Pair Plasmas: Implications for Particle Acceleration

While relativistic magnetic reconnection in pair plasmas has emerged in recent years as a candidate for the origin of radiation from extreme astrophysical environments, the corresponding subrelativistic pair-plasma regime has remained less explored, leaving open the question of how relativistic physics affects reconnection. In this paper, we investigate the differences between these regimes by contrasting two-dimensional particle-in-cell simulations of reconnection in pair plasmas with relativistic magnetization (σ ≫ 1) and subrelativistic magnetization (σ < 1). By utilizing unprecedentedly large domain sizes and outflow boundary conditions, we demonstrate that lowering the magnetization results in a change in the reconnection geometry from a plasmoid chain to a Petschek geometry, where laminar exhausts bounded by slow-mode shocks emanate from a single diffusion region. We attribute this change to the reduced plasmoid production rate in the low-σ case: When the secondary tearing rate is sufficiently low, plasmoids are too few in number to prevent the system from relaxing into a stable Petschek configuration. This geometric change also affects particle energization: We show that while high-σ plasmoid chains generate power-law energy spectra, low-σ Petschek exhausts merely heat incoming plasma and yield negligible nonthermal acceleration. These results have implications for predicting the global current sheet geometry and the resulting energy spectra in a variety of systems.

Plasma astrophysics↗
DOE OSTI2026

Constraints on the Correlation of IceCube Neutrinos with a Tracer of Nearby Large-scale Structure

The IceCube Neutrino Observatory has observed extragalactic astrophysical neutrinos with an apparently isotropic distribution. Only a small fraction of the observed astrophysical neutrinos can be explained by known sources. Neutrino production is thought to occur in energetic environments that are ultimately powered by the gravitational collapse of dense regions of the large-scale mass distribution in the universe. Whatever their identity, neutrino sources likely trace this large-scale mass distribution. The clustering of neutrinos with a tracer of the large-scale structure may provide insight into the distribution of neutrino sources with respect to redshift and the identity of neutrino sources. We implement a two-point angular cross correlation of the Northern sky track events with an infrared galaxy catalog derived from the Wide-field Infrared Survey Explorer (WISE) and Two Micron All Sky Survey (2MASS) source catalogs, which trace the nearby large-scale structure. No statistically significant correlation is found between the neutrinos and this infrared galaxy catalog. We find that ≤54% of the diffuse muon neutrino flux can be attributed to sources correlated with the galaxy catalog with 90% confidence. Additionally, when assuming that the neutrino source comoving density evolves following a power law in redshift, dN s /dV ∝ (1 + z) k , we find that sources with negative evolution, in particular k < −1.75, are disfavored at the 90% confidence level.

Abbasi, R. [Loyola University, Chicago, IL (United↗
DOE OSTI2026

The Role of Magnetic Reconnection in Energizing Protons and Heavier Ions at the Heliospheric Current Sheet

During near-Sun crossings of the heliospheric current sheet (HCS), Parker Solar Probe (PSP) observed populations of high-energy protons and heavier ions, indicating possible energization by magnetic reconnection up to 10 s–100 s keV nucleon −1 . Here we study ion acceleration by magnetic reconnection at the HCS. To estimate ion energization, we solve the Parker transport equation coupled to a large-scale 2D MHD reconnection simulation. We find that multiple ion species develop power-law distributions with both spectral index and high-energy cutoff E max consistent with in situ data. By accounting for the injection physics determined by kinetic simulations, we confirm that the charge-to-mass ratio scales as E max ∝ (Q/M) α with α ∼ 0.8–1.1, approximately consistent with PSP measurements in the broader range α ∼ 0.6–1.7. In the limit where ions are injected at the same energy per nucleon, α can be as low as ∼0.3. These findings further support the role of magnetic reconnection in producing high-energy heavy ions at the HCS.

Murtas, Giulia [West Virginia Univ., Morgantown, W↗
DOE OSTI2026

Magnetohydrodynamic Wave Mode Conversion in a Pseudostreamer Structure

The interaction of a fast magnetoacoustic wave with a magnetic null point in a pseudostreamer magnetic configuration is modeled. The incident fast wave is excited by an impulsive velocity driver located in the corona outside the Alfvén-acoustic equipartition layer. The fast wave front experiences refraction and wraps around the null point, while the front segment approaching the null point along the magnetic bisector passes through it. In the equipartition layer surrounding the null point, part of the fast wave energy is converted into slow magnetoacoustic waves. Slow waves propagate outward from the null point along the magnetic separatrices. The fronts of slow waves exhibit a characteristic chevron shape. The fast-to-slow mode conversion efficiency increases with the width of the initial perturbation, following a power-law dependence with an index of approximately 1.06. For initial wave amplitudes below about 10% of the local Alfvén speed, the dependence of the conversion efficiency on the initial amplitude is weak. After interacting with the null point, the amplitudes of the fast waves decrease significantly owing to energy conversion and the spreading of wave energy over a broad spatial domain. In contrast, the energy of slow waves remains concentrated in the vicinity of the magnetic separatrices, resulting in less pronounced geometric attenuation of these waves. These findings are consistent with observational reports of the interaction of a coronal fast wave with a pseudostreamer.

Zhong, Yu [Univ. of Warwick, Coventry (United King↗
DOE OSTI2026

Luminosity Function of Quasars at 1.0 < z < 3.5 from SDSS and DESI

We present a study of the evolution of type 1 quasars at 1.0 < z < 3.5, covering the peak epoch of quasar activity. The quasar evolution has been extensively explored by a variety of previous works and the derived quasar luminosity functions (QLFs) are not very consistent with each other, presumably due to the complexities introduced by different quasar selection techniques and associated completeness corrections. We use a new strategy to construct QLFs based on a library of all known quasars. We focus on a wide region of ∼1700 deg 2 and a deep field of ∼265 deg 2 that have rich spectroscopic data primarily from the Sloan Digital Sky Survey and DESI. We then apply traditional color cuts in the rest-frame UV–optical to select quasar candidates and use the quasar library to identify them. Our final sample consists of 62,426 quasars at 1.0 < z < 3.5, with a high completeness (∼96%) and a high purity (∼93%) in the color selection. Simple color cuts can potentially minimize selection biases for the study of quasar evolution. We derive binned QLFs and characterize them using a double power-law model. Sample incompleteness and contamination are considered as part of the uncertainties in the calculation. Compared to previous results, our QLFs are slightly higher at the faint end, and also higher at the bright end at 2.5 < z < 3.5. The QLFs suggest that the quasar evolution at 1.0 < z < 2.5 can be well described by the pure luminosity evolution model, while at 2.5 < z < 3.5, it can be described by either the pure luminosity evolution or the pure density evolution model.

Luminosity function↗
DOE OSTI2023

Fantastic Fits with fantasy of Active Galactic Nuclei Spectra: Exploring the Fe II Emission near the Hα Line

In this study, a refined approach for multicomponent fitting of active galactic nuclei (AGNs) spectra is presented utilizing the newly developed Python code Fully Automated pythoN Tool for AGN Spectra analYsis (fantasy). AGN spectra are modeled by simultaneously considering the underlying broken power-law continuum, predefined emission line lists, and an Fe II model, which is here extended to cover the wavelength range 3700–11000 Å. The Fe II model, founded solely on atomic data, effectively describes the extensive emission of the complex iron ion in the vicinity of the Hγ and Hβ lines, as well as near the Hα line, which was previously rarely studied. The proposed spectral fitting approach is tested on a sample of high-quality AGN spectra from the Sloan Digital Sky Survey Data Release 17. The results indicate that when Fe II emission is present near Hβ, it is also detected redward from Hα, potentially contaminating the broad Hα line wings and thus affecting the measurements of its flux and width. The production of Fe II emission is found to be strongly correlated with Eddington luminosity and appears to be controlled by a similar mechanism as the hydrogen Balmer lines. The study highlights the benefits of fitting type 1 AGN spectra with the fantasy code, pointing that it may be used as a robust tool for analyzing a large number of AGN spectra in the coming spectral surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Stochastic Approach to Reconstruct Gamma-Ray-burst Light Curves

Gamma-ray bursts (GRBs), as they are observed at high redshift ( z = 9.4), are vital to cosmological studies and investigating Population III stars. To tackle these studies, we need correlations among relevant GRB variables with the requirement of small uncertainties on their variables. Thus, we must have good coverage of GRB light curves (LCs). However, gaps in the LC hinder the precise determination of GRB properties and are often unavoidable. Therefore, extensive categorization of GRB LCs remains a hurdle. We address LC gaps using a stochastic reconstruction, wherein we fit two preexisting models (the Willingale model; W07; and a broken power law; BPL) to the observed LC, then use the distribution of flux residuals from the original data to generate data to fill in the temporal gaps. We also demonstrate a model-independent LC reconstruction via Gaussian processes. At 10% noise, the uncertainty of the end time of the plateau, its correspondent flux, and the temporal decay index after the plateau decreases by 33.3%, 35.03%, and 43.32% on average for the W07, and by 33.3%, 30.78%, 43.9% for the BPL, respectively. The uncertainty of the slope of the plateau decreases by 14.76% in the BPL. After using the Gaussian process technique, we see similar trends of a decrease in uncertainty for all model parameters for both the W07 and BPL models. These improvements are essential for the application of GRBs as standard candles in cosmology, for the investigation of theoretical models, and for inferring the redshift of GRBs with future machine-learning analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

IceCat-1: The IceCube Event Catalog of Alert Tracks

We present a catalog of likely astrophysical neutrino track-like events from the IceCube Neutrino Observatory. IceCube began reporting likely astrophysical neutrinos in 2016, and this system was updated in 2019. The catalog presented here includes events that were reported in real time since 2019, as well as events identified in archival data samples starting from 2011. We report 275 neutrino events from two selection channels as the first entries in the catalog, the IceCube Event Catalog of Alert Tracks, which will see ongoing extensions with additional alerts. The Gold and Bronze alert channels respectively provide neutrino candidates with a 50% and 30% probability of being astrophysical, on average assuming an astrophysical neutrino power-law energy spectral index of 2.19. For each neutrino alert, we provide the reconstructed energy, direction, false-alarm rate, probability of being astrophysical in origin, and likelihood contours describing the spatial uncertainty in the alert's reconstructed location. We also investigate a directional correlation of these neutrino events with gamma-ray and X-ray catalogs, including 4FGL, 3HWC, TeVCat, and Swift-BAT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

GRB 221009A: The B.O.A.T. Burst that Shines in Gamma Rays

We present a complete analysis of Fermi Large Area Telescope (LAT) data of GRB 221009A, the brightest gamma-ray burst (GRB) ever detected. The burst emission above 30 MeV detected by the LAT preceded, by 1 s, the low-energy (<10 MeV) pulse that triggered the Fermi Gamma-Ray Burst Monitor (GBM), as has been observed in other GRBs. The prompt phase of GRB 221009A lasted a few hundred seconds. It was so bright that we identify a bad time interval of 64 s caused by the extremely high flux of hard X-rays and soft gamma rays, during which the event reconstruction efficiency was poor and the dead time fraction quite high. The late-time emission decayed as a power law, but the extrapolation of the late-time emission during the first 450 s suggests that the afterglow started during the prompt emission. We also found that high-energy events observed by the LAT are incompatible with synchrotron origin, and, during the prompt emission, are more likely related to an extra component identified as synchrotron self-Compton (SSC). A remarkable 400 GeV photon, detected by the LAT 33 ks after the GBM trigger and directionally consistent with the location of GRB 221009A, is hard to explain as a product of SSC or TeV electromagnetic cascades, and the process responsible for its origin is uncertain. Because of its proximity and energetic nature, GRB 221009A is an extremely rare event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Hard Synchrotron Spectra from Magnetically Dominated Plasma Turbulence

Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically dominated (${\sigma }_{0}\gg 1$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $\sim {\sigma }_{0}{\gamma }_{\mathrm{th}0}$ (here, ${\gamma }_{\mathrm{th}0}$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are preferentially oriented in the plane perpendicular to the field. This energy-dependent anisotropy leads to a synchrotron spectral flux $\nu {F}_{\nu }\propto {\nu }^{s}$ that is much harder than for isotropic particles. Remarkably, for ${\sigma }_{0}\gg 1$ we find that the angle-integrated spectral slope in the slow cooling regime is $s\sim 0.5-0.7$ for a wide range of turbulence fluctuations, $0.25\lesssim \delta {B}_{\mathrm{rms}0}^{2}/{B}_{0}^{2}\lesssim 4$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $s\gtrsim 0.7$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Further, our findings may help explain the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of pulsar wind nebulae.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Bounds on Velocity-dependent Dark Matter–Proton Scattering from Milky Way Satellite Abundance

We use the latest measurements of the Milky Way satellite population from the Dark Energy Survey and Pan-STARRS1 to infer the most stringent astrophysical bound to date on velocity-dependent interactions between dark matter particles and protons. We model the momentum-transfer cross section as a power law of the relative particle velocity v with a free normalizing amplitude, σ MT = σ 0 v n , to broadly capture the interactions arising within the nonrelativistic effective theory of dark matter–proton scattering. The scattering leads to a momentum and heat transfer between the baryon and dark matter fluids in the early universe, ultimately erasing structure on small physical scales and reducing the abundance of low-mass halos that host dwarf galaxies today. From the consistency of observations with the cold collisionless dark matter paradigm, using a new method that relies on the most robust predictions of the linear perturbation theory, we infer an upper limit on σ0 of 1.4 × 10 -23 , 2.1 × 10 -19 , and 1.0 × 10 -12 cm 2 , for interaction models with n = 2, 4, and 6, respectively, for a dark matter particle mass of 10 MeV. These results improve observational limits on dark matter–proton scattering by orders of magnitude and thus provide an important guide for viable sub-GeV dark matter candidates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Find the Gap: Black Hole Population Analysis with an Astrophysically Motivated Mass Function

We introduce a novel black hole mass function that realistically models the physics of pair-instability supernovae with a minimal number of parameters. Applying this to all events in the LIGO-Virgo Gravitational-Wave Transient Catalog 2 (GWTC-2), we detect a peak at M BHMG = 46$^{+17}_{-6}$ M ⊙ . Repeating the analysis without the black holes from the event GW190521, we find this feature at M BHMG =54± 6 M ⊙ . These results establish the edge of the anticipated "black hole mass gap" at a value compatible with the expectation from standard stellar structure theory. The mass gap manifests itself as a discontinuity in the mass function and is populated by a distinct, less-abundant population of higher-mass black holes. We find that the primary black hole population scales with power-law index -1.95 ± 0.51 (-1.97 ± 0.44) with (without) GW190521, consistent with models of star formation. Using Bayesian techniques, we establish that our mass function fits a new catalog of black hole masses approximately as well as pre-existing phenomenological mass functions. We also remark on the implications of these results for constraining or discovering new phenomena in nuclear and particle physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the Jet–Ejecta Interaction in 3D GRMHD Simulations of a Binary Neutron Star Merger Aftermath

Short γ-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta's isotropic equivalent mass along the pole at the time of the BH formation is ≲10 -4 M ⊙ , setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the postmerger remnant disk contains a strong large-scale vertical magnetic field, ≳1015 G; and (iii) if the jets are weak (≲10 50 erg), the ejecta's isotropic equivalent mass along the pole must be small (≲10 -2 M ⊙ ). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization (≲1), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope ≳3), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope ~1.5).

Jets↗
DOE OSTI2022

Ion and Electron Acceleration in Fully Kinetic Plasma Turbulence

Turbulence is often invoked to explain the origin of nonthermal particles in space and astrophysical plasmas. By means of 3D fully kinetic particle-in-cell simulations, we demonstrate that turbulence in low-β plasmas (β is the ratio of plasma pressure to magnetic pressure) accelerates ions and electrons into a nonthermal energy distribution with a power-law energy range. The ion spectrum is harder than the electron one, and both distributions get steeper for higher β. We show that the energization of electrons is accompanied by a significant energy-dependent pitch-angle anisotropy, with most electrons moving parallel to the local magnetic field, while ions stay roughly isotropic. We demonstrate that particle injection from the thermal pool occurs in regions of high current density. Parallel electric fields associated with magnetic reconnection are responsible for the initial energy gain of electrons, whereas perpendicular electric fields control the overall energization of ions. Our findings have important implications for the origin of nonthermal particles in space and astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the Statistics of Elsasser Increments in Solar Wind and Magnetohydrodynamic Turbulence

In this Letter we investigate the dependency with scale of the empirical probability distribution functions (PDF) of Elsasser increments using large sets of WIND data (collected between 1995 and 2017) near 1 au. The empirical PDF are compared to the ones obtained from high-resolution numerical simulations of steadily driven, homogeneous reduced MHD turbulence on a 2048 3 rectangular mesh. A large statistical sample of Alfvénic increments is obtained by using conditional analysis based on the solar wind average properties. The PDF tails obtained from observations and numerical simulations are found to have exponential behavior in the inertial range, with an exponential decrement that satisfies power laws of the form α l ∝ l –μ , where l is the scale size, with μ between 0.17 and 0.25 for observations and 0.43 for simulations. PDF tails were extrapolated assuming their exponential behavior extends to arbitrarily large increments in order to determine structure function scaling laws at very high orders. Our results point to potentially universal scaling laws governing the PDF of Elsasser increments and to an alternative approach to investigate high-order statistics in solar wind observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Radiative Reconnection-powered TeV Flares from the Black Hole Magnetosphere in M87

Abstract Active galactic nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region near the nucleus. However, the emission mechanism for these flares is not well understood. Recent high-resolution general-relativistic magnetohydrodynamic simulations show the occurrence of episodic magnetic reconnection events that can power flares near the black hole event horizon. In this work, we analyze the radiative properties of the reconnecting current layer under the extreme plasma conditions applicable to the black hole in M87 from first principles. We show that abundant pair production is expected in the vicinity of the reconnection layer, to the extent that the produced secondary pair plasma dominates the reconnection dynamics. Using analytic estimates backed by two-dimensional particle-in-cell simulations, we demonstrate that in the presence of strong synchrotron cooling, reconnection can produce a hard power-law distribution of pair plasma imprinted in the outgoing synchrotron (up to a few tens of MeV) and the inverse-Compton signal (up to TeV). We produce synthetic radiation spectra from our simulations, which can be directly compared with the results of future multiwavelength observations of M87* flares.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Limit on Supernova Emission in the Brightest Gamma-Ray Burst, GRB 221009A

We present photometric and spectroscopic observations of the extraordinary gamma-ray burst (GRB) 221009A in search of an associated supernova. Some past GRBs have shown bumps in the optical light curve that coincide with the emergence of supernova spectral features, but we do not detect any significant light-curve features in GRB 221009A, nor do we detect any clear sign of supernova spectral features. Using two well-studied GRB-associated supernovae (SN 2013dx, M r,max =-19.54;$ SN 2016jca, M r,max =-19.04$) at a similar redshift as GRB 221009A (z = 0.151), we modeled how the emergence of a supernova would affect the light curve. If we assume the GRB afterglow to decay at the same rate as the X-ray data, the combination of afterglow and a supernova component is fainter than the observed GRB brightness. For the case where we assume the best-fit power law to the optical data as the GRB afterglow component, a supernova contribution should have created a clear bump in the light curve, assuming only extinction from the Milky Way. If we assume a higher extinction of E(B - V) = 1.74 mag (as has been suggested elsewhere), the supernova contribution would have been hard to detect, with a limit on the associated supernova of M r,max ≈ -$19.54. We do not observe any clear supernova features in our spectra, which were taken around the time of expected maximum light. The lack of a bright supernova associated with GRB 221009A may indicate that the energy from the explosion is mostly concentrated in the jet, leaving a lower energy budget available for the supernova.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Chocolate Chip Cookie Model: Dust-to-metal Ratio of H ii Regions

Using a sample of face-on star-forming galaxies selected from the Sloan Digital Sky Survey, we statistically derive the typical optical depth $\tau$ cl of individual H II regions based on the "Chocolate Chip Cookie" model of Lu et al. (2022). By binning galaxies into stellar mass and gas-phase metallicity bins and interpreting $\tau$ cl as the dust-to-gas ratio (DGR) of H II regions, we further investigate the correlations between DGR and stellar mass and between DGR and gas-phase metallicity, respectively. We find that DGR increases monotonically with the stellar mass of galaxies. At a given stellar mass, DGR shows a linear correlation with the gas-phase metallicity, which implies a constant dust-to-metal ratio (DTM) of galaxies at a given stellar mass. These results adequately indicate that the DTM of galaxies is simply a function of their stellar masses. In terms of gas-phase metallicity, because of the mass–metallicity relation, DTM increases with increasing metallicity with a power-law index 1.45 in the low-metallicity region, while it remains constant at the high-metallicity end.

79 ASTRONOMY AND ASTROPHYSICS↗