Molecular absorption of Cherenkov light at CTAO

This paper analyzes how molecular absorption and Rayleigh scattering affect Cherenkov light transmission at the CTAO observatory, using simulations to propose an atmospheric monitoring and calibration strategy that mitigates systematic uncertainties for high-energy gamma-ray astronomy.

The Gist

Imagine the Cherenkov Telescope Array (CTAO) as a giant, ultra-sensitive pair of eyes looking up at the sky, trying to catch invisible "ghosts" called gamma rays. These ghosts don't have eyes of their own, so when they crash into Earth's atmosphere, they create a brief, beautiful flash of blue light (Cherenkov light), like a sonic boom made of light. The telescopes catch this flash to figure out where the ghost came from and how powerful it was.

However, the atmosphere isn't just empty space; it's a thick, shifting soup of air, water, and gases. Sometimes, this soup gets a little "dirty" or "thick," making it harder for the blue flash to reach the telescopes.

This paper is essentially a weather report for the telescope's vision, specifically focusing on two invisible ingredients in the air: Ozone and Nitrogen Oxides.

Here is the breakdown of what the scientists found, using some everyday analogies:

1. The Problem: The "Foggy Window" Effect

Think of the atmosphere as a giant window through which the telescopes look.

• Rayleigh Scattering: This is like the window being slightly dusty. It scatters light in all directions. We know this happens all the time, and we can predict it pretty well.
• Molecular Absorption (The Real Issue): This is like someone smearing a dark, sticky substance (like grease or ink) on the window. This substance eats the light before it reaches the camera. The paper focuses on Ozone (a gas that usually protects us from the sun) and Nitrogen Oxides.

2. The Culprits: The "Traveling Gas Clouds"

Ozone and Nitrogen Oxides are weird gases. Unlike oxygen or nitrogen, which are mixed evenly throughout the air (like sugar dissolved perfectly in a cup of tea), these gases are clumpy. They form pockets or clouds that move around.

• The "Stratosphere-to-Troposphere Transport" (STT) Event: Imagine a giant, invisible elevator in the sky. Sometimes, this elevator drops a heavy bucket of concentrated ozone from the upper atmosphere (the stratosphere) down into the lower atmosphere (the troposphere) where the telescopes are.
• The Result: Suddenly, the "window" gets much darker. The telescopes are looking through a thicker layer of ozone than usual.

3. The Investigation: Checking the "Menu"

The scientists looked at data from two telescope locations:

• CTAO-North: Located in the Canary Islands (Spain).
• CTAO-South: Located in the Atacama Desert (Chile).

They used super-computers to simulate what happens when these "ozone buckets" drop. They asked: If the air gets 20% darker with ozone, how much does the telescope's view get messed up?

4. The Findings: It's Not a Big Deal... Mostly

Here is the good news and the "watch out" warning:

• The "Big Picture" is Fine: For most observations, especially high-energy gamma rays, the extra ozone is like looking through a slightly tinted pair of sunglasses. The telescopes can still see clearly enough. The error is small (less than 1-2%).
• The "Low Energy" Warning: This is where it gets tricky. Low-energy gamma rays create very faint, dim flashes of light. If you add a layer of ozone, it's like trying to see a candle through a dirty window. The light gets dimmed significantly.
  • The Impact: During these "ozone drop" events, the telescopes might miss about 3% to 5% of the low-energy events they should have seen. It's like a security camera failing to spot a small, quiet intruder because the lens got smudged.
• Nitrogen Oxides: The scientists checked these too, but they found them to be like a speck of dust on the lens. They barely matter at all.

5. The Solution: A "Smart Cleaning Schedule"

The paper concludes that we don't need to panic, but we do need a plan.

