Mass distribution of neutron stars in binary systems

Technische Zusammenfassung: Massenverteilung von Neutronensternen in Binärsystemen

Problemstellung
Beobachtungen deuten darauf hin, dass die Massenverteilung bekannter Neutronensterne (NS) in Binärsystemen ein bimodales Muster aufweist, mit Peaks bei etwa 1,4 M$_\odot$ und $\sim$1,8 M$_\odot$. Während der Peak bei 1,4 M$_\odot$ häufig auf Standard-Kernkollaps-Supernovae oder die Entwicklung durch gemeinsame Hüllen zurückgeführt wird, bleibt der Ursprung des Peaks bei $\sim$1,8 M$_\odot$ umstritten. Bisherige theoretische Modelle, die das NS-Massenwachstum durch Akkretion zu erklären versuchten, waren oft durch die Eddington-Grenze begrenzt oder berücksichtigten Super-Eddington-Akkretionsregime nicht vollständig. Darüber hinaus haben Unsicherheiten bezüglich der Effizienz des Massentransfers, des Propeller-Effekts und der Entwicklung des Magnetfelds genaue Vorhersagen der finalen Massenverteilung akkretierender NS behindert. Diese Studie untersucht, ob ein Modell der Super-Eddington-Akkretion die beobachtete bimodale Massenverteilung von NS in Binärsystemen erfolgreich reproduzieren kann.

Methodik
Die Autoren wählten einen hybriden Ansatz der Populationszusammenstellung, der zwei verschiedene Codes kombiniert:

1. Binary Star Evolution (BSE): Diente zur Simulation der Entstehung und initialen Evolution von $10^7$ primordialen Binärsystemen. Dieser Code generierte eine Population von NS-Binärsystemen mit Hauptreihen-Kompanen unter Einbeziehung von Anfangsmassenfunktionen, Massenverhältnissen, Exzentrizitäten und natalen Kick-Geschwindigkeiten. Systeme wurden gefiltert, um diejenigen zu identifizieren, die einen stabilen Massentransfer ermöglichen, wobei solche ausgeschlossen wurden, die zu instabiler gemeinsamer Hüllenentwicklung oder extremen Massentransferraten ($>10^4 \dot{M}_{Edd}$) neigen.
2. Modules for Experiments in Stellar Astrophysics (MESA): Ausgewählte NS-Binärsysteme aus dem BSE-Ausgang wurden detaillierten Evolutionsberechnungen mit MESA (Version r24.08.1) unterzogen. Diese Phase konzentrierte sich auf den Akkretionsprozess unter stabilen Bedingungen der Roche-Lappen-Überfüllung.

Wichtige physikalische Modelle und Annahmen:

• Super-Eddington-Akkretion: Im Gegensatz zu früheren Modellen, die durch die Eddington-Grenze begrenzt waren, integriert diese Studie die Physik der Super-Eddington-Akkretion, wie sie von Lipunov (1987) beschrieben wird. Das Modell unterscheidet zwischen sub-kritischen und super-kritischen Regimen basierend auf der Beziehung zwischen dem Alfvén-Radius ($R_a$) und dem Sphärisierungsradius ($R_s$).
• Akkretionseffizienz: Der Anteil des transferierten Masses, der vom NS akkretiert wird ($f_{mt,a}$), wird durch den Akkretionszustand bestimmt. Im super-kritischen Regime wird überschüssiges Material durch Strahlungsdruck in der Nähe des NS weggeblasen und trägt dabei spezifischen Bahnimpuls ab. Für die zurückbehaltene Masse wird eine Akkretionseffizienz von 0,5 angenommen.
• Entwicklung des Magnetfelds: Das NS-Magnetfeld nimmt als Funktion der akkretierten Masse ($\Delta M$) gemäß einem exponentiellen Gesetz ab und reduziert sich von einem Anfangswert von $5 \times 10^{12}$ G auf ein Minimum von $5 \times 10^8$ G.
• Spin-Entwicklung: Das Modell verfolgt die NS-Spinperiode durch drei Phasen: Radio (keine Akkretion), Propeller (Auswurf von Materie) und Akkretion (Spin-up). Der Übergang zwischen diesen Phasen hängt von den relativen Positionen des Lichtzylinders, des Korotationsradius und des Magnetosphärenradius ab.
• Verlust von Drehimpuls: Die Simulation berücksichtigt Gravitationswellenstrahlung, magnetische Bremsung und den Verlust von Drehimpuls durch ausgeworfene Masse und die Bildung einer zirkumbinären Scheibe.

