Relativistic magnetohydrodynamics in the early Universe

Technische Samenvatting: Relativistische Magnetohydrodynamica in het Vroege Universum

Probleemstelling
De studie van primordiale magnetische velden, fluïdumperturbaties en de productie van gravitatiegolven in het vroege universum steunt vaak op magnetohydrodynamische (MHD) simulaties. Historisch gezien opereerden deze simulaties onder de aanname van subrelativistische bulkvloeistofbeweging ($u^2 \ll 1$) en gingen ze frequent uit van de limiet waarbij de Lorentzfactor $\gamma^2 \approx 1$. Hoewel deze benadering geldt voor materiedominante fluïda, introduceert het ongeëxamineerde fouten in stralingsgedomineerde tijdperken of wanneer de toestandsvergelingsparameter $c_s^2$ van de orde eenheid is (bijv. $c_s^2 = 1/3$). Eerdere formuleringen hebben vaak de tijdderivaten van de Lorentzfactor ($\partial_\tau \gamma^2$) verwaarloosd, wat leidde tot ontbrekende niet-lineaire termen van de orde $U^2/\ell$ (waarbij $U$ en $\ell$ karakteristieke snelheid- en lengteschalen zijn). Deze omissies beïnvloeden de behoudswetten van energie en impuls, de generatie van vorticiteit en de voortplantingssnelheden van MHD-golven, wat potentieel kan leiden tot superluminale Alfvén-snelheden indien niet gecorrigeerd.

Methodologie
De auteurs afleiden de behoudswetten voor MHD in een homogene, isotrope en expanderende Friedmann–Litemître–Robertson–Walker (FLRW) achtergrond. Het werk verloopt via verschillende theoretische stappen:

1. Geometrisch Kader: Het FLRW-metriek wordt vastgesteld met behulp van conforme tijd $\tau$ en comoving coördinaten, wat een generieke $\alpha$-tijddefinitie mogelijk maakt om verschillende schalingseigenschappen van fluïdumvariabelen te verkennen.
2. Perfect Fluïdum Dynamica: De auteurs leiden de bewegingsvergelijkingen af voor een perfect fluïdum met een constante toestandsvergelding $p = c_s^2 \rho$. Zij presenteren deze vergelijkingen in twee vormen:
   • Behoudsvorm: Evolutie van de stress-energie-tensorcomponenten $\tilde{T}^{0\mu}$.
   • Niet-behoudsvorm: Evolutie van de primitieve variabelen (comoving energiedichtheid $\tilde{\rho}$ en peculiaire snelheid $u$).
     Cruciaal is dat de afleiding termen bevat die $\partial_\tau \gamma^2$ en $\partial_\tau u^2$ omvatten, die vaak worden weggelaten in de subrelativistische limiet.
3. Imperfecte Fluïda: Het kader wordt uitgebreid om eerst-orde imperfecte fluïdumeffecten (viscositeit en warmtegeleiding) te bevatten met behulp van de Klassieke Irreversibele Thermodynamica (CIT), waarbij covariante expressies voor Navier-Stokes viscositeit en Fouriers thermische geleidbaarheid worden geboden.
4. Elektromagnetisme: Maxwells vergelijkingen worden herzien in de expanderende achtergrond, waarbij comoving elektromagnetische velden ($\tilde{E}, \tilde{B}$) en een covariante gegeneraliseerde wet van Ohm worden geïntroduceerd. De verplaatsingsstroom wordt geanalyseerd om de ideale MHD-limiet in het hoog-geleidende regime van het vroege universum te rechtvaardigen.
5. MHD Koppeling: De fluïdumvergelijkingen worden gekoppeld aan de Maxwell-vergelijkingen via de Lorentz-kracht. De auteurs analyseren lineaire perturbaties om geluids-, Alfvén- en magnetosonische golven te bestuderen.
6. Correcties: Een "Boris-correctie" wordt aangepast aan de relativistische MHD-context om te waarborgen dat Alfvén-snelheden subluminaal blijven, zelfs wanneer de verplaatsingsstroom wordt genegeerd.

Belangrijkste Bijdragen en Resultaten

• Gecorrigeerde Subrelativistische Vergelijkingen: Het artikel toont aan dat voorgaande subrelativistische MHD-vergelijkingen (bijv. die gebruikt in de Pencil Code) fouten bevatten in de coëfficiënten van convectieve termen $(u \cdot \nabla) \ln \tilde{\rho}$ en snelheidstermen. Specifiek introduceert de tijdsafgeleide van $\gamma^2$ een correctiefactor van $1/2$ in bepaalde termen voor stralingsgedomineerde fluïda ($c_s^2 = 1/3$) en wijzigt de coëfficiënten voor generieke $c_s^2$. De gecorrigeerde niet-behoudsvorm vergelijkingen worden gepresenteerd in Eq. (1.2) en (6.31).
• Vorticiteitsproductie: De auteurs tonen aan dat vorticiteit ($\omega = \nabla \times u$) geen topologische invariant is in het vroege universum. Zelfs in de afwezigheid van externe forcering en barocline termen, kan vorticiteit worden gegenereerd uit een initieel rotatievrij snelheidveld door relativistische correcties die $c_s^2 \neq 0$ en de expansie van het universum omvatten.
• Schaling en Conforme Invariantie: Het artikel beoordeelt diverse schalingen van fluïdumvariabelen (comoving, super-comoving) om Hubble-frictietermen te minimaliseren. Het identificeert dat voor subrelativistische stromen met $c_s^2 \neq 1/3$, een specifieke schaling ($\beta = 3(1+c_s^2)$) de energievergelijking conform afvlak is maakt, terwijl de impulsvergelijking een Hubble-frictieterm behoudt.
• Relativistische Alfvén-snelheid en Boris-correctie: De analyse onthult dat het negeren van de verplaatsingsstroom in ideale MHD kan leiden tot superluminale Alfvén-snelheden ($v_A > 1$) wanneer $B_0^2/\rho_0$ groot is. De auteurs passen de Boris-correctie aan de relativistische MHD-vergelijkingen aan, waarbij de impulsvergelijking wordt gewijzigd om te garanderen dat $v_A \to v_A \gamma_A \leq 1$, waarbij $\gamma_A = (1+v_A^2)^{-1}$.
• Transportcoëfficiënten: Schattingen voor schuifviscositeit, bulkviscositeit en thermische geleidbaarheid in het primordiale plasma worden herzien, wat bevestigt dat het plasma op grote schaal effectief een perfect fluïdum is, hoewel dissipatieve effecten noodzakelijk zijn voor realistische turbulentiemodellering.

Significantie
Dit werk biedt de eerste volledig relativistische beschrijving van MHD-vergelijkingen in een expanderende achtergrond die ook de correcties voor de subrelativistische limiet bevat die voorheen over het hoofd werden gezien. De auteurs stellen dat deze correcties niet louter hogere-orde relativistische effecten zijn, maar van dezelfde orde als convectieve derivaten in niet-lineaire fluïdumdynamica. Bijgevolg beargumenteert het artikel dat toekomstige simulaties van magnetohydrodynamica in het vroege universum, met name die betrokken bij stralingsgedomineerde fluïda of hoogenergetische fasen, deze gecorrigeerde termen moeten opnemen om niet-lineaire dynamica, turbulentie en vorticiteitsgeneratie accuraat te modelleren. De afgeleide vergelijkingen dienen als een theoretische fundering voor de komende ontwikkelingen in numerieke codes zoals de Pencil Code en CosmoLattice.