Isochrones in primordial magnetic field evolution

Resumo Técnico: Isócronas na Evolução de Campos Magnéticos Primordiais

Enunciado do Problema
A evolução dos campos magnéticos primordiais (PMFs) durante a era de radiação do Universo é governada pela decadência turbulenta e por um cascateamento inverso, que aumenta a escala de comprimento característica do campo enquanto diminui sua intensidade. Uma ferramenta de diagnóstico central nesta área é a relação entre a velocidade de Alfvén ($v_A$) e a escala integral magnética ($\xi_M$). Modelos teóricos sugerem que, em tempos suficientemente tardios, a evolução dos PMFs segue "isócronas" definidas pela relação $v_A = \xi_M/t^*$, onde $t^*$ é o tempo decorrido desde a geração.

No entanto, existe uma discrepância entre essa expectativa teórica e as simulações numéricas. Diferentes mecanismos de magnetogênese produzem variadas escalas de comprimento iniciais e intensidades de campo, levando a condições iniciais que não se situam nas isócronas assintóticas esperadas no tempo nominal de início da simulação ($t=0$). Além disso, sabe-se que o tempo real de decadência da turbulência magnetohidrodinâmica (MHD) excede o tempo de Alfvén nominal por um fator dependente do número de Lundquist ($Lu$). O artigo aborda o desafio de reconciliar diversas condições iniciais com isócronas universais e determinar o "tempo próprio" necessário para alinhar os dados de simulação de tempo inicial com as previsões teóricas de tempo tardio.

Metodologia
Os autores empregam simulações numéricas bidimensionais de turbulência MHD em decaimento utilizando o Pencil Code. O estudo foca no caso 2D, onde a grandeza conservada é a anastrofia ($I_A = \langle A_z^2 \rangle$), em vez da helicidade magnética (3D helicoidal) ou do integral de Hosking (3D não-helicoidal).

• Condições Iniciais: As simulações variam a posição inicial do pico do espectro de energia magnética ($k_B$) e a força inicial do campo ($B_0$). Três configurações iniciais principais (Pontos A, B e C) são testadas para explorar diferentes combinações de $v_A$ inicial e $\xi_M$.
• Campos de Velocidade: A maioria das execuções começa com velocidade inicial zero ($u_0=0$). Um subconjunto das execuções introduz um campo de velocidade inicial finito ($u_0 \neq 0$) com uma amplitude relativa de $u_0/B_0 \approx 0.3$ para testar seu efeito na trajetória evolutiva.
• Formas Espectrais: O estudo examina o impacto da inclinação espectral inicial, comparando inclinações de intervalo inercial padrão ($\alpha_{ine} = -2$) contra casos com intervalos inerciais muito curtos ou ausentes ($\alpha_{ine} = -20$).
• Resolução e Parâmetros: As simulações são realizadas com resoluções variando de $1024^2$ a $8192^2$ pontos de malha, cobrindo números de Lundquist de $Lu \approx 70$ a $Lu \approx 1100$. O número de Prandtl magnético é fixado em $Pr_M = 1$.

Principais Contribuições e Resultados

1. Definição de Tempo Próprio e Isócronas:
   O artigo demonstra que, embora as condições iniciais possam ser arbitrárias, elas podem ser mapeadas em isócronas universais através do reajuste do eixo de tempo nominal. Os autores definem um "tempo próprio" ($t_{proper}$) para cada simulação. Isso é alcançado ajustando o tempo de Alfvén instantâneo ($t_A = \xi_M/v_A$) ao tempo nominal $t$ usando uma relação linear $t_{fit}^A = mt + t_0$.

   • A inclinação $m$ fornece o fator de decadência $C_M = 1/m$, que quantifica o quanto o tempo de decadência é mais longo em comparação ao tempo de Alfvén. O estudo encontra que $C_M$ varia de aproximadamente 10 a 24, aumentando com $Lu$ e saturando próximo a 50 para $Lu$ muito alto.
   • O intercepto $t_0$ (escalonado por $C_M$) fornece um tempo de deslocamento ($t_{offset}$) que, quando adicionado ao tempo nominal, alinha a trajetória da simulação com a isócrona teórica.

2. Trajetórias Universais:
   Ao aplicar essas correções de tempo, os autores mostram que trajetórias evolutivas de simulações com valores vastamente diferentes de $v_A$ e $\xi_M$ colapsam em isócronas paralelas. Isso confirma que os parâmetros do campo magnético residem em trajetórias universais mesmo em tempos iniciais, desde que o eixo do tempo seja devidamente deslocado.

3. Impacto da Velocidade Inicial:
   A presença de um campo de velocidade inicial de força realista ($u_0/B_0 \approx 0.3$) é encontrada como produtora de uma trajetória evolutiva mais "reta" no diagrama $v_A$ vs. $\xi_M$. Em contraste, execuções começando com velocidade zero exibem uma fase de ajuste inicial onde o campo de velocidade se constrói, causando uma leve curvatura ascendente na trajetória antes de se estabelecer no comportamento de decadência assintótica.

4. Conservação da Anastrofia:
   O estudo quantifica a decadência da anastrofia devido à resistividade finita. O expoente de decadência $s$ (onde $I_A \propto t^{-s}$) é encontrado como diminuindo com o aumento do número de Lundquist, seguindo aproximadamente $s \approx 3.1 Lu^{-2/3}$. Para $Lu$ alto, a anastrofia é quase conservada ao longo de escalas de tempo cosmológicas.

5. Evolução Espectral:
   As simulações confirmam que o espectro de energia se desloca para números de onda menores (cascateamento inverso) de maneira invariante de forma. O estudo também observa que, mesmo com uma inclinação espectral inicial incompatível (por exemplo, carecendo de um intervalo inercial), o sistema desenvolve rapidamente o intervalo inercial esperado, e o ponto $(\xi_M, v_A)$ eventualmente se alinha com a trajetória de uma execução padrão, embora com um valor de anastrofia conservada ligeiramente diferente.

Significância e Alegações
O artigo alega resolver o conflito aparente entre diversos cenários de magnetogênese inicial e o conceito de isócronas universais. Os autores argumentam que a restrição sobre as condições iniciais permitidas não é física, mas sim uma consequência de uma definição de tempo arbitrária. Ao determinar o deslocamento de tempo próprio e o fator de decadência $C_M$, qualquer configuração de campo inicial pode ser mapeada para a trajetória evolutiva universal.

O trabalho enfatiza que comparar diferentes mecanismos de magnetogênese requer o reajuste do tempo nominal de modo que o valor inicial de $v_A t$ concorde com a escala inicial $\xi_M$. Os autores concluem que, embora existam artefatos relacionados à resolução finita e à fase de ajuste inicial (em execuções de velocidade zero), a física subjacente da turbulência MHD em decaimento leva a isócronas previsíveis e universais que são independentes do mecanismo de geração específico, desde que o eixo do tempo seja corretamente calibrado. O estudo não pretende identificar definitivamente a origem dos campos magnéticos extragalácticos, mas fornece o arcabouço teórico necessário para restringir tais origens frente a dados observacionais.