Matter power spectrum induced by primordial magnetic fields: from the linear to the non-linear regime

Basic Information

• Paper ID: 2410.02676
• Title: Matter power spectrum induced by primordial magnetic fields: from the linear to the non-linear regime
• Authors: Pranjal Ralegankar, Enrico Garaldi, Matteo Viel
• Classification: astro-ph.CO (Cosmology and Non-Galactic Astrophysics)
• Publication Date: October 2024 (Submitted to JCAP)
• Paper Link: https://arxiv.org/abs/2410.02676

Abstract

This paper presents the first study using dedicated magnetohydrodynamic (MHD) simulations to investigate the effects of primordial magnetic fields (PMFs) on the matter power spectrum from the linear to the non-linear regime. While linear theory predicts that PMFs enhance the matter power spectrum at small scales, at sufficiently small scales, PMF-induced baryon perturbations backreact on the magnetic field, invalidating the linear approximation. The study reveals that in the non-linear regime, the dimensionless baryon power spectrum saturates to values of $\mathcal{O}(1)$, contrasting with previous theoretical expectations based on the magnetic Jeans scale. Furthermore, the results indicate that previous studies overestimated the total matter power spectrum by several orders of magnitude near the non-linear transition.

Research Background and Motivation

Problem Background

1. Significance of Primordial Magnetic Fields: Magnetic fields are ubiquitous in the universe, observed from galaxies and galaxy clusters to cosmic voids. PMFs may originate from early universe physics such as inflation or phase transitions, providing a unique window into early universe physics.
2. Impact on Structure Formation: PMFs induce large baryon density perturbations at small scales through the Lorentz force, affecting the matter power spectrum. The absence of this enhancement effect can be used to constrain PMF strength.
3. Theoretical Limitations: Existing constraints rely primarily on linear theory, but at small scales, PMF-induced perturbations become sufficiently large to backreact on the magnetic field itself, rendering linear theory inapplicable.

Research Motivation

• First application of complete non-linear MHD equations to study PMF effects on density perturbations
• Verification and correction of theoretical expectations based on the magnetic Jeans scale
• Provision of numerical foundations for more accurate PMF constraints

Core Contributions

1. First MHD Simulations: Conducted the first coupled gravitational-MHD simulations starting from the recombination epoch to study time-dependent effects of PMFs on the matter power spectrum
2. Non-linear Behavior Discovery: Found that the baryon power spectrum saturates to a constant $\mathcal{O}(1)$ value in the non-linear regime rather than being suppressed, challenging conventional understanding of the magnetic Jeans scale
3. Batchelor Spectrum Study: First calculation of the matter power spectrum for PMFs with Batchelor spectrum, extending the parameter space of PMF research
4. Semi-analytical Fitting: Provided semi-analytical fitting formulas combining linear theory and numerical results for use in future structure formation studies
5. Constraint Revision: Discovered that previous studies overestimated the matter power spectrum by several orders of magnitude near the non-linear transition, necessitating reassessment of PMF constraints

Methodology Details

Problem Definition

Investigation of how PMFs affect the evolution of baryon and dark matter density perturbations in the post-recombination universe, with particular focus on the transition from linear to non-linear regimes.

Theoretical Framework

Governing Equations

Post-recombination, the baryon fluid can be treated as an ideal conductor, with magnetic field evolution governed by ideal MHD equations:

$\frac{\partial \vec{B}}{\partial t} = \frac{1}{a}\nabla \times (\vec{v}_b \times \vec{B})$

Baryon motion is affected by PMF Lorentz forces:

$\frac{\partial \vec{v}_b}{\partial t} + H\vec{v}_b + \frac{(\vec{v}_b \cdot \nabla)\vec{v}_b}{a} + \frac{c_b^2}{a}\nabla\delta_b = \frac{(\nabla \times \vec{B}) \times \vec{B}}{4\pi a^5 \rho_b} - \frac{\nabla\phi}{a}$

Linear Theory Solutions

