[論文レビュー] Relativistic astrospheres of gamma-ray binaries: modeling of non-thermal processes

A. M. Bykov, A. E. Petrov+1 more|arXiv (Cornell University)|Mar 12, 2026

Astrophysics and Cosmic Phenomena被引用数 0

ひとこと要約

The paper conducts 2D and 3D relativistic MHD simulations of pulsar wind nebulae in gamma-ray binaries to study how strong, Gauss-range stellar wind magnetic fields shape the nebula, cocoon, and relativistic clumps, and to assess PeV particle acceleration and non-thermal emission.

ABSTRACT

A long standing problem in high energy astrophysics is the nature of galactic accelerators of particles with energies above PeV. Such objects are sources of galactic cosmic rays and can produce PeV-regime photons observed by ground-based observatories. Among very likely accelerators are astrospheres of pulsars in gamma-ray binaries. These binaries have long been observed as bright sources of TeV gamma-rays. Recently, 2D relativistic magnetohydrodynamic (rMHD) simulations have shown that the astrospheres can accelerate particles to energies well above PeV, provided that they harbor a Gauss-range magnetic field. Such a strong field is necessary in the region of two colliding winds: the relativistic outflow of the pulsar or accreting black hole and the wind of its stellar companion, a massive early-type star. Here, the wind collision region is explored as the site of PeV protons acceleration. The local structure of colliding flows is illustrated using rMHD simulations of a powerful pulsar wind in 2D and 3D models. The relativistic outflow of a pulsar or black hole, evolving inside the strongly magnetized stellar wind, have an elongated shape and surrounded by a kind of magnetic cocoon providing favorable conditions for acceleration of ultra high energy ions. The simulated spectra of particles, accelerated by intermittent relativistic turbulence in these systems, have piece-wise power-law shape and extend well above PeV energies for powerful outflows. The model indicated that gamma-ray binaries harboring a powerful relativistic outflow, produced either by a pulsar or accreting black hole, can be bright sources of synchrotron MeV-regime photons and multi-PeV regime gamma-rays, as recently detected from galactic microquasars like Cyg X-3. The Gauss-range magnetic field of a massive star wind strongly influences the non-thermal emission of gamma-ray binaries with relativistic companions.

研究の動機と目的

• Motivate and quantify how strong magnetic fields in massive-star winds affect particle acceleration and non-thermal emission in gamma-ray binaries.
• Characterize the local wind collision region structure and its evolution under Gauss-range stellar wind fields.
• Assess the validity of planar (2D) models for capturing key physics compared to full 3D simulations.
• Identify conditions under which relativistic clumps form and contribute to ultra-high-energy particle acceleration.

提案手法

• Perform self-consistent relativistic MHD–PIC modeling of the wind collision region in gamma-ray binaries with prescribed strong stellar wind magnetic fields.
• Compare planar 2D and full 3D rMHD simulations across a range of B_sw values (sub-Gauss to Gauss) to study nebula morphology and magnetic cocoon formation.
• Analyze the formation and role of relativistic clumps in particle acceleration using Monte-Carlo modeling (Appendix B).
• Estimate how magnetic pressure and ram pressure balance sets nebula shape through equation B_sw ≈ sqrt(2 E_dot / (u_pwn r^2)).
• Investigate the impact of wind anisotropy and field-velocity alignment on shock structure and magnetization.
• Explore how clumps with Lorentz factors Γ>3 can upscatter pre-accelerated particles to PeV energies.

実験結果

リサーチクエスチョン

• RQ1 Does planar (2D) modeling adequately reproduce the key features of 3D relativistic astrospheres in gamma-ray binaries?
• RQ2 How does the stellar wind magnetic field strength (B_sw) influence the nebula morphology, magnetic cocoon formation, and acceleration sites for PeV-scale protons?
• RQ3 Can relativistic clumps in the nebula provide efficient re-acceleration to energies above PeV, and under which conditions do they form?
• RQ4 What are the observable implications for non-thermal emission (synchrotron MeV photons and multi-PeV gamma-rays) in systems with powerful relativistic outflows?

主な発見

Model	α	σ0	n_sw, cm^-3	B_sw, G	ψ, deg
1	45°	0.03	300	0.1	90
2	45°	0.03	300	0.1	90
3	45°	0.1	300	0.1	90
4	45°	1	300	0.1	90
5	85°	0.03	300	0.1	90
6	85°	0.03	300	0.1	90
7	45°	0.03	3000	0.4	45
8	45°	0.03	30000	0.4	45
9	45°	0.03	30000	0.1	45
10	45°	0.03	30000	0.1	45
11	45°	10	30000	0.4	45

• The 2D and 3D nebulae structures are similar across tested models, especially for B_sw in the 2–3 G range, indicating planar models can capture essential physics in these regimes.
• A dense magnetic cocoon forms around the nebula, enhancing confinement and enabling PeV-scale proton acceleration within short timescales (< 10^5 s, much shorter than orbital periods).
• Increasing B_sw elongates the nebula along the magnetic field and strengthens the magnetic cocoon, while also reducing the overall nebular volume and shifting the termination shock closer to the pulsar.
• Relativistic clumps (Γ ≳ 2–6 in different models) arise in 2D and 3D simulations, with sub-AU to per-cent-of-AU scales; these clumps can provide efficient re-acceleration, producing spectra with a second PeV-energy hump when Γ≳3.
• A Monte-Carlo treatment shows that embedding clumps with Γ>3 in colliding flows can yield a second high-energy spectral component beyond PeV, significantly hardening the particle spectrum.
• The clumps persist for hours and can be larger than 1 AU in some configurations, enabling substantial re-acceleration of pre-accelerated particles in Gauss-range fields.
• 3D runs exhibit more uniform pressure and less turbulence than 2D, but still produce fast-moving clumps; their properties depend on grid resolution and magnetization σ0 of the pulsar wind.