Optical Super-orbital Modulation of SMC X-1: Disk Precession and a Revised Pulsar Mass

Technische Samenvatting: Optische Super-orbitale Modulatie van SMC X-1: Schijfprecessie en een Herzien Pulsarmassa

Probleemstelling
De observationele bepaling van de ondergrens voor de massa van neutronensterren is cruciaal voor het beperken van de fysica van kerninstortings-supernova's (CCSNe). De hoog-massieve röntgen-dubbelster SMC X-1 is een cruciaal object in deze context, omdat eerdere massaschattingen de pulsarmassa in de buurt van, of potentieel onder, de theoretische ondergrens van $\approx 1.1 M_\odot$ plaatsten. Deze eerdere schattingen leunden op het bepalen van de massa-verhouding van de dubbelster en de orbitale inclinatie, voornamelijk via spectroscopische radiale-snelheidsmetingen (RV) van de donorster en modellering van de optische ellipsoïdale modulatie. Echter, SMC X-1 vertoont complexe optische super-orbitale modulaties (variaties gesynchroniseerd met de $\sim 40\text{--}60$ dagen durende röntgen-super-orbitale cyclus) die niet volledig werden meegenomen in eerdere statische modellen. Waarnemingen vanaf de grond misten de fotometrische precisie en temporele bemonstering om deze variaties binnen individuele orbitale cycli op te lossen, wat leidde tot potentiële systematische fouten in de afgeleide parameters van de dubbelster.

Methodologie
De auteurs analyseerden optische lichtkrommes van SMC X-1 met hoge cadans, verkregen door de Transiting Exoplanet Survey Satellite (TESS) over zes sectoren, gecombineerd met langetermijn-röntgenmonitoringdata van de Monitor of All-sky X-ray Image (MAXI).

1. Observationele Analyse: De studie identificeerde systematische variaties in de optische orbitale lichtkrommes die gesynchroniseerd waren met de röntgen-super-orbitale modulatie. Specifiek varieerden de diepte van het minimum bij inferieure conjunctie en de asymmetrie van de dubbel-gepiekte lichtkromme afhankelijk van de röntgen-toestand (opkomend, afnemend, aan, of uit).
2. Modelontwikkeling: Een gemodificeerde code voor ellipsoïdale modulatie werd ontwikkeld in Python (met gebruik van NumPy, SciPy, Astropy en CuPy voor parallelisatie). In tegenstelling tot statische modellen, integreert deze code een precesserende accretieschijf die de volgende aspecten verandert:
   • De geometrie van röntgen-irradiatie op het oppervlak van de donorster.
   • De geometrie van optisch/ultraviolette (OUV) reprocessing op het oppervlak van de schijf.
   • De wederzijdse occultatie tussen de donor, de schijf en de pulsar.
     Het model berekent de temperatuur en zwarte-straling van oppervlakte-elementen (SE's) op zowel de donor als de schijf, rekening houdend met zwaartekrachtsverdonkering, röntgenverwarming en OUV-irradiatie.
3. Parameterschatting: Markov Chain Monte Carlo (MCMC) bemonstering werd uitgevoerd om het model aan te passen aan de TESS-optische en MAXI-röntgenlichtkrommes. De analyse testte twee scenario's voor de massa-verhouding ($q = M_{NS}/M_s$): één gebaseerd op waargenomen RV's ($\log q = -1.17$) en een tweede waarbij $q$ met 20% werd verhoogd ($\log q = -1.09$) om de hypothese te testen dat de waargenomen donor-RV onderschat is.

Belangrijkste Resultaten

• Correlaties: De studie bevestigde sterke correlaties tussen de röntgen-super-orbitale flux en optische kenmerken: het minimum bij inferieure conjunctie is ondiep tijdens röntgen-"uit"-toestanden en diep tijdens "aan"-toestanden; de asymmetrie van de dubbele piek correleert met de afgeleide van de röntgenflux (bijvoorbeeld, de linkerpiek is hoger tijdens de opkomende röntgenfase).
• Modelvalidatie: Het gemodificeerde ellipsoïdale model slaagde erin de waargenomen optische en röntgenlichtkrommes te reproduceren, wat observationeel bewijs levert dat de super-orbitale modulatie wordt aangedreven door een precesserende accretieschijf. De analyse suggereert dat de schijfprecessie retrograde is.
• Bias in Radiale Snelheid: De intense röntgen-irradiatie van de donorster verschuift het centrum van optische emissie weg van het zwaartepunt. Het model voorspelt dat dit effect de waargenomen amplitude van de radiale snelheid met ongeveer 20% onderschat (de schijnbare RV-kromme lijkt "afgeknipt" in vergelijking met een echte sinusgolf).
• Herziene Massaschatting:
  • Met de standaard massa-verhouding ($\log q = -1.17$) is de afgeleide pulsarmassa $M_{NS} = 1.11 \pm 0.07 M_\odot$, consistent met eerdere studies.
  • Bij correctie voor de 20% RV-onderschatting (door $q$ te verhogen naar $\log q = -1.09$), is de herziene pulsarmassa $M_{NS} = 1.35 \pm 0.10 M_\odot$.
• Discrepanties: Het model vereiste een röntgenluminositeit ($L_X \approx 4 \times 10^{39}$ erg s$^{-1}$) die aanzienlijk hoger is dan typische waargenomen waarden om de intensiteit van de optische modulatie te reproduceren. De auteurs suggereren dat dit kan wijzen op anisotrope röntgenemissie (bundeling richting de donor) of een misalignement van spin en baan, hoewel het laatste moeilijk te verenigen is met de circulaire baan van het systeem, tenzij massa-accretie tijdens een gemeenschappelijke-hulsfase inefficiënt was.

Betekenis
Het artikel beweert dat de systematische variaties in de optische lichtkromme, aangedreven door de precesserende schijf, eerder niet op een manier werden gemodelleerd die de massabepaling beïnvloedde. Door een model te ontwikkelen dat rekening houdt met de veranderende irradiatie-geometrie, tonen de auteurs aan dat de radiale snelheid van de donor waarschijnlijk onderschat is door door irradiatie veroorzaakte verschuivingen in het emissiecentrum. Correctie voor dit effect levert een herziene pulsarmassa op van $\approx 1.35 M_\odot$. Deze waarde plaatst de SMC X-1-pulsar ruim binnen het typische bereik voor neutronensterren, waardoor de spanning met theoretische ondergrenzen voor kerninstortings-supernova's wordt opgeheven, en suggereert dat eerdere massaschattingen bevooroordeeld waren door niet-gecorrigeerde irradiatie-effecten.