Constraining the Pulsar Beaming Fraction with TeV-Selected Galactic Pulsar Wind Nebulae and unidentified TeV Sources

This study estimates pulsar beaming fractions across radio, X-ray, and gamma-ray bands using TeV-selected pulsar wind nebulae and unidentified sources, revealing survey-dependent variations likely caused by selection effects and pulsar age differences, which can be reconciled within a unified framework employing a time-dependent opening angle.

The Gist

The Big Question: How Many Pulsars Are We Missing?

Imagine the Milky Way galaxy is a giant, dark city at night. Scattered throughout this city are millions of pulsars. You can think of a pulsar as a cosmic lighthouse. It spins incredibly fast and shoots out powerful beams of light (radio waves, X-rays, and gamma rays) from its poles.

However, there's a catch: You only see the lighthouse if you are standing in the path of its beam. If the beam is pointing away from Earth, the pulsar is invisible to us, even though it's right there.

Astronomers call the fraction of the sky covered by these beams the "beaming fraction."

• If the beam is a wide floodlight, the fraction is high (we see many).
• If the beam is a narrow laser pointer, the fraction is low (we miss many).

Knowing this fraction is crucial. If we only count the pulsars we see, we might think there are only 4,000 of them. But if the beams are narrow, the real number could be 20,000 or more. This matters because pulsars are the "parents" of double neutron stars, which crash together to create gravitational waves (the ripples in space-time that LIGO detects). If we underestimate the number of pulsars, we might be wrong about how often these cosmic crashes happen.

The Problem: Old Estimates Were Guessing

For years, scientists tried to guess this "beaming fraction" by building complex computer models of how pulsars work. But these models are like trying to guess the shape of a hidden object by feeling it in the dark—they depend heavily on assumptions that might be wrong.

The New Trick: Looking at the "Glow" Instead of the "Light"

This paper proposes a clever new way to solve the puzzle. Instead of looking at the pulsar's beam directly, the authors looked at the Pulsar Wind Nebula (PWN).

The Analogy:
Imagine a lighthouse standing in the middle of a foggy lake.

1. The Beam: The light shooting out of the lighthouse (the pulsar beam). You only see it if you are in the right spot.
2. The Fog: The lighthouse also blows a giant, glowing fog bank around itself (the PWN). This fog glows in high-energy "TeV" light.

The Key Insight:
The fog (PWN) glows in all directions. It doesn't matter which way the lighthouse beam is pointing; the fog is visible to everyone.

• If we see the fog AND the beam, we know the beam is pointing at us.
• If we see the fog but no beam, we know the beam is pointing away from us.

By counting how many "fogs" have visible beams versus how many don't, the scientists can calculate exactly how narrow the beams are, without needing to guess the physics of the lighthouse itself.

The Experiment: Three Different Cameras

The team used data from three different "cameras" (telescopes) that look at the sky in TeV light:

1. H.E.S.S. (High Energy Stereoscopic System)
2. HAWC (High Altitude Water Cherenkov Observatory)
3. LHAASO (Large High Altitude Air Shower Observatory)

They counted the "fogs" (PWNe) and the "missing beams" (Unidentified sources) for each camera.

The Surprising Result: It Depends on the Camera

Here is where it gets interesting. The three cameras didn't agree on the answer!

• H.E.S.S. (The "Zoom Lens"): This camera has very sharp vision but a narrow field of view. It found that the beams are wider (about 23–36% of the sky).
  • Why? It's really good at spotting young, bright, compact "fogs" near the center of the galaxy. It missed the old, faint, spread-out ones.
• HAWC & LHAASO (The "Wide-Angle Lens"): These cameras see a huge chunk of the sky but have slightly fuzzier vision. They found that the beams are much narrower (only about 4–13% of the sky).
  • Why? These cameras are great at spotting old, giant, spread-out "fogs" that have been expanding for thousands of years. They found a lot of "missing beams" in these older systems.

The Metaphor:
Imagine you are trying to count how many people are wearing red hats in a stadium.

• H.E.S.S. is like a person with binoculars looking at the front row. They see the bright red hats clearly and think, "Wow, lots of red hats!"
• HAWC/LHAASO is like a person with a wide-angle camera looking at the whole stadium, including the back rows. They see that in the back rows, the red hats are faded and hard to spot, so they think, "Actually, there are very few red hats."

