Isochrones in primordial magnetic field evolution

Technisches Resümee: Isochronen in der Entwicklung primordialer Magnetfelder

Problemstellung
Die Entwicklung primordialer Magnetfelder (PMFs) während der strahlungsdominierten Ära des Universums wird durch turbulente Zerfallsprozesse und einen inversen Kaskadeneffekt bestimmt, welcher die charakteristische Längenskala des Feldes vergrößert, während die Feldstärke abnimmt. Ein zentrales Diagnosewerkzeug in diesem Bereich ist die Beziehung zwischen der Alfvén-Geschwindigkeit ($v_A$) und der magnetischen Integralskala ($\xi_M$). Theoretische Modelle legen nahe, dass die Entwicklung der PMFs zu ausreichend späten Zeiten den „Isochronen“ folgt, die durch die Relation $v_A = \xi_M/t^*$ definiert sind, wobei $t^*$ die seit der Entstehung vergangene Zeit darstellt.

Es besteht jedoch eine Diskrepanz zwischen dieser theoretischen Erwartung und numerischen Simulationen. Unterschiedliche Mechanismen der Magnetogenese erzeugen variierende initiale Längenskalen und Feldstärken, was zu Anfangsbedingungen führt, die zum nominellen Simulationsstartzeitpunkt ($t=0$) nicht auf den erwarteten asymptotischen Isochronen liegen. Zudem ist bekannt, dass die tatsächliche Zerfallszeit der magnetohydrodynamischen (MHD) Turbulenz die nominelle Alfvén-Zeit um einen Faktor übersteigt, der von der Lundquist-Zahl ($Lu$) abhängt. Die vorliegende Arbeit adressiert die Herausforderung, diverse Anfangsbedingungen mit universellen Isochronen in Einklang zu bringen und die „korrekte Zeit“ (proper time) zu bestimmen, die erforderlich ist, um die Daten der frühen Simulationsphasen mit den theoretischen Vorhersagen für späte Zeiten abzugleichen.

Methodik
Die Autoren verwenden zweidimensionale numerische Simulationen von zerfallender MHD-Turbulenz unter Verwendung des Pencil Code. Die Studie konzentriert sich auf den 2D-Fall, bei dem die erhaltene Größe die Anastrophie ($I_A = \langle A_z^2 \rangle$) ist, anstatt der magnetischen Helizität (3D helikal) oder des Hosking-Integrals (3D nicht-helikal).

• Anfangsbedingungen: Die Simulationen variieren die initiale Position des magnetischen Energiepeaks ($k_B$) sowie die initiale Feldstärke ($B_0$). Drei primäre initiale Konfigurationen (Punkte A, B und C) werden getestet, um verschiedene Kombinationen von initialem $v_A$ und $\xi_M$ zu untersuchen.
• Geschichtsfelder: Die meisten Durchläufe starten mit einer initialen Geschwindigkeit von Null ($u_0=0$). Eine Teilmenge der Durchläufe führt ein endliches initiales Geschichtsfeld ($u_0 \neq 0$) mit einer relativen Amplitude von $u_0/B_0 \approx 0,3$ ein, um dessen Effekt auf den Evolutionspfad zu testen.
• Spektrale Formen: Die Studie untersucht den Einfluss der initialen Spektralsteigung, indem sie Standard-Inertialbereichssteigungen ($\alpha_{ine} = -2$) mit Fällen vergleicht, in denen der Inertialbereich sehr kurz oder gar nicht vorhanden ist ($\alpha_{ine} = -20$).
• Auflösung und Parameter: Die Simulationen werden mit Auflösungen von $1024^2$ bis $8192^2$ Gitterpunkten durchgeführt, was Lundquist-Zahlen von $Lu \approx 70$ bis $Lu \approx 1100$ abdeckt. Die magnetische Prandtl-Zahl ist auf $Pr_M = 1$ fixiert.

Wesentliche Beiträge und Ergebnisse

