Search

We describe a novel proposal for inflationary magnetogenesis by identifying the non-Abelian sector of Spectator Chromo Natural Inflation (SCNI) with the $\rm{SU(2)}_{\rm L}$ sector of the Standard Model. This mechanism relies on the recently discovered attractor of SCNI in the strong backreaction regime, where the gauge fields do not decay on super-horizon scales and their backreaction leads to a stable new trajectory for the rolling axion field. The large super-horizon gauge fields are partly transformed after the electroweak phase transition into electromagnetic fields. The strength and correlation length of the resulting helical magnetic fields depend on the inflationary Hubble scale and the details of the SCNI sector. For suitable parameter choices we show that the strength of the resulting magnetic fields having correlation lengths around $1\, {\rm {Mpc}}$ are consistent with the required intergalactic magnetic fields for explaining the spectra of high energy $\gamma$ rays from distant blazars., Comment: 29 pages, 8 figures, 1 table

The Faraday rotation effect, quantified by the Rotation Measure (RM), is a powerful probe of the large-scale magnetization of the Universe - tracing magnetic fields not only on galaxy and galaxy cluster scales but also in the intergalactic Medium (IGM; referred to as $\mathrm{RM}_{\text{IGM}}$). The redshift dependence of the latter has extensively been explored with observations. It has also been shown that this relation can help to distinguish between different large-scale magnetization scenarios. We study the evolution of this $\mathrm{RM}_{\text{IGM}}$ for different primordial magnetogenesis scenarios to search for the imprints of primordial magnetic fields (PMFs; magnetic fields originating in the early Universe) on the redshift-dependence of $\mathrm{RM}_{\text{IGM}}$. We use cosmological magnetohydrodynamic (MHD) simulations for evolving PMFs during large-scale structure formation, coupled to the light cone analysis to produce a realistic statistical sample of mock $\mathrm{RM}_{\text{IGM}}$ images. We study the predicted behavior for the cosmic evolution of $\mathrm{RM}_{\text{IGM}}$ for different correlation lengths of PMFs, and provide fitting functions for their dependence on redshifts. We compare these mock RM trends with the recent analysis of the the LOw-Frequency ARray (LOFAR) RM Grid and find that large-scale-correlated PMFs should have (comoving) strengths $\lesssim 0.75$ nanoGauss, if originated during inflation with the scale invariant spectrum and (comoving) correlation length $\sim 19$ cMpc/h or $ \lesssim 30$ nanoGauss if they originated during phase-transition epochs with the comoving correlation length $\sim 1$ cMpc/h. Our findings agree with previous observations and confirm the results of semi-analytical studies, showing that upper limits on the PMF strength decrease as their coherence scales increase., Comment: 16 pages, 9 figures, comments welcome!

At high energies, the dynamics of a plasma with charged fermions can be described in terms of chiral magnetohydrodynamics. Using direct numerical simulations, we demonstrate that chiral magnetic waves (CMWs) can produce a chiral asymmetry $\mu_5 = \mu_\mathrm{L} - \mu_\mathrm{R}$ from a spatially fluctuating (inhomogeneous) chemical potential $\mu = \mu_\mathrm{L} + \mu_\mathrm{R}$, where $\mu_\mathrm{L}$ and $\mu_\mathrm{R}$ are the chemical potentials of left- and right-handed electrically charged fermions, respectively. If the frequency of the CMW is less than or comparable to the characteristic growth rate of the chiral dynamo instability, the magnetic field can be amplified on small spatial scales. The growth rate of this small-scale chiral dynamo instability is determined by the spatial maximum value of $\mu_5$ fluctuations. Therefore, the magnetic field amplification occurs during periods when $\mu_5$ reaches temporal maxima during the CMW. If the small-scale chiral dynamo instability leads to a magnetic field strength that exceeds a critical value, which depends on the resistivity and the initial value of $\mu$, magnetically dominated turbulence is produced. Turbulence gives rise to a large-scale dynamo instability, which we find to be caused by the magnetic alpha effect. Our results have consequences for the dynamics of certain high-energy plasmas, such as the early Universe., Comment: 17 pages, 11 figures, PRD (in press)

We consider the effects of backreaction on axion-SU(2) dynamics during inflation. We use the linear evolution equations for the gauge field modes and compute their backreaction on the background quantities numerically using the Hartree approximation. We show that the spectator chromo-natural inflation attractor is unstable when back-reaction becomes important. Working within the constraints of the linear mode equations, we find a new dynamical attractor solution for the axion field and the vacuum expectation value of the gauge field, where the latter has an opposite sign with respect to the chromo-natural inflation solution. Our findings are of particular interest to the phenomenology of axion-SU(2) inflation, as they demonstrate the instability of the usual trajectory due to large backreaction effects. The viable parameter space of the model becomes significantly altered, provided future non-Abelian lattice simulations confirm the existence of the new dynamical attractor. In addition, the backreaction effects lead to characteristic oscillatory features in the primordial gravitational wave background that are potentially detectable with upcoming gravitational wave detectors., Comment: 25 pages, 13 figures, 2 tables

We study the gravitational wave (GW) spectrum produced by acoustic waves in the early universe, such as would be produced by a first order phase transition, focusing on the low-frequency side of the peak. We confirm with numerical simulations the Sound Shell model prediction of a steep rise with wave number $k$ of $k^9$ to a peak whose magnitude grows at a rate $(H/k_\text{p})H$, where $H$ is the Hubble rate and $k_\text{p}$ the peak wave number, set by the peak wave number of the fluid velocity power spectrum. We also show that hitherto neglected terms give a shallower part with amplitude $(H/k_\text{p})^2$ in the range $H \lesssim k \lesssim k_\text{p}$, which in the limit of small $H/k$ rises as $k$. This linear rise has been seen in other modelling and also in direct numerical simulations. The relative amplitude between the linearly rising part and the peak therefore depends on the peak wave number of the velocity spectrum and the lifetime of the source, which in an expanding background is bounded above by the Hubble time $H^{-1}$. For slow phase transitions, which have the lowest peak wave number and the loudest signals, the acoustic GW peak appears as a localized enhancement of the spectrum, with a rise to the peak less steep than $k^9$. The shape of the peak, absent in vortical turbulence, may help to lift degeneracies in phase transition parameter estimation at future GW observatories., Comment: 21 pages, 10 figures, Accepted for publication in JCAP

