CTAO Simulations for Potential PeVatron Candidates

This study utilizes CTAO simulations with Gammapy to evaluate the observatory's ability to detect hadronic gamma-ray signatures from four PeVatron candidates, concluding that while 100 hours of observation are needed to distinguish proton cut-off energies up to 600 TeV, CTAO can likely exclude Cassiopeia A, RX J1713.7-3946, and HESS J1731-347 as PeVatron sources, leaving HAWC J2227+610 inconclusive.

The Gist

Imagine the universe is a giant, chaotic construction site. For decades, scientists have been trying to figure out where the most powerful "bricks" in the universe come from. These bricks are Cosmic Rays—tiny, super-fast particles that zoom through space at nearly the speed of light.

The big mystery is: What machine is powerful enough to build these bricks?

Scientists suspect that Supernova Remnants (SNRs)—the glowing, expanding clouds left behind after a massive star explodes—are the factories. Specifically, they think these factories might be "PeVatrons," machines capable of accelerating particles to PeV (Peta-electronvolt) energies. That's a million billion electron volts! It's like the difference between a bicycle and a rocket ship.

However, we can't see these factories directly. The particles are charged, so they get bounced around by magnetic fields in space, like a pinball losing its path. To find the factories, we have to look for the smoke: high-energy gamma rays. These are neutral messengers that travel in a straight line from the source to our telescopes.

The Problem: A Blurry Camera

Until now, our telescopes have been like old, grainy cameras. They can see the smoke, but the picture is fuzzy. We can't tell if the smoke is coming from a PeVatron (a proton factory) or just a regular electron accelerator. It's like trying to tell the difference between a Ferrari and a bicycle just by looking at a blurry photo of their exhaust fumes.

The Solution: The CTAO Super-Camera

Enter the Cherenkov Telescope Array Observatory (CTAO). Think of this as the next-generation, high-definition, 4K camera for the universe. It's a massive array of telescopes in both the Northern and Southern hemispheres, designed to see gamma rays with incredible clarity and sensitivity.

This paper is a simulation study. The authors didn't wait for the telescope to be fully built; they used powerful computers to ask: "If we point this super-camera at four specific cosmic factories, what will we see?"

The Four Suspects

The team picked four famous "suspects" (Supernova Remnants) that might be PeVatrons:

1. RX J1713.7-3946
2. Cassiopeia A
3. HESS J1731-347
4. HAWC J2227+610

They ran simulations to see if CTAO could finally solve the mystery for each one.

The Investigation: Two Main Tests

1. The "Cut-Off" Test (Finding the Speed Limit)

Every accelerator has a speed limit. If a factory is a true PeVatron, it can push particles up to 1 PeV. If it's a weaker factory, the particles stop accelerating at a lower energy (a "cut-off").

• The Analogy: Imagine trying to guess the top speed of a race car. If you only have a blurry photo, you can't tell if it's going 100 mph or 200 mph. But if you have a high-speed camera (CTAO), you can see exactly where the car starts to slow down.
• The Result: The team found that CTAO needs to stare at a source for at least 100 hours to clearly see this "speed limit."
  • For three of the suspects (RX J1713, Cassiopeia A, and HESS J1731), the simulation showed they stop accelerating way before they reach PeV speeds. They are like fast cars, but not rocket cars. CTAO can confidently say: "Not a PeVatron."
  • For the fourth suspect (HAWC J2227+610), the picture is still a bit fuzzy. It might be a PeVatron, but the data isn't quite enough to be 100% sure yet. It remains a "maybe."

2. The "PeVatron Test" (The Final Verdict)

The team used a special mathematical tool called the PeVatron Test Statistic (PTS). Think of this as a lie detector test for the data.

• They asked the data: "Are you coming from a machine that goes to 1 PeV, or a weaker one?"
• The Verdict: For the first three suspects, the "lie detector" screamed "NO!" with a confidence level so high (5 sigma) that it's practically a scientific fact. They are not PeVatrons.
• For HAWC J2227+610, the lie detector said, "I'm not sure. The data is too quiet to tell."

The Big Takeaway

This paper is essentially a blueprint for the future. It tells us:

