Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ErGe”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The Emission and Suppression of Line Features in Luminous Transients

Featureless optical and ultraviolet (UV) spectra are a puzzling signature to emerge from recent observations of luminous fast blue optical transients (LFBOTs) and some tidal disruption events (TDEs). We describe the landscape of source and gas properties that are expected to form H, He I , and He II emission lines, and map spectral types to the parameter space of luminosity and system radius. Using one-dimensional radiative transfer calculations, we show that high source luminosities (L > 10 44 erg s −1 ) and compact ejecta radii (r < 10 14 cm) produce featureless spectra, due to the high temperature and ionization state of the emitting medium. Intermediate luminosities and moderately compact systems can generate He II -dominated spectra, while lower luminosities and more extended atmospheres result in conspicuous H and He I emission. Large expansion velocities (v ≥ 0.1c) can further broaden lines such that they blend into the continuum. Featureless UV spectra may require even more extreme ionization environments or velocities in order to suppress the many intrinsically strong metal lines at those wavelengths. Applying this framework to understand the absence of features observed in LFBOTs and featureless TDEs, we find that, under the optically thick, quasi-thermal conditions considered here, nonhomologous, compact outflows are likely necessary for featurelessness to persist in optical and UV spectra.

Aspegren, Olivia [University of California, Berkel↗
DOE OSTI2026

Characterizing the γ-Ray Emission from Low-luminosity AGN

A majority of the active galactic nuclei (AGN) in the local Universe are classified as low-luminosity AGN (LLAGN), having bolometric luminosities ≲10 42 erg s −1 . Although high-energy γ-ray emission is predicted from both the jets and disks of LLAGN, to date only four have been detected by the Fermi Large Area Telescope (Fermi-LAT). In this work, we therefore conduct a comprehensive study of all the LLAGN from the Palomar spectroscopic survey of bright, northern galaxies, including both subthreshold and detected γ-ray sources, using 14.4 yr of LAT data. Our analysis results in a new detection of one LLAGN, as well as a detection of the subthreshold population using a stacking technique. We find that the signal from the subthreshold sample is consistent with being dominated by star formation activity, although a contribution from compact jets or a mixed contribution from jetted and nonjetted systems is also feasible. On the other hand, the individually detected LLAGN are likely dominated by jet emission. We perform detailed spectral modeling for a subset of these sources and find that the γ-ray signal can be explained by synchrotron self-Compton radiation if the inner jet emission region is weakly magnetized, with its total energy density being strongly particle dominated and only slowly moving. With this work we also publicly release our Python-based stacking library for analyzing subthreshold source populations with the LAT based on a proven technique used in numerous studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A New Member of the Fast and Furious Family: A Relativistic and Time-variable UV Outflow in a Luminous Quasar

We report the fastest quasar outflow first detected in the ultraviolet, via variable C IV and Si IV absorption at outflow velocities −77,000 km s −1 to at least −90,000 km s −1 , in the radio-quiet quasar SDSS J231854.31+243954.2 (J2318). J2318 is a weak-lined quasar in the rest-frame ultraviolet, but Gemini GNIRS spectroscopy reveals an Hα redshift of z = 2.6781 ± 0.0004. A 20 yr photometric time series shows peak-to-peak variability of 0.5 mag in the g band. The C IV outflow strengthened monotonically over three epochs spanning ∼2.2 rest-frame years. The existence of such a high-velocity outflow implies that models of quasar outflows must be able to either accelerate gas to 0.3c while still preserving C IV and Si IV ions, or enable the formation of C IV and Si IV ions in gas, which has been accelerated to 0.3c. Virial estimates reveal a black-hole mass of 1.65 × 10 9 M ⊙ , which leads to an Eddington luminosity and Eddington ratio of 2.4 × 10 47 erg s −1 and 0.45, respectively. Using very conservative assumptions, the UV-absorbing outflow alone has an estimated mass loss of >0.82 M ⊙ yr −1 and a kinetic luminosity ratio L kin /L bol ≥ 0.75%. The lower limit is just above the threshold usually cited for significant feedback on the host galaxy. Comparison to PDS 456, the only other known quasar with a UV-absorbing outflow at 0.3c, suggests that the true $\dot{M}$ and Lkin/Lbol could be up to 2 orders of magnitude larger.

Seaton, Lucas M. [York University] (ORCID:00090006↗
DOE OSTI2026

Neutrino and Electromagnetic Signatures from Superluminous Supernovae: A Case Study for SN 2017egm

Superluminous supernovae (SLSNe) are rare transients that are ∼10–100 times more luminous than ordinary stellar explosions, reaching peak optical luminosities ∼10 44 –10 45 erg s −1 . The energy source powering SLSNe remains uncertain. In this work, we explore the multiwavelength and multimessenger signatures of the scenario in which SLSNe are powered by a newly born millisecond magnetar. We model the dynamical evolution and emission from the coupled system comprised of the magnetar, wind, nebula, and supernova ejecta, consistently evaluating the pair multiplicity of the wind and nebula regions, and the bulk wind Lorentz factor governing the e + − e − injection spectra in the nebula. We compute the thermal and nonthermal electromagnetic signatures, neutrino signatures, and investigate their detection prospects. For SN 2017egm, the nearest observed SLSNe, our prediction for high-energy gamma rays matches the recent detection by the Fermi Large Area Telescope. For neutrinos, using SN 2017egm a canonical SLSNe, we find that in the era of the Vera C. Rubin Observatory, a stacking analysis with upcoming neutrino observatories can lead to 3σ detection significance of neutrino events from a population of SLSNe within a decade of operation.

Mukhopadhyay, Mainak [Fermi National Accelerator L↗
DOE OSTI2021

Placing High-redshift Quasars in Perspective: A Catalog of Spectroscopic Properties from the Gemini Near Infrared Spectrograph–Distant Quasar Survey

Here, we present spectroscopic measurements for 226 sources from the Gemini Near Infrared Spectrograph–Distant Quasar Survey (GNIRS-DQS). Being the largest uniform, homogeneous survey of its kind, it represents a flux-limited sample ($m_i$ ≲ 19.0 mag, H ≲ 16.5 mag) of Sloan Digital Sky Survey (SDSS) quasars at 1.5 ≲ $\mathcal{z}$ ≲ 3.5 with a monochromatic luminosity ($\lambda {L}_{\lambda }$) at 5100 Å in the range of 10 44 –10 46 erg s –1 . A combination of the GNIRS and SDSS spectra covers principal quasar diagnostic features, chiefly the C IV λ1549, Mg ii λλ2798, 2803, Hβ λ4861, and [O III ] λλ4959, 5007 emission lines, in each source. The spectral inventory will be utilized primarily to develop prescriptions for obtaining more accurate and precise redshifts, black hole masses, and accretion rates for all quasars. Additionally, the measurements will facilitate an understanding of the dependence of rest-frame ultraviolet–optical spectral properties of quasars on redshift, luminosity, and Eddington ratio, and test whether the physical properties of the quasar central engine evolve over cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

BASS. XXVI. DR2 Host Galaxy Stellar Velocity Dispersions

We present new central stellar velocity dispersions for 484 Sy 1.9 and Sy 2 from the second data release of the Swift/BAT AGN Spectroscopic Survey (BASS DR2). This constitutes the largest study of velocity dispersion measurements in X-ray-selected obscured active galactic nuclei (AGN) with 956 independent measurements of the Ca II H and K λ3969, 3934 and Mg I λ5175 region (3880–5550 Å) and the calcium triplet region (8350–8730 Å) from 642 spectra mainly from VLT/X-Shooter or Palomar/DoubleSpec. Our sample spans velocity dispersions of 40–360 km s 1 , corresponding to 4–5 orders of magnitude in black hole mass (M BH = 10 5.5-9.6 M ⊙ ), bolometric luminosity (L bol ~ 10 42–46 erg s -1 ), and Eddington ratio (L/L Edd ~ 10 -5 to 2). For 281 AGN, our data and analysis provide the first published central velocity dispersions, including six AGN with low-mass black holes (M BH = 10 5.5-6.5 M ⊙ ), discovered thanks to high spectral resolution observations (σ inst ~ 25 km s -1 ). The survey represents a significant advance with a nearly complete census of velocity dispersions of hard X-ray–selected obscured AGN with measurements for 99% of nearby AGN (z < 0.1) outside the Galactic plane (|b| > 10°). The BASS AGN have much higher velocity dispersions than the more numerous optically selected narrow-line AGN (i.e., ~150 versus ~100 km s -1 ) but are not biased toward the highest velocity dispersions of massive ellipticals (i.e., >250 km s -1 ). Despite sufficient spectral resolution to resolve the velocity dispersions associated with the bulges of small black holes (~10 4–5 M ⊙ ), we do not find a significant population of super-Eddington AGN. Using estimates of the black hole sphere of influence from velocity dispersion, direct stellar and gas black hole mass measurements could be obtained with existing facilities for more than ~100 BASS AGN.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Probing the Diffuse Ly α Emission on Cosmological Scales: Ly α Emission Intensity Mapping Using the Complete SDSS-IV eBOSS

Based on Sloan Digital Sky Survey Data Release 16, we have detected the large-scale structure of Lyα emission in the universe at redshifts z = 2–3.5 by cross-correlating quasar positions and Lyα emission imprinted in the residual spectra of luminous red galaxies. We apply an analytical model to fit the corresponding Lyα surface brightness profile and multipoles of the redshift-space quasar–Lyα emission cross-correlation function. The model suggests an average cosmic Lyα luminosity density of $6.6^{+3.3}_{-3.1}$ x 10 40 erg s -1 cMpc -3 , a ~2σ detection with a median value about 8–9 times those estimated from deep narrowband surveys of Lyα emitters at similar redshifts. Although the low signal-to-noise ratio prevents us from a significant detection of the Lyα forest–Lyα emission cross-correlation, the measurement is consistent with the prediction of our best-fit model from quasar–Lyα emission cross-correlation within current uncertainties. We rule out the scenario where the Lyα photons mainly originate from quasars. We find that Lyα emission from star-forming galaxies, including contributions from that concentrated around the galaxy centers and that in diffuse Lyα-emitting halos, is able to explain the bulk of the Lyα luminosity density inferred from our measurements. Ongoing and future surveys can further improve the measurements and advance our understanding of the cosmic Lyα emission field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

IceCube Search for Neutrinos Coincident with Compact Binary Mergers from LIGO-Virgo’s First Gravitational-wave Transient Catalog

Using the IceCube Neutrino Observatory, we search for high-energy neutrino emission coincident with compact binary mergers observed by the LIGO and Virgo gravitational-wave (GW) detectors during their first and second observing runs. In this paper, we present results from two searches targeting emission coincident with the sky localization of each GW event within a 1000 s time window centered around the reported merger time. One search uses a model-independent unbinned maximum-likelihood analysis, which uses neutrino data from IceCube to search for pointlike neutrino sources consistent with the sky localization of GW events. The other uses the Low-Latency Algorithm for Multi-messenger Astrophysics, which incorporates astrophysical priors through a Bayesian framework and includes LIGO-Virgo detector characteristics to determine the association between the GW source and the neutrinos. No significant neutrino coincidence is seen by either search during the first two observing runs of the LIGO-Virgo detectors. We set upper limits on the time-integrated neutrino emission within the 1000 s window for each of the 11 GW events. These limits range from 0.02 to 0.7 GeV cm -2 . We also set limits on the total isotropic equivalent energy, E iso , emitted in high-energy neutrinos by each GW event. These limits range from 1.7 × 10 51 to 1.8 × 10 55 erg. We conclude with an outlook for LIGO-Virgo observing run O3, during which both analyses are running in real time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evidence for X-Ray Emission in Excess to the Jet-afterglow Decay 3.5 yr after the Binary Neutron Star Merger GW 170817: A New Emission Component

For the first ~3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n < 0.01 cm -3 . We report on observational evidence for an excess of X-ray emission at δt > 900 days after the merger. With L x ≈ 5 × 10 38 erg s -1 at 1234 days, the recently detected X-ray emission represents a ≥ 3.2σ (Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance ≥ 3.1σ, depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥ 3.7σ. A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented

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

Bridging the Bondi and Event Horizon Scales: 3D GRMHD Simulations Reveal X-shaped Radio Galaxy Morphology

X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up ≲10% of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole (SMBH), SMBH merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic magnetohydrodynamic (GRMHD) simulation for the first time follows magnetized rotating gas from outside the SMBH sphere of influence of radius RB to the SMBH of gravitational radius Rg at the largest scale separation, RB/Rg = 103, to date. Initially, our axisymmetric system of constant-density hot gas contains a weak vertical magnetic field and rotates in the equatorial plane of a rapidly spinning SMBH. We seed the gas with small-scale 2% level pressure perturbations. Infalling gas forms an accretion disk, and the SMBH launches relativistically magnetized collimated jets reaching well outside RB. Under the pressure of the infalling gas, the jets intermittently turn on and off, erratically wobble, and inflate pairs of cavities in different directions, resembling an X-shaped jet morphology. Synthetic X-ray images reveal multiple pairs of jet-powered shocks and cavities. Large-scale magnetic flux accumulates on the SMBH, becomes dynamically important, and leads to a magnetically arrested disk state. The SMBH accretes at 2% of the Bondi rate ($\dot{M}\simeq 2.4\times {10}^{-3}{M}_{\odot }\,{\mathrm{yr}}^{-1}$ for M87*) and launches twin jets at η = 150% efficiency. These jets are powerful enough (P jets ≃ 2 × 10 44 erg s -1 ) to escape along the SMBH spin axis and end the short-lived intermittent jet state, whose transient nature can account for the rarity of XRGs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

XMM and NuSTAR Observations of an Optically Quiescent Quasar

Optically quiescent quasars (OQQs) represent a recently systematized class of infrared-luminous active galactic nuclei (AGNs) that have galaxy-like optical continua. They may represent an interesting, brief phase in the AGN life cycle, e.g., either cocooned within high-covering-factor media or indicative of recent triggering, though their nature remains unclear. Here, we present the first targeted simultaneous X-ray observations of an OQQ, our previously identified prototype, SDSS J075139.06+402811.2 at z = 0.587. The source is significantly detected over 0.5–16 keV with XMM-Newton and NuSTAR, unambiguously confirming the presence of current accretion activity. Spectral modeling yields an intrinsic luminosity L 2–10 keV ≈ 4.4 × 10 43 erg s -1 , well within the AGN regime, but underluminous relative to its infrared power. It is lightly obscured, with log N H [cm -2 ] ≈ 22.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Three-dimensional Dynamics of Strongly Twisted Magnetar Magnetospheres: Kinking Flux Tubes and Global Eruptions

Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement of ψ crit ≳ π . In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲10 43 erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

FAST Discovery of a Fast Neutral Hydrogen Outflow

Abstract In this letter, we report the discovery of a fast neutral hydrogen outflow in SDSS J145239.38+062738.0, a merging radio galaxy containing an optical type I active galactic nucleus (AGN). This discovery was made through observations conducted by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) using redshifted 21 cm absorption. The outflow exhibits a blueshifted velocity likely up to ∼−1000 km s −1 with respect to the systemic velocity of the host galaxy with an absorption strength of ∼−0.6 mJy beam −1 corresponding to an optical depth of 0.002 atv= −500 km s −1 . The mass outflow rate ranges between 2.8 × 10 −2 and 3.6M ⊙ yr −1 , implying an energy outflow rate ranging between 4.2 × 10 39 and 9.7 × 10 40 erg s −1 , assuming 100 K s< 1000 K. Plausible drivers of the outflow include the starbursts, AGN radiation, and radio jet, the last of which is considered the most likely culprit according to the kinematics. By analyzing the properties of the outflow, AGN, and jet, we find that if the Hioutflow is driven by the AGN radiation, the AGN radiation does not seem powerful enough to provide negative feedback, whereas the radio jet shows the potential to provide negative feedback. Our observations contribute another example of a fast outflow detected in neutral hydrogen and demonstrate the capability of FAST in detecting such outflows.

Astronomy & Astrophysics↗
DOE OSTI2025

The Kangaroo’s First Hop: The Early Fast Cooling Phase of EP250108a/SN 2025kg

Fast X-ray transients are a rare and poorly understood population of events. Previously difficult to detect in real time, the launch of the Einstein Probe with its Wide-field X-ray Telescope has led to a rapid expansionof the sample and allowed the exploration of these transients across the electromagnetic spectrum. EP250108a is a recently detected example linked to an optical counterpart, SN 2025kg, or “the kangaroo.” Together with a companion Letter we present our observing campaign and analysis of this event. In this letter, we focus on the early evolution of the optical counterpart over the first 6 days, including our measurement of the redshift of z = 0.17641. We compare to other supernovae and fast transients showing similar features, finding significant similarities with SN 2006aj and SN 2020bvc, and show that the source is well modelled by a rapidly expanding cooling blackbody. We show the observed X-ray and radio properties are consistent with a collapsar-powered jet that is low energy (≲10 51 erg) and/or fails to break out of the dense material surrounding it. While we examine the possibility that the optical emission emerges from the shock produced as the supernova ejecta expand into a dense shell of circumstellar material, due to our X-ray and radio inferences, we favour a model where it arises from a shocked cocoon resulting from a trapped jet. This makes SN 2025 one of the few examples of this currently observationally rare event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

EP 250108a/SN 2025kg: Observations of the Most Nearby Broad-line Type Ic Supernova Following an Einstein Probe Fast X-Ray Transient

With a small sample of fast X-ray transients (FXTs) with multiwavelength counterparts discovered to date, their progenitors and connections to γ-ray bursts (GRBs) and supernovae (SNe) remain ambiguous. Here, we present photometric and spectroscopic observations of SN 2025kg, the SN counterpart to the FXT EP 250108a. At z = 0.17641, this is the closest known SN discovered following an Einstein Probe (EP) FXT. We show that SN 2025kg’s optical spectra reveal the hallmark features of a broad-lined Type Ic SN. Its light-curve evolution and expansion velocities are comparable to those of GRB-SNe, including SN 1998bw, and two past FXT-SNe. We present JWST/NIRSpec spectroscopy taken around SN 2025kg’s maximum light, and find weak absorption due to He I 1.0830 μm and 2.0581 μm and a broad, unidentified emission feature at ∼4–4.5 μm. Further, we observe broadened Hα in optical data at 42.5 days that is not detected at other epochs, indicating interaction with H-rich material. From its light curve, we derive a 56 Ni mass of 0.2–0.6 M ⊙ . Together with our companion Letter, our broadband data are consistent with a trapped or low-energy (≲10 51 erg) jet-driven explosion from a collapsar with a zero-age main-sequence mass of 15–30 M ⊙ . Finally, we show that the sample of EP FXT-SNe supports past estimates that low-luminosity jets seen through FXTs are more common than successful (GRB) jets, and that similar FXT-like signatures are likely present in at least a few percent of the brightest Type Ic-BL SNe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Origin of Pulsed Radio Emission from Magnetars

Extended periods of radio pulsations have been observed for six magnetars, displaying characteristics different from those of ordinary pulsars. In this Letter, we argue that radio emission is generated in a closed, twisted magnetic flux bundle originating near the magnetic pole and extending beyond 100 km from the magnetar. The electron–positron flow in the twisted bundle has to carry electric current and, at the same time, experiences a strong drag from the radiation field of the magnetar. This combination forces the plasma into a “radiatively locked” state with a sustained two-stream instability, generating radio emission. We demonstrate this mechanism using novel first-principles simulations that follow the plasma behavior by solving the relativistic Vlasov equation with the discontinuous Galerkin method. First, using one-dimensional simulations, we demonstrate how radiative drag induces the two-stream instability, sustaining turbulent electric fields. When extended to two dimensions, the system produces electromagnetic waves, including superluminal modes capable of escaping the magnetosphere. We measure their frequency and emitted power and incorporate the local simulation results into a global magnetospheric model. The model explains key features of the observed radio emission from magnetars: its appearance after an X-ray outburst, wide pulse profiles, luminosities ∼10 30 erg s −1 , and a broad range of frequencies extending up to ∼100 GHz.

Astronomical simulations↗
DOE OSTI2026

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS↗