Inclination Bias in Techniques Used to Identify Be Star Candidates

Technische Samenvatting: Inclinatiebias in technieken gebruikt voor de identificatie van Be-ster kandidaat-objecten

Probleemstelling
Klassieke Be-sterren zijn snel roterende hoofdreekssterren die gekenmerkt worden door waterstof-Balmer-emissie (met name H$\alpha$) als gevolg van een circumstellaire decretieschijf. Een hardnekkig observationeel raadsel is het schijnbare gebrek aan hoog-inclinatie ($i > 80^\circ$) Be-sterren in monsters afgeleid van diverse surveys. Ho evenredig kan dit duiden op een fysieke afwezigheid van dergelijke sterren of gecorreleerde rotatieassen in clusters, maar het kan ook voortkomen uit selectiebiases die inherent zijn aan de methoden die worden gebruikt om Be-ster kandidaat-objecten te identificeren. Dit artikel onderzoekt of spectroscopische en fotometrische technieken voor het detecteren van Be-ster kandidaat-objecten systematische biases introduceren met betrekking tot de inclinatiehoek ($i$) van het systeem, waarbij $i$ wordt gedefinieerd als de hoek tussen de stellaire rotatieas en de gezichtslijn van de waarnemer.

Methodologie
Om deze biases te kwantificeren, genereerden de auteurs een synthetische steekproef van 20.000 Be-sterren met behulp van radiatieve transfercodes (Bedisk en Beray). De steekproef werd als volgt geconstrueerd:

• Stellaire Eigenschappen: Centrale B-type sterren werden getrokken uit een Salpeter initiële massafunctie, met massa's variërend van 3 tot 14 $M_\odot$. Stellaire parameters (straal, luminositeit, $T_{\text{eff}}$, $\log g$) werden overgenomen van Ekström et al. (2012).
• Schijfmodellen: Aan elke ster werd een circumstellaire schijf toegewezen met parameters getrokken uit specifieke distributies: equatoriale dichtheid ($\rho_0$) uit een log-normale distributie, machtswet-index ($n$) uit een normale distributie, en schijfradius ($r_d$) en inclinatie ($i$) getrokken uit een uniforme distributie. De schijfdichtheid volgt een geparametriseerd model dat afhankelijk is van de schaalhoogte en radiale afname.
• Synthetische Observabelen: Spectrale energieverdelingen (SED's) en H$\alpha$-profielen werden berekend voor kijk-inclinaties variërend van $0^\circ$ tot $90^\circ$ (inclusclusief specifieke hoog-inclinatie stappen van $83^\circ$, $85^\circ$ en $87^\circ$).
• Steekproeffiltering: Synthetische sterren met H$\alpha$-profielen die ononderscheidbaar zijn van zuiver fotosferische sterren (een verschil van $<2\%$) werden verworpen om ervoor te zorgen dat de steekproef detecteerbare Be-sterren vertegenwoordigde.
• Bias-testen: De synthetische steekproef werd onderworpen aan drie specifieke identificatiemethoden:
  1. Spectroscopisch: Het piekzoekalgoritme van Hou et al. (2016) toegepast op LAMOST-resolutie spectra ($R \approx 1.800$).
  2. Fotometrisch (Methode A): De kleur-kleur drempelmethode van Iqbal & Keller (2013) gebruikmakend van $R - \text{H}\alpha$ versus $V - I$.
  3. Fotometrisch (Methode B): De kleur-magnitude drempelmethode van Milone et al. (2018) gebruikmakend van $F814W$ versus $F656N - F814W$.

Belangrijkste Resultaten
De studie onthult duidelijke inclinatiebiases voor elke methode:

• Spectroscopische Methode (Hou et al. 2016): Deze methode vertoont een substantiële bias tegen hoge inclinaties ($i > 80^\circ$). Ongeveer 24% van de synthetische Be-sterren werd verworpen als kandidaat. Het afwijzingspercentage stijgt scherp voor $i \ge 80^\circ$, waarbij Be-sterren meer dan twee keer zo vaak worden verworpen vergeleken met sterren bij lagere inclinaties. Dit wordt toegeschreven aan schilabsorptie (shell absorption) in edge-on systemen, wat de emissiepieken die nodig zijn voor het piekzoekalgoritme (condities C1, C2 en C3) onderdrukt. Omgekeerd vertoont de methode een bias ten gunste van lage inclinaties ($i \lesssim 50^\circ$).
• Fotometrische Methode A (Iqbal & Keller 2013): Deze methode is bevooroordeeld tegen zeer hoge inclinaties ($i \approx 90^\circ$), met een afwijzingspercentage van $\approx 52\%$ bij $90^\circ$. Het is echter minder bevooroordeeld tegen lage inclinaties vergeleken met de spectroscopische methode. De bias tegen hoge inclinaties wordt gedreven door diepe absorptiedalen in de H$\alpha$-lijn bij edge-on weergaven, wat de $R - \text{H}\alpha$ magnitude vergroot en sterren onder de detectiedrempel duwt.
• Fotometrische Methode B (Milone et al. 2018): Deze methode vertoont een sterke bias tegen lage inclinaties ($i < 50^\circ$), met een afwijzingspercentage van $\approx 58\%$ in totaal. In tegenstelling tot de andere methoden is deze relatief onbevooroordeeld bij zeer hoge inclinaties ($i \approx 90^\circ$). De bias tegen lage inclinaties ontstaat omdat systemen met een lage inclinatie vaak niet de sterke H$\alpha$-emissie hebben die nodig is om hen significant naar links van de drempellijn in het kleur-magnitude diagram te verschuiven.
• Gecombineerd Effect: Beide fotometrische methoden resulteren in een overschot aan gedetecteerde kandidaten bij matige inclinaties ($50^\circ < i < 80^\circ$). De waarschijnlijkheidsdichtheidsfuncties afgeleid van de geobserveerde steekproeven verschillen significant van de verwachte willekeurige $\sin i$ distributie.

Betekenis en Conclusies
Het artikel concludeert dat de schijnbare afwezigheid van hoog-inclinatie Be-sterren in observationele steekproeven niet als een fysieke realiteit kan worden aangenomen zonder correctie voor methodologische biases.

• De spectroscopische methode van Hou et al. (2016) is zwaar bevooroordeeld tegen hoge inclinaties, wat potentieel het tekort aan edge-on Be-sterren in LAMOST-gebaseerde catalogi verklaart.
• Fotometrische methoden introduceren verschillende biases; Iqbal & Keller (2013) is bevooroordeeld tegen hoge inclinaties, terwijl Milone et al. (2018) bevooroordeeld is tegen lage inclinaties.
• Deze biases moeten worden meegewogen wanneer Be-sterren worden gebruikt als proeven voor gecorreleerde rotatieassen in open clusters of bij het interpreteren van de statistische distributie van Be-ster fracties.
• De auteurs suggereren dat het combineren van complementaire methoden (bijv. spectroscopisch en fotometrisch) de totale inclinatiebias potentieel kan verminderen.
• Tot slot merkt het artikel op dat een aanzienlijk deel van de model Be-sterren (24% tot 58%) wordt verworpen door deze standaard selectiecriteria, wat impliceert dat observationele bepalingen van Be-ster fracties de werkelijke populatie aanzienlijk kunnen onderschatten als deze biases niet worden gecorrigeerd.