Conversion of electromagnetic energy into magnetohydrodynamic energy occurs when the electric conductivity changes from negligible to finite values. This process is relevant during the epoch of reheating of the early Universe at the end of inflation and before the emergence of the radiation-dominated era. We find that conversion into kinetic and thermal energies is primarily the result of electric energy dissipation and that the magnetic energy plays only a secondary role in this process. This means that, since electric energy dominates over magnetic energy during inflation and reheating, significant amounts of electric energy can be converted into magnetohydrodynamic energy when conductivity emerges early enough, before the relevant length scales become stable., Comment: 22 pages, 16 figures, 3 tables, published in Entropy, special issue on Energy Transfer and Dissipation in Plasma Turbulence

In the standard model of particle physics, the chiral anomaly can occur in relativistic plasmas and plays a role in the early Universe, protoneutron stars, heavy-ion collisions, and quantum materials. It gives rise to a magnetic instability if the number densities of left- and right-handed electrically charged fermions are unequal. Using direct numerical simulations, we show this can result just from spatial fluctuations of the chemical potential, causing a chiral dynamo instability, magnetically driven turbulence, and ultimately a large-scale magnetic field through the magnetic alpha effect., Comment: 5 pages, 4 figures

In the primordial plasma, at temperatures above the scale of electroweak symmetry breaking, the presence of chiral asymmetries is expected to induce the development of helical hypermagnetic fields through the phenomenon of chiral plasma instability. It results in magnetohydrodynamic turbulence due to the high conductivity and low viscosity and sources gravitational waves that survive in the universe today as a stochastic polarized gravitational wave background. In this article, we show that this scenario only relies on Standard Model physics, and therefore the observable signatures, namely the relic magnetic field and gravitational background, are linked to a single parameter controlling the initial chiral asymmetry. We estimate the magnetic field and gravitational wave spectra, and validate these estimates with 3D numerical simulations., Comment: 8 pages, 4 figures, 1 table; v2 updated to match published version, results unchanged and PRD accepted; v3: fixed broken reference link

In decaying magnetohydrodynamic (MHD) turbulence with a strong magnetic field, the spectral magnetic energy density is known to increase with time at small wavenumbers $k$, provided the spectrum at low $k$ is sufficiently steep. This process is called inverse cascading and occurs for an initial Batchelor spectrum, where the magnetic energy per linear wavenumber interval increases like $k^4$. For an initial Saffman spectrum that is proportional to $k^2$, however, inverse cascading has not been found in the past. We study here the case of an intermediate $k^3$ spectrum, which may be relevant for magnetogenesis in the early Universe during the electroweak epoch. This case is not well understood in view of the standard Taylor expansion of the magnetic energy spectrum for small $k$. Using high resolution MHD simulations, we show that also in this case there is inverse cascading with a strength just as expected from the conservation of the Hosking integral, which governs the decay of an initial Batchelor spectrum. Even for shallower $k^\alpha$ spectra with spectral index $\alpha>3/2$, our simulations suggest a spectral increase at small $k$ with time $t$ proportional to $t^{4\alpha/9-2/3}$. The critical spectral index of $\alpha=3/2$ is related to the slope of the spectral envelope in the Hosking phenomenology. Our simulations with $2048^3$ mesh points now suggest inverse cascading even for an initial Saffman spectrum., Comment: 18 pages, 9 figures, 3 tables, submitted to J. Plasma Physics, improved evidence for inverse cascading for nonhelical $k^2$ spectrum

We study the evolution of magnetic fields coupled with chiral fermion asymmetry in the framework of chiral magnetohydrodynamics with zero initial total chirality. The initial magnetic field has a turbulent spectrum peaking at a certain characteristic scale and is fully helical with positive helicity. The initial chiral chemical potential is spatially uniform and negative. We consider two opposite cases where the ratio of the length scale of the chiral plasma instability (CPI) to the characteristic scale of the turbulence is smaller and larger than unity. These initial conditions might be realized in cosmological models such as certain types of axion inflation. The magnetic field and chiral chemical potential evolve with inverse cascading in such a way that the magnetic helicity and chirality cancel each other at all times. The CPI time scale is found to determine mainly the time when the magnetic helicity spectrum attains negative values at high wave numbers. The turnover time of the energy-carrying eddies, on the other hand, determines the time when the peak of the spectrum starts to shift to smaller wave numbers via an inverse cascade. The onset of helicity decay is determined by the time when the chiral magnetic effect becomes efficient at the peak of the initial magnetic energy spectrum. When spin flipping is important, the chiral chemical potential vanishes and the magnetic helicity becomes constant, which leads to a faster increase of the correlation length, as expected from magnetic helicity conservation. This also happens when the initial total chirality is imbalanced. Our findings have important implications for baryogenesis after axion inflation., Comment: 24 pages, 21 figures, 6 tables, corrected Fig.17, Phys. Rev. D, in press

In plasmas composed of massless electrically charged fermions, chirality can be interchanged with magnetic helicity while preserving the total chirality through the quantum chiral anomaly. The decay of turbulent energy in plasmas such as those in the early Universe and compact stars is usually controlled by certain conservation laws. In the case of zero total chirality, when the magnetic helicity density balances with the appropriately scaled chiral chemical potential to zero, the total chirality no longer determines the decay. We propose that in such a case, an adaptation to the Hosking integral, which is conserved in nonhelical magnetically dominated turbulence, controls the decay in turbulence with helicity balanced by chiral fermions. We show, using a high resolution numerical simulation, that this is indeed the case. The magnetic energy density decays and the correlation length increases with time just like in nonhelical turbulence with vanishing chiral chemical potential. But here, the magnetic helicity density is nearly maximum and shows a scaling with time $t$ proportional to $t^{-2/3}$. This is unrelated to the $t^{-2/3}$ decay of magnetic {\it energy} in fully helical magnetic turbulence. The modulus of the chiral chemical potential decays in the same fashion. This is much slower than the exponential decay previously expected in theories of asymmetric baryon production from the hypermagnetic helicity decay after axion inflation., Comment: 6 pages, 4 figures, corrected Fig.2

We study the propagation of cosmological gravitational wave (GW) backgrounds from the early radiation era until the present day in modified theories of gravity. Comparing to general relativity (GR), we study the effects that modified gravity parameters, such as the GW friction $\alpha_{\rm M}$ and the tensor speed excess $\alpha_{\rm T}$, have on the present-day GW spectrum. We use both the WKB estimate, which provides an analytical description but fails at superhorizon scales, and numerical simulations that allow us to go beyond the WKB approximation. We show that a constant $\alpha_{\rm T}$ makes relatively insignificant changes to the GR solution, especially taking into account the constraints on its value from GW observations by the LIGO--Virgo collaboration, while $\alpha_{\rm M}$ can introduce modifications to the spectral slopes of the GW energy spectrum in the low-frequency regime depending on the considered time evolution of $\alpha_{\rm M}$. The latter effect is additional to the damping or growth occurring equally at all scales that can be predicted by the WKB approximation. In light of the recent observations by pulsar timing array (PTA) collaborations, and the potential observations by future detectors such as SKA, LISA, DECIGO, BBO, or ET, we show that, in most of the cases, constraints cannot be placed on the effects of $\alpha_{\rm M}$ and the initial GW energy density $\mathcal{E}_{\rm GW}^*$ separately, but only on the combined effects of the two, unless the signal is observed at different frequency ranges. In particular, we provide some constraints on the combined effects from the reported PTA observations., Comment: 36 pages, 11 figures, 2 tables

