Disk-Regulated Mass Transfer Between Rotating Non-Degenerate Stars: Insights from Be and sdOB Binaries

Technische Zusammenfassung: Scheibenregulierte Massentransfer zwischen rotierenden nicht-degenerierten Sternen

Problemstellung
Der Massentransfer zwischen nicht-degenerierten Sternen ist ein kritischer, jedoch noch unzureichend verstandener Prozess in der Binärentwicklung. Eine zentrale Unsicherheit liegt in der Akkretionseffizienz – dem Bruchteil der übertragenen Masse, den der Akkretor beibelt. Standardmäßige Binärpopulation-Synthesen (BPS) und detaillierte Sternentwicklungs-Simulationen (z. B. unter Verwendung von MESA) verwenden typischerweise ein „rotationsbegrenztes“ Akkretionsschema. In diesem Rahmen wird die Akkretion eingeschränkt, sobald der Akkretor sich seiner kritischen Breakup-Rotationsrate nähert, was häufig zu sehr niedrigen Akkretionseffizienzen führt.

Bestehende Beobachtungsdaten von Be-Stern-Binärsystemen, insbesondere Be+subdwarf O/B-Typ (sdOB) Systemen, deuten jedoch auf signifikant höhere Akkretionseffizienzen hin (oft $\gtrsim 30\%$ bis $>50\%$). Bestehende rotationslimitierte Modelle unterschätzen systematisch die Massen von Be-Sternen in diesen Systemen und können den beobachteten Parameterraum von Be+sdOB-Binärsystemen nicht reproduzieren. Vorherige Studien haben zudem gezeigt, dass weder rotationslimitierte noch thermisch limitierte Preskriptionen allein die Effizienzanforderungen erfüllen können, die aus beobachteten Be+sdOB-Populationen abgeleitet wurden.

Methodik
Die Autoren schlagen eine physikalisch motivierte „Disk-Stern-Kopplung“ zur Massenakkretions-Preskription vor, um diese Diskrepanzen zu adressen. Der zentrale physikalische Mechanismus postuliert, dass eine Akkretions- oder Dekreationsscheibe den Transport des Drehimpulses zum Akkretor reguliert, wodurch die Massenakkretion fortgesetzt werden kann, selbst wenn der Stern nahe seiner kritischen Rotationsgrenze rotiert.

1. Drehimpulsregulation: Basierend auf den stationären dünnen Scheibenlösungen von Popham & Narayan (1991) modellieren die Autoren den spezifischen Drehimpuls, der an den Stern akkretiert wird ($j_{acc}$), als Funktion der Sternrotationsrate ($\Omega/\Omega_{crit}$).

   • Bei niedrigen Rotationsraten entspricht $j_{acc}$ in etwa dem Kepler'schen Wert an der Sternoberfläche.
   • Wenn der Stern sich einer Schwellenwert-Rotationsrate nähert (angenommen als $\Omega/\Omega_{crit} \approx 0,9$), ermöglicht der Scherstress der Scheibe, dass der Drehimpuls vom Stern nach außen in die Scheibe transportiert wird. Infolgedessen sinkt $j_{acc}$ drastisch ab, wird potenziell negativ, was eine kontinuierliche Massenakkretion ermöglicht, ohne dass der Stern seine kritische Rotation überschreitet.
   • Die Autoren implementieren eine geglättete Funktionsform für $j_{acc}$, um die numerische Stabilität in 1D-Simulationen zu gewährleisten und Diskontinuitäten zu vermeiden, die durch einen strikten Schwellenwert entstehen würden.

2. Erhaltungssätze und Effizienzberechnung: Das Modell nimmt an, dass ein Überschuss an Drehimpuls, der aus dem Stern extrahiert wird, durch die Scheibe transportiert und über Massenverlust am äußeren Rand der Scheibe ($R_{out}$) entfernt wird. Durch Anwendung der Massen- und Drehimpulserhaltung wird die Akkretionseffizienz ($\eta = \dot{M}_{acc}/\dot{M}_{tr}$) konsistent berechnet. Die Effizienz hängt vom Verhältnis des spezifischen Drehimpulses am Scheibenrand zum spezifischen Drehimpuls am Akkretionsradius ab.

3. Simulationsrahmen: Die Autoren haben diese Preskription innerhalb des POSYDON-Frameworks implementiert und dabei MESA (Modules for Experiments in Stellar Astrophysics) verwendet, um Binärmodelle zu entwickeln.

   • Sie deaktivierten das implizite Wind-Schema, das zur Vermeidung von superkritischer Rotation dient, und verließen sich stattdessen auf explizite rotationsverstärkte Winde sowie die reduzierte Akkretionsrate im superkritischen Regime, um eine subkritische Rotation wiederherzustellen.
   • Ein Grid von Binärmodellen wurde mit Sonnenmetallizität simuliert, wobei die Primärmassen ($2-9 M_\odot$), Massenverhältnisse ($0,4-0,95$) und Orbitperioden variiert wurden.
   • Die Studie untersuchte auch den Einfluss von Core-Overshooting-Behandlungen (Vergleich zwischen Standard-POSYDON-Einstellungen und einer reduzierten Overshooting-Preskription) auf die finalen Komponentenmassen.

Wesentliche Beiträge und Ergebnisse

• Analytische Abschätzungen: Die Autoren haben analytische Ausdrücke für die Akkretionseffizienz hergeleitet. Sie zeigen, dass die Effizienz im Disk-Akkretions-Regime von der Sternrotation, dem Massenverhältnis und dem Verhältnis des Sternradius zum Roche-Lobe-Radius abhängt. Im Direct-Impact-Regime hängt sie von der Rotation und dem Massenverhältnis ab. Die gesamten massengewichteten Akkretionseffizienzen liegen nahe bei den Werten, die in der Nähe der Schwellenwert-Rotationsrate erwartet werden, wo der spezifische Akkretionsdrehimpuls stark abfällt.
• Vergleich der Binärentwicklung:
  • Rotationslimitiertes Modell: In diesem Standardmodell erreichen Akkretoren ihre kritische Rotation nach der Akkretion nur eines kleinen Bruchteils der übertragenen Masse ($\sim 0,2 M_\odot$ in einem Testfall), was zu niedrigen Effizienzen führt ($<10\%$ für die meisten Systeme). Dies führt dazu, dass Be-Sterne systematisch weniger massereich sind als beobachtet.
  • Disk-Stern-Kopplungsmodell: Dieses Modell erlaubt es Akkretoren, signifikant mehr Masse zu akkretieren (z. B. $\sim 1,6 M_\odot$ in der schnellen Phase des Testfalls), wobei die Effizienzen $\sim 40-60\%$ erreichen. Die Scheibe extrahiert den Drehimpuls, was es dem Stern ermöglicht, sich auszudehnen und zu akkretieren, während er nahe der kritischen Rotationsschwelle bleibt.
• Reproduktion von Beobachtungen:
  • Bei Anwendung auf Be+sdOB-Systeme scheitert das rotationslimitierte Modell daran, Be-Sterne mit Massen zu erzeugen, die mit Beobachtungen übereinstimmen, insbesondere bei Systemen mit weiten Orbits (Fall B).
  • Das Disk-Stern-Kopplungsmodell erzeugt Be-Sterne, deren Massen einen signifikanten Teil des beobachteten Bereichs abdecken.
  • Die Autoren stellten fest, dass die finalen Komponentenmassen sensitiv gegenüber dem angenommenen Overshooting sind. Ein Modell, das die Disk-Stern-Kopplungs-Preskription mit reduziertem Overshooting kombiniert, liefert Komponentenmassen für Be+sdOB-Systeme, die in besserer Übereinstimmung mit Beobachtungen stehen als die Standardmodelle mit hohem Overshooting.
• Orbitalperioden: Während das Disk-Stern-Kopplungsmodell die Massenerhaltung erhöht, sind die resultierenden Orbitalperioden vergleichbar mit oder etwas länger als die des rotationslimitierten Modells. Dies liegt daran, dass der konservativere Massentransfer den gesamten Massenverlust aus dem System reduziert, was zu extremeren Massenverhältnissen führt, die die Orbitexpansion beschleunigen und den Drehimpulsenverlust durch den Scheibenrand teilweise kompensieren.

Bedeutung und Ansprüche
Das Paper behauptet, dass die Diskrepanz zwischen beobachteten Be+sdOB-Binärsystemeigenschaften und aktuellen Simulationen durch die Verwendung eines diskregulierten Massentransfermodells gelöst werden kann. Die Autoren argumentieren:

1. Rotationslimitierte Akkretionspreskriptionen unterschätzen systematisch die Massen von Be-Sternen in Be+sdOB-Systemen.
2. Ein Disk-Stern-Kopplungsmechanismus, der eine Regulierung des Drehimpulstransports ermöglicht, erlaubt höhere Akkretionseffizienzen, die mit Beobachtungsdaten konsistent sind (speziell mit der Anforderung von Effizienzen $>50\%$ in vielen Systemen).
3. Die Kombination dieser neuen Akkretionspreskription mit reduziertem Core-Overshooting bietet eine physikalisch konsistentere Erklärung für die beobachteten Massen und Orbitalperioden von Be+sdOB-Binärsystemen.

Die Autoren räumen ein, dass die tatsächliche Effizienz von unsicheren Faktoren wie der Scheibenstruktur, Magnetfeldern und den spezifischen Bedingungen der Scheibenbildung und des Massenverlusts (z. B. via L2-Punkte oder Hot Spots) abhängt. Sie kommen jedoch zu dem Schluss, dass das Disk-Stern-Kopplungsmodell einen notwendigen physikalischen Mechanismus darstellt, um die Existenz massereicher Be-Sterne in weiten Binärsystemen zu erklären, in denen der tidal Spin-Down ineffizient ist.