Relativistic jets from millisecond proto-magnetars

arXiv:2601.07918 · astro-ph.HE, gr-qc · Submitted 2026-01-12 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Today's paper: "Relativistic jets from millisecond proto-magnetars".

Jocelyn: This study presents three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations of neutrino-heated winds from rapidly rotating,

Vera: First, who's behind it and why it matters.

Title and authors: Vera: So, what did the researchers actually find when they ran these simulations? They summarized their findings by looking at how the outflows are structured and how they partition energy between different components.

Jocelyn: I see that the main finding revolves around a multidimensional effect: centrifugal forces at lower latitudes create a dense, sub-relativistic outflow, which in turn acts to confine and direct any less baryon-loaded outflows emerging from higher latitudes.

Subrahmanyan: That structure is key because it leads to the formation of a structured bipolar jet with its peak magnetization occurring along the spin axis up to about sigma ∼ thirty to one hundred which they suggest is sufficient for reaching bulk Lorentz factors of Γ∞ around one hundred on larger scales.

Vera: That's a very specific detail about the magnetization levels, Subrahmanyan. It shows how the magnetic field gets concentrated where it needs to be for a jet to form effectively.

Jocelyn: I also found that this entire process of collimation happens intrinsically at small radii, even without needing an external envelope of ejecta to provide that confinement, which is a neat point.

Subrahmanyan: That internal collimation mechanism, where the heavy equatorial wind channels the magnetic flux toward the rotation axis, is a crucial piece of physics because it bypasses one potential obstacle for jet formation.

Vera: It seems they are showing that these proto-magnetars naturally produce an anisotropic structure that matches what we expect from short gamma-ray bursts in terms of opening angles.

Jocelyn: And the results show the energy partition, indicating that the jet releases a significant amount of energy, with a rate estimated around E˙ ≈ E˙EM ∼ one thousand forty-nine to one thousand fifty erg s−one in material with Lorentz factors from ten down to ten zero.

Subrahmanyan: That energy release estimate suggests that the resulting outflow is energetic enough and clean enough in principle to power a short-duration gamma-ray burst, which really supports the idea of a magnetar powering the prompt emission phase.

The paper's summary: Vera: Now that we know what they found, I think it’s interesting to consider how this work pushes things forward in terms of future research or refinement. What are the suggested improvements for these types of simulations?

Jocelyn: The paper suggests moving beyond just looking at the steady-state properties and incorporating more dynamic aspects, specifically exploring the full time-dependent evolution during the cooling and spin-down phases of these proto-magnetars.

Subrahmanyan: That’s a necessary step because a proto-magnetar isn't static; its magnetic field configuration is constantly changing as it cools and spins down, so a time-dependent model would give us a more complete picture of the outflow over time.

Vera: I also noticed they mention the need to incorporate self-consistent remnant magnetic field configurations into future studies, which suggests that the next step is ensuring these simulations aren't just using simplified starting points for those fields.

Jocelyn: And looking ahead, incorporating a model for the interaction between these magnetar winds and any surrounding accretion disks could provide even richer data on how energy is actually transferred in a real astrophysical environment.

Subrahmanyan: From my side, I think the next big challenge is integrating these results with other components of the explosion physics, like modeling how this jet interacts with the surrounding medium to see if it produces the expected afterglow signatures.

Vera: So, it seems like they're setting up a path for linking these initial conditions to later observational evidence by focusing on those more complex, time-dependent interactions.

Jocelyn: It’s interesting because they are trying to bridge the gap between the initial formation of the object and the observable high-energy phenomena we see later.

The paper's improvements: Vera: So, to wrap up our discussion on "Relativistic jets from millisecond proto-magnetars," it seems this paper provides strong support for a scenario where these rapidly rotating, strongly magnetized objects can indeed launch relativistic jets within seconds of formation.

Jocelyn: I think the main implication is that this offers a physical mechanism explaining the coexistence of luminous transients and gamma-ray bursts, specifically by showing how these proto-magnetars naturally produce both ultra-relativistic polar outflows and sub-relativistic equatorial components.

Subrahmanyan: Precisely; it shows that the physics governing rapid rotation and strong magnetic fields creates an intrinsic anisotropy that satisfies the requirements for ultra-relativistic motion needed for short GRBs, linking them firmly to compact object explosions.

Vera: It’s a solid piece of work that connects the detailed simulation results about magnetization and energy partition directly to the observed characteristics of these events.

Jocelyn: It gives us a clearer picture of how we might interpret multi-messenger signals by understanding this internal structure, which is really useful for pulsar surveys and gravitational wave follow-up.

