Isochrones in primordial magnetic field evolution

Resumen Técnico: Isócronas en la Evolución de los Campos Magnéticos Primordiales

Planteamiento del Problema
La evolución de los campos magnéticos primordiales (PMF) durante la era de radiación dominante del Universo está gobernada por el decaimiento turbulento y una cascada inversa, que aumenta la escala de longitud característica del campo mientras disminuye su intensidad. Una herramienta de diagnóstico central en este campo es la relación entre la velocidad de Alfvén ($v_A$) y la escala integral magnética ($\xi_M$). Los modelos teóricos sugieren que, a tiempos suficientemente tardíos, la evolución de los PMF sigue "isócronas" definidas por la relación $v_A = \xi_M/t^*$, donde $t^*$ es el tiempo transcurrido desde la generación.

Sin embargo, existe una discrepancia entre esta expectativa teórica y las simulaciones numéricas. Diferentes mecanismos de magnetogénesis producen diversas escalas de longitud inicial y fuerzas de campo, lo que conduce a condiciones iniciales que no se encuentran en las isócronas asintóticas esperadas en el tiempo de inicio nominal de la simulación ($t=0$). Además, se sabe que el tiempo real de decaimiento de la turbulencia magnetohidrodinámica (MHD) excede el tiempo de Alfvén nominal por un factor dependiente del número de Lundquist ($Lu$). El artículo aborda el desafío de reconciliar diversas condiciones iniciales con isócronas universales y determinar el "tiempo propio" requerido para alinear los datos de las simulaciones de tiempo temprano con las predicciones teóricas de tiempo tardío.

Metodología
Los autores emplean simulaciones numéricas bidimensionales de turbulencia MHD en decaimiento utilizando el Código Pencil. El estudio se centra en el caso 2D, donde la cantidad conservada es la anastrofía ($I_A = \langle A_z^2 \rangle$), en lugar de la helicidad magnética (3D helicoidal) o el integral de Hosking (3D no helicoidal).

• Condiciones Iniciales: Las simulaciones varían la posición inicial del pico del espectro de energía magnética ($k_B$) y la fuerza de campo inicial ($B_0$). Se prueban tres configuraciones iniciales principales (Puntos A, B y C) para explorar diferentes combinaciones de $v_A$ inicial y $\xi_M$.
• Campos de Velocidad: La mayoría de las ejecuciones comienzan con una velocidad inicial cero ($u_0=0$). Un subconjunto de las ejecuciones introduce un campo de velocidad inicial finito ($u_0 \neq 0$) con una amplitud relativa de $u_0/B_0 \approx 0.3$ para probar su efecto en la trayectoria evolutiva.
• Formas Espectrales: El estudio examina el impacto de la pendiente espectral inicial, comparando pendientes de rango inercial estándar ($\alpha_{ine} = -2$) con casos de rangos inerciales muy cortos o ausentes ($\alpha_{ine} = -20$).
• Resolución y Parámetros: Las simulaciones se realizan con resoluciones que van desde $1024^2$ hasta $8192^2$ puntos de malla, cubriendo números de Lundquist desde $Lu \approx 70$ hasta $Lu \approx 1100$. El número de Prandtl magnético se fija en $Pr_M = 1$.

Contribuciones Clave y Resultados

1. Definición de Tiempo Propio e Isócronas:
   El artículo demuestra que, si bien las condiciones iniciales pueden ser arbitrarias, estas pueden mapearse en isócronas universales mediante el restablecimiento del eje de tiempo nominal. Los autores definen un "tiempo propio" ($t_{proper}$) para cada simulación. Esto se logra ajustando el tiempo de Alfvén instantáneo ($t_A = \xi_M/v_A$) al tiempo nominal $t$ utilizando una relación lineal $t_{fit}^A = mt + t_0$.

   • La pendiente $m$ produce el factor de decaimiento $C_M = 1/m$, que cuantifica cuánto más largo es el tiempo de decaimiento en comparación con el tiempo de Alfvén. El estudio encuentra que $C_M$ oscila aproximadamente entre 10 y 24, aumentando con $Lu$ y saturándose cerca de 50 para valores de $Lu$ muy altos.
   • La intersección $t_0$ (escalada por $C_M$) proporciona un tiempo de desfase ($t_{offset}$) que, al sumarse al tiempo nominal, alinea la trayectoria de la simulación con la isócrona teórica.

2. Trayectorias Universales:
   Al aplicar estas correcciones de tiempo, los autores muestran que las trayectorias evolutivas de simulaciones con valores de $v_A$ inicial y $\xi_M$ vastamente diferentes colapsan sobre isócronas paralelas. Esto confirma que los parámetros del campo magnético se encuentran en trayectorias universales incluso en tiempos tempranos, siempre que el eje de tiempo se ajuste correctamente.

3. Impacto de la Velocidad Inicial:
   Se observa que la presencia de un campo de velocidad inicial de fuerza realista ($u_0/B_0 \approx 0.3$) produce una trayectoria evolutiva más "recta" en el diagrama $v_A$ vs. $\xi_M$. En contraste, las ejecuciones que comienzan con velocidad cero exhiben una fase de ajuste inicial donde el campo de velocidad se construye, causando una ligera curvatura ascendente en la trayectoria antes de establecerse en el comportamiento de decaimiento asintótico.

4. Conservación de la Anastrofía:
   El estudio cuantifica el decaimiento de la anastrofía debido a la resistividad finita. Se encuentra que el exponente de decaimiento $s$ (donde $I_A \propto t^{-s}$) disminuye con el aumento del número de Lundquist, siguiendo aproximadamente $s \approx 3.1 Lu^{-2/3}$. Para un $Lu$ alto, la anastrofía se conserva casi íntegramente durante las escalas de tiempo cosmológicas.

5. Evolución Espectral:
   Las simulaciones confirman que el espectro de energía se desplaza hacia números de onda más bajos (cascada inversa) de una manera invariante en su forma. El estudio también señala que, incluso con un desajuste en la pendiente espectral inicial (por ejemplo, la falta de un rango inercial), el sistema desarrolla rápidamente el rango inercial esperado, y el punto $(\xi_M, v_A)$ eventualmente se alinea con la trayectoria de una ejecución estándar, aunque con un valor de anastrofía conservada ligeramente diferente.

Significado y Reivindicaciones
El artículo afirma resolver el conflicto aparente entre los diversos escenarios de magnetogénesis y el concepto de isócronas universales. Los autores argumentan que la restricción sobre las condiciones iniciales permisibles no es física, sino una consecuencia de una definición de tiempo arbitraria. Al determinar el desfase del tiempo propio y el factor de decaimiento $C_M$, cualquier configuración de campo inicial puede mapearse a la trayectoria evolutiva universal.

El trabajo enfatiza que comparar diferentes mecanismos de magnetogénesis requiere restablecer el tiempo nominal de tal manera que el valor inicial de $v_A t$ coincida con la escala inicial $\xi_M$. Los autores concluyen que, si bien existen artefactos relacionados con la resolución finita y la fase de ajuste inicial (en ejecuciones con velocidad cero), la física subyacente del decaimiento de la turbulencia MHD conduce a isócronas universales predecibles que son independientes del mecanismo de generación específico, siempre que el eje de tiempo se calibre correctamente. El estudio no pretende identificar definitivamente el origen de los campos magnéticos extragalácticos, sino que proporciona el marco teórico necesario para constreñir tales orígenes frente a los datos observacionales.