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 calculations. If the turbulence is assumed to have an energy-carrying scale that is about a hundredth of the Hubble scale 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 GW production varies significantly with how the turbulence is driven. Detectability of turbulence at the electroweak scale by the planned Laser Interferometer Space Antenna requires anywhere from 0.1 to 10 percent of the thermal plasma energy density to be in plasma motions or magnetic fields, depending on the model of the driving process. At frequencies below the peak value, our results predict a new universal form of the spectrum with more power than previously thought. The enhanced low-frequency tail makes turbulence at significantly higher energy scales detectable.

Hydromagnetic turbulence produced during phase transitions in the early universe can be a powerful source of stochastic gravitational waves (GWs). GWs can be modelled by the linearised spatial part of the Einstein equations sourced by the Reynolds and Maxwell stresses. We have implemented two different GW solvers into the {\sc Pencil Code}---a code which uses a third order timestep and sixth order finite differences. Using direct numerical integration of the GW equations, we study the appearance of a numerical degradation of the GW amplitude at the highest wavenumbers, which depends on the length of the timestep---even when the Courant--Friedrichs--Lewy condition is ten times below the stability limit. This degradation leads to a numerical error, which is found to scale with the third power of the timestep. A similar degradation is not seen in the magnetic and velocity fields. To mitigate numerical degradation effects, we alternatively use the exact solution of the GW equations under the assumption that the source is constant between subsequent timesteps. This allows us to use a much longer timestep, which cuts the computational cost by a factor of about ten.

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.

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.

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.

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.