Subrahmanyan: I just want to say that the work lays a foundation for understanding the immediate aftermath of these events, setting up a path where future detailed studies can probe the full evolution of these magnetars.

Vera: That's exactly what it is, Subrahmanyan; a solid starting point for understanding this fascinating area of astrophysics.

Jocelyn: We're really looking forward to seeing how this framework helps us interpret the next set of observational data we get from our telescopes and detectors.

Conclusion: Vera: So, to wrap up our discussion on "Relativistic jets from millisecond proto-magnetars," this paper confirms that these objects can launch relativistic jets within seconds of forming, which is really exciting news for short GRBs.

Jocelyn: I think the main implication is that we have a new physical mechanism for prompt emission in those gamma-ray bursts, and it’s tied directly to the physics of compact object collapse.

Subrahmanyan: Indeed, it suggests that millisecond proto-magnetars provide a natural source for ultra-relativistic outflows, which helps us explain why we see so much energy in these explosions.

Vera: It’s really neat how the simulation showed that the structure—the dense equatorial wind confining the jet—actually matches what we infer from observational constraints on jet opening angles.

Jocelyn: That connection between the internal dynamics of a proto-magnetar and observable signatures is exactly why this research is so important for us in pulsar surveys.

Subrahmanyan: I think the link to nucleosynthesis, where this outflow dictates the chemical fingerprint of a kilonova, is where things get really interesting for cosmology.

Vera: It’s true; understanding that energy partition gives us a way to predict what kind of heavy elements we should expect from these magnetar-powered transients.

Jocelyn: And knowing how much energy is actually channeled into the jet material helps us narrow down the possible progenitor scenarios in our data analysis.

Subrahmanyan: Ultimately, this paper solidifies a pathway where we can move from just observing these events to understanding the underlying physics of the engine itself.

Vera: That’s a great summary of what makes this study so compelling, Subrahmanyan; it really ties together the dynamics and the potential observational consequences.

Jocelyn: I think we should keep an eye on how this framework helps us interpret future data from gravitational waves and electromagnetic telescopes in the context of these compact object explosions.

Subrahmanyan: Absolutely, and I'm eager to see how these simulation results inform our next theoretical models on compact object evolution under extreme conditions.

Dhruv K. Desai, Luciano Combi, Daniel M. Siegel, Brian D. Metzger

Institute of Physics, University of Greifswald · Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada · Department of Physics, University of Guelph, Guelph, Ontario N1G 2W1, Canada · Department of Physics and Columbia Astrophysics Laboratory, Columbia University · Center for Computational Astrophysics, Flatiron Institute

astro-ph.HE, gr-qc

Submitted: 2026-01-12

Updated: 2026-09-30

Comments: 16 pages, 6 figures. v2: accepted version, published in ApJL

Journal ref: ApJL 1009, L49 (2026)

DOI: 10.3847/2041-8213/aea3e9

License: http://creativecommons.org/licenses/by-sa/4.0/

Importance score: 86/100

The gist: This study presents three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations of neutrino-heated winds from rapidly rotating, strongly magnetized millisecond proto-magnetars to

Key concepts

Proto-Magnetar: These are hot, dense remnants formed in events like neutron star mergers or supernovae. They can possess extremely high rotational energy and magnetic fields (up to 10^15 G), making them potential central engines for powerful explosions like gamma-ray bursts.
GRMHD Simulation: This is a complex computer model that solves the equations governing gravity, fluid motion, and electromagnetism in a relativistic setting. It tracks how matter and magnetic fields interact dynamically over time to see if they can produce jets.
Jet Collimation: The paper found that jets narrow themselves intrinsically at very small radii (less than 100 km) due to the dense, fast wind near the equator. This equatorial outflow acts like a natural funnel, channeling energy and magnetic flux toward the rotation axis to form a focused polar jet.
Peak Magnetization ($\sigma$): This measures how strong the magnetic field is relative to other energy components in a jet. The simulation found that the polar jet reaches peak magnetization values of 30 to 100, which is necessary for the outflow to achieve ultra-relativistic speeds required for short GRBs.

Terminology

Core — 8× in this paper · explained in 2 episodes
Boltzmann Equation — 2× in this paper · explained in 2 episodes
Anisotropy — 1× in this paper · explained in 3 episodes
Initial Data — 1× in this paper · explained in 2 episodes
alpha — 1× in this paper · explained in 1 episode

Summary

