How many VHE gamma-ray binaries with young pulsars can be observed?

This paper employs population synthesis calculations to estimate the number of observable Galactic VHE gamma-ray binaries containing young pulsars, accounting for the anisotropic interaction between pulsar winds and massive star winds that influences particle acceleration, radiation, and photon absorption.

The Gist

The Cosmic Detective Story: Why Are We Missing So Many "Gamma-Ray Stars"?

Imagine the Milky Way galaxy as a massive, bustling city. In this city, there are "cosmic power plants" called Gamma-Ray Binaries. These are special duos: a super-dense, spinning dead star (a pulsar) and a giant, fiery living star (a massive star).

When these two dance around each other, they crash their powerful "winds" (streams of particles) together. This collision acts like a giant cosmic particle accelerator, shooting out high-energy gamma rays. We call these the "gamma-ray loud" systems.

The Mystery:
Astronomers have only found about a dozen of these power plants so far. But when the authors of this paper ran a computer simulation of how stars are born and die in our galaxy, the math said there should be hundreds of them.

So, where are the missing hundreds? Are they hiding? Are they broken? Or are we just looking in the wrong direction?

This paper argues that the missing stars aren't broken; they are just hiding behind a curtain of magnetic fields.

---

The Three Reasons We Can't See Them All

The authors used a "population synthesis" model (basically, a super-advanced video game that simulates the life of millions of star systems) to figure this out. Here are the three main reasons why most of these systems are invisible to us:

1. The "Spotlight" Effect (Anisotropy)

The Analogy: Imagine a lighthouse. If it shines a light in all directions, everyone on the coast sees it. But if it's a laser pointer, only the person standing exactly in the beam sees the light. Everyone else sees darkness.

The Science:
Usually, we think these star systems shoot gamma rays out in a sphere, like a balloon expanding. But this paper suggests that because the massive star has a strong magnetic field (like a giant invisible magnet), it forces the gamma rays to shoot out in a narrow, laser-like beam.

• The Result: If you aren't standing exactly in the path of that laser beam, the system looks completely dark to you. Since the beams are narrow and point in random directions, we are likely missing most of them because we aren't standing in the right spot.

2. The "Dance Floor" Problem (Orbital Geometry)

