Isochrones in primordial magnetic field evolution

Résumé Technique : Isochrones dans l'évolution des champs magnétiques primordiaux

Énoncé du Problème
L'évolution des champs magnétiques primordiaux (PMF) durant l'ère de la domination du rayonnement de l'Univers est régie par une décroissance turbulente et une cascade inverse, qui augmentent l'échelle de longueur caractéristique du champ tout en diminuant sa force. Un outil de diagnostic central dans ce domaine est la relation entre la vitesse d'Alfvén ($v_A$) et l'échelle d'intégration magnétique ($\xi_M$). Les modèles théoriques suggèrent qu'à des époques suffisamment tardives, l'évolution des PMF suit des « isochrones » définies par la relation $v_A = \xi_M/t^*$, où $t^*$ est le temps écoulé depuis la génération.

Cependant, une divergence existe entre cette attente théorique et les simulations numériques. Différents mécanismes de magnétogenèse produisent des échelles de longueur initiales et des forces de champ variables, conduisant à des conditions initiales qui ne se situent pas sur les isochrones asymptotiques attendues au temps nominal de début de simulation ($t=0$). De plus, le temps de décroissance réel de la turbulence magnétohydrodynamique (MHD) est connu pour excéder le temps d'Alfvén nominal d'un facteur dépendant du nombre de Lundquist ($Lu$). L'article traite du défi consistant à réconcilier diverses conditions initiales avec des isochrones universelles et à déterminer le « temps propre » requis pour aligner les données de simulation précoces avec les prédictions théoriques tardives.

Méthodologie
Les auteurs utilisent des simulations numériques bidimensionnelles de turbulence MHD décroissante à l'aide du code Pencil. L'étude se concentre sur le cas 2D, où la quantité conservée est l'anastrophie ($I_A = \langle A_z^2 \rangle$), plutôt que l'hélicité magnétique (hélicité 3D) ou l'intégrale de Hosking (non-hélicité 3D).

• Conditions Initiales : Les simulations font varier la position initiale du pic du spectre d'énergie magnétique ($k_B$) et la force initiale du champ ($B_0$). Trois configurations initiales principales (Points A, B et C) sont testées pour explorer différentes combinaisons de $v_A$ et $\xi_M$ initiaux.
• Champs de Vélocité : La plupart des calculs commencent avec une vitesse initiale nulle ($u_0=0$). Un sous-ensemble de calculs introduit un champ de vitesse initial fini ($u_0 \neq 0$) avec une amplitude relative $u_0/B_0 \approx 0,3$ pour tester son effet sur la trajectoire évolutive.
• Formes Spectrales : L'étude examine l'impact de la pente spectrale initiale, comparant les pentes de l'intervalle inertiel standard ($\alpha_{ine} = -2$) à des cas présentant des intervalles inertiels très courts ou absents ($\alpha_{ine} = -20$).
• Résolution et Paramètres : Les simulations sont effectuées avec des résolutions allant de $1024^2$ à $8192^2$ points de maillage, couvrant des nombres de Lundquist de $Lu \approx 70$ à $Lu \approx 1100$. Le nombre de Prandtl magnétique est fixé à $Pr_M = 1$.

Contributions Clés et Résultats

1. Définition du Temps Propre et des Isochrones :
   L'article démontre que bien que les conditions initiales puissent être arbitraires, elles peuvent être projetées sur des isochrones universelles en réinitialisant l'axe temporel nominal. Les auteurs définissent un « temps propre » ($t_{proper}$) pour chaque simulation. Ceci est réalisé en ajustant le temps d'Alfvén instantané ($t_A = \xi_M/v_A$) au temps nominal $t$ à l'aide d'une relation linéaire $t_{fit}^A = mt + t_0$.

   • La pente $m$ donne le facteur de décroissance $C_M = 1/m$, qui quantifie à quel point le temps de décroissance est plus long que le temps d'Alfvén. L'étude trouve que $C_M$ varie d'environ 10 à 24, augmentant avec $Lu$ et saturant près de 50 pour des $Lu$ très élevés.
   • L'ordonnée à l'origine $t_0$ (mise à l'échelle par $C_M$) fournit un temps de décalage ($t_{offset}$) qui, lorsqu'il est ajouté au temps nominal, aligne la trajectoire de la simulation avec l'isochrone théorique.

2. Trajectoires Universelles :
   En appliquant ces corrections temporelles, les auteurs montrent que les trajectoires évolutives de simulations ayant des valeurs de $v_A$ et $\xi_M$ initiales très différentes convergent vers des isochrones parallèles. Cela confirme que les paramètres du champ magnétique se situent sur des trajectoires universelles même aux temps précoces, à condition que l'axe temporel soit correctement décalé.

3. Impact de la Vélocité Initiale :
   La présence d'un champ de vélocité initial d'une intensité réaliste ($u_0/B_0 \approx 0,3$) produit une trajectoire évolutive plus « droite » dans le diagramme $v_A$ vs $\xi_M$. En revanche, les calculs commençant avec une vitesse nulle présentent une phase d'ajustement initiale où le champ de vitesse se construit, provoquant une légère courbure ascendante de la trajectoire avant de se stabiliser dans le comportement de décroissance asymptotique.

4. Conservation de l'Anastrophie :
   L'étude quantifie la décroissance de l'anastrophie due à la résistivité finie. L'exposant de décroissance $s$ (où $I_A \propto t^{-s}$) diminue avec l'augmentation du nombre de Lundquist, suivant approximativement $s \approx 3,1 Lu^{-2/3}$. Pour des $Lu$ élevés, l'anastrophie est presque conservée sur des échelles de temps cosmologiques.

5. Évolution Spectrale :
   Les simulations confirment que le spectre d'énergie se déplace vers des nombres d'onde plus bas (cascade inverse) de manière invariante en forme. L'étude note également que même avec une pente spectrale initiale discordante (par exemple, dépourvue d'un intervalle inertiel), le système développe rapidement l'intervalle inertiel attendu, et le point $(\xi_M, v_A)$ finit par s'aligner sur la trajectoire d'un calcul standard, bien qu'avec une valeur d'anastrophie conservée légèrement différente.

Signification et Revendications
L'article prétend résoudre le conflit apparent entre les divers scénarios de magnétogénèse et le concept d'isochrones universelles. Les auteurs soutiennent que la restriction sur les conditions initiales permises n'est pas physique mais résulte d'une définition arbitraire du temps. En déterminant le décalage de temps propre et le facteur de décroissance $C_M$, toute configuration de champ initiale peut être projetée sur la trajectoire évolutive universelle.

Le travail souligne que la comparaison de différents mécanismes de magnétogénèse nécessite de réinitialiser le temps nominal de sorte que la valeur initiale de $v_A t$ concorde avec l'échelle initiale $\xi_M$. Les auteurs concluent que bien que des artefacts liés à la résolution finie et à la phase d'ajustement initiale (dans les calculs à vitesse nulle) existent, la physique sous-jacente de la turbulence MHD décroissante conduit à des isochrones prévisibles et universelles, indépendantes du mécanisme de génération spécifique, pourvu que l'axe temporel soit correctement calibré. L'étude ne prétend pas identifier de manière définitive l'origine des champs magnétiques extragalactiques, mais fournit le cadre théorique nécessaire pour contraindre de telles origines par rapport aux données observationnelles.