This study presents three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations of neutrino-heated winds from rapidly rotating, strongly magnetized millisecond proto-magnetars to determine if these objects can self-consistently launch relativistic jets within the first few seconds after their formation. The research is significant because it addresses a key question regarding the central engines of gamma-ray bursts (GRBs), suggesting that millisecond proto-magnetars can indeed power ultra-relativistic outflows, thereby providing a physical mechanism for the prompt emission phase of short GRBs and linking these events to the broader class of compact-object explosions.

Proto-Magnetar Formation and Context

Proto-neutron stars (PNSs) are hot, dense remnants formed in core-collapse supernovae, neutron star mergers, and white dwarf accretion-induced collapse. These objects can contain significant rotational energy, which dynamo processes may amplify into magnetic fields of the order ≳ 1015 G. These millisecond proto-magnetars are linked to r-process heavy element production and are invoked as power sources behind superluminous supernovae (SLSNe) and gamma-ray bursts (GRBs). The paper focuses on whether these objects can launch ultra-relativistic jets at early times, as needed for short GRBs.

Simulation Methodology

The researchers employ a three-dimensional general-relativistic MHD code based on the Einstein Toolkit. Key numerical features include:

Solving ideal GRMHD equations with a finite-volume scheme using WENO-Z reconstruction and the approximate HLLE Riemann solver for hydrodynamic variables.
Evolving the magnetic field through the electromagnetic 4-potential in a generalized Lorenz gauge using an upwind constraint transport scheme to maintain solenoidal constraint to machine precision.
Using a one-moment (M0) approximation of the general-relativistic Boltzmann equation for neutrino transport, adapted from D. Radice et al. (2016), evolved along null coordinate radial rays on a semi-spherical grid with 1300 × 25 × 50 grid points.
Using the SFHo nuclear EOS and extending it to low rest-mass densities, and imposing density floors and magnetization ceilings to maintain numerical stability in fast, highly magnetized outflows.

Wind Structure and Angular Stratification

The simulations reveal an essential multidimensional effect: centrifugal forces strongly enhance mass loss near the rotational equator, producing a dense, sub-relativistic outflow (v ∼ 0.1c). This equatorial wind naturally confines and collimates less baryon-loaded outflows emerging from higher latitudes. This leads to the formation of a structured bipolar jet with a peak magnetization along the pole up to σ ∼ 30 − 100, sufficient to reach bulk Lorentz factors Γ∞ ∼ 100 on larger scales.

Jet Collimation Mechanism

A key finding is that jet collimation arises intrinsically at small radii (≲ 100 km), without requiring confinement by an external envelope of supernova or merger ejecta. The heavy equatorial wind acts as an effective confining medium, channeling magnetic flux and energy toward the rotation axis. This mechanism produces a narrow polar outflow whose energetics and opening angles are broadly consistent with those inferred for short gamma-ray bursts, with a jet half-opening angle of ≃ 8◦.

Implications for GRB Power

The results demonstrate that sufficiently rapidly rotating (P ∼ 1 ms), highly magnetized proto-magnetars (Bp ≳ 1015 G) naturally produce a strongly anisotropic outflow structure. The polar jet reaches peak magnetization σ ∼ 30 − 100, satisfying the requirement for ultra-relativistic motion. The resulting energy partition shows that the proto-magnetar jet releases E˙ ≈ E˙EM ∼ 1049 − 1050 erg s−1 in material with Γ∞ ∼ 10 − 100, in principle sufficiently luminous and clean to power a short-duration gamma-ray burst. This supports a scenario where a brief protomagnetar phase can power the prompt emission of short GRBs within seconds of remnant formation.

Conclusion and Future Work

The simulations support the conclusion that millisecond proto-magnetars can launch relativistic jets within seconds of formation. The resulting partition of energy between ultra-relativistic and sub-relativistic components provides a natural explanation for the coexistence of GRB emission and luminous, magnetarpowered transients in compact-object explosions. Future work will focus on incorporating self-consistent remnant magnetic field configurations, full time-dependent evolution during cooling and spin-down, and the interaction between magnetar winds and surrounding accretion disks.

Key Results Summary:

(Note: The paper enumerates results implicitly through the description of findings rather than a numbered list for all points.)

Improvements for AI systems

As a fastidious researcher, I have analyzed this manuscript, Relativistic jets from millisecond proto-magnetars, and identified several areas where AI systems could be significantly enhanced or specialized using these findings.

Here are the specific improvements and the resulting capabilities of an improved AI system:

)1. Enhanced Physics Simulation & Modeling

The paper presents complex 3D General Relativistic Magnetohydrodynamic (GRMHD) simulations incorporating neutrino transport (M0 approximation).