Primordial magnetic fields (PMFs) are possible candidates for explaining the observed magnetic fields in galaxy clusters. Two competing scenarios of primordial magnetogenesis have been discussed in the literature: inflationary and phase-transitional. We study the amplification of both large- and small-scale correlated magnetic fields, corresponding to inflation- and phase transition-generated PMFs, in a massive galaxy cluster. We employ high-resolution magnetohydrodynamic cosmological zoom-in simulations to resolve the turbulent motions in the intracluster medium. We find that the turbulent amplification is more efficient for the large-scale inflationary models, while the phase transition-generated seed fields show moderate growth. The differences between the models are imprinted on the spectral characteristics of the field (such as the amplitude and the shape of the magnetic power spectrum) and therefore on the final correlation length. We find a one order of magnitude difference between the final strengths of the inflation- and phase transition-generated magnetic fields, and a factor of 1.5 difference between their final coherence scales. Thus, the final configuration of the magnetic field retains information about the PMF generation scenarios. Our findings have implications for future extragalactic Faraday rotation surveys with the possibility of distinguishing between different magnetogenesis scenarios., Comment: 21 pages, 13 figures, comments welcome

The magnetic fields in galaxy clusters and probably also in the interstellar medium are believed to be generated by a small-scale dynamo. Theoretically, during its kinematic stage, it is characterized by a Kazantsev spectrum, which peaks at the resistive scale. It is only slightly shallower than the Saffman spectrum that is expected for random and causally connected magnetic fields. Causally disconnected fields have the even steeper Batchelor spectrum. Here we show that all three spectra are present in the small-scale dynamo. During the kinematic stage, the Batchelor spectrum occurs on scales larger than the energy-carrying scale of the turbulence, and the Kazantsev spectrum on smaller scales within the inertial range of the turbulence -- even for a magnetic Prandtl number of unity. In the saturated state, the dynamo develops a Saffman spectrum on large scales, suggestive of the build-up of long-range correlations. At large magnetic Prandtl numbers, elongated structures are seen in synthetic synchrotron emission maps showing the parity-even E polarization. We also observe a significant excess in the E polarization over the parity-odd B polarization at subresistive scales, and a deficiency at larger scales. This finding is at odds with the observed excess in the Galactic microwave foreground emission, which is believed to be associated with larger scales. The E and B polarizations may be highly non-Gaussian and skewed in the kinematic regime of the dynamo. For dust emission, however, the polarized emission is always nearly Gaussian, and the excess in the E polarization is much weaker., Comment: 14 pages, 25 figures, 5 tables, MNRAS, in press

The Saffman helicity invariant of Hosking and Schekochihin (2021, PRX 11, 041005), which we here call the Hosking integral, has emerged as an important quantity that may govern the decay properties of magnetically dominated nonhelical turbulence. Using a range of different computational methods, we confirm that this quantity is indeed gauge-invariant and nearly perfectly conserved in the limit of large Lundquist numbers. For direct numerical simulations with ordinary viscosity and magnetic diffusivity operators, we find that the solution develops in a nearly self-similar fashion. In a diagram quantifying the instantaneous decay coefficients of magnetic energy and integral scale, we find that the solution evolves along a line that is indeed suggestive of the governing role of the Hosking integral. The solution settles near a line in this diagram that is expected for a self-similar evolution of the magnetic energy spectrum. The solution will settle in a slightly different position when the magnetic diffusivity decreases with time, which would be compatible with the decay being governed by the reconnection time scale rather than the Alfv\'en time., Comment: 19 pages, 7 figures, 2 tables; accepted by JPP

Hydrodynamic and magnetohydrodynamic turbulence in the early Universe can drive gravitational waves (GWs) and imprint their spectrum onto that of GWs, which might still be observable today. We study the production of the GW background from freely decaying magnetohydrodynamic turbulence from helical and nonhelical initial magnetic fields. To understand the produced GW spectra, we develop a simple model on the basis of the evolution of the magnetic stress tensor. We find that the GW spectra obtained in this model reproduce those obtained in numerical simulations if we consider the detailed time evolution of the low-frequency tail of the stress spectrum from numerical simulations. We also show that the shapes of the produced GW frequency spectra are different for helical and nonhelical cases for the same initial magnetic energy spectra. Such differences can help distinguish helical and nonhelical initial magnetic fields from a polarized background of GWs -- especially when the expected circular polarization cannot be detected directly., Comment: 12 pages, 12 figures, 3 tables, Accepted for publication in Phys. Rev. D

We revisit the big bang nucleosynthesis (BBN) limits on primordial magnetic fields and/or turbulent motions accounting for the decaying nature of turbulent sources between the time of generation and BBN. This leads to larger estimates for the gravitational wave (GW) signal than previously expected. We address the detection prospects through space-based interferometers (for GWs generated around the electroweak energy scale) as well as pulsar timing arrays and astrometric missions (for GWs generated around the quantum chromodynamics energy scale)., Comment: 9 pages; 6 figures; discussions and references added; typos fixed; conclusions unchanged

We study the leading-order nonlinear gravitational waves (GWs) produced by an electromagnetic (EM) stress in reheating magnetogenesis scenarios. Both nonhelical and helical magnetic fields are considered. By numerically solving the linear and leading-order nonlinear GW equations, we find that the GW energy from the latter is usually larger. We compare their differences in terms of the GW spectrum and parameterize the GW energy difference due to the nonlinear term, $\Delta\mathcal{E}_{\rm GW}$, in terms of EM energy $\mathcal{E}_{\rm EM}$ as $\Delta\mathcal{E}_{\rm GW}=(\tilde p\mathcal{E}_{\rm EM}/k_*)^3$, where $k_*$ is the characteristic wave number, $\tilde p=0.84$ and $0.88$ are found in the nonhelical and helical cases, respectively, with reheating around the QCD energy scale, while $\tilde p=0.45$ is found at the electroweak energy scale. We also compare the polarization spectrum of the linear and nonlinear cases and find that adding the nonlinear term usually yields a decrease in the polarization that is proportional to the EM energy density. We parameterize the fractional polarization suppression as $|\Delta \mathcal{P}_{\rm GW}/\mathcal{P}_{\rm GW}|=\tilde r \mathcal{E}_{\rm EM}/k_*$ and find $\tilde r = 1.2 \times 10^{-1}$, $7.2 \times 10^{-4}$, and $3.2 \times 10^{-2}$ for the helical cases with reheating temperatures $T_{\rm r} = 300 {\rm TeV}$, $8 {\rm GeV}$, and $120 {\rm MeV}$, respectively. Prospects of observation by pulsar timing arrays, space-based interferometers, and other novel detection proposals are also discussed., Comment: 17 pages, 7 figures, 3 tables, to be submitted to PRD