In the linear limit, the magnetic field remains constant, and density perturbations can be expressed as:

$\delta_b^{\text{lin}} = -\xi_b(k,a)\frac{3M_{\text{Pl}}^2}{\rho_{m0}}S_B(k)$

where $S_B = \nabla \cdot [(\nabla \times \vec{B}) \times \vec{B}]/(4\pi\rho_{b0})$ is the magnetic field source term.

Magnetic Damping Scale

When baryon velocity perturbations become non-linear ($v_b/l \sim aH$), the damping scale can be estimated as:

$\lambda_D \sim \frac{v_A}{aH} \approx 0.1\text{Mpc}\left(\frac{\sqrt{\langle B^2\rangle}}{\text{nG}}\right)$

Numerical Simulation Setup

Initial Conditions

• Simulations begin from the recombination epoch ($z_{\text{in}} = 1090$)
• Magnetic field power spectrum: $P_B(k) = Ak^{n_B}e^{-k^2\lambda_{Di}^2}$
• Focus on spectral indices $n_B = -2$ and $n_B = 2$
• Initial density perturbations set to zero, neglecting inflationary initial conditions

Simulation Parameters

MHD simulations performed using the AREPO code, including:

• Dark matter particles and adaptive mesh gas
• Hybrid tree-PM-direct summation gravitational algorithm
• Godunov method for solving Riemann problems
• IllustrisTNG galaxy formation model

Experimental Setup

Simulation Configuration

Seven different simulations (A-G) were conducted with main parameters including:

• Magnetic field strength: $B_{1\text{Mpc}} = 10^{-3} - 0.8$ nG
• Spectral indices: $n_B = -2.9, -2.0, 2.0$
• Box sizes: 2-20 Mpc/h
• Particle numbers: $128^3 - 256^3$

Evaluation Metrics

• Dimensionless power spectrum: $\Delta(k) = k^3P(k)/(2\pi^2)$
• Evolution of magnetic field, baryon, and total matter power spectra
• Comparison with linear theory predictions

Experimental Results

Main Findings

1. Three Evolutionary Stages

Simulations identified three distinct stages when PMF strength exceeds ~0.05 nG:

1. Turbulent Stage ($z > 100$): PMFs induce turbulent motions, with baryon density perturbations reaching $\mathcal{O}(1)$ values at scales smaller than the turbulent scale
2. Gravity-Dominated Stage ($100 > z > 30$): Gravity overcomes Lorentz force, with perturbations growing at all scales
3. Halo Formation Stage ($z < 30$): Matter collapses to form halos, with magnetic fields strengthened through dynamo mechanisms

2. Non-linear Behavior

• Large-scale Verification: Simulations perfectly verify linear theory predictions at large scales
• Small-scale Saturation: Baryon power spectrum saturates to $\mathcal{O}(1)$ values at $k > \lambda_D^{-1}$ rather than the expected suppression
• Dark Matter Following: Dark matter perturbations follow baryon behavior through gravitational coupling

3. Spectral Index Dependence

• $n_B = -2$: Large-scale behavior $\Delta \propto k^{2n_B+10} = k^6$
• $n_B = 2$: Large-scale behavior $\Delta \propto k^7$, independent of $n_B$ and $B_{1\text{Mpc}}$

Quantitative Results

Magnetic Damping Scale Fitting

Correction coefficients obtained through simulation fitting:

• $\kappa_D(n_B = -2) = 0.8$
• $\kappa_D(n_B = 2) = 0.9$

Semi-analytical Fitting Parameters

For cases with $n_B < -1.5$: $\Delta_{\text{fit}} = \Delta^{\text{lin}}(k,a)[1 + (k\lambda_J)^p]^{-(2n_B+10)/p}$

Fitting parameters ($n_B = -2$):

• Baryons: $\kappa_J^b = 1.57$, $p_b = 1.95$
• Dark matter: $\kappa_J^{\text{DM}} = 0.774$, $p_{\text{DM}} = 1.86$
• Total matter: $\kappa_J^m = 1.17$, $p_m = 1.94$