1. CTAO is a game-changer: It will be 100 times more precise than our current tools, allowing us to see the "tail" of the energy spectrum where the real secrets hide.
2. Patience is key: To solve these mysteries, we need to stare at the sky for a long time (100+ hours).
3. We are clearing the deck: We can likely rule out several famous candidates as PeVatrons, which helps scientists focus their search on the real ones.
4. The hunt continues: HAWC J2227+610 is still the most promising candidate, but we need more data (and maybe help from other telescopes) to confirm if it's the ultimate cosmic particle accelerator.

In short, this study is like a detective saying, "We have a new, super-sharp magnifying glass. We've checked four suspects, and three are innocent. One is still on the list, but we need to keep watching to catch them in the act."

Technical Summary

Here is a detailed technical summary of the paper "CTAO Simulations for Potential PeVatron Candidates."

1. Problem Statement

Cosmic Rays (CRs) below the "knee" of the energy spectrum (∼PeV) are believed to originate from Galactic sources, specifically Supernova Remnants (SNRs) acting as particle accelerators (PeVatrons). However, directly observing CRs is impossible due to magnetic deflection. Instead, scientists rely on neutral messengers, specifically high-energy $\gamma$-rays produced by hadronic interactions (proton-proton collisions).

Current ground-based Imaging Atmospheric Cherenkov Telescopes (IACTs) face significant limitations:

• Sensitivity: They lack the sensitivity to detect faint fluxes at very high energies (VHE).
• Energy Range: Their energy coverage often cuts off before the proton spectrum reaches the PeV scale, making it difficult to distinguish between leptonic (electron-driven) and hadronic (proton-driven) emission mechanisms.
• Ambiguity: For several candidate SNRs (RX J1713.7-3946, Cassiopeia A, HESS J1731-347, and HAWC J2227+610), existing data is insufficient to confirm if they are true PeVatrons (capable of accelerating protons to $\ge$ 1 PeV).

The paper aims to investigate the capabilities of the Cherenkov Telescope Array Observatory (CTAO) to resolve these ambiguities, specifically by simulating observations to determine if CTAO can distinguish hadronic emission, measure proton cut-off energies, and confirm or exclude PeVatron candidates.

2. Methodology

The study employs a simulation-based approach using the Gammapy Python package and the Naima radiative modeling library.

• Source Selection: Four SNR candidates were selected based on previous multi-wavelength (MWL) analyses showing significant hadronic contributions:
  1. RX J1713.7-3946
  2. Cassiopeia A
  3. HESS J1731-347
  4. HAWC J2227+610 (G106.3+2.7)
• Radiative Modeling:
  • The study assumes a hadronic origin where protons follow an Exponential Cut-off Power-Law (ECPL) distribution: $f_p(E) \propto (E/E_0)^{\alpha_p} \exp(-E/E_{p,cut})^\beta$.
  • Parameters (spectral index, cut-off energy, amplitude, magnetic field, density) were derived from fitting existing MWL data using Naima.
  • Model Selection: Among various radiative models (synchrotron, bremsstrahlung, inverse Compton, pion-decay), the pion-decay model was chosen for CTAO flux simulation as it provided the best fit (lowest log-likelihood) and required fewer parameters.
• Simulation Setup:
  • Instrument: CTAO (North and South sites) using the Prod5 Instrument Response Functions (IRFs).
  • Observation Strategy: Simulations were run for 50, 100, and 200 hours of observation time.
  • Data Generation: Gammapy Maker class generated datasets based on the spectral models. The simulated high-energy CTAO data was combined with real low-energy archival data to create a "Combined Dataset" (CD) for spectral fitting.
• Analysis Metrics:
  • Flux Precision: Comparison of error bars between real data and CTAO simulations.
  • Cut-off Differentiation: Using the Test Statistic ($\Delta$TS) to determine if CTAO can distinguish between different proton cut-off energies (e.g., 300 TeV vs. 600 TeV vs. 1 PeV).
  • PeVatron Test Statistic (PTS): A likelihood ratio test defined to check if the proton cut-off energy ($E_{cut,p}$) is consistent with 1 PeV.
    • $S_{PTS} < -5$: Source is not a PeVatron (5$\sigma$ confidence).
    • $S_{PTS} \ge 5$: Source is a PeVatron.
    • $|S_{PTS}| < 5$: Inconclusive.