Both are right about what they see, but they are looking at different groups of people.

The Solution: The Beam Shrinks with Age

The authors realized the discrepancy wasn't a mistake; it was a clue. They used a computer simulation to show that pulsar beams get narrower as the pulsar gets older.

• Young Pulsars (H.E.S.S. sample): The beams are wide open (like a floodlight).
• Old Pulsars (HAWC/LHAASO sample): The beams have shrunk down to a narrow laser pointer.

This explains why H.E.S.S. (seeing young ones) got a high number, and HAWC (seeing old ones) got a low number.

Why This Matters

1. Fixing the Census: We now have a better way to estimate the true number of pulsars in the galaxy. It's not just a guess; it's based on the "fog" they leave behind.
2. Gravitational Waves: This helps us predict how often neutron stars crash into each other. If beams are narrow, there are more hidden pulsars, which means more potential crashes waiting to happen.
3. Future Telescopes: The paper mentions a future telescope called CTAO. It will be the ultimate camera—sharp enough to see the details and wide enough to see the whole stadium. It will finally settle the debate and give us the perfect count.

In a Nutshell

The universe is full of cosmic lighthouses. We can't see them all because their beams are narrow. By looking at the glowing "fog" (nebula) they create, we can count the ones we miss. We discovered that these beams start wide when the lighthouses are young and shrink as they age. This helps us understand the true population of these mysterious stars and predict the cosmic collisions that shake the fabric of space-time.

Technical Summary

1. Problem Statement

The pulsar beaming fraction ($f_b$) is a critical parameter defined as the fraction of the full solid angle swept by a pulsar's emission beam during one rotation. It quantifies the probability that a randomly oriented observer's line of sight intersects the beam. Accurate determination of $f_b$ is essential for:

• Reconstructing the intrinsic Galactic pulsar population from observed samples (correcting for "missing" pulsars).
• Constraining theoretical models of pulsar emission geometry (e.g., polar cap vs. outer gap models).
• Estimating the merger rate of double neutron stars (DNS) and the corresponding gravitational-wave detection rates.

Current estimates of $f_b$ rely heavily on Monte Carlo (MC) simulations that assume specific beam geometry models (e.g., polar cap for radio, outer gap for $\gamma$-rays). These model-dependent approaches yield values ranging from $\sim 0.1$ to $0.5$ and often show discrepancies between radio and $\gamma$-ray populations. Furthermore, recent predictions for DNS merger rates based on these fractions conflict with the lack of confirmed DNS merger detections by LIGO/Virgo, suggesting systematic uncertainties in $f_b$ may be larger than previously thought.

2. Methodology

The authors propose a geometry-independent method to estimate $f_b$ using TeV-selected Pulsar Wind Nebulae (PWNe) and Unidentified (Unid) TeV sources.

• Core Assumption:
  1. Isotropy: TeV emission from PWNe is approximately isotropic (unlike the beamed emission of the pulsar itself) because it arises from inverse-Compton scattering of relativistic electrons distributed throughout the nebula.
  2. Sample Definition:
     • PWNe Sample: TeV sources where the associated pulsar beam does intersect our line of sight (detected in radio, $\gamma$-ray, or X-ray).
     • Unid Sample: TeV sources in the Galactic plane where no pulsar is detected. These are interpreted as PWNe powered by pulsars whose beams do not intersect our line of sight.
• Calculation:
  The beaming fraction is estimated as the ratio of the number of detected pulsars ($N_{\nu, \text{PWN}}$) to the total number of TeV-selected sources (PWNe + Unid):
  $$f_{b,\nu} = \frac{N_{\nu, \text{PWN}}}{N_{\text{PWN}} + N_{\text{Unid}}}$$
  This is calculated independently for three major TeV surveys: H.E.S.S., HAWC, and LHAASO, and for different observational bands (Radio, $\gamma$-ray, X-ray).
• Data Sources:
  • TeV sources from TeVCat (Wakely & Horan 2008).
  • Pulsar catalogs: ATNF (radio), Fermi-LAT ($\gamma$-ray), and X-ray catalogs.
  • The study excludes sources with ambiguous associations (e.g., $\gamma$-ray binaries) and treats candidates uniformly as PWNe to maximize sample size.