• Current Status: The telescopes are currently calibrated assuming the air is "average."
• The Proposal: Since these ozone clouds happen frequently (especially in certain seasons), the scientists suggest a monitoring system.
  • Imagine a weather app that doesn't just tell you if it's raining, but tells you, "Hey, there's a thick ozone cloud moving over the telescope right now."
  • If the software knows this, it can adjust its calculations. It's like telling the camera, "The lens is dirty right now, so let's boost the brightness settings to compensate."

Summary

Think of the CTAO as a high-end camera trying to take photos of fireworks at night. Usually, the air is clear. But sometimes, a sudden gust of wind blows a cloud of smoke (ozone) right in front of the lens.

• For bright fireworks (high energy): You barely notice the smoke.
• For dim fireworks (low energy): The smoke makes them disappear.

This paper says: "We have a way to track when that smoke cloud is there. If we use that info to adjust our camera settings, we won't miss any of the dim fireworks, ensuring our data is perfect."

This ensures that when scientists study the universe, they aren't accidentally blaming a smudged lens for a missing star.

Technical Summary

1. Problem Statement

The Cherenkov Telescope Array Observatory (CTAO) is designed to detect very-high-energy gamma rays by observing the Cherenkov light produced by extensive air showers (EAS) in the atmosphere. The accuracy of CTAO's scientific performance relies heavily on precise atmospheric modeling. While Rayleigh scattering (by air molecules) and particulate scattering (by aerosols) are well-understood extinction mechanisms, molecular absorption by non-well-mixed gases (specifically ozone and nitrogen oxides) introduces variable systematic uncertainties.

The core problem addressed is the potential impact of Stratosphere-to-Troposphere Transport (STT) events and seasonal variations in ozone ($O_3$) and nitrogen oxides ($NO_x$) on the transmission of Cherenkov light. These events can cause rapid, significant changes in atmospheric optical depth, potentially biasing energy reconstruction and effective area calculations, particularly for low-energy gamma-ray observations ($<150$ GeV). The study aims to quantify these effects to determine if specific calibration strategies are required to meet CTAO's systematic uncertainty budget (energy scale uncertainty $<10\%$).

2. Methodology

The study employs a multi-step approach combining atmospheric data analysis, optical depth modeling, and Monte Carlo simulations:

• Data Sources:

  • Ozone: Retrieved from the ECMWF ERA-5 reanalysis dataset (chosen for its high spatial resolution of $0.25^\circ$ and integration into the CTAO pipeline) and the EAC4 dataset.
  • Nitrogen Oxides ($NO_x$): Retrieved from the EAC4 dataset, which offers more advanced chemical modeling than ERA-5.
  • Validation: Cross-validated with the GDAS analysis dataset (ds083.2).
  • Timeframe: Analysis covers data from 2019–2024 for both CTAO-North (La Palma, Spain) and CTAO-South (Atacama, Chile).

• Profile Generation:

  • Molecular Absorption Profiles (MAPs): Mass mixing ratios were converted to number densities using atmospheric density profiles (from ERA-5). Absorption cross-sections were sourced from the HITRAN database (for $NO_2$) and a combination of HITRAN and Ref. [39] (for $O_3$).
  • Molecular Extinction Profiles (MEPs): MAPs were combined with Rayleigh scattering profiles using the Beer-Lambert law to create total extinction profiles.
  • Scenarios: The study modeled "average" seasonal profiles (Summer/Winter) and specific "extreme" STT events (e.g., a deep STT event at CTAO-South in June 2020 and a shallow, intense event at CTAO-North in February 2021).

• Simulation & Analysis:

  • Fast Simulation: Used the testeff tool (part of sim_telarray) to convolve atmospheric absorption profiles with the Cherenkov emission spectrum and telescope optical throughput. This provided first-order estimates of efficiency drops.
  • Full Simulation: Generated 300 million gamma-ray showers using CORSIKA (PROD6 configuration) at a $20^\circ$ zenith angle. These were processed through sim_telarray using the specific MEPs to analyze image intensities and trigger effective areas for Large-Sized Telescopes (LSTs), Medium-Sized Telescopes (MSTs), and Small-Sized Telescopes (SSTs).