The Analogy: Imagine two dancers. One is spinning wildly (the pulsar), and the other is wearing a giant, thick, swirling skirt (the massive star's disk).

• If the spinning dancer stays far away, they never touch the skirt.
• If they get close, they crash into the skirt.

The Science:
Many of these systems have very long, oval-shaped orbits. The pulsar spends most of its time far away, where the collision is weak and the gamma rays are dim. It only gets close enough to create a bright "flare" of gamma rays for a short time when it swings near the massive star.

• The Result: If we look at the system when the pulsar is far away, we see nothing. We only catch a glimpse when they are close. If we aren't watching at the exact right moment, we miss the show.

3. The "Traffic Jam" (Absorption)

The Analogy: Imagine trying to shout a message across a crowded, foggy room. If you shout, the fog might swallow your voice before it reaches the other side.

The Science:
The massive star is incredibly bright and hot. It floods the area with light. When the gamma rays try to escape, they sometimes crash into this starlight and get absorbed or turned into something else before they can reach Earth.

• The Result: Even if the system is shining, the "fog" of the star's own light might block the gamma rays from reaching our telescopes.

---

The Big Picture: Why Does This Matter?

You might ask, "Why do we care about finding a few more invisible stars?"

The answer lies in Cosmic Rays.

• The Problem: Earth is constantly bombarded by high-energy particles (cosmic rays) coming from space. We know they exist, but we don't know exactly where they come from.
• The Theory: These "hidden" gamma-ray binaries are likely the factories that create the most powerful cosmic rays in the galaxy (up to "PeV" energies, which is a quadrillion electron volts).
• The Conclusion: If there are only a dozen visible factories, they can't produce enough cosmic rays to explain what we see on Earth. But if there are hundreds of hidden factories (as the paper predicts), they could easily explain the cosmic rain hitting our planet.

The Takeaway

The authors are saying: "Don't worry, the math works out. There are likely hundreds of these cosmic power plants in our galaxy. We just can't see them because they are like laser pointers in the dark, or they are only visible for a split second during their dance, or their light is blocked by fog."

As our telescopes get better and we look at the sky for longer periods, we expect to start finding more of these hidden gems, solving the mystery of where our galaxy's most energetic particles come from.

Technical Summary

1. Problem Statement

The origin of Galactic PeV (Peta-electronvolt) cosmic rays remains an open question. While supernova remnants are the primary candidates for cosmic rays up to PeV energies, their gamma-ray spectra typically cut off around 100 TeV. Conversely, a class of sources known as Gamma-Ray Loud Binaries (GRLBs), specifically those containing a young massive star (OB or Be-type) and a relativistic companion (pulsar or black hole), show potential to accelerate particles well above 100 TeV, potentially up to PeV energies.

However, there is a significant discrepancy between the theoretical population of such systems and the observed number. Current gamma-ray observatories have detected only about a dozen GRLBs, whereas population synthesis models suggest there should be hundreds of such systems in the Galaxy. The authors aim to resolve this "deficiency" by investigating the physical conditions (magnetic fields, orbital geometry, and anisotropy) that might render the majority of these predicted sources unobservable in Very High Energy (VHE) gamma-rays.

2. Methodology

The authors employed a multi-faceted approach combining population synthesis, magnetohydrodynamic (MHD) simulations, and particle propagation modeling:

• Population Synthesis:

  • Used a modified version of the BSE (Binary Star Evolution) code to simulate the evolution of massive binary systems in the Milky Way.
  • Assumed a constant star formation rate of $2 M_\odot \text{yr}^{-1}$ and solar metallicity.
  • Modeled the formation of Be+PSR (Be-star + Pulsar) and OB+PSR binaries through two channels: Stable Mass Transfer (SMT) and Common Envelope (CE) evolution.
  • Calculated distributions for orbital periods, eccentricities, companion masses, and pulsar spin-down luminosities ($\dot{E}$).
  • Included a module for the spin evolution of magnetized neutron stars to determine if they remain in the "ejector" (pulsar) phase.

• MHD Simulations:

  • Performed relativistic MHD (rMHD) simulations using the PLUTO code to model the interaction zone between the pulsar wind and the massive star's wind.
  • Investigated the structure of the Pulsar Wind Nebula (PWN) bubble in the presence of a strong, Gauss-range magnetic field carried by the stellar wind ($B_{sw}$).

• Particle Acceleration and Propagation:

  • Modeled Fermi-type I acceleration in the wind collision zone.
  • Simulated the propagation of accelerated leptons and protons using rMHD-PIC (Particle-in-Cell) techniques.
  • Analyzed radiative losses (synchrotron and inverse Compton) and the resulting anisotropy of the emission.

3. Key Contributions and Physical Mechanisms

A. The "Deficiency" of Observed Systems

The population synthesis predicts:

• ~100 Be+PSR binaries with spin-down luminosity $\dot{E} > 10^{33} \text{ erg s}^{-1}$.
• ~200 Be+PSR binaries with periods $< 10$ days.
• ~100 Be+PSR binaries with $\dot{E} > 10^{35} \text{ erg s}^{-1}$.
  In contrast, only a handful (e.g., PSR B1259-63, PSR J2032+4127, LSI 61 303) are observed. The paper argues this is not due to a lack of formation, but rather observational bias caused by physical mechanisms.

B. The Role of Strong Magnetic Fields

The authors posit that the stellar winds of massive stars carry strong magnetic fields (Gauss range, $B \sim 1-10$ G) at orbital distances of $\sim 1$ AU.

• Acceleration: These fields are crucial for confining particles, allowing acceleration to PeV energies for protons and multi-TeV for leptons.
• Anisotropy: The strong magnetic field prevents the PWN bubble from expanding isotropically. Instead, it elongates the bubble along the magnetic field lines, creating a "magnetic jet" or tube.

C. Anisotropic Emission and Viewing Angles

This is the central finding of the paper. The VHE emission is highly anisotropic due to two factors:

1. Geometric Anisotropy: The PWN bubble is elongated along the stellar wind's magnetic field ($B_{sw}$). Particles emit gamma-rays primarily along their velocity vectors, which are aligned with $B_{sw}$. Consequently, the VHE flux is concentrated in a narrow cone (opening angle $\chi_0 \sim 5^\circ - 11^\circ$).
2. Orbital Geometry:
   • Be-disks: The equatorial disks of Be-stars are often inclined relative to the orbital plane. The pulsar only passes through the dense, magnetized disk (where acceleration is most efficient) for a fraction of the orbit.
   • Eccentricity: Highly eccentric orbits mean the pulsar spends most of its time far from the star, where magnetic fields are weaker and acceleration is inefficient.

Result: A system may be a powerful PeV accelerator, but if the observer's line of sight does not align with the magnetic jet axis during the active phase, the source appears faint or invisible in VHE gamma-rays.

D. Absorption Effects

The paper also notes that $\gamma-\gamma$ absorption (interaction with stellar photons) is highly geometry-dependent. The anisotropic emission pattern can either enhance or mitigate absorption depending on the viewing angle relative to the star and the magnetic field.

4. Results

• Predicted Population: The model predicts a Galactic population of ~100 Be+PSR binaries capable of accelerating particles to PeV energies with $\dot{E} > 10^{33} \text{ erg s}^{-1}$.
• Flux Estimates: For a typical system (e.g., parameters similar to PSR J2032+4127), the VHE flux can be as high as $F > 10^{-12} \text{ erg cm}^{-2} \text{ s}^{-1}$, but only if the observer is within the narrow emission cone ($\chi_0 \approx 11^\circ$ for 1 TeV particles).
• Observability: Due to the narrow emission cone and the specific orbital phases required (periastron passage through the disk), the "duty cycle" for detection is extremely low. Most of the predicted ~100 systems are likely missed by current surveys because they are viewed from "off-axis" angles.
• Cosmic Ray Contribution: While individual weaker sources ($\dot{E} \sim 10^{33} \text{ erg s}^{-1}$) are faint, a population of ~100 such sources could collectively contribute significantly to the Galactic PeV cosmic ray flux, potentially explaining the "knee" in the cosmic ray spectrum better than a few rare, extremely powerful sources.

5. Significance

• Resolving the Population Discrepancy: The paper provides a robust physical explanation for why the observed number of GRLBs is so low compared to theoretical predictions: strong anisotropy induced by stellar magnetic fields.
• PeVatron Candidates: It reinforces the hypothesis that GRLBs are viable Galactic PeVatrons (sources of PeV cosmic rays), provided they possess strong magnetic fields in the wind collision zone.
• Future Observations: The findings suggest that future VHE observatories (like CTA or LHAASO) need to monitor these systems over long periods and consider orbital phase and geometry to detect the "hidden" population. It also implies that the lack of isotropic PeV cosmic ray anisotropy ($<1\%$) supports a model with many distributed, weaker sources rather than a few dominant ones.
• Magnetic Field Constraints: The study highlights the necessity of strong (Gauss-range) magnetic fields in massive star winds, a parameter that can be tested via future spectropolarimetric observations (e.g., VLT measurements of LS 5039).

In conclusion, the authors argue that the "missing" VHE gamma-ray binaries are not missing in existence but are hidden in plain sight due to the highly anisotropic nature of their emission, dictated by the interaction of pulsar winds with magnetized stellar winds.