Šum jsou data: KIC 8462852

Problém: co je šum?

Standardní astronomický pipeline pro zpracování Kepler fotometrie prochází kroky: kalibrace instrumentálních efektů, odstranění trendu, identifikace outlierů, vyhlazení, modelování dipů. Výsledkem je PDCSAP (Pre-search Data Conditioning SAP) — „vyčištěná“ světelná křivka. Vše, co pipeline odstraní, je implicitně klasifikováno jako šum.

Tento přístup je optimalizovaný pro hledání tranzitů — periodických, tvarově definovaných signálů. Pro KIC 8462852, kde jsou dipy aperiodické a morfologicky variabilní, je situace jiná. Vzniká otázka: co když pipeline odstraňuje nejen instrumentální artefakty, ale i nízkoamplitudové astrofyzikální signály?

Otázku nelze zodpovědět spekulací. Lze ji zodpovědět šumovým budgetem: systematicky identifikovat každý známý zdroj šumu, kvantifikovat jeho příspěvek k celkové varianci, odstranit ho, a analyzovat, co zbývá.

Dataset a metodologie

Data: Kepler SAP (Simple Aperture Photometry) pro KIC 8462852, 18 kvartálů (Q0–Q17), 65 916 platných bodů, kadence 29,4 min. Používáme výhradně SAP, nikoli PDCSAP, aby každý krok čištění byl transparentní a kontrolovatelný.

Pracujeme také s diagnostickými kanály obsaženými v FITS souborech: SAP_FLUX_ERR (fotonový šum), SAP_BKG (pozadí), POS_CORR1/2 (korekce pozice — pointing jitter), MOM_CENTR1/2 (centroid pozice).

Pro každý kvartál nezávisle provádíme tři kroky. Prvním krokem je polynomiální detrend (stupeň 2) — odstraňuje pomalý instrumentální drift způsobený změnami teploty, stárnutím detektoru a změnami orientace. Druhým krokem je regrese proti pointing jitteru — rezidua po detrendu korelujeme s POS_CORR1, POS_CORR2 a jejich kvadratickými a křížovými členy (5 regresorů). Třetím krokem je regrese proti background variaci — zbylá rezidua korelujeme s normalizovaným SAP_BKG a jeho kvadratickým členem.

Fotonový šum (SAP_FLUX_ERR) má medián 62 ppm a tvoří méně než 0,1 % celkové variance. Je zanedbatelný a neprovádíme proti němu žádnou korekci — je to ireducibilní fyzikální limit měření.

Šumový budget

Celková variance SAP reziduí po polynomiálním detrendu: 5 281 ppm RMS.

Pointing jitter (POS_CORR korelace): 1 810 ppm, 11,7 % variance. Toto je dominantní instrumentální zdroj — mechanické chvění teleskopu na oběžné dráze.

Background variace (SAP_BKG korelace): 869 ppm, 2,7 % variance. Rozptýlené světlo a plošné pozadí.

Fotonový šum: 62 ppm, pod 0,1 % variance. Zanedbatelný.

Celkem instrumentální: 14,4 % variance.

Čistý reziduál po odstranění všech identifikovaných instrumentálních zdrojů: 4 885 ppm, 85,5 % variance. Tato zbytková variance nelze vysvětlit regresí proti dostupným instrumentálním diagnostikám. To neznamená, že je nutně astrofyzikální — mohou existovat instrumentální mechanismy, které dostupné regresory nezachycují (common-mode artefakty, teplotní cykly, focus drift). Rozlišení vyžaduje kontrolní vzorek hvězd ze stejného CCD modulu.

85,5 % celkové variance SAP reziduí nelze vysvětlit regresí proti identifikovaným instrumentálním diagnostikám.

Analýza čistého reziduálu

Definujeme šest oken pro srovnávací analýzu. Každé okno má délku 90–190 dní. Pre-dipová okna jsou definována s bezpečnostním bufferem — končí 20 nebo 50 dní před minimem dipu, aby nedocházelo k leakage z dipu samotného do analyzovaného okna.

Kontrolní okna: klidné období BJD 400–600 (referenční, 9 340 bodů), klidné období 2 BJD 900–1 090 (sekundární kontrola, 8 811 bodů).