Primordial magnetic fields could explain the large-scale magnetic fields present in the Universe. Inflation and phase transitions in the early Universe could give rise to such fields with unique characteristics. We investigate the magneto-hydrodynamic evolution of these magnetogenesis scenarios with cosmological simulations. We evolve inflation-generated magnetic fields either as (i) uniform (homogeneous) or as (ii) scale-invariant stochastic fields, and phase transition-generated ones either as (iii) helical or as (iv) non-helical fields from the radiation-dominated epoch. We find that the final distribution of magnetic fields in the simulated cosmic web shows a dependence on the initial strength and the topology of the seed field. Thus, the observed field configuration retains information on the initial conditions at the moment of the field generation. If detected, primordial magnetic field observations would open a new window for indirect probes of the early universe. The differences between the competing models are revealed on the scale of galaxy clusters, bridges, as well as filaments and voids. The distinctive spectral evolution of different seed fields produces imprints on the correlation length today. We discuss how the differences between rotation measures from highly ionized regions can potentially be probed with forthcoming surveys., Comment: 28 pages, 17 figures, accepted for publication in ApJ

We study the dynamics of magnetic fields in chiral magnetohydrodynamics, which takes into account the effects of an additional electric current related to the chiral magnetic effect in high energy plasmas. We perform direct numerical simulations, considering weak seed magnetic fields and inhomogeneities of the chiral chemical potential mu_5 with a zero mean. We demonstrate that a small-scale chiral dynamo can occur in such plasmas if fluctuations of mu_5 are correlated on length scales that are much larger than the scale on which the dynamo growth rate reaches its maximum. Magnetic fluctuations grow by many orders of magnitude due to the small-scale chiral dynamo instability. Once the nonlinear backreaction of the generated magnetic field on fluctuations of mu_5 sets in, the ratio of these scales decreases and the dynamo saturates. When magnetic fluctuations grow sufficiently to drive turbulence via the Lorentz force before reaching maximum field strength, an additional mean-field dynamo phase is identified. The mean magnetic field grows on a scale that is larger than the integral scale of turbulence after the amplification of the fluctuating component saturates. The growth rate of the mean magnetic field is caused by a magnetic alpha effect that is proportional to the current helicity. With the onset of turbulence, the power spectrum of mu_5 develops a universal k^(-1) scaling independently of its initial shape, while the magnetic energy spectrum approaches a k^(-3) scaling., Comment: 22 pages, 22 figures

In relativistic magnetized plasmas, asymmetry in the number densities of left- and right-handed fermions, i.e., a non-zero chiral chemical potential mu_5, leads to an electric current along the magnetic field. This causes a chiral dynamo instability for a uniform mu_5, but our simulations reveal a dynamo even for fluctuating mu_5 with zero mean. It produces magnetically-dominated turbulence and generates mean magnetic fields via the magnetic alpha effect. Eventually, a universal scale-invariant k^(-1) spectrum of mu_5 and a k^(-3) magnetic spectrum are formed independently of the initial condition., Comment: 6 pages, 4 figures, published in Phys. Rev. Lett

Using numerical simulations of helical inflationary magnetogenesis in a low reheating temperature scenario, we show that the magnetic energy spectrum is strongly peaked at a particular wavenumber that depends on the reheating temperature. Gravitational waves (GWs) are produced at frequencies between 3 nHz and 50 mHz for reheating temperatures between 150 MeV and 3x10^5 GeV, respectively. At and below the peak frequency, the stress spectrum is always found to be that of white noise. This implies a linear increase of GW energy per logarithmic wavenumber interval, instead of a cubic one, as previously thought. Both in the helical and nonhelical cases, the GW spectrum is followed by a sharp drop for frequencies above the respective peak frequency. In this magnetogenesis scenario, the presence of a helical term extends the peak of the GW spectrum and therefore also the position of the aforementioned drop toward larger frequencies compared to the case without helicity. This might make a difference in it being detectable with space interferometers. The efficiency of GW production is found to be almost the same as in the nonhelical case, and independent of the reheating temperature, provided the electromagnetic energy at the end of reheating is fixed to be a certain fraction of the radiation energy density. Also, contrary to the case without helicity, the electric energy is now less than the magnetic energy during reheating. The fractional circular polarization is found to be nearly hundred per cent in a certain range below the peak frequency range., Comment: 16 pages, 8 figures, 4 tables, submitted to ApJ

We use direct numerical simulations of decaying primordial hydromagnetic turbulence with helicity to compute the resulting gravitational wave (GW) production and its degree of circular polarization. We find a clear dependence of the polarization of the resulting GWs on the fractional helicity of the turbulent source and we show that driven magnetic fields produce GWs more efficiently than magnetic fields that are initially present, leading to larger spectral amplitudes. The helicity does not have a huge impact on the maximum spectral amplitude in any of the two types of turbulence considered. However, the GW spectrum at wave numbers away from the peak becomes smaller for larger values of the magnetic fractional helicity. The degree of circular polarization approaches zero at frequencies below the peak, and reaches its maximum at the peak. At higher frequencies, it stays finite if the magnetic field is initially present, and it approaches zero if it is driven. We predict that the spectral peak of the GW signal can be detected by LISA if the turbulent energy density is at least $\sim 3\%$ of the radiation energy density, and the characteristic scale is a hundredth of the horizon at the electroweak scale. We show that the resulting GW polarization is unlikely to be detectable by the anisotropies induced by our proper motion in the dipole response function of LISA. Such signals can, however, be detectable by cross-correlating data from the LISA-Taiji network for turbulent energy densities of $\sim 5\%$, and fractional helicity of 0.5 to 1. Second-generation space-based GW detectors, such as BBO and DECIGO, would allow for the detection of a larger range of the GW spectrum and smaller amplitudes of the magnetic field., Comment: 49 pages, 18 figures, 1 table, accepted in JCAP