1. Definition von korrekter Zeit und Isochronen:
   Die Arbeit zeigt, dass die Anfangsbedingungen zwar beliebig sein können, aber auf universelle Isochronen abgebildet werden können, indem die nominelle Zeitachse zurückgesetzt wird. Die Autoren definieren eine „korrekte Zeit“ ($t_{proper}$) für jede Simulation. Dies wird erreicht, indem die instantane Alfvén-Zeit ($t_A = \xi_M/v_A$) mittels einer linearen Relation $t_{fit}^A = mt + t_0$ an die nominelle Zeit $t$ angepasst wird.

   • Die Steigung $m$ liefert den Zerfallsfaktor $C_M = 1/m$, der quantifiziert, um wie viel länger die Zerfallszeit im Vergleich zur Alfvén-Zeit ist. Die Studie stellt fest, dass $C_M$ etwa zwischen 10 und 24 liegt, mit $Lu$ ansteigt und bei sehr hohen $Lu$-Werten gegen 50 konvergiert.
   • Der Achsenabschnitt $t_0$ (skaliert durch $C_M$) liefert einen Zeitversatz ($t_{offset}$), der, wenn er zur nominellen Zeit addiert wird, die Trajektorie der Simulation mit der theoretischen Isochrone in Einklang bringt.

2. Universelle Pfade:
   Durch Anwendung dieser Zeitkorrekturen zeigen die Autoren, dass die Evolutionspfade aus Simulationen mit völlig unterschiedlichen Werten für $v_A$ und $\xi_M$ auf parallele Isochronen kollabieren. Dies bestätigt, dass die magnetischen Feldparameter selbst in frühen Phasen auf universellen Pfaden liegen, sofern die Zeitachse korrekt verschoben wird.

3. Einfluss der initialen Geschwindigkeit:
   Das Vorhandensein eines initialen Geschichtsfeldes realistischer Stärke ($u_0/B_0 \approx 0,3$) führt zu einem „geraderen“ Evolutionspfad im $v_A$ vs. $\xi_M$-Diagramm. Im Gegensatz dazu zeigen Durchläufe, die mit einer Geschwindigkeit von Null starten, eine initiale Anpassungsphase, in der sich das Geschichtsfeld aufbaut, was zu einer leichten Aufwärtskrümmung des Pfades führt, bevor er sich in das asymptotische Zerfallsverhalten einpendelt.

4. Anastrophie-Erhaltung:
   Die Studie quantifiziert den Zerfall der Anastrophie aufgrund endlicher Resistivität. Der Zerfallsexponent $s$ (wobei $I_A \propto t^{-s}$) nimmt mit steigender Lundquist-Zahl ab und folgt näherungsweise $s \approx 3,1 Lu^{-2/3}$. Für hohe $Lu$-Werte ist die Anastrophie über kosmologische Zeitskalen hinweg nahezu erhalten.

5. Spektrale Evolution:
   Die Simulationen bestätigen, dass sich das Energiespektrum in einer forminvarianten Weise zu niedrigeren Wellenzahlen verschiebt (inverser Kaskadeneffekt). Die Studie stellt zudem fest, dass das System selbst bei einer nicht übereinstimmenden initialen Spektralsteigung (z. B. ohne Inertialbereich) schnell den erwarteten Inertialbereich entwickelt und der $(\xi_M, v_A)$-Punkt schließlich mit dem Pfad eines Standard-Durchlaufs übereinstimmt, wenn auch mit einem leicht unterschiedlichen Wert der erhaltenen Anastrophie.

Bedeutung und Behauptungen
Die Arbeit behauptet, den scheinbaren Konflikt zwischen diversen primordialen Magnetogenese-Szenarien und dem Konzept universeller Isochronen gelöst zu haben. Die Autoren argumentieren, dass die Beschränkung auf zulässige Anfangsbedingungen nicht physikalischer Natur ist, sondern eine Folge einer willkürlichen Zeitdefinition. Durch die Bestimmung des korrekten Zeitversatzes und des Zerfallsfaktors $C_M$ kann jede initiale Magnetfeldkonfiguration auf den universellen Evolutionspfad abgebildet werden.

Die Arbeit betont, dass der Vergleich verschiedener Magnetogenese-Mechanismen ein Zurücksetzen der nominellen Zeit erfordert, sodass der initiale Wert von $v_A t$ mit der initialen Skala $\xi_M$ übereinstimmt. Die Autoren kommen zu dem Schluss, dass trotz Artefakten im Zusammenhang mit endlicher Auflösung und der initialen Anpassungsphase (bei Null-Geschwindigkeits-Durchläufen) die zugrunde liegende Physik des zerfallenden MHD-Turbulenzprozesses zu vorhersagbaren, universellen Isochronen führt, die unabhängig vom spezifischen Erzeugungsmechanismus sind, sofern die Zeitachse korrekt kalibriert wird. Die Studie beansprucht nicht, den Ursprung extragalaktischer Magnetfelder definitiv identifiziert zu haben, liefert jedoch den notwendigen theoretischen Rahmen, um solche Ursprünge gegen Beobachtungsdaten einzugrenzen.