Extending TESS flare frequency distributions with CHEOPS: Power-law versus lognormal

Probleemstelling

Stervlekkens (stellar flares) zijn intense uitbarstingen van straling veroorzaakt door magnetische reconnectie, vooral frequent op M-dwergsterren. Het begrijpen van de frequentie van deze vlekken is cruciaal voor het beoordelen van de bewoonbaarheid van exoplaneten, aangezien energieke vlekken atmosferen kunnen eroderen of juist prebiotische chemie kunnen stimuleren.

Traditioneel worden Flare Frequency Distributions (FFD's) gemodelleerd als een machtswet ($dN/dE \propto E^{-\alpha}$). Echter, recente studies suggereren dat deze aanname onjuist kan zijn en dat een lognormale verdeling een betere beschrijving biedt, wat impliceert dat vlekken een ander generatiemechanisme hebben (bijv. multiplicatieve superpositie van kleinere gebeurtenissen in plaats van zelf-georganiseerde criticaliteit).

De kernproblemen die deze studie aanpakt zijn:

1. Observatiebias: Lage-energie vlekken worden vaak gemist door instrumentele beperkingen, wat de FFD kan vervormen (schijnbare afvlakking bij lage energieën).
2. Energie-metriek: Er is een verschil tussen het gebruik van Equivalent Duration (ED, direct gemeten uit fotometrie) en bolometrische energie ($E_{bol}$, fysiek relevant voor habitabiliteit). De conversie van ED naar $E_{bol}$ hangt af van de sterrenluminositeit, die zelf lognormaal verdeeld is. Dit kan een oorspronkelijke machtswet in ED verstoren tot een schijnbare lognormale verdeling in $E_{bol}$.
3. Complexiteit: Veel vlekken hebben complexe structuren (meerdere pieken). Het behandelen hiervan als één enkele gebeurtenis versus het ontleden in componenten beïnvloedt de afgeleide FFD.

Methodologie

De auteurs hebben een uitgebreide analyse uitgevoerd op een steekproef van 110 M-dwergsterren (types M0 tot M5) met data van twee complementaire ruimtemissies: TESS (Transiting Exoplanet Survey Satellite) en CHEOPS (CHaracterising ExOPlanet Satellite).

Kernstappen in de analyse:

1. Data-Verwerking:
   • TESS: Lange tijdbaselines (tot 93 sectoren), lage cadentie (20s of 120s).
   • CHEOPS: Korte tijdbaselines (~52 dagen totaal), ultra-hoge cadentie (tot 3s) en hoge fotometrische precisie.
   • Detrending met Savitzky-Golay filters om rotatiemodulatie te verwijderen.
2. Vlek Detectie en Decompositie:
   • Gebruik van het pakket AltaiPony voor detectie.
   • Decompositie: Complexe vlekken met meerdere pieken werden ontbonden in individuele componenten met behulp van het Llamaradas Estelares model en de Akaike Information Criterion (AIC) om het beste aantal componenten te bepalen.
3. Correctie voor Bias (Injection-Recovery):
   • Kunstmatige vlekken werden geïnjecteerd in de lichtkrommen om de detectiekans en energieherstelratio te kwantificeren.
   • Dit corrigeerde voor de onderdetectie van zwakke vlekken en de systematische onderschatting van energie door detrending.
4. Energie Calibratie:
   • Conversie van ED naar bolometrische energie ($E_{bol}$) met behulp van een zwartlichaamspectrum van 10.000 K en instrument-specifieke correctiefactoren ($\epsilon$).
5. Statistische Modellering:
   • Constructie van FFD's voor zowel ED als $E_{bol}$.
   • Vergelijking van drie modellen: Machtswet, Lognormaal, en Afgeknotte Machtswet (truncated power law).
   • Gebruik van Vuong's likelihood ratio test voor niet-geneste modellen, aangevuld met stabilised Kolmogorov-Smirnov (sKS) tests en exceedance tests om de staart van de verdeling te analyseren.

Belangrijkste Bijdragen

• Combinatie van Missies: Voor het eerst worden TESS en CHEOPS data gecombineerd om een FFD te construeren die zes ordes van grootte in bolometrische energie bestrijkt ($10^{30}$ tot $>10^{35}$ erg).
• Decompositie van Complexe Vlekken: Systematische ontleding van complexe vlekken, waarbij CHEOPS door zijn hogere cadentie bijna twee keer zo veel complexe vlekken detecteert als TESS.
• Bias-Correctie: Toepassing van een robuuste injection-recovery methode om de waargenomen FFD te corrigeren voor instrumentele onvolledigheid, wat essentieel is voor het onderscheiden van fysieke effecten van observatie-artefacten.

Resultaten

1. ED-gebaseerde FFD's:
   • FFD's gebaseerd op Equivalent Duration volgen nauwkeurig een machtswet.
   • Na correctie voor detectiebias wordt de helling ($\alpha$) steiler (bijv. van 1.67 naar 1.95 voor CHEOPS), wat aangeeft dat er meer zwakke vlekken zijn dan oorspronkelijk waargenomen.
2. Bolometrische Energie ($E_{bol}$) FFD's:
   • FFD's gebaseerd op $E_{bol}$ wijken significant af van een pure machtswet.
   • Ze worden beter beschreven door een lognormale verdeling of, nog beter, een afgeknotte machtswet (truncated power law).
   • De afwijking wordt toegeschreven aan de convolutie van de machtswet in ED met de lognormale verdeling van de sterrenluminositeiten in de steekproef.
3. De "Break" bij Hoge Energie:
   • Er is een duidelijke afwijking van de machtswet bij ongeveer $1.8 \times 10^{35}$ erg.
   • De exceedance test en sKS-test tonen aan dat deze afwijking waarschijnlijk het gevolg is van beperkte bemonstering (we hebben niet lang genoeg geobserveerd om de zeldzaamste, meest energieke superflares te zien) en niet noodzakelijk een intrinsieke fysieke truncatie.
4. Invloed op Habitabiliteit:
   • De schijnbare afvlakking bij lage energieën (die eerder als bewijs voor lognormaal gedrag werd gezien) is een artefact van observatiebias. Na correctie volgt de verdeling een machtswet tot de truncatiegrens.
   • De frequentie van superflares ($>10^{34}$ erg) is lager dan wat een pure machtswet-extrapolatie zou voorspellen, wat de zorgen over atmosferische erosie op exoplaneten iets kan temperen.

Significantie en Conclusie

De studie concludeert dat de keuze van de energie-metriek en de correctie voor observatiebias cruciaal zijn voor het interpreteren van vlekstatistieken.

• Fysica: De onderliggende vlekgeneratie lijkt een machtswet te volgen in termen van ED, maar de conversie naar bolometrische energie en de beperkte bemonstering van extreme gebeurtenissen leiden tot schijnbare afwijkingen.
• Toekomstige Missies: Huidige instrumenten (TESS/CHEOPS) kunnen niet betrouwbaar vlekken boven $10^{35}$ erg bemonsteren. De toekomstige PLATO-missie, met zijn langere baselines en hoge precisie, zal nodig zijn om beide regimes (lage en hoge energie) volledig te karakteriseren en te bepalen of de truncatie fysiek of statistisch van aard is.
• Toepassing: De gecorrigeerde FFD's bieden realistischere input voor simulaties van exoplanetatmosferen en voor het ontwerpen van toekomstige missies om vlekken te detecteren.

Kortom, de schijnbare noodzaak voor lognormale modellen in FFD's is grotendeels te wijten aan observatiebias en de convolutie met sterrenluminositeit; na correctie blijft een machtswet (met een mogelijke truncatie bij zeer hoge energieën door bemonsteringslimieten) de beste beschrijving.