3. Key Contributions

• Quantification of STT Impact: The paper provides the first detailed assessment of how specific STT events (both deep and shallow) alter the optical depth in the CTAO wavelength range (280–750 nm).
• Tool Development: Development of software to generate MAPs and MEPs from public atmospheric data, formatted for compatibility with CTAO simulation tools.
• Differentiation of Gases: A clear distinction in the impact of Ozone vs. Nitrogen Oxides. The study confirms Ozone is the dominant absorber in the relevant UV/visible range, while $NO_x$ effects are negligible.
• Calibration Strategy Proposal: A data-driven recommendation for atmospheric monitoring and potential calibration procedures to mitigate systematic biases.

4. Key Results

Atmospheric Variability

• Ozone: Exhibits strong seasonal cycles and dynamic transport events.
  • CTAO-South: STT events can double the ozone mixing ratio, persisting for days.
  • CTAO-North: Shallow STT events can increase ozone mixing ratios by a factor of six, though they often remain confined to lower altitudes (below 500 hPa).
• Nitrogen Oxides: $NO_2$ contributes to extinction but to a much lesser degree than ozone. $NO$ is negligible at night due to oxidation to $NO_2$.

Impact on Light Transmission (Fast Simulations)

• Ozone Extinction: During STT events, ozone-related extinction can exceed Rayleigh scattering in the Hartley band (200–310 nm).
• Efficiency Drop:
  • CTAO-South (STT event): LST detection efficiency dropped by ~0.7%.
  • CTAO-North (STT event): LST detection efficiency dropped by ~2.1%.
  • $NO_2$ Impact: Negligible ($\le 0.02\%$).
• Temperature Dependence: While absorption cross-sections vary with temperature (up to 8% difference in pure ozone extinction between 193 K and 293 K), the effect on total telescope efficiency (including Rayleigh scattering and optics) is minimal ($\le 0.3\%$).

Impact on Performance (Full Simulations)

• Image Intensity:
  • Images recorded during STT events are dimmer.
  • For low-energy showers (fewer than 1,000 photoelectrons), intensity reductions reach 2–3% at CTAO-North and ~1% at CTAO-South.
• Trigger Effective Area:
  • CTAO-South: Reduction is significant only at low energies ($\le 80$ GeV), affecting primarily LSTs.
  • CTAO-North: Reduction is more pronounced, reaching up to 5% for energies below 40 GeV due to the shallow nature of the event affecting the lower atmosphere where low-energy showers develop.
  • At energies $>1$ TeV, the impact is within statistical uncertainty limits.

5. Significance and Recommendations

• Systematic Uncertainty Budget: The CTAO requires the atmospheric contribution to the energy scale uncertainty to be $\le 8\%$. Preliminary estimates suggest controlling molecular absorption effects to the ~1% level is sufficient.
• Current Status: No formal calibration for ozone is currently mandated. However, the study demonstrates that uncorrected STT events can introduce biases (2–5% in effective area/intensity) that approach or exceed the tolerance for low-energy science goals.
• Recommendations:
  1. Monitoring: Implement regular monitoring of ozone mixing ratio profiles, especially during seasons prone to STT events (Dec–May for North; Jun–Jan for South).
  2. Calibration: Develop a strategy to apply ozone-specific corrections to Instrument Response Functions (IRFs) when STT events are detected, particularly for low-energy observations.
  3. Integration: The tools developed in this study can be integrated into the CTAO Data Processing and Preservation System (DPPS) to automate the generation of tailored atmospheric profiles.

Conclusion: While molecular absorption by ozone is generally a minor factor compared to Rayleigh scattering, dynamic events like STT can cause non-negligible, energy-dependent biases. Proactive monitoring and potential calibration are recommended to ensure CTAO meets its high-precision scientific requirements, especially for low-energy gamma-ray astronomy. Nitrogen oxide variability is deemed negligible for CTAO operations.