We present three-dimensional direct numerical simulations of the production of magnetic fields and gravitational waves (GWs) in the early Universe during a low energy scale matter-dominated post-inflationary reheating era, and during the early subsequent radiative era, which is strongly turbulent. The parameters of the model are determined such that it avoids a number of known physical problems and produces magnetic energy densities between 0.2% and 2% of the critical energy density at the end of reheating. During the subsequent development of a turbulent magnetohydrodynamic cascade, magnetic fields and GWs develop a spectrum that extends to higher frequencies in the millihertz (nanohertz) range for models with reheating temperatures of around 100 GeV (150 MeV) at the beginning of the radiation-dominated era. However, even though the turbulent cascade is fully developed, the GW spectrum shows a sharp drop for frequencies above the peak value. This suggests that the turbulence is less efficient in driving GWs than previously thought. The peaks of the resulting GW spectra may well be in the range accessible to space interferometers, pulsar timing arrays, and other facilities., Comment: 20 pages, 15 figures, 4 tables, Astrophys. J. 920, 26 (2021)

We consider a generic dispersive massive gravity theory and numerically study its resulting modified energy and strain spectra of tensor gravitational waves (GWs) sourced by (i) fully developed turbulence during the electroweak phase transition (EWPT) and (ii) forced hydromagnetic turbulence during the QCD phase transition (QCDPT). The GW spectra are then computed in both spatial and temporal Fourier domains. We find, from the spatial spectra, that the slope modifications are weakly dependent on the eddy size at QCDPT, and, from the temporal spectra, that the modifications are pronounced in the $1$--$10{\rm nHz}$ range -- the sensitivity range of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) -- for a graviton mass $m_{\rm g}$ in the range $2\times10^{-23}{\rm eV}\lesssim m_{\rm g}c^2\lesssim7\times10^{-22}{\rm eV}$., Comment: 15 pages, 3 figures, 2 tables, comments welcome

We study the gravitational wave (GW) signal sourced by primordial turbulence that is assumed to be present at cosmological phase transitions like the electroweak and quantum chromodynamics phase transitions. We consider various models of primordial turbulence, such as those with and without helicity, purely hydrodynamical turbulence induced by fluid motions, and magnetohydrodynamic turbulence whose energy can be dominated either by kinetic or magnetic energy, depending on the nature of the turbulence. We also study circularly polarized GWs generated by parity violating sources such as helical turbulence. Our ultimate goal is to determine the efficiency of GW production through different classes of turbulence. We find that the GW energy and strain tend to be large for acoustic or irrotational turbulence, even though its tensor mode amplitude is relatively small at most wave numbers. Only at very small wave numbers is the spectral tensor mode significant, which might explain the efficient GW production in that case., Comment: 25 pages, 8 figures

We perform numerical simulations of gravitational waves (GWs) induced by hydrodynamic and hydromagnetic turbulent sources that might have been present at cosmological quantum chromodynamic (QCD) phase transitions. For turbulent energies of about 4% of the radiation energy density, the typical scale of such motions may have been a sizable fraction of the Hubble scale at that time. The resulting GWs are found to have an energy fraction of about $10^{-9}$ of the critical energy density in the nHz range today and may already have been observed by the NANOGrav collaboration. This is further made possible by our findings of shallower spectra proportional to the square root of the frequency for nonhelical hydromagnetic turbulence. This implies more power at low frequencies than for the steeper spectra previously anticipated. The behavior toward higher frequencies depends strongly on the nature of the turbulence. For vortical hydrodynamic and hydromagnetic turbulence, there is a sharp drop of spectral GW energy by up to five orders of magnitude in the presence of helicity, and somewhat less in the absence of helicity. For acoustic hydrodynamic turbulence, the sharp drop is replaced by a power law decay, albeit with a rather steep slope. Our study supports earlier findings of a quadratic scaling of the GW energy with the magnetic energy of the turbulence and inverse quadratic scaling with the peak frequency, which leads to larger GW energies under QCD conditions., Comment: 14 pages, 8 figures, 3 tables; page numbers, figure captions, references updated and added; discussions expanded; results unchanged and PRD accepted

Relic gravitational waves (GWs) can be produced by primordial magnetic fields. However, not much is known about the resulting GW amplitudes and their dependence on the details of the generation mechanism. Here we treat magnetic field generation through the chiral magnetic effect (CME) as a generic mechanism and explore its dependence on the speed of generation (the product of magnetic diffusivity and characteristic wavenumber) and the speed characterizing the maximum magnetic field strength expected from the CME. When the latter exceeds the former (regime I), the regime applicable to the early universe, we obtain an inverse cascade with moderate GW energy that scales with the third power of the magnetic energy. When the generation speed exceeds the CME limit (regime II), the GW energy continues to increase without a corresponding increase of magnetic energy. In the early kinematic phase, the GW energy spectrum (per linear wavenumber interval) has opposite slopes in both regimes and is characterized by an inertial range spectrum in regime I and a white noise spectrum in regime II. The occurrence of these two slopes is shown to be a generic consequence of a nearly monochromatic exponential growth of the magnetic field. The resulting GW energy is found to be proportional to the fifth power of the limiting CME speed and the first power of the generation speed., Comment: 15 pages, 15 figures, 4 tables, ApJ

We perform direct numerical simulations to compute the net circular polarization of gravitational waves from helical (chiral) turbulent sources in the early Universe for a variety of initial conditions, including driven (stationary) and decaying turbulence. We investigate the resulting gravitational wave signal assuming different turbulent geneses such as magnetically or kinetically driven cases. Under realistic physical conditions in the early Universe we compute numerically the wave number-dependent polarization degree of the gravitational waves. We find that the spectral polarization degree strongly depends on the initial conditions. The peak of the spectral polarization degree occurs at twice the typical wavenumber of the source, as expected, and for fully helical decaying turbulence, it reaches its maximum of nearly 100\% {\it only} at the peak. We determine the temporal evolution of the turbulent sources as well as the resulting gravitational waves, showing that the dominant contribution to their spectral energy density happens shortly after the activation of the source. Only through an artificially prolonged decay of the turbulence can further increase of the gravitational wave amplitude be achieved. We estimate the detection prospects for the net polarization, arguing that its detection contains {\it clean} information (including the generation mechanisms, time, and strength) about the sources of possible parity violations in the early Universe., Comment: 9 pages, 5 figures; discussions added; references updated and added; conclusions unchanged; Phys. Rev. Research in press