Comparison with Previous Studies

Figure 8 shows that the semi-analytical fitting from this study compared to literature predictions:

• Results from Shaw & Lewis (2012) and Cruz et al. (2023) overestimate the matter power spectrum by several orders of magnitude near the non-linear transition
• Peak power spectrum differences reach $10^3-10^6$ times

Related Work

PMF Constraint Methods

1. CMB Anisotropies: Most model-independent constraints, with upper limits at nG scales
2. Large-scale Structure: Including Lyman-α forest, dwarf galaxy abundance, stellar mass density, etc.
3. Reionization: PMFs affect early star formation and reionization history
4. 21cm Signal: Promising probe for future observations

Theoretical Development

• Linear Theory: Foundational framework established by Kim, Olinto & Rosner (1996) and others
• Magnetic Jeans Scale: Concept analogous to thermal Jeans scale, predicting perturbation suppression at this scale
• Numerical Simulations: This paper represents the first complete MHD+gravitational simulation study

Conclusions and Discussion

Main Conclusions

1. Linear Theory Verification: Perfect verification of linear theory predictions at large scales
2. Challenge to Magnetic Jeans Picture: Baryon power spectrum saturates rather than being suppressed at small scales
3. Correction of Previous Constraints: Discovery that literature overestimated power spectrum in non-linear regime
4. Extended Parameter Space: First study of Batchelor spectrum effects

Limitations

1. Resolution Constraints: Small-scale turbulence may not be fully resolved
2. Pre-recombination Evolution: Neglects evolution before recombination
3. Parameter Coverage: Only two spectral indices studied
4. Precision Estimates: Results in non-linear regime have precision of $\mathcal{O}(1)$

Future Directions

1. Astrophysical Applications: Application of fitting functions to galaxy formation simulations
2. Observational Constraints: Reassessment of PMF constraints based on structure formation
3. Extended Studies: Including helical magnetic fields, additional spectral indices, etc.
4. High-resolution Simulations: Better resolution of turbulent cascade processes

In-depth Evaluation

Strengths

1. Pioneering Work: First complete MHD+gravitational simulations, filling an important gap
2. Technical Rigor: Advanced AREPO code with complete physical processes considered
3. Important Discoveries: Challenges conventional understanding of magnetic Jeans scale with significant theoretical implications
4. Practical Value: Semi-analytical fitting provides important reference for future research
5. Robust Results: Multiple simulations verify result robustness

Limitations

1. Parameter Restrictions: Only two spectral indices studied, limited coverage
2. Resolution Constraints: Acknowledges potential incomplete resolution of small-scale turbulence
3. Initial Conditions: Neglect of pre-recombination evolution may affect small-scale result accuracy
4. Physical Simplifications: Some subgrid physics processes may affect results

Impact

1. Theoretical Contribution: Corrects understanding of PMF non-linear effects
2. Observational Impact: May necessitate reassessment of existing PMF constraints
3. Methodological Value: Provides exemplar for future MHD cosmological simulations
4. Cross-disciplinary Significance: Contributes to both magnetized plasma physics and cosmology

Application Scenarios

1. PMF Constraints: Reassessment of structure formation-based constraints
2. Galaxy Formation: Provides initial conditions for galaxy formation simulations including PMF effects
3. 21cm Cosmology: Provides theoretical foundation for future 21cm observation predictions
4. Early Universe Physics: Exploration of structure formation in magnetized early universe

References

The paper cites 84 important references covering PMF theory, observational constraints, numerical methods, and related topics. Key references include:

• Kim, Olinto & Rosner (1996): Foundational work on PMF-induced density perturbations
• Shaw & Lewis (2012): Important PMF constraint study
• Jedamzik & Saveliev (2019): Stringent CMB constraints
• Springel (2010): Core AREPO code reference

This paper represents a milestone in primordial magnetic field cosmology, not only verifying linear theory predictions but more importantly discovering new physical behavior in the non-linear regime, establishing an important foundation for future development in this field.