3. Key Contributions

• Quantification of CTAO Precision: Demonstrated that CTAO can reduce flux uncertainties by a factor of $\sim$100 compared to current instruments, particularly in the multi-TeV range.
• Cut-off Energy Sensitivity: Established the energy limit at which CTAO can effectively distinguish between different proton cut-off scenarios.
• PeVatron Exclusion/Confirmation: Applied the PTS metric to rigorously test the PeVatron hypothesis for four specific candidates using only simulated CTAO data.
• Observation Time Requirements: Determined the minimum exposure time required to resolve spectral features (cut-offs) effectively.

4. Results

A. Flux Simulation and Precision

• CTAO simulations show significantly smaller error bars than current IACTs.
• Table 4 indicates flux error reductions:
  • RX J1713.7-3946: $3.01 \times 10^{-12} \to 2.76 \times 10^{-13}$erg cm$^{-2}$s$^{-1}$.
  • HAWC J2227+610: $6.31 \times 10^{-13} \to 1.95 \times 10^{-13}$erg cm$^{-2}$s$^{-1}$.
• Combining real low-energy data with simulated CTAO high-energy data improved the determination of proton cut-off energies ($E_{p,cut}$). For HESS J1731-347, the cut-off value shifted by 118.89%, highlighting the necessity of high-energy data for accurate modeling.

B. Cut-off Differentiation Study (HAWC J2227+610)

• Simulations tested proton cut-offs at 300 TeV, 600 TeV, and 1 PeV.
• Observation Time: 50 hours was insufficient to distinguish cut-offs. 100 hours was the minimum required to begin separating fluxes, while 200 hours provided the best statistical significance.
• Detection Limit: CTAO can effectively distinguish fluxes with different cut-offs up to $\sim$600 TeV. Above this energy, flux separation diminishes due to sensitivity limits, making it difficult to distinguish between a 600 TeV cut-off and a 1 PeV cut-off without extremely long observation times.

C. PeVatron Test Statistic (PTS) Results

Using the PTS metric with 50 hours of simulated observation:

• RX J1713.7-3946: $S_{PTS} = -48$ (Excluded as PeVatron).
• Cassiopeia A: $S_{PTS} = -26$ (Excluded as PeVatron).
• HESS J1731-347: $S_{PTS} = -10$ (Excluded as PeVatron).
• HAWC J2227+610: $S_{PTS} = -2$ (Inconclusive).
  • The negative values for the first three sources indicate a rejection of the PeVatron hypothesis with $>5\sigma$ significance.
  • HAWC J2227+610 remains inconclusive because its true cut-off is likely near or above 1 PeV, and CTAO alone (without MWL data) cannot definitively distinguish the tail of the spectrum from a 1 PeV cut-off.

5. Significance and Conclusions

• Confirmation of CTAO's Potential: The study confirms that CTAO will revolutionize the study of Galactic PeVatrons by providing high-precision flux measurements in the multi-TeV to hundreds of TeV range.
• Resolution of Ambiguity: CTAO can definitively exclude several SNRs (Cassiopeia A, RX J1713.7-3946, HESS J1731-347) as PeVatron candidates, resolving long-standing debates about their acceleration capabilities.
• Limitations and Future Work:
  • For the most promising candidate (HAWC J2227+610), CTAO alone may not be sufficient to confirm the PeVatron nature due to source confusion and the need for data beyond 600 TeV.
  • Recommendations: Future studies must combine CTAO data with Multi-Wavelength (MWL) observations and data from other future observatories (e.g., ASTRI Mini-array, SWGO) to disentangle leptonic and hadronic components and extend the energy reach.
• Observational Strategy: A minimum of 100 hours of observation is recommended for candidates to effectively constrain spectral cut-offs.

In summary, this paper provides a robust framework for using CTAO simulations to validate PeVatron candidates, demonstrating that while CTAO will likely exclude current "false positives," it will require complementary data to definitively confirm the first Galactic PeVatron.