Wichtige Ergebnisse

1. Mechanismus des Massenzuwachses: Die Studie zeigt, dass das NS-Massenwachstum stark von der orbitalen Periode des Binärsystems und der Masse des Gebersterns abhängt. Im Super-Eddington-Modell können NS über längere Zeiträume in einem Akkretionsgleichgewichtszustand verbleiben, was eine signifikante Massenanreicherung ermöglicht (bis zu $\sim$1,0 M$_\odot$ Wachstum von einer Geburtsmasse von 1,4 M$_\odot$). Dies steht im Gegensatz zum Modell von Li et al. (2021) (Eddington-begrenzt), bei dem NS häufig in eine verlängerte Propeller-Phase eintreten, die das Massewachstum begrenzt.
2. Reproduktion der Bimodalität: Die Ergebnisse der Populationszusammenstellung reproduzieren erfolgreich die beobachtete bimodale Verteilung:
   • Der Peak bei $\sim$1,8 M$_\odot$: Dieser Peak entsteht primär aus NS-Binärsystemen mit Gebermassen $M_d < \sim$1,6 M$_\odot$ und orbitalen Perioden $P_{orb} < 20$ Tage. In diesen Systemen ermöglicht die Massentransferrate eine effiziente, quasi-gleichgewichtige sub-Eddington-Akkretion über lange Zeitskalen.
   • Der Peak bei 1,4 M$_\odot$: Dieser Peak wird NS-Systemen zugeschrieben, die eine Entwicklung durch gemeinsame Hüllen durchlaufen (oft unter Beteiligung massereicherer Geber oder spezifischer orbitaler Konfigurationen), wobei das Massewachstum minimal ist.
   • Schwanz hoher Massen: Systeme mit Gebern mittlerer Masse ($1,6 \le M_d \le 3,0$ M$_\odot$) und kurzen orbitalen Perioden können ausreichend Material akkretieren, um 2,1 M$_\odot$ zu überschreiten.
3. Verteilung der Spinperioden: Das Modell sagt NS-Spinperioden voraus, die für Röntgen-Binärsysteme mit geringer Masse (LMXBs) zwischen 1 ms und 15 ms konzentriert sind. Die Ergebnisse sind etwas größer als einige Beobachtungsdaten, was die Autoren auf die Annahme eines minimalen Magnetfelds von $5 \times 10^8$ G zurückführen, die verhindert, dass die Akkretionsscheibe zu nahe an die NS-Oberfläche reicht.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass die Einbeziehung der Physik der Super-Eddington-Akkretion entscheidend ist, um den Peak hoher Massen in der NS-Massenverteilung zu erklären. Indem es NS ermöglicht, auch dann effizient zu akkretieren, wenn die Massentransferraten die Eddington-Grenze überschreiten, löst das Modell die Diskrepanz zwischen theoretischen Vorhersagen geringen Massenzuwachses und der beobachteten Population massereicher NS. Die Autoren behaupten, dass ihr Modell eine kohärente physikalische Erklärung für die bimodale Verteilung liefert, ohne exotische Entstehungskanäle für den Peak hoher Massen zu benötigen, und diesen stattdessen auf die Standard-Binärentwicklung mit Super-Eddington-Akkretion zurückführt.

Einschränkungen und Unsicherheiten
Die Autoren erkennen mehrere Unsicherheiten in ihrem Modell an:

• Akkretionseffizienz: Die angenommene Effizienz von 0,5 und die Behandlung des Massenverlusts in der Nähe des NS sind vereinfacht.
• Abklingrate des Magnetfelds: Die Abklingrate und die minimale Feldstärke sind Parameter, die die Spin-Entwicklung und die Akkretionseffizienz erheblich beeinflussen.
• Zustandsgleichung (EoS): Das Trägheitsmoment ist auf $10^{45}$ g cm$^2$ festgelegt, was nahe am unteren Limit liegt. Variationen der NS-Zustandsgleichung könnten das Trägheitsmoment verändern und sich auf die Zeitskalen des Spin-ups sowie den Übergang in das Propeller-Regime auswirken.
• Ausgeschlossene Kanäle: Die Studie quantifiziert nicht den Beitrag von NS, die durch kollapsinduzierte Akkretion (AIC) entstanden sind oder mit bimodalen Massenverteilungen geboren wurden, sondern konzentriert sich stattdessen auf NS, die durch Standard-Kernkollaps und anschließende Binärentwicklung entstanden sind.