Motivated by a scenario of magnetogenesis in which a homogeneous magnetic field is generated during inflation, we study the magnetohydrodynamic evolution of the primordial plasma motions for two kinds of initial conditions -- (i) a spatially homogeneous field with an unlimited correlation length, and (ii) a zero flux scale-invariant statistically homogeneous magnetic field. In both cases, we apply, for a short initial time interval, monochromatic forcing at a certain wave number so that the correlation length is finite, but much smaller than the typical length scale of turbulence. In particular, we investigate the decay of nonhelical and helical hydromagnetic turbulence. We show that, in the presence of a homogeneous magnetic field, the decay of helical and nonhelical small-scale fields can occur rapidly. This is a special property of a system with a perfectly homogeneous magnetic field, which is sometimes considered as a local approximation to a slowly varying background field. It can never change and acts as an imposed magnetic field. This is in a sharp contrast to the case of a statistically homogeneous magnetic field, where we recover familiar decay properties: a much slower decay of magnetic energy and a faster growth of the correlation length, especially in the case with magnetic felicity. The result suggests that a homogeneous magnetic field, if generated during inflation, should persist under the influence of small-scale fields and could be the origin of the large-scale magnetic field in the Universe., Comment: 16 pages, 7 figures, 1 table. Journal version

The large-scale magnetic field of the Milky Way is thought to be created by an alpha-Omega dynamo, which implies that it should have opposite handedness North and South of the Galactic midplane. Here we attempt to detect a variation in handedness using polarization data from the Wilkinson Microwave Anisotropy Probe. Previous analyzes of the parity-even and parity-odd parts of linear polarization of the global dust and synchrotron emission have focused on quadratic correlations in spectral space of, and between, these two components. Here, by contrast, we analyze the parity-odd polarization itself and show that it has, on average, opposite signs in Northern and Southern Galactic hemispheres. Comparison with a Galactic mean-field dynamo model shows broad qualitative agreement and reveals that the sign of the observed hemispheric dependence of the azimuthally averaged parity-odd polarization is not determined by the sign of alpha, but by the sense of differential rotation., Comment: 7 pages, 4 figures, 1 table, submitted to ApJL

For velocity and magnetic fields, the turbulent pressure and, more generally, the squared fields such as the components of the turbulent stress tensor, play important roles in astrophysics. For both one and three dimensions, we derive the equations relating the energy spectra of the fields to the spectra of their squares. We solve the resulting integrals numerically and show that for turbulent energy spectra of Kolmogorov type, the spectral slope of the stress spectrum is also of Kolmogorov type. For shallower turbulence spectra, the slope of the stress spectrum quickly approaches that of white noise, regardless of how blue the spectrum of the field is. For fully helical fields, the stress spectrum is elevated by about a factor of two in the subinertial range, while that in the inertial range remains unchanged. We discuss possible implications for understanding the spectrum of primordial gravitational waves from causally generated magnetic fields during cosmological phase transitions in the early universe. We also discuss potential diagnostic applications to the interstellar medium, where polarization and scintillation measurements characterize the square of the magnetic field., Comment: 7 pages, 9 figures, ApJ, in press

We perform direct numerical simulations of magnetohydrodynamic turbulence in the early universe and numerically compute the resulting stochastic background of gravitational waves and relic magnetic fields. These simulations do not make the simplifying assumptions of earlier analytic work. If the turbulence is assumed to have an energy-carrying scale that is about a hundredth of the Hubble radius at the time of generation, as expected in a first-order phase transition, the peak of gravitational wave power will be in the mHz frequency range for a signal produced at the electroweak scale. The efficiency of gravitational wave (GW) production varies significantly with how the turbulence is driven. Detectability of turbulence at the electroweak scale by the planned Laser Interferometer Space Antenna (LISA) requires anywhere from 0.1% to 10% of the thermal plasma energy density to be in plasma motions or magnetic fields, depending on the model of the driving process. Our results predict a new universal form below the spectral peak frequency that is shallower than previously thought. This implies larger values of the GW energy spectra in the low-frequency range. This extends the range where turbulence is detectable with LISA to lower frequencies, corresponding to higher energy scales than the assumed energy-carrying scale., Comment: 11 pages, 1 appendix, 7 + 1 (appendix) figures, 2 tables, published in PRD, updated figures and corrected normalization error

Blazar observations point toward the possible presence of magnetic fields over intergalactic scales of the order of up to $\sim1\,$Mpc, with strengths of at least $\sim10^{-16}\,$G. Understanding the origin of these large-scale magnetic fields is a challenge for modern astrophysics. Here we discuss the cosmological scenario, focussing on the following questions: (i) How and when was this magnetic field generated? (ii) How does it evolve during the expansion of the universe? (iii) Are the amplitude and statistical properties of this field such that they can explain the strengths and correlation lengths of observed magnetic fields? We also discuss the possibility of observing primordial turbulence through direct detection of stochastic gravitational waves in the mHz range accessible to LISA., Comment: Invited talk given at FM8 "New Insights in Extragalactic Magnetic Fields", XXX IAU GA, Vienna, Austria, August 29, 2018; 4 pages; 2 figures; accepted for publication

The chiral magnetic effect (CME) is a quantum relativistic effect that describes the appearance of an additional electric current along a magnetic field. It is caused by an asymmetry between the number densities of left- and right-handed fermions, which can be maintained at high energies when the chirality flipping rate can be neglected, for example in the early Universe. The inclusion of the CME in the Maxwell equations leads to a modified set of MHD equations. We discuss how the CME is implemented in the PENCIL CODE. The CME plays a key role in the evolution of magnetic fields since it results in a dynamo effect associated with an additional term in the induction equation. This term is formally similar to the $\alpha$ effect in classical mean-field MHD. However, the chiral dynamo can operate without turbulence and is associated with small spatial scales that can be, in the case of the early Universe, orders of magnitude below the Hubble radius. A chiral $\alpha_\mu$ effect has also been identified in mean-field theory. It occurs in the presence of turbulence but is not related to kinetic helicity. Depending on the plasma parameters, chiral dynamo instabilities can amplify magnetic fields over many orders of magnitude. These instabilities can affect the propagation of MHD waves, which is demonstrated by our DNS. We also study the coupling between the evolution of the chiral chemical potential and the ordinary chemical potential, which is proportional to the sum of the number densities of left- and right-handed fermions. An important consequence of this coupling is the emergence of chiral magnetic waves (CMWs). We confirm numerically that linear CMWs and MHD waves are not interacting. Our simulations suggest that the chemical potential has only a minor effect on the non-linear evolution of the chiral dynamo., Comment: 22 pages, 10 figures, in press at the GAFD special issue "Physics and Algorithms of the Pencil Code"