• Validation via MC Simulation:
  To explain discrepancies between surveys, the authors perform a population synthesis simulation (following Shi & Ng 2024). They generate a synthetic pulsar population with evolving spin periods, magnetic fields, and inclination angles, applying survey-specific selection effects (sensitivity, angular resolution) to reproduce the observed $f_b$ values.

3. Key Contributions

• Novel Tracer: Utilizes TeV-selected PWNe as a tracer for the total pulsar population (both beamed and un-beamed), bypassing the need for explicit assumptions about the pulsar's emission geometry.
• Survey-Specific Analysis: Instead of combining all TeV data, the study analyzes H.E.S.S., HAWC, and LHAASO separately to isolate survey-dependent selection effects (angular resolution, energy range, sky coverage).
• Unified Framework: Demonstrates that the disparate $f_b$ values derived from different surveys can be reconciled within a single physical framework by introducing a time-dependent opening angle ($\rho(t)$).

4. Results

• Survey Discrepancy:
  • H.E.S.S.: Yields a higher beaming fraction, $f_b \sim 0.23 - 0.36$ (for $\dot{E} > 10^{36}$ erg/s).
  • HAWC/LHAASO: Yields significantly lower values, $f_b \sim 0.04 - 0.13$ (for $\dot{E} > 10^{34.5}$ erg/s).
  • The discrepancy is a factor of 2–4, which is statistically significant even after accounting for Poisson noise.
• Band Independence: Within a single survey, $f_b$ is consistent across radio, $\gamma$-ray, and X-ray bands. This suggests that for young pulsars associated with TeV PWNe, the emission geometries of different bands are correlated.
• Physical Interpretation of Discrepancy:
  • Angular Resolution & Energy Range: H.E.S.S. has high angular resolution but a narrow field of view, making it sensitive to compact, young, high-energy PWNe. HAWC and LHAASO have wide fields of view and are sensitive to extended, older, lower-energy sources (TeV halos).
  • Selection Bias: HAWC/LHAASO samples are biased toward older pulsars ($\tau_c \gtrsim 10^5$ yr) with extended nebulae. H.E.S.S. samples are biased toward younger pulsars.
• MC Simulation Findings:
  • To reproduce the low $f_b$ of HAWC/LHAASO, the simulation requires a smaller beam opening angle and a larger maximum age ($t_{max}$) compared to H.E.S.S.
  • A time-dependent opening angle model, $\rho(t) = \frac{90^\circ}{1+(t/t_c)^\beta}$, successfully reproduces the observed fractions for all surveys.
  • Best-fit parameters: The opening angle decreases with age. For radio, $t_c \approx 5.75 \times 10^4$ yr; for $\gamma$-ray, $t_c \approx 7.01 \times 10^4$ yr. This implies that older pulsars (dominating HAWC/LHAASO samples) have narrower beams, leading to a lower probability of detection.

5. Significance

• Resolution of Model Conflicts: The study resolves the apparent conflict between different $f_b$ estimates by showing they are not contradictory but rather reflect different evolutionary stages of the pulsar population selected by different instruments.
• Evolution of Beam Geometry: It provides observational evidence that the pulsar beam opening angle shrinks over time (on a timescale of $\sim 10^5$ years). This supports physical models where the emission region (e.g., outer gap or current sheet) evolves or where magnetic axis alignment occurs over time.
• Implications for Gravitational Waves: The finding that older pulsars have narrower beams suggests that the intrinsic population of DNS progenitors might be larger than previously estimated if based solely on young pulsar samples. This could help reconcile the predicted vs. observed DNS merger rates.
• Future Outlook: The authors highlight that the upcoming Cherenkov Telescope Array (CTAO), with its superior sensitivity and angular resolution, will be crucial for detecting faint, unresolved PWNe and refining these beaming fraction estimates, potentially reducing the current systematic uncertainties.

In conclusion, this paper establishes that the pulsar beaming fraction is not a static constant but evolves with the age of the pulsar, and that survey characteristics significantly bias the inferred values. A unified model incorporating time-dependent beam narrowing reconciles these observations.