Pre-dipová okna: 90 dní končící 20 dní před D800 (3 395 bodů), 100 dní končících 50 dní před D800 (3 976 bodů), 90 dní končících 20 dní před D1520 (3 441 bodů), 100 dní končících 50 dní před D1520 (4 376 bodů).

Výsledky, normalizované na klidné období (400–600):

Kontrola (900–1090): 1,51× — přirozená baseline variabilita hvězdy.
Pre-D800, buffer 50 dní: 3,55×.
Pre-D800, buffer 20 dní: 7,91×.
Pre-D1520, buffer 50 dní: 7,50×.
Pre-D1520, buffer 20 dní: 22,83×.

Gradient je konzistentní: čím blíže k dipu, tím vyšší variance. S bufferem 50 dní (nejkonzervativnější odhad, žádný možný dip leakage) je pre-dipová variance stále 3,6–7,5× nad referencí, zatímco kontrolní okno nepřesahuje 1,5×.

Pre-dipová variance v čistém reziduálu je 3,6–7,5× nad referencí i s konzervativním 50denním bufferem, zatímco kontrolní okno nepřesahuje 1,5×.

Autokorelace

Autokorelační intenzita (Σ ACF²) čistého reziduálu ukazuje odlišný vzorec než variance. Pre-D800 s bufferem 20 dní: 0,91× (srovnatelné s klidným obdobím). Pre-D1520 s bufferem 20 dní: 1,19×.

To znamená: zvýšená variance, ale ne zvýšená korelační struktura. Signál má vyšší amplitudu, ale nikoliv jinou časovou organizaci. Toto je konzistentní s interpretací, že jde o nízkoamplitudové dimming eventy (mikro-dipy, prekurzory), nikoli o koherentní změnu režimu.

Zvýšená variance bez zvýšené autokorelační struktury — konzistentní s mikro-dipy/prekurzory, nikoli s koherentní změnou režimu.

Interpretace a instrumentální alternativa

Toto je ústřední bod celé analýzy.

Zvýšená variance v před-dipových oknech přežívá odstranění tří identifikovaných instrumentálních zdrojů (pointing jitter, background variace, per-quarter polynomiální drift). To je nutná, ale nikoliv dostatečná podmínka pro astrofyzikální interpretaci.

Instrumentální původ nelze vyloučit bez kontrolního vzorku. Nicméně zvýšená variance přežívá odstranění identifikovaných regresorů a je reprodukovatelná pro oba dipy, což klade na instrumentální alternativu dvě podmínky. Za prvé, identifikované instrumentální regresory vysvětlují 14,4 % celkové variance — i pokud existují další neidentifikované instrumentální zdroje, musely by specificky korelovat s před-dipovými okny a nikoliv s kontrolními okny, aby vysvětlily pozorovaný kontrast 3,6–7,5× vs 1,5×. Za druhé, efekt je reprodukovatelný pro oba nezávislé dipy (D800 i D1520) a přežívá konzervativní 50denní buffer.

Zbývají však alternativy, které nelze vyloučit bez dalších dat. Kepler má dlouhodobé driftové struktury, quarter-dependent systematiky a common-mode artefakty (teplotní cykly, focus drift, rolling band noise), které nejsou plně zachyceny regresí proti POS_CORR a SAP_BKG. Tyto systematiky nemají kauzální vazbu na budoucí chování hvězdy, ale mají vlastní časovou strukturu, která se může statisticky potkat s před-dipovým obdobím.

Definitivní rozlišení vyžaduje kontrolní vzorek: 30–50 klidných hvězd ze stejného CCD modulu v každém kvartálu. Pokud tyto hvězdy vykazují srovnatelnou before-event variabilitu v náhodných časových oknech, efekt je instrumentální. Pokud ne, instrumentální původ je vyloučen. Tento test dosud nebyl proveden v kontextu pre-dipové reziduální analýzy.

Co to znamená

Pozitivní nález (podmíněný): V čistém reziduálu — po odstranění pointing jitteru, background variace a instrumentálního trendu — je variance systematicky zvýšená v období desítek dní před oběma hlavními dipy. Efekt je reprodukovatelný (D800 i D1520), přežívá konzervativní buffery (50 dní) a je výrazně nad kontrolní variabilitou (3,6–7,5× vs 1,5×). Interpretujeme to jako kandidátní signál nízkoamplitudové fotometrické aktivity předcházející hlubokým dipům. Definitivní potvrzení vyžaduje kontrolní vzorek ze stejného CCD modulu.