Gradient- and curl-type or E- and B-type polarizations have been routinely analyzed to study the physics contributing to the cosmic microwave background polarization and galactic foregrounds. They characterize the parity-even and parity-odd properties of the underlying physical mechanisms, for example hydromagnetic turbulence in the case of dust polarization. Here we study spectral correlation functions characterizing the parity-even and parity-odd parts of linear polarization for homogeneous and inhomogeneous turbulence to show that only the inhomogeneous helical case can give rise to a parity-odd polarization signal. We also study nonhelical turbulence and suggest that a strong nonvanishing (here negative) skewness of the E polarization is responsible for an enhanced ratio of the EE to the BB (quadratic) correlation in both helical and nonhelical cases. This could explain the enhanced EE/BB ratio observed recently for dust polarization. We close with a preliminary assessment of using linear polarization of the Sun to characterize its helical turbulence without being subjected to the pi ambiguity that magnetic inversion techniques have to address., Comment: 11 pages, 13 figures, 4 tables, ApJ, in press

We investigate the statistical properties of isotropic, stochastic, Gaussian distributed, helical magnetic fields characterized by different shapes of the energy spectra at large length scales and study the associated realizability condition. We discuss smoothed magnetic fields that are commonly used when the primordial magnetic field is constrained by observational data. We are particularly interested in scale-invariant magnetic fields that can be generated during the inflationary stage by quantum fluctuations. We determine the correlation length of such magnetic fields and relate it to the infrared cutoff of perturbations produced during inflation. We show that this scale determines the observational signatures of the inflationary magnetic fields on the cosmic microwave background. At smaller scales, the scale-invariant spectrum changes with time. It becomes a steeper weak-turbulence spectrum at progressively larger scales. We show numerically that the critical length scale where this happens is the turbulent-diffusive scale, which increases with the square root of time., Comment: 24 pages, 9 figures. Discussion expanded, references added, results unchanged. Accepted for publication in JCAP

An asymmetry in the number density of left- and right-handed fermions is known to give rise to a new term in the induction equation that can result in a dynamo instability. At high temperatures, when a chiral asymmetry can survive for long enough, this chiral dynamo instability can amplify magnetic fields efficiently, which in turn drive turbulence via the Lorentz force. While it has been demonstrated in numerical simulations that this chiral magnetically driven turbulence exists and strongly affects the dynamics of the magnetic field, the details of this process remain unclear. The goal of this paper is to analyse the energetics of chiral magnetically driven turbulence and its effect on the generation and dynamics of the magnetic field using direct numerical simulations. We study these effects for different initial conditions, including a variation of the initial chiral chemical potential and the magnetic Prandtl number, Pm. In particular, we determine the ratio of kinetic to magnetic energy, Upsilon, in chiral magnetically driven turbulence. Within the parameter space explored in this study, Upsilon reaches a value of approximately 0.064-0.074 -independently of the initial chiral asymmetry and for Pm=1. Our simulations suggest, that Upsilon decreases as a power law when increasing Pm by decreasing the viscosity. While the exact scaling depends on the details of the fitting criteria and the Reynolds number regime, an approximate result of Upsilon(Pm)=0.1 Pm^{-0.4} is reported. Using the findings from our numerical simulations, we analyse the energetics of chiral magnetically driven turbulence in the early Universe., Comment: 22 pages, 11 figures, in press at the GAFD special issue "Recent Developments in Natural Dynamos"

Using direct numerical simulations (DNS), we study laminar and turbulent dynamos in chiral magnetohydrodynamics (MHD) with an extended set of equations that accounts for an additional contribution to the electric current due to the chiral magnetic effect (CME). This quantum phenomenon originates from an asymmetry between left- and right-handed relativistic fermions in the presence of a magnetic field and gives rise to a chiral dynamo. We show that the magnetic field evolution proceeds in three stages: (1) a small-scale chiral dynamo instability; (2) production of chiral magnetically driven turbulence and excitation of a large-scale dynamo instability due to a new chiral effect (alpha_mu effect); and (3) saturation of magnetic helicity and magnetic field growth controlled by a conservation law for the total chirality. The $\alpha_\mu$ effect becomes dominant at large fluid and magnetic Reynolds numbers and is not related to kinetic helicity. The growth rate of the large-scale magnetic field and its characteristic scale measured in the numerical simulations agree well with theoretical predictions based on mean-field theory. The previously discussed two-stage chiral magnetic scenario did not include stage (2) during which the characteristic scale of magnetic field variations can increase by many orders of magnitude. Based on the findings from numerical simulations, the relevance of the CME and the chiral effects revealed in the relativistic plasma of the early universe and of proto-neutron stars are discussed., Comment: 24 pages, 21 figures

We present new simulations of decaying hydromagnetic turbulence for a relativistic equation of state relevant to the early universe. We compare helical and nonhelical cases either with kinetically or magnetically dominated initial fields. Both kinetic and magnetic initial helicities lead to maximally helical magnetic fields after some time, but with different temporal decay laws. Both are relevant to the early universe, although no mechanisms have yet been identified that produce magnetic helicity with strengths comparable to the big bang nucleosynthesis limit at scales comparable to the Hubble horizon at the electroweak phase transition. Nonhelical magnetically dominated fields could still produce picoGauss magnetic fields under most optimistic conditions. Only helical magnetic fields can potentially have nanoGauss strengths at scales up to 30 kpc today., Comment: 17 pages, 13 figures; submitted to PRD

We show that in decaying hydromagnetic turbulence with initial kinetic helicity, a weak magnetic field eventually becomes fully helical. The sign of magnetic helicity is opposite to that of the kinetic helicity - regardless of whether or not the initial magnetic field was helical. The magnetic field undergoes inverse cascading with the magnetic energy decaying approximately like t^{-1/2}. This is even slower than in the fully helical case, where it decays like t^{-2/3}. In this parameter range, the product of magnetic energy and correlation length raised to a certain power slightly larger than unity, is approximately constant. This scaling of magnetic energy persists over long time scales. At very late times and for domain sizes large enough to accommodate the growing spatial scales, we expect a cross-over to the t^{-2/3} decay law that is commonly observed for fully helical magnetic fields. Regardless of the presence or absence of initial kinetic helicity, the magnetic field experiences exponential growth during the first few turnover times, which is suggestive of small-scale dynamo action. Our results have applications to a wide range of experimental dynamos and astrophysical time-dependent plasmas, including primordial turbulence in the early universe., Comment: 15 pages, 14 figures, submitted to Phys. Rev. Fluids