Improvement: Develop a specialized AI model trained on the output data from these GRMHD simulations (density, velocity profiles, magnetization maps across different spin/field parameter spaces).
Improved AI Capability: The system could perform inverse modeling or parameter inference. Given observational constraints (e.g., inferred jet half-opening angles, observed energy partitions), the AI could rapidly predict the required initial conditions (spin period P and magnetic field strength Bp) of a proto-magnetar that would produce those specific jet structure. This moves beyond simple forward modeling to constraint discovery.

)2. Multi-Messenger Data Interpretation & Correlation

The paper explicitly links the resulting outflow partitioning (ultra-relativistic polar vs. sub-relativistic equatorial components) to observed electromagnetic counterparts (GRB energies, supernova/kilonova ejecta).

Improvement: Build a deep learning model specifically designed to map the internal structure of a simulated magnetar wind (e.g., its angular stratification and energy partition shown in Figure 3 and 4) directly onto observable multi-messenger signatures.
Improved AI Capability: The system could ingest raw data from gravitational wave events (like GW170817), electromagnetic transients (GRB light curves), and kilonova spectra, and use the learned mapping to instantly determine if the observed energy partitioning is consistent with a jet launched by a millisecond proto-magnetar versus other compact object engines.

)3. Nucleosynthesis Yield Prediction

The paper notes that the outflow structure dictates nucleosynthetic yields (e.g., neutron-rich ejecta powering kilonova emission).

Improvement: Train a physics-informed neural network (PINN) using the derived relationships in Equations (B5) and (B6) to predict the elemental abundance patterns of matter ejected by these specific magnetar outflows, accounting for magnetic field geometry.
Improved AI Capability: The system could predict the chemical fingerprint of a kilonova based on whether it was powered by a magnetar jet versus standard core-collapse material, allowing astrophysicists to distinguish progenitor types based on the resulting heavy element synthesis signature.

)4. High-Dimensional Parameter Space Exploration

The study explores a high-dimensional parameter space (BP, P, Rν¯e).

Improvement: Implement an advanced Bayesian optimization or reinforcement learning agent specifically tailored for navigating the complex constraints of GRMHD simulations to find regions of high physical interest (e.g., the region where jet magnetization σj is maximized).
Improved AI Capability: The AI could autonomously search vast, computationally expensive simulation spaces to identify novel physical regimes—such as the optimal combination of rotation and magnetic field strength—that maximize jet Lorentz factors while minimizing baryon loading, effectively discovering new physical constraints on compact object evolution that human intuition might miss.

)5. Automated Literature Synthesis and Hypothesis Generation

The paper synthesizes results from various previous models (1D vs. 3D, different initial conditions).

Improvement: Create a Large Language Model (LLM) fine-tuned exclusively on the text, equations, and results of this paper and its cited literature to perform highly specialized cross-referencing.
Improved AI Capability: When presented with a new simulation result or an observational anomaly, the system could instantly synthesize relevant conclusions from this paper regarding jet collimation mechanisms, equatorial wind enhancement, or energy partitioning ratios, providing immediate, context-specific scientific commentary and identifying potential gaps in current theoretical models.

Abstract

Rapidly rotating, strongly magnetized neutron stars (``millisecond proto-magnetars'') formed in stellar core-collapse, neutron star mergers, and white dwarf accretion-induced collapse have long been proposed as central engines of gamma-ray bursts (GRB) and accompanying supernovae/kilonovae. However, during the first few seconds after birth, neutrino heating drives baryon-rich winds from the neutron star surface, potentially limiting the magnetization and achievable Lorentz factors of the outflow and casting doubt on whether proto-magnetars can launch ultra-relativistic jets at early times, as needed to power short-duration GRB. We present 3D general-relativistic magnetohydrodynamic simulations of neutrino-heated proto-magnetar winds that incorporate M0 neutrino transport. While the global wind properties broadly agree with previous analytic estimates calibrated to one-dimensional models, our simulations reveal essential multidimensional effects. For rapidly rotating models with spin periods P = 1 ms, centrifugal forces strongly enhance mass loss near the rotational equator, producing a dense, sub-relativistic outflow (0.1c). This equatorial wind naturally confines and collimates less baryon-loaded outflows emerging from higher latitudes, leading to the formation of a structured bipolar jet with a peak magnetization up to 30-100 along the pole, sufficient to reach bulk Lorentz factors 100 on larger scales. The resulting angular stratification of the outflow energy into ultra-relativistic polar and sub-relativistic equatorial components is broadly consistent with the observed partition between beaming-corrected GRB energies and supernova/kilonova ejecta. Our results demonstrate that millisecond proto-magnetars can launch relativistic jets within seconds of formation and highlight their potential role in powering the diverse electromagnetic counterparts of compact-object explosions.