Limity: Netvrdíme, že jde o „fázový přechod“ nebo „strukturální reorganizaci.“ Autokorelační analýza neukazuje změnu korelační struktury, pouze zvýšenou amplitudu. To je konzistentní s jednoduššší interpretací: hvězda začíná vykazovat drobné dimming eventy dříve, než přijde hlavní dip. Netvrdíme, že můžeme predikovat dipy. Netvrdíme, že PDCSAP pipeline je špatný — je optimalizovaný pro jiný účel.

Co jsme na cestě vyloučili: Počáteční analýza na PDCSAP datech ukazovala variance ratio 87–91×. Robustnostní testy odhalily, že většina tohoto signálu pocházela z PDC pipeline artefaktů a dip leakage. Po systematickém čištění zůstává 3,6–7,5× — řádově menší, ale stále statisticky významný a kauzálně interpretovatelný. Transparentnost tohoto procesu je záměrná. Ukazuje, jak snadno může nedostatečná kontrola systematik vést k přehnané interpretaci — a proč je šumový budget nutný krok, nikoli volitelný doplněk.

Metodologický rozdíl

Standardní pipeline pro KIC 8462852: kalibrace → detrending → outlier removal → vyhlazení → modelování dipů → vše ostatní = noise floor.

Náš pipeline: kalibrace → identifikace a kvantifikace každého instrumentálního zdroje → transparentní odstranění → analýza toho, co zbývá → šumový budget jako součást výsledku.

Standardní přístup je optimalizovaný pro detekci tranzitů — periodických, tvarově definovaných signálů. Pro aperiodické jevy typu KIC 8462852 může systematicky odstraňovat nízkoamplitudovou variabilitu, která nese informaci o stavu systému mezi hlavními událostmi.

Výsledek noise-preserving přístupu: tam, kde standardní pipeline vidí „čistou křivku s dipy a šumem,“ náš přístup identifikuje instrumentální příspěvky (14,4 %), nevysvětlený reziduál (85,5 %) a kandidátní pre-dipovou aktivitu, jejíž potvrzení vyžaduje kontrolní vzorek.

Testovatelné predikce

Pokud je pre-dipová zvýšená variabilita reálná, pak by měla být pozorovatelná i v jiných fotometrických kampaních KIC 8462852: LCO pozemní data (2017–), TESS data (pokud jsou dostupná). Měla by být achromatická nebo chromaticky konzistentní s prachovým modelem (Boyajian et al. 2018 ukazují diferenciální reddening). Měla by být nepřítomná u kontrolních hvězd podobného spektrálního typu ve stejném Kepler poli.

Šumový budget by měl být replikován nezávislou skupinou s odlišným pipeline — pokud 85,5 % nevysvětlené variance je robustní výsledek, měl by být reprodukovatelný.

Co tvrdíme a co ne

Tvrdíme: Kompletní šumový budget SAP fotometrie KIC 8462852 ukazuje, že 85,5 % variance reziduí nelze vysvětlit regresí proti identifikovaným instrumentálním diagnostikám. V čistém reziduálu po odstranění těchto zdrojů je variance v před-dipových oknech systematicky zvýšená oproti kontrolním klidným obdobím, a to i s konzervativním 50denním bufferem. Tento efekt je kandidátním signálem astrofyzikálního původu. Definitivní rozlišení mezi astrofyzikálním a dosud neidentifikovaným instrumentálním původem vyžaduje kontrolní vzorek hvězd ze stejného CCD modulu.

Netvrdíme: Nemáme prediktivní model. Netvrdíme fázový přechod ani strukturální reorganizaci. Netvrdíme, že naše čištění je kompletní — existují instrumentální zdroje (common-mode artefakty, focus drift, rolling band noise), které regrese proti POS_CORR a SAP_BKG nezachytí. Netvrdíme, že efekt je unikátní pro KIC 8462852 — k tomu by bylo třeba systematické analýzy většího vzorku.

Kde se potkáváme s konvenční fyzikou

KIC 8462852 zůstává nevysvětlená. Prachové oblaky, komety, cirkumstelární materiál -- žádné vysvětlení nepřežilo všechna data (Boyajian et al. 2016). Konvenční přístup: analyzuj dipy, zbytek je šum. Co když šum obsahuje informaci, kterou dipy nemají?