The presence of asymmetry between fermions of opposite handedness in plasmas of relativistic particles can lead to exponential growth of a helical magnetic field via a small-scale chiral dynamo instability known as the chiral magnetic effect. Here, we show, using dimensional arguments and numerical simulations, that this process produces through the Lorentz force chiral magnetically driven turbulence. A k^{-2} magnetic energy spectrum emerges via inverse transfer over a certain range of wavenumbers k. The total chirality (magnetic helicity plus normalized chiral chemical potential) is conserved in this system. Therefore, as the helical magnetic field grows, most of the total chirality gets transferred into magnetic helicity until the chiral magnetic effect terminates. Quantitative results for height, slope, and extent of the spectrum are obtained. Consequences of this effect for cosmic magnetic fields are discussed., Comment: 7 pages, 4 figures, 1 table, published in ApJL

We consider a scale-invariant helical magnetic field generated during inflation. We show that, if the mean magnetic helicity density of such a field is measured, it can be used to determine a lower bound on the duration of inflation. Upper bounds can be used to derive constraints on the minimal duration of inflation if one assumes that the magnetic field generated during inflation is helical. Using three-dimensional simulations, we show that an initially scale-invariant field develops, which is similar both with and without magnetic helicity. In the fully helical case, however, the magnetic field appears to have a more pronounced folded structure., Comment: 19 pages, 7 figures, JCAP accepted; results unchanged, discussions added, references added and updated, typos fixed

We perform numerical simulations of decaying hydrodynamic and magnetohydrodynamic turbulence. We classify our time-dependent solutions by their evolutionary tracks in parametric plots between instantaneous scaling exponents. We find distinct classes of solutions evolving along specific trajectories toward points on a line of self-similar solutions. These trajectories are determined by the underlying physics governing individual cases, while the infrared slope of the initial conditions plays only a limited role. In the helical case, even for a scale-invariant initial spectrum (inversely proportional to wavenumber k), the solution evolves along the same trajectory as for a Batchelor spectrum (proportional to k^4., Comment: 5 pages, 4 figures, with 3 pages supplemental material, published in PRL

We study the evolution of primordial magnetic fields in an expanding cosmic plasma. For this purpose we present a comprehensive theoretical model to consider the evolution of MHD turbulence that can be used over a wide range of physical conditions, including cosmological and astrophysical applications. We model different types of decaying cosmic MHD turbulence in the expanding universe and characterize the large-scale magnetic fields in such a medium. Direct numerical simulations of freely decaying MHD turbulence are performed for different magnetogenesis scenarios: magnetic fields generated during cosmic inflation as well as electroweak and QCD phase transitions in the early universe. Magnetic fields and fluid motions are strongly coupled due to the high Reynolds number in the early universe. Hence, we abandon the simple adiabatic dilution model to estimate magnetic field amplitudes in the expanding universe and include turbulent mixing effects on the large-scale magnetic field evolution. Numerical simulations have been carried out for non-helical and helical magnetic field configurations. The numerical results show the possibility of inverse transfer of energy in magnetically dominated non-helical MHD turbulence. On the other hand, decay properties of helical turbulence depend on whether the turbulent magnetic field is in a weakly or a fully helical state. Our results show that primordial magnetic fields can be considered as a seed for the observed large-scale magnetic fields in galaxies and clusters. Bounds on the magnetic field strength are obtained and are consistent with the upper and lower limits set by observations of extragalactic magnetic fields., Comment: 20 pages, 6 figures; Invited talk given at "Turbulent Mixing and Beyond" Workshop, 4-9 August 2014, ICTP, Italy; v. 3 matches version published in Physica Scripta

In the presence of magnetic helicity, inverse transfer from small to large scales is well known in magnetohydrodynamic (MHD) turbulence and has applications in astrophysics, cosmology, and fusion plasmas. Using high resolution direct numerical simulations of magnetically dominated self-similarly decaying MHD turbulence, we report a similar inverse transfer even in the absence of magnetic helicity. We compute for the first time spectral energy transfer rates to show that this inverse transfer is about half as strong as with helicity, but in both cases the magnetic gain at large scales results from velocity at similar scales interacting with smaller-scale magnetic fields. This suggests that both inverse transfers are a consequence of a universal mechanisms for magnetically dominated turbulence. Possible explanations include inverse cascading of the mean squared vector potential associated with local near two-dimensionality and the shallower $k^2$ subinertial range spectrum of kinetic energy forcing the magnetic field with a $k^4$ subinertial range to attain larger-scale coherence. The inertial range shows a clear $k^{-2}$ spectrum and is the first example of fully isotropic magnetically dominated MHD turbulence exhibiting weak turbulence scaling., Comment: 4 pages, 4 figures, submitted to PRL, plus 9 pages of supplemental material with 12 more figures

We review our current understanding of galactic dynamo theory, paying particular attention to numerical simulations both of the mean-field equations and the original three-dimensional equations relevant to describing the magnetic field evolution for a turbulent flow. We emphasize the theoretical difficulties in explaining non-axisymmetric magnetic fields in galaxies and discuss the observational basis for such results in terms of rotation measure analysis. Next, we discuss nonlinear theory, the role of magnetic helicity conservation and magnetic helicity fluxes. This leads to the possibility that galactic magnetic fields may be bi-helical, with opposite signs of helicity and large and small length scales. We discuss their observational signatures and close by discussing the possibilities of explaining the origin of primordial magnetic fields., Comment: 28 pages, 15 figure, to appear in Lecture Notes in Physics "Magnetic fields in diffuse media", Eds. E. de Gouveia Dal Pino and A. Lazarian

Using one-dimensional models, we show that a helical magnetic field with an appropriate sign of helicity can compensate the Faraday depolarization resulting from the superposition of Faraday-rotated polarization planes from a spatially extended source. For radio emission from a helical magnetic field, the polarization as a function of the square of the wavelength becomes asymmetric with respect to zero. Mathematically speaking, the resulting emission occurs then either at observable or at unobservable (imaginary) wavelengths. We demonstrate that rotation measure (RM) synthesis allows for the reconstruction of the underlying Faraday dispersion function in the former case, but not in the latter. The presence of positive magnetic helicity can thus be detected by observing positive RM in highly polarized regions in the sky and negative RM in weakly polarized regions. Conversely, negative magnetic helicity can be detected by observing negative RM in highly polarized regions and positive RM in weakly polarized regions. The simultaneous presence of two magnetic constituents with opposite signs of helicity is shown to possess signatures that can be quantified through polarization peaks at specific wavelengths and the gradient of the phase of the Faraday dispersion function. Similar polarization peaks can tentatively also be identified for the bi-helical magnetic fields that are generated self-consistently by a dynamo from helically forced turbulence, even though the magnetic energy spectrum is then continuous. Finally, we discuss the possibility of detecting magnetic fields with helical and non-helical properties in external galaxies using the Square Kilometre Array., Comment: 12 pages, 12 figures, ApJ, in press (with 3-D turbulence results now included)