Analytic insight into the physics of SASI II. Spiral instability of the prograde mode in a rotating stellar core

Technische Zusammenfassung: Analytischer Einblick in die Physik der SASI II. Spirale Instabilität des prograden Modus in einem rotierenden Sternkern

Problemstellung
Während des Kernkollapses eines rotierenden massereichen Sterns treibt die Standing Accretion Shock Instability (SASI) nicht-achsensymmetrische Bewegungen an, die spezifische Signaturen in Neutrino- und Gravitationswellensignalen hinterlassen. Während frühere Studien etablierten, dass Rotation das Wachstum prograder Spiralmodi begünstigt, blieb der physikalische Mechanismus hinter dieser Destabilisierung ungeklärt, insbesondere in Regimen, in denen die Rotationsraten langsamer als die SASI-Frequenz sind. Zudem beruhten bestehende perturbatieve Analysen oft auf zylindrischen oder sphärisch-äquatorialen Approximationen, welche die poloidale Struktur des Flusses vernachlässigten, was deren Genauigkeit bei der Interpretation numerischer Simulationen rotierender Vorläufersterne einschränkte. Zusätzlich fehlte für den Einfluss der Rotation auf die Oszillationsfrequenzen dominanter SASI-Modi ein robuster analytischer Rahmen.

Methodik
Diese Studie erweitert die perturbative Analyse gestauter Akkretionsschocks in sphärischer Geometrie um die Effekte von Rotation und differentieller Rotation. Der Ansatz umfasst:

1. Modellierung des stationären Flusses: Es wird ein vereinfachtes Modell der Äquatorialebene konstruiert, das einen idealen Gaszustand mit einem adiabatischen Index $\gamma = 4/3$ annimmt, der von relativistischem Elektrondegenerationsdruck dominiert wird. Der Fluss beinhaltet eine Kühlfunktion, die von Dichte und Druck abhängt, wobei das Gravitationspotenzial in der Newtonschen Näherung behandelt wird.
2. Lineare Perturbationsanalyse: Die linearen Perturbationen werden durch Spherical Harmonics $Y^m_\ell$ mit poloidaler Zahl $\ell$ und Azimutalzahl $m$ charakterisiert. Die Analyse nutzt einen Satz gekoppelter Differentialgleichungen für den Massenfluss, die Entropie, den spezifischen Drehimpuls und die Vortizität.
3. Erhaltungseigenschaften: Um die Schwierigkeiten bei der Lösung der poloidalen Struktur in einem rotierenden Vorläuferstern zu umgehen, nutzt die Studie eine Erhaltungseigenschaft der radialen Vortizität. Dies ermöglicht es, das System zu reduzieren und numerisch mittels einer Shooting-Methode zu lösen.
4. Adiabatische Näherung: Ein adiabatisches Modell wird entwickelt, um den Beitrag adiabatischer Prozesse zur Destabilisierung zu isolieren. Dies reduziert die Ordnung des Differentialsystems und ermöglicht eine Integralgleichungsformulierung für die Eigenfrequenzen, indem der post-shock akustische Bereich als selbst erzwungener Oszillator behandelt wird.
5. Analytische Näherung: Die Studie leitet explizite analytische Ausdrücke für die Eigenfrequenzen der dominanten SASI-Modi ($\ell=m=1$ und $\ell=m=2$) durch Anpassung der numerischen Ergebnisse über einen Bereich von Schockradien und Rotationsraten her.

Zentrale Beiträge und Ergebnisse

• Analytische Eigenfrequenzformeln: Die Arbeit liefert analytische Approximationen für den Realteil der Eigenfrequenzen ($\omega_r$) der dominanten Spiralmodi ($\ell=m=1$ und $\ell=m=2$) mit einer Genauigkeit von etwa 10 %, wenn der Schockradius $r_{sh}$ etwa das 1,5-fache des Radius der maximalen Verzögerung ($r_\nabla$) überschreitet. Die Frequenz wird als Funktion des Schockradius und des spezifischen Drehimpulses $L$ ausgedrückt:
  $$\omega_{fit}^r(\ell, m) \sim \frac{|v_{sh}|}{r_{sh}} F_{\ell,m}\left(\frac{r_{sh}}{r_\nabla}\right) + \frac{mL}{r_{sh}r_\nabla} G_{\ell,m}\left(\frac{r_{sh}}{r_\nabla}\right)$$
  wobei $F$ und $G$ dimensionslose Funktionen sind, die aus numerischen Fits abgeleitet wurden.
• Mechanismus der Destabilisierung: Die Studie erklärt, dass die Destabilisierung prograder Moden durch differentielle Rotation durch Phasenmischung zwischen der advektiven Anregung und der akustischen Struktur getrieben wird. Differentielle Rotation erhöht die radiale Wellenlänge der Vortizitäts-Perturbationen, die mit dem prograden Modus assoziiert sind. Dies führt zu einer besseren Phasenanpassung mit der großen radialen Skala der akustischen Struktur, wodurch die Kopplungseffizienz erhöht wird.
• Rolle der Rotationsrate:
  • In Regimen moderater Rotation (in denen kein Korotationsradius existiert) ist der nachteilige Effekt auf die Phasenmischung bei kleinen Radien signifikant aufgrund des steilen nach innen gerichteten Anstiegs der Rotationsfrequenz ($\propto 1/r^2$).
  • In Regimen stärkerer Rotation mit einem Korotationsradius identifiziert die stationäre Phasenapproximation die dominante advektive-akustische Kopplung zwischen der Korotationszone und dem Schock.
• Modus-Dominanz: Die Analyse bestätigt, dass die Wachstumsrate des prograden Modus näherungsweise linear mit dem spezifischen Drehimpuls ansteigt. Mit zunehmender Rotation kann der fundamentale Modus $\ell=m=2$ gegenüber $\ell=m=1$ dominant werden. Das Verhältnis ihrer Frequenzen hängt von der Rotationsrate ab und kann im Vergleich zum nicht-rotierenden Fall um das bis zu 2,3-fache steigen.
• Unabhängigkeit von Kühlungsdetails: Es wurde festgestellt, dass die Oszillationsfrequenz des dominanten Modus schwach von der detaillierten Formulierung der Neutrinoemission (Parameter der Kühlfunktion) abhängt, sofern der Schockradius ausreichend groß ist ($r_{sh} \gtrsim 1.5 r_\nabla$). Die Eigenschaften werden am besten erfasst, wenn der Schockradius in Einheiten von $r_\nabla$ statt des Radius des Proto-Neutronensterns gemessen wird.

Bedeutung
Die Arbeit beansprucht, die ersten analytischen Approximationen für SASI-Eigenfrequenzen bereitzustellen, die den Effekt der Rotation explizit einschließen. Durch die Erweiterung des „selbst erzwungenen Oszillator“-Rahmens auf rotierende Flüsse bietet sie eine physikalische Erklärung für die lineare Verstärkung der prograden Moden-Wachstumsraten durch differentielle Rotation – ein Phänomen, das zuvor numerisch beobachtet, aber analytisch nicht verstanden wurde. Die abgeleiteten Formeln ermöglichen die potenzielle Rückgewinnung stellarer Eigenschaften (wie den Schockradius und den spezifischen Drehimpuls) aus beobachteten Neutrino- und Gravitationswellenfrequenzen. Die Studie hebt zudem die Einschränkungen früherer zylindrischer oder äquatorialer Approximationen hervor, indem sie die Bedeutung der vollständigen poloidalen Struktur für die genaue Interpretation von SASI-Frequenzen in rotierenden Vorläufersternen demonstriert.