The tensor modes reentering the Hubble radius when the plasma is dominated by a stiff fluid lead to a spectral energy density whose blue slope depends on the total post-inflationary sound speed. This result gets however corrected by a secondary (gauge-dependent) term coming from the curvature inhomogeneities that reenter all along the same stage of expansion. In comparison with the first-order result, the secondary contribution is shown to be always suppressed inside the sound horizon and its effect on the total spectral energy density of the relic gravitons is therefore negligible for all phenomenological purposes. It is also suggested that the effective anisotropic stress of the curvature inhomogeneities can be obtained from the functional derivative of the second-order action of curvature inhomogeneities with respect to the background metric.

In the same way as the energy density associated with the tensor modes of the geometry modifies the evolution of the curvature perturbations, the scalar modes may also indirectly affect the cosmic backgrounds of relic gravitons by inducing higher-order corrections that are only superficially gauge-invariant. This spurious gauge-invariance gets manifest when the effective anisotropic stresses, computed in different coordinate systems, are preliminarily expressed in a form that only depends on the curvature inhomogeneities defined on comoving orthogonal hypersurfaces and on their corresponding time derivatives. Using this observation we demonstrate in general terms that the higher-order contributions derived in diverse coordinate systems coincide when the wavelengths are smaller than the sound horizon defining the evolution of the curvature inhomogeneities but they lead to sharply different results in the opposite limit. A similar drawback arises when the energy density of the relic gravitons is derived from competing energy-momentum pseudo-tensors and should be consistently taken into account in the related phenomenological discussions.

The diffuse backgrounds of relic gravitons with frequencies ranging between the aHz band and the GHz region encode the ultimate information on the primeval evolution of the plasma and on the underlying theory of gravity well before the electroweak epoch. While the temperature and polarization anisotropies of the microwave background radiation probe the low-frequency tail of the graviton spectra, during the next score year the pulsar timing arrays and the wide-band interferometers (both terrestrial and hopefully space-borne) will explore a much larger frequency window encompassing the nHz domain and the audio band. The salient theoretical aspects of the relic gravitons are reviewed in a cross-disciplinary perspective touching upon various unsettled questions of particle physics, cosmology and astrophysics.

When the different wavelengths of the scalar and tensor modes of the geometry are all assigned on the same space-like hypersurface the maximally amplified frequencies of the spectrum remain smaller than the Planck mass only if the duration of a stage of accelerated expansion and the corresponding tensor to scalar ratio are severely constrained. All the different wavelengths can be initialized on the same space-like hypersurface at the onset of inflation but this strategy and the related conclusions are plausible only for classical inhomogeneities. We argue that a whole class of potential constraints is easily evaded provided the different wavelengths of the quantum fields are assigned as soon as they cross the corresponding Planckian hypersurfaces. In this case the Cauchy data for the mode functions depend on the wavenumber so that larger wavelengths start evolving earlier while shorter wavelengths are assigned later. Within this strategy the duration of a conventional inflationary phase and the corresponding tensor to scalar ratio are not constrained but the large-scale power spectra inherit specific large-scale corrections that remain however unobservable.

The energy density, the pressure and the anisotropic stress of the relic gravitons do not have a unique gauge-invariant and frame-invariant expression since the equivalence principle ultimately forbids the localization of the energy-momentum of the gravitational field. All the strategies proposed so far suggest compatible descriptions inside the effective horizon but they lead to sharply different answers when the wavelengths of the gravitons exceed the Hubble radius by remaining shorter than the typical extension of causally connected patches. We present a coherent discussion of the energy-momentum pseudo-tensors of the relic gravitons in the different kinematical regimes with the aim of scrutinizing the mutual consistency of the competing suggestions. The various proposals are systematically compared by deducing and analyzing the explicit form of the observables in realistic physical situations. General lessons are drawn on the most plausible parametrization of the energy-momentum pseudo-tensor of the relic gravitons.

The coherence of the relic gravitons is investigated within a quantum mechanical perspective. After introducing the notion and the properties of the generalized Glauber correlators valid in the tensor case, the degrees of first- and second-order coherence are evaluated both inside and beyond the effective horizon. The inclusive approach (encompassing the polarizations of the gravitons) is contrasted with the exclusive approximation where the total intensity is calculated either from a single polarization or even from a single mode of the field. While the relic gravitons reentering the effective horizon after the end of a quasi-de Sitter stage of expansion are first-order coherent, the Hanbury Brown-Twiss correlations always exhibit a super-Poissonian statistics with different quantitative features that depend on the properties of their initial states and on the average over the tensor polarizations.

If the inflaton and the quintessence fields are identified, the background geometry evolves through a stiff epoch undershooting the expansion rate of a radiation-dominated plasma. For some classes of inflationary potentials this scenario is at odds with the current observational evidence since the corresponding tensor-to-scalar ratio is too large. Quintessential inflation is analyzed when the gravitational action is supplemented by a contribution quadratic in the Einstein-Hilbert term. In the Palatini formulation the addition such a term does not affect the scalar modes during the inflationary phase and throughout the course of the subsequent stiff epoch but it suppresses the tensor power spectrum and the tensor-to-scalar ratio. While in the Palatini formulation the power-law potentials leading to a quintessential inflationary dynamics are again viable, the high-frequency spike of the relic graviton spectrum is squeezed and the whole signal is suppressed at least when the higher-order contributions appearing in the action are explicitly decoupled from the inflaton.

The degree of second-order coherence of the relic gravitons produced from the vacuum is super-Poissonian and larger than in the case of a chaotic source characterized by a Bose-Einstein distribution. If the initial state does not minimize the tensor Hamiltonian and has a dispersion smaller than its averaged multiplicity, the overall statistics is by definition sub-Poissonian. Depending on the nature of the sub-Poissonian initial state, the final degree of second-order coherence of the quanta produced by stimulated emission may diminish (possibly even below the characteristic value of a chaotic source) but it always remains larger than one (i.e. super-Poissonian). When the initial statistics is Poissonian (like in the case of a coherent state or for a mixed state weighted by a Poisson distribution) the degree of second-order coherence of the produced gravitons is still super-Poissonian. Even though the quantum origin of the relic gravitons inside the Hubble radius can be effectively disambiguated by looking at the corresponding Hanbury Brown-Twiss correlations, the final distributions caused by different initial states maintain their super-Poissonian character which cannot be altered.

The polarizations of the tensor modes of the geometry evolving in cosmological backgrounds are treated as the components of a bispinor whose dynamics follows from an appropriate gauge-invariant action. This novel framework bears a close analogy with the (optical) Jones calculus and leads to a compact classification of the various interactions able to polarize the relic gravitons.

We investigate the impact of the post-inflationary thermal histories on the cosmic graviton spectrum caused by the inflationary variation of their refractive index. Depending on the frequency band, the spectral energy distribution can be mildly red, blue or even violet. Wide portions of the parameter space lead to potentially relevant signals both in the audio range (probed by the advanced generation of terrestrial interferometers) and in the mHz band (where space-borne detectors could be operational within the incoming score year). The description of the refractive index in conformally related frames is clarified.

If the refractive index of the tensor modes increases during a conventional inflationary stage of expansion the relic graviton spectrum is tilted towards high frequencies. Two apparently diverse parametrizations of this effect are shown to be related by a rescaling of the four-dimensional metric through a conformal factor that involves the refractive index itself. Non-monotonic spectra with a maximum in the MHz region correspond to a limited variation of the refractive index terminating well before the end of inflation. After exploring a general approach encompassing the ones proposed so far, we estimate the required sensitivity for the direct detection of the predicted gravitational radiation and demonstrate that the allowed regions of the parameter space are within reach for some of the planned detectors operating either in the audio band (like Ligo/Virgo and Kagra) or in the mHz band (like Lisa, Bbo and Decigo).

We show that the spectral energy distribution of relic gravitons mildly increases for frequencies smaller than the $\mu$Hz and then flattens out whenever the refractive index of the tensor modes is dynamical during a quasi-de Sitter stage of expansion. For a conventional thermal history the high-frequency plateau ranges between the mHz and the audio band but it is supplemented by a spike in the GHz region if a stiff post-inflationary phase precedes the standard radiation-dominated epoch. Even though the slope is blue at intermediate frequencies, it may become violet in the MHz window. For a variety of post-inflationary histories, including the conventional one, a dynamical index of refraction leads to a potentially detectable spectral energy density in the kHz and in the mHz regions while all the relevant phenomenological constraints are concurrently satisfied.

Prior to photon decoupling magnetic random fields of comoving intensity in the nano-Gauss range distort the temperature and the polarization anisotropies of the microwave background, potentially induce a peculiar B-mode power spectrum and may even generate a frequency-dependent circularly polarized V-mode. We critically analyze the theoretical foundations and the recent achievements of an interesting trialogue involving plasma physics, general relativity and astrophysics.

During a de Sitter stage of expansion the spectator fields of different spin are constrained by the critical density bound and by further requirements determined by their specific physical nature. The evolution of spectator electric fields in conformally flat background geometries is occasionally concocted by postulating the existence of ad hoc currents but this apparently innocuous trick violates the second law of thermodynamics. Such a problem occurs, in particular, for those configurations (customarily employed for the analysis of the Schwinger effect in four-dimensional de Sitter backgrounds) leading to an electric energy density which is practically unaffected by the expansion of the underlying geometry. The obtained results are compared with more mundane situations where Joule heating develops in the early stages of a quasi-de Sitter phase.

The dynamics of bouncing universes is characterized by violating certain coordinate invariant restrictions on the total energy-momentum tensor, customarily referred to as energy conditions. Although there could be epochs where the null energy condition is locally violated, it may perhaps be enforced in an averaged sense. Explicit examples of this possibility are investigated in different frameworks.

The degrees of quantum coherence of cosmological perturbations of different spins are computed in the large-scale limit and compared with the standard results holding for a single mode of the electromagnetic field in an optical cavity. The degree second-order coherence of curvature inhomogeneities (and, more generally, of the scalar modes of the geometry) reproduces faithfully the optical limit. For the vector and tensor fluctuations the numerical values of the normalized degrees of second-order coherence in the zero-time delay limit are always larger than unity (which is the Poisson benchmark value) but differ from the corresponding expressions obtainable in the framework of the single-mode approximation. General lessons are drawn on the quantum coherence of large-scale cosmological fluctuations.

Owing to the analogy with the ordinary photons in the visible range of the electromagnetic spectrum, the Glauber theory is generalized to address the quantum coherence of the gauge field fluctuations parametrically amplified during an inflationary stage of expansion. The first and second degrees of quantum coherence of relic photons are then computed beyond the effective horizon defined by the evolution of the susceptibility. In the zero-delay limit the Hanbury Brown-Twiss correlations exhibit a super-Poissonian statistics which is however different from the conventional results of the single-mode approximation customarily employed, in quantum optics, to classify the coherence properties of visible light. While in the case of large-scale curvature perturbations the degrees of quantum coherence coincide with the naive expectation of the single-mode approximation, the net degree of second-order coherence computed for the relic photons diminishes thanks to the effect of the polarizations. We suggest that the Hanbury Brown-twiss correlations are probably the only tool to assess the quantum or classical origin of the large-scale magnetic fluctuations and of the corresponding curvature perturbations.

The quantum mechanical generation of hypermagnetic and hyperlectric fields in four-dimensional conformally flat background geometries rests on the simultaneous continuity of the effective horizon and of the extrinsic curvature across the inflationary boundary. The junction conditions for the gauge fields are derived in general terms and corroborated by explicit examples with particular attention to the limit of a sudden (but nonetheless continuous) transition of the effective horizon. After reducing the dynamics to a pair of integral equations related by duality transformations, we compute the power spectra and deduce a novel class of logarithmic corrections which turn out to be, however, numerically insignificant and overwhelmed by the conductivity effects once the gauge modes reenter the effective horizon. In this perspective the magnetogenesis requirements and the role of the postinflationary conductivity are clarified and reappraised. As long as the total duration of the inflationary phase is nearly minimal, quasi-flat hypermagnetic power spectra are comparatively more common than in the case of vacuum initial data.

The maximally gyrotropic configurations of the hypermagnetic field at the electroweak epoch can induce a stochastic background of relic gravitational waves with comoving frequencies ranging from the $\mu$Hz to the kHz. Using two complementary approaches we construct a physical template family for the emission of the gravitational radiation produced by the hypermagnetic knots. The current constraints and the presumed sensitivities of the advanced wide-band interferometers (both terrestrial and space-borne) are combined to infer that the lack of observations at intermediate frequencies may invalidate the premise of baryogenesis models based (directly or indirectly) on the presence of gyrotropic configurations of the hypermagnetic field at the electroweak epoch. Over the intermediate frequency range the spectral energy density of the gravitational waves emitted by the hypermagnetic knots at the electroweak scale can exceed the inflationary signal even by nine orders of magnitude without affecting the standard bounds applicable on the stochastic backgrounds of gravitational radiation. The signal of hypermagnetic knots can be disambiguated, at least in principle, since the the produced gravitational waves are polarized.

Even if the Higgs field does not affect the evolution of the background geometry, its massive inhomogeneities induce large-scale gauge fields whose energy density depends on the slow-roll parameters, on the effective scalar mass and, last but not least, on the dimensionless coupling to the space-time curvature. Since the non-Abelian gauge modes are screened, the non-minimal coupling to gravity predominantly affects the evolution of the hypercharge and electromagnetic fields. While in the case of minimal coupling the obtained constraints are immaterial, as soon as the coupling increases beyond one fourth the produced fields become overcritical. We chart the whole parameter space of this qualitatively new set of bounds. Whenever the limits on the curvature coupling are enforced, the magnetic field may still be partially relevant for large-scale magnetogenesis and exceed $10^{-20}$ G for the benchmark scale of the protogalactic collapse.

Bouncing solutions are obtained from a generally covariant action characterized by a potential which is a nonlocal functional of the dilaton field at two separated space-time points. Gradient instabilities are shown to arise in this context but they are argued to be nongeneric. After performing a gauge-invariant and frame-invariant derivation of the evolution equations of the fluctuations, a heuristic criterion for the avoidance of pathological instabilities is proposed and corroborated by a number of explicit examples that turn out to be compatible with a quasi-flat spectrum of curvature inhomogeneities for typical wavelengths larger than the Hubble radius.

The tensor to scalar ratio is affected by the evolution of the large-scale gauge fields potentially amplified during an inflationary stage of expansion. After deriving the exact evolution equations for the scalar and tensor modes of the geometry in the presence of dynamical gauge fields, it is shown that the tensor to scalar ratio is bounded from below by the dominance of the adiabatic contribution and it cannot be smaller than one thousands whenever the magnetogenesis is driven by a single inflaton field.

The curvature inhomogeneities are systematically scrutinized in the framework of the Glauber approach. The amplified quantum fluctuations of the scalar and tensor modes of the geometry are shown to be first-order coherent while the interference of the corresponding intensities is larger than in the case of Bose-Einstein correlations. After showing that the degree of second-order coherence does not suffice to characterize unambiguously the curvature inhomogeneities, we argue that direct analyses of the degrees of third and fourth-order coherence are necessary to discriminate between different correlated states and to infer more reliably the statistical properties of the large-scale fluctuations. We speculate that the moments of the multiplicity distributions of the relic phonons might be observationally accessible thanks to new generations of instruments able to count the single photons of the Cosmic Microwave Background in the THz region.

The evolution equations of anomalous magnetohydrodynamics are derived in the extreme relativistic regime and contrasted with the treatment of hydromagnetic nonlinearities pioneered by Lichnerowicz in the absence of anomalous currents. In particular we explore the situation where the conventional vector currents are complemented by the axial-vector currents arising either from the pseudo Nambu-Goldstone bosons of a spontaneously broken symmetry or because of finite fermionic density effects. After expanding the generally covariant equations in inverse powers of the conductivity, the relativistic analog of the magnetic diffusivity equation is derived in the presence of vortical and magnetic currents. While the anomalous contributions are generally suppressed by the diffusivity, they are shown to disappear in the perfectly conducting limit. When the flow is irrotational, boost-invariant and with vanishing four-acceleration the corresponding evolution equations are explicitly integrated so that the various physical regimes can be directly verified.

Stochastic backgrounds of relic gravitons of cosmological origin extend from frequencies of the order of the aHz up to the GHz range. Since the temperature and polarization anisotropies constrain the low frequency normalization of the spectra, in the concordance paradigm the strain amplitude corresponding to the frequency window of wide-band interferometers turns out to be, approximately, nine orders of magnitude smaller than the astounding signal recently reported and attributed to a binary black hole merger. The backgrounds of relic gravitons expected from the early Universe are compared with the stochastic foregrounds stemming from the estimated multiplicity of the astrophysical sources. It is suggested that while the astrophysical foregrounds are likely to dominate between few Hz and 10 kHz, relic gravitons with frequencies exceeding 100 kHz represent a potentially uncontaminated signal for the next generation of high-frequency detectors currently under scrutiny.

The viscous inhomogeneities of a relativistic plasma determine a further class of entropic modes whose amplitude must be sufficiently small since curvature perturbations are observed to be predominantly adiabatic and Gaussian over large scales. When the viscous coefficients only depend on the energy density of the fluid the corresponding curvature fluctuations are shown to be almost adiabatic. After addressing the problem in a gauge-invariant perturbative expansion, the same analysis is repeated at a non-perturbative level by investigating the nonlinear curvature inhomogeneities induced by the spatial variation of the viscous coefficients. It is demonstrated that the quasiadiabatic modes are suppressed in comparison with a bona fide adiabatic solution. Because of its anomalously large tensor to scalar ratio the quasiadiabatic mode cannot be a substitute for the conventional adiabatic paradigm so that, ultimately, the present findings seems to exclude the possibility of a successful accelerated dynamics solely based on relativistic viscous fluids. If the dominant adiabatic mode is not affected by the viscosity of the background a sufficiently small fraction of entropic fluctuations of viscous origin cannot be a priori ruled out.

Provided the quantum fluctuations are amplified in the presence of a classical gauge field configuration the resulting curvature perturbations exhibit a mild statistical anisotropy which should be sufficiently weak not to conflict with current observational data. The curvature power spectra induced by weakly anisotropic initial states are computed here for the first time when the electric and the magnetic gauge couplings evolve at different rates as it happens, for instance, in the relativistic theory of van der Waals interactions. After recovering the results valid for coincident gauge couplings, the constraints imposed by the isotropy and the homogeneity of the initial states are discussed. The obtained bounds turn out to be more stringent than naively expected and cannot be ignored when discussing the underlying magnetogenesis scenarios.

We investigate the production of the hypermagnetic gyrotropy when the electric and magnetic gauge couplings evolve at different rates, as it happens in the the relativistic theory of the Van der Waals forces. If a pseudo-scalar interaction breaks the duality symmetry of the corresponding equations, the gyrotropic configurations of the hypermagnetic fields can be amplified from the vacuum during an inflationary stage of expansion. After charting the parameter space of the model in terms of the rates of evolution of the magnetic and electric gauge couplings, we identify the regions where the gyrotropy is sufficiently intense to seed the baryon asymmetry of the Universe at the electroweak epoch while the backreaction constraints, the strong coupling bounds and the other astrophysical limits are concurrently satisfied.

The equations of anomalous magnetohydrodynamics describe an Abelian plasma where conduction and chiral currents are simultaneously present and constrained by the second law of thermodynamics. At high frequencies the magnetic currents play the leading role and the spectrum is dominated by two-fluid effects. The system behaves instead as a single fluid in the low-frequency regime where the vortical currents induce potentially large hypermagnetic fields. After deriving the physical solutions of the generalized Appleton-Hartree equation, the corresponding dispersion relations are scrutinized and compared with the results valid for cold plasmas. Hypermagnetic knots and fluid vortices can be concurrently present at very low frequencies and suggest a qualitatively different dynamics of the hydromagnetic nonlinearities.

When the gauge fields have derivative couplings to scalars, like in the case of the relativistic theory of Van der Waals (or Casimir-Polder) interactions, conformal invariance is broken but the magnetic and electric susceptibilities are not bound to coincide. We analyze the formation of large-scale magnetic fields in slow-roll inflation and find that they are generated at the level of a few hundredths of a nG and over typical length scales between few Mpc and $100$ Mpc. Using a new time parametrization that reduces to conformal time but only for coincident susceptibilities, the gauge action is quantized while the evolution equations of the corresponding mode functions are more easily solvable. The power spectra depend on the normalized rates of variation of the two susceptibilities (or of the corresponding gauge couplings) and on the absolute value of their ratio at the beginning of inflation. We pin down explicit regions in the parameter space where all the physical requirements (i.e. the backreaction constraints, the magnetogenesis bounds and the naturalness of the initial conditions of the scenario) are jointly satisfied. Weakly coupled initial data are favoured if the gauge couplings are of the same order at the end of inflation. Duality is systematically used to simplify the analysis of the wide parameter space of the model.

The non-perturbative curvature inhomogeneities induced by relativistic viscous fluids are not conserved in the large-scale limit. However when the bulk viscosity is a function of the total energy density of the plasma (or of the trace of the extrinsic curvature) the relevant evolution equations develop a further symmetry preventing the non-linear growth of curvature perturbations. In this situation the fully inhomogeneous evolution can be solved to leading order in the gradient expansion. Over large-scales both the acceleration and the curvature inhomogeneities are determined by the bulk viscosity coefficients. Conversely the shear viscosity does not affect the evolution of the curvature and does not produce any acceleration. The curvature modes analyzed here do not depend on the choice of time hypersurfaces and are invariant for infinitesimal coordinate transformations in the perturbative regime.

The dynamical evolution of the refractive index of the tensor modes of the geometry produces a specific class of power spectra characterized by a blue (i.e. slightly increasing) slope which is directly determined by the competition of the slow-roll parameter and of the rate of variation of the refractive index. Throughout the conventional stages of the inflationary and post-inflationary evolution, the microwave background anisotropies measurements, the pulsar timing limits and the big-bang nucleosythesis constraints set stringent bounds on the refractive index and on its rate of variation. Within the physically allowed region of the parameter space the cosmic background of relic gravitons leads to a potentially large signal for the ground based detectors (in their advanced version) and for the proposed space-borne interferometers. Conversely, the lack of direct detection of the signal will set a qualitatively new bound on the dynamical variation of the refractive index.

In a recent Letter [Phys. Rev. Lett. 114 091301 (2105)] the cause of the acceleration of the present Universe has been identified with the shear viscosity of an imperfect relativistic fluid even in the absence of any bulk viscous contribution. The gist of this comment is that the shear viscosity, if anything, can only lead to an accelerated expansion over sufficiently small scales well inside the Hubble radius.

When the predecoupling plasma is thermodynamically reversible its fluctuations are classified in terms of the adiabatic and entropic modes. A different category of physical solutions, so far unexplored, arises when the inhomogeneities of the viscosity coefficients induce computable curvature perturbations. The viscous modes are explicitly illustrated and compared with the conventional isocurvature solutions.

Cosmological singularities are often discussed by means of a gradient expansion that can also describe, during a quasi-de Sitter phase, the progressive suppression of curvature inhomogeneities. While the inflationary event horizon is being formed the two mentioned regimes coexist and a uniform expansion can be conceived and applied to the evolution of spatial gradients across the protoinflationary boundary. It is argued that conventional arguments addressing the preinflationary initial conditions are necessary but generally not sufficient to guarantee a homogeneous onset of the conventional inflationary stage.

In conformally flat background geometries the long wavelength gravitons can be described in the fluid approximation and they induce scalar fluctuations both during inflation and in the subsequent radiation-dominated epoch. While this effect is minute and suppressed for a de Sitter stage of expansion, the fluctuations of the energy-momentum pseudo-tensor of the graviton fluid lead to curvature perturbations that increase with time all along the post-inflationary evolution. An explicit calculation of these effects is presented for a standard thermal history and it is shown that the growth of the curvature perturbations caused by the long wavelength modes is approximately compensated by the slope of the power spectra of the energy density, pressure and anisotropic stress of the relic gravitons.

Provided the consistency relations are not violated, the recent Bicep2 observations pin down the absolute normalization, the spectral slope and the maximal frequency of the cosmic graviton background produced during inflation. The properly normalized spectra are hereby computed from the lowest frequencies (of the order of the present Hubble rate) up to the highest frequency range in the GHz region. Deviations from the conventional paradigm cannot be excluded and are examined by allowing for different physical possibilities including, in particular, a running of the tensor spectral index, an explicit breaking of the consistency relations and a spike in the high-frequency tail of the spectrum coming either from a post-inflationary phase dominated by a stiff fluid of from the contribution of waterfall fields in a hybrid inflationary context. The direct determinations of the tensor to scalar ratio at low frequencies, if confirmed by the forthcoming observations, will also affect and constrain the high-frequencies uncertainties. The limits on the cosmic graviton backgrounds coming from wide-band interferometers (such as Ligo/Virgo, Lisa and Bbo/Decigo) together with a more accurate scrutiny of the tensor B mode polarization at low frequencies will set direct bounds on the post-inflationary evolution and on other unconventional completions of the standard lore.

The quasinormal mode of a gravitating and magnetized fluid in a spatially flat, isotropic and homogeneous cosmological background is derived in the presence of the fluid sources of anisotropic stress and of the entropic fluctuations of the plasma. The obtained gauge-invariant description involves a system of two coupled differential equations whose physical content is analyzed in all the most relevant situations. The Cauchy problem of large-scale curvature perturbations during the radiation dominated stage of expansion can be neatly formulated and its general solution is shown to depend on five initial data assigned when the relevant physical wavelengths are larger than the particle horizon. The consequences of this approach are explored.

This paper investigates the Faraday effect as a different source of B mode polarization. The E mode polarization is Faraday rotated provided a stochastic large-scale magnetic field is present prior to photon decoupling. In the first part of the paper we discuss the case where the tensor modes of the geometry are absent and we argue that the B mode recently detected by the Bicep2 collaboration cannot be explained by a large-scale magnetic field rotating, through the Faraday effect, the well established E mode polarization. In this case, the observed temperature autocorrelations would be excessively distorted by the magnetic field. In the second part of the paper the formation of Faraday rotation is treated as a stationary, random and Markovian process with the aim of generalizing a set of scaling laws originally derived in the absence of the tensor modes of the geometry. We show that the scalar, vector and tensor modes of the brightness perturbations can all be Faraday rotated even if the vector and tensor parts of the effect have been neglected, so far, by focussing the attention on the scalar aspects of the problem. The mixing between the power spectra of the E mode and B mode polarizations involves a unitary transformation depending nonlinearly on the Faraday rotation rate. The present approach is suitable for a general scrutiny of the polarization observables and of their frequency dependence.

If we posit the validity of the consistency relations, the tensor spectral index and the relative amplitude of the scalar and tensor power spectra are both fixed by a single slow roll parameter. The physics of the protoinflationary transition can break explicitly the consistency relations causing a reduction of the inflationary curvature scale in comparison with the conventional lore. After a critical scrutiny, we argue that the inflationary curvature scale, the total number of inflationary efolds and, ultimately, the excursion of the inflaton across its Planckian boundary are all characterized by a computable theoretical error. While these considerations ease some of the tensions between the Bicep2 data and the other satellite observations, they also demand an improved understanding of the protoinflationary transition whose physical features may be assessed, in the future, through a complete analysis of the spectral properties of the B mode autocorrelations.

As repeatedly speculated in the past, the linear polarization of the Cosmic Microwave Background can be rotated via the Faraday effect. An economic explanation of the recent Bicep2 observations, not relying on long-wavelength tensor modes of the geometry, would stipulate that the detected B mode comes exclusively from a Faraday rotated E mode polarization. We show hereunder that this interpretation is ruled out by the existing upper limits on the B mode polarization obtained by independent experiments at observational frequencies much lower than the operating frequency of the Bicep2 experiment. We then derive the fraction of the observed B mode polarization ascribable to the Faraday effect and suggest a dedicated experimental strategy for its detection.

The formation of the microwave background polarization anisotropies is investigated when the stochastic Faraday rate is stationary, random and Markovian. The scaling properties of the polarization power spectra and of their nonlinear combinations are scrutinized as a function of the comoving frequency. It is argued that each frequency channel of a given experiment measuring simultaneously the E-mode and the B-mode spectra can be analyzed in this framework with the aim of testing the physical origin of the polarization in a model-independent perspective.

Anisotropic inflationary background geometries are analyzed in the context of an extended gauge action where the electric and magnetic susceptibilities are not bound to coincide and depend on the inflaton field. After deriving various classes of solutions with electric and magnetic hairs, we discuss the problem of the initial boundary conditions of the shear parameter and consider a globally neutral plasma as a possible relic of a preinflationary stage of expansion. While electric hairs are washed out by the finite value of the protoinflationary conductivity, magnetic hairs can persist and introduce a tiny amount of shear causing a different inflationary rate of expansion along orthogonal spatial directions. The plasma interactions are a necessary criterion to discriminate between physical and unphysical initial conditions but they are not strictly sufficient to warrant the stability of a given magnetic solution.

We investigate what can be said about the interaction of scalar fields with Abelian gauge fields during a quasi-de Sitter phase of expansion and under the assumption that the electric and the magnetic susceptibilities do not coincide. The duality symmetry, transforming the magnetic susceptibility into the inverse of the electric susceptibility, exchanges the magnetic and electric power spectra. The mismatch between the two susceptibilities determines an effective refractive index affecting the evolution of the canonical fields. The constraints imposed by the duration of the inflationary phase and by the magnetogenesis requirements pin down the rate of variation of the susceptibilities that is consistent with the observations of the magnetic field strength over astrophysical and cosmological scales but avoids back-reaction problems. The parameter space of this magnetogenesis scenario is wider than in the case when the susceptibilities are equal, as it happens when the inflaton or some other spectator field is solely coupled to the standard gauge kinetic term.

Anomalous symmetries induce currents which can be parallel rather than orthogonal to the hypermagnetic field. Building on the analogy with charged liquids at high magnetic Reynolds numbers, the persistence of anomalous currents is scrutinized for parametrically large conductivities when the plasma approximation is accurate. Different examples in globally neutral systems suggest that the magnetic configurations minimizing the energy density with the constraint that the helicity be conserved coincide, in the perfectly conducting limit, with the ones obtainable in ideal magnetohydrodynamics where the anomalous currents are neglected. It is argued that this is the rationale for the ability of extending to anomalous magnetohydrodynamics the hydromagnetic solutions characterized by finite gyrotropy. The generally covariant aspects of the problem are addressed with particular attention to conformally flat geometries which are potentially relevant for the description of the electroweak plasma prior to the phase transition.

We explore what can be said on the effective temperature and sound speed of a statistical ensemble of fluid phonons present at the onset of a conventional inflationary phase. The phonons are the actual normal modes of the gravitating and irrotational fluid that dominates the protoinflationary dynamics. The bounds on the tensor to scalar ratio result in a class of novel constraints involving the slow roll parameter, the sound speed of the phonons and the temperature of the plasma prior to the onset of inflation. If the current size of the Hubble radius coincides with the inflationary event horizon redshifted down to the present epoch, the sound speed of the phonons can be assessed from independent measurements of the tensor to scalar ratio and of the tensor spectral index.

The temperature and polarization anisotropies of the Cosmic Microwave Background are analyzed under the hypothesis that the same inflationary seed accounting for protogalactic magnetism also affects the Einstein-Boltzmann hierarchy whose initial conditions are assigned for typical correlation scales larger than the Hubble radius after matter-radiation equality but before decoupling. Since the primordial gauge spectrum depends on a combination of pivotal parameters of the concordance model, the angular power spectra of the temperature and of the polarization are computed, for the first time, in the presence of a putative large-scale magnetic field of inflationary origin and without supplementary hypotheses.

This analysis aims at exploring what can be said about the growth rate of magnetized inhomogeneities under two concurrent hypotheses: a phase of quasi-de Sitter dynamics driven by a single inflaton field and the simultaneous presence of a spectator field coupled to gravity and to the gauge sector. Instead of invoking ad hoc correlations between the various components, the system of scalar inhomogeneities is diagonalized in terms of two gauge-invariant quasi-normal modes whose weighted sum gives the curvature perturbations on comoving orthogonal hypersurfaces. The predominance of the conventional adiabatic scalar mode implies that the growth rate of magnetized inhomogeneities must not exceed 2.2 in Hubble units if the conventional inflationary phase is to last about 70 efolds and for a range of slow roll parameters between 0.1 and 0.001. Longer and shorter durations of the quasi-de Sitter stage lead, respectively, either to tighter or to looser bounds which are anyway more constraining than the standard backreaction demands imposed on the gauge sector. Since a critical growth rate of order 2 leads to a quasi-flat magnetic energy spectrum, the upper bounds on the growth rate imply a lower bound on the magnetic spectral index. The advantages of the uniform curvature gauge are emphasized and specifically exploited throughout the treatment of the multicomponent system characterizing this class of problems.

Primary graviton spectra, produced via stimulated emission from an initial Bose-Einstein distribution, are enhanced for typical scales larger than the redshifted thermal wavelength. A mixed state of phonons induces a secondary graviton spectrum which is hereunder computed in terms of three parameters (i.e. the number of phonon species, the tensor-to-scalar ratio and the thermal wavelengths of the mixture). The primary and secondary graviton spectra are shown to be sensitive, respectively, to the first-order and second-order correlation properties of the initial quantum mixture so that the semiclassical theory is argued to be generally inadequate in this context. For particular values of the parameters the secondary contribution may turn out to be comparable with the primary spectrum over large-scales.

As soon as the geometry expands quasi-exponentially the plasma sources are customarily tuned to zero as if the duration of the inflationary phase was immaterial for the gauge field fluctuations at large-scale. The serendipitous disappearance of the plasma (or even the partial neglect of its effects) depends on the symmetries of the system which are, in four space-time dimensions and in the simplest Abelian case, the invariance under Weyl rescaling and the electromagnetic duality symmetry. The quantum, thermal and conducting initial conditions of inflationary magnetogenesis are classified and discussed with the aim of determining when plasma effects can be effectively disregarded. The speculative implications of a non-degenerate monopole plasma for the conservation of the large-scale electric flux are briefly examined.

Quantum and thermal fluctuations of an irrotational fluid are studied across the transition regime connecting a protoinflationary phase of decelerated expansion to an accelerated epoch driven by a single inflaton field. The protoinflationary inhomogeneities are suppressed when the transition to the slow roll phase occurs sharply over space-like hypersurfaces of constant energy density. If the transition is delayed, the interaction of the quasi-normal modes related, asymptotically, to fluid phonons and inflaton quanta leads to an enhancement of curvature perturbations. It is shown that the dynamics of the fluctuations across the protoinflationary boundaries is determined by the monotonicity properties of the pump fields controlling the energy transfer between the background geometry and the quasi-normal modes of the fluctuations. After corroborating the analytical arguments with explicit numerical examples, general lessons are drawn on the classification of the protoinflationary transition.

We consider a globally neutral Lorentzian plasma as a possible remnant of a preinflationary stage of expansion and pose the problem of the suitable initial conditions for the evolution of the large-scale electromagnetic inhomogeneities. During the protoinflationary regime the Weyl invariance of the Ohmic current guarantees that the comoving conductivity is approximately constant. The subsequent breaking of Weyl invariance by the masses of the charge carriers drives the conductivity to zero. The newly derived conducting initial conditions for the amplification of large-scale magnetic fields are contrasted with the conventional vacuum initial conditions. It is shown, in a specific class of examples, that when the number of inflationary efolds is close to minimal the effects of the conducting initial conditions cannot be neglected.

The adiabatic inhomogeneities of the scalar curvature lead to a compressible flow affecting the dynamics of the hydromagnetic nonlinearities. The influence of the plasma on the evolution of a putative magnetic field is explored with the aim of obtaining an effective description valid for sufficiently large scales. The bulk velocity of the plasma, computed in the framework of the LambdaCDM scenario, feeds back into the evolution of the magnetic power spectra leading to a (nonlocal) master equation valid in Fourier space and similar to the ones discussed in the context of wave turbulence. Conversely, in physical space, the magnetic power spectra obey a Schroedinger-like equation whose effective potential depends on the large-scale curvature perturbations. Explicit solutions are presented both in physical space and in Fourier space. It is argued that curvature inhomogeneities, compatible with the WMAP 7yr data, shift to lower wavenumbers the magnetic diffusivity scale.

After electron-positron annihilation and prior to photon decoupling the magnetic Reynolds number is approximately twenty orders of magnitude larger than its kinetic counterpart which is, in turn, smaller than one. In this globally neutral system the large-scale inhomogeneities are provided by the spatial fluctuations of the scalar curvature. Owing to the analogy with the description of Markovian conducting fluids in the presence of acoustic fluctuations, the evolution equations of a putative magnetic field are averaged over the large-scale flow determined by curvature perturbations. General lessons are drawn on the typical diffusion scale of magnetic inhomogeneities. It is speculated that Reynolds numbers prior to electron-positron annihilation can be related to the entropy contained in the Hubble volume during the various stages of the evolution of the conducting plasma.

Exact solutions of the perturbations equations in the magnetized LambdaCDM scenario are presented. They apply during the dark ages and, more specifically, after the baryons are freed from the drag of the photons. The magnetized growth rate of matter perturbations is compared with the growth index obtained in the concordance paradigm and under the assumption that dark energy does not cluster for a redshift window ranging from the epoch of reionization to the stage of dark-energy dominance. The constraints derived from this analysis are shown to be qualitatively complementary and quantitatively competitive with the bounds stemming from the analysis of the distortion patterns induced by the magnetized adiabatic mode on the temperature and polarization anisotropies of the Cosmic Microwave Background.

The growth rate of matter perturbations is computed in a magnetized environment for the LambdaCDM and wCDM paradigms. It is argued that the baryons do not necessarily follow into the dark matter potential wells after they are released from the drag of the photons. The baryonic evolution equations inherit a forcing term whose explicit form depends on the plasma description and can be deduced, for instance, in the resistive magnetohydrodynamical approximation. After deriving an analytical expression for the growth rate applicable when dark energy does not cluster, the effects of relativistic corrections and of the inhomogeneities associated with the other species of the plasma are taken into account numerically. The spectral amplitudes and slopes of the stochastic magnetic background are selected to avoid appreciable distortions in the measured temperature and polarization anisotropies of the Cosmic Microwave Background. The growth of structures in the current paradigms of structure formation represents a complementary probe of large-scale magnetism in the same way as the shape of the growth factor and the associated indices can be used, in the conventional lore, to discriminate between competing scenarios of dark energy or even to distinguish different models of gravity.

The evolution equations of the vorticities of the electrons, ions and photons in a pre-decoupling plasma are derived, in a fully inhomogeneous geometry, by combining the general relativistic gradient expansion and the drift approximation within the Adler-Misner-Deser decomposition. The vorticity transfer between the different species is discussed in this novel framework and a set of general conservation laws, connecting the vorticities of the three-component plasma with the magnetic field intensity, is derived. After demonstrating that a source of large-scale vorticity resides in the spatial gradients of the geometry and of the electromagnetic sources, the total vorticity is estimated to lowest order in the spatial gradients and by enforcing the validity of the momentum constraint. By acknowledging the current bounds on the tensor to scalar ratio in the (minimal) tensor extension of the $\Lambda$CDM paradigm the maximal comoving magnetic field induced by the total vorticity turns out to be, at most, of the order of $10^{-37}$ G over the typical comoving scales ranging between 1 and 10 Mpc. While the obtained results seem to be irrelevant for seeding a reasonable galactic dynamo action, they demonstrate how the proposed fully inhomogeneous treatment can be used for the systematic scrutiny of pre-decoupling plasmas beyond the conventional perturbative expansions.

The properties of magnetized plasmas are always investigated under the hypothesis that the relativistic inhomogeneities stemming from the fluid sources and from the geometry itself are sufficiently small to allow for a perturbative description prior to photon decoupling. The latter assumption is hereby relaxed and pre-decoupling plasmas are described within a suitable expansion where the inhomogeneities are treated to a given order in the spatial gradients. It is argued that the (general relativistic) gradient expansion shares the same features of the drift approximation, customarily employed in the description of cold plasmas, so that the two schemes are physically complementary in the large-scale limit and for the low-frequency branch of the spectrum of plasma modes. The two-fluid description, as well as the magnetohydrodynamical reduction, are derived and studied in the presence of the spatial gradients of the geometry. Various solutions of the coupled system of evolution equations in the anti-Newtonian regime and in the quasi-isotropic approximation are presented. The relation of this analysis to the so-called separate Universe paradigm is outlined. The evolution of the magnetized curvature perturbations in the nonlinear regime is addressed for the magnetized adiabatic mode in the plasma frame.

The quantum theory of optical coherence is applied to the scrutiny of the statistical properties of the relic inflaton quanta. After adapting the description of the quantized scalar and tensor modes of the geometry to the analysis of intensity correlations, the normalized degrees of first-order and second-order coherence are computed in the concordance paradigm and are shown to encode faithfully the statistical properties of the initial quantum state. The strongly bunched curvature phonons are not only super-Poissonian but also super-chaotic. Testable inequalities are derived in the limit of large angular scales and can be physically interpreted in the light of the tenets of Hanbury Brown-Twiss interferometry. The quantum mechanical results are compared and contrasted with different situations including the one where intensity correlations are the result of a classical stochastic process. The survival of second-order correlations (not necessarily related to the purity of the initial quantum state) is addressed by defining a generalized ensemble where super-Poissonian statistics is an intrinsic property of the density matrix and turns out to be associated with finite volume effects which are expected to vanish in the thermodynamic limit.

The statistics of the curvature quanta generated during a stage of inflationary expansion is used to derive a count response model for the large-scale phonons determining, in the concordance lore, the warmer and the cooler spots of the large-scale temperature inhomogeneities. The multiplicity distributions for the counting statistics are shown to be generically overdispersed in comparison with conventional Poissonian regressions. The generalized count response model deduced hereunder accommodates an excess of correlations in the regime of high multiplicities and prompts dedicated analyses with forthcoming data collected by instruments of high angular resolution and high sensitivity to temperature variations per pixel.

The secondary spectra of the gravitons induced by a waterfall-like field are computed and the general bounds on the spectral energy density of the tensor modes of the geometry are translated into explicit constraints on the amplitude and slope of the waterfall spectrum. The obtained results are compared with the primary gravitational wave spectra of the concordance model and of its neighboring extensions as well as with the direct Ligo/Virgo bounds on stochastic backgrounds of relic gravitons. Space-borne interferometers (such as Lisa, Bbo, Decigo) seem to be less relevant but their potential implications are briefly outlined.

The multiplicity distributions produced by the variation of time-dependent gravitational fields in a conformally flat background geometry belong to the same class of infinitely divisible distributions found, for fixed centre of mass energies and symmetric (pseudo)rapidity intervals, in charged multiplicities produced in $pp$, $p\overline{p}$ and in heavy ion collisions. Apparently unrelated multiplicity distributions are classified in terms of the (positive) discrete representations of the $SU(1,1)$ group. The gravitational analogy suggest a global high-energy asymptote for the distributions measured in $pp$ and $p\overline{p}$ collisions. Second-order cross correlations between positively and negatively charged distributions represent a relevant diagnostic for a closer scrutiny of the multiparticle final state.

When the cosmic microwave photons scatter electrons just prior to the decoupling of matter and radiation, magnetic fields do contribute to the Stokes matrix as well as to the scalar, vector and tensor components of the transport equations for the brightness perturbations. The magnetized electron-photon scattering is hereby discussed in general terms by including, for the first time, the contribution of magnetic fields with arbitrary direction and in the presence of the scalar, vector and tensor modes of the geometry. The propagation of relic vectors and relic gravitons is discussed for a varying magnetic field orientation and for different photon directions. The source terms of the transport equations in the presence of the relativistic fluctuations of the geometry are also explicitly averaged over the magnetic field orientations and the problem of a consistent account of the small-scale and large-scale magnetic field is briefly outlined.

A primordial degree of circular polarization of the Cosmic Microwave Background is not observationally excluded. The hypothesis of primordial dichroism can be quantitatively falsified if the plasma is magnetized prior to photon decoupling since the initial V-mode polarization affects the evolution of the temperature fluctuations as well as the equations for the linear polarization. The observed values of the temperature and polarization angular power spectra are used to infer constraints on the amplitude and on the spectral slope of the primordial V-mode. Prior to photon decoupling magnetic fields play the role of polarimeters insofar as they unveil the circular dichroism by coupling the V-mode power spectrum to the remaining brightness perturbations. Conversely, for angular scales ranging between 4 deg and 10 deg the joined bounds on the magnitude of circular polarization and on the magnetic field intensity suggest that direct limits on the V-mode power spectrum in the range of 0.01 mK could directly rule out pre-decoupling magnetic fields in the range of 10-100 nG. The frequency dependence of the signal is located, for the present purposes, in the GHz range.

The tensor contribution to the $V$-mode polarization induced by a magnetized plasma at last scattering vanishes exactly. Conversely a polarized background of relic gravitons cannot generate a $V$-mode polarization. The reported results suggest that, in the magnetized $\Lambda$CDM paradigm, the dominant source of circular dichroism stems from the large-scale fluctuations of the spatial curvature.

When the last electron-photon scattering takes place in a magnetized environment, the degree of circular polarization of the outgoing radiation depends upon the magnetic field strength. After deriving the scattering matrix of the process, the generalized radiative transfer equations are deduced in the presence of the relativistic fluctuations of the geometry and for all the four brightness perturbations. The new system of equations is solved under the assumption that the incident radiation is not polarized. The induced V-mode polarization is analyzed both analytically and numerically. The corresponding angular power spectra are calculated and compared with the measured (or purported) values of the linear polarizations (i.e. E-mode and B-mode) as they arise in the concordance model and in its neighboring extensions. Possible connections between the V-mode polarization of the Cosmic Microwave background and the topological properties of the magnetic flux lines prior to equality are outlined and briefly explored in analogy with the physics of magnetized sun spots.

The V-mode polarization of the Cosmic Microwave Background is discussed in a weakly magnetized plasma. The VV and VT angular power spectra are computed for adiabatic initial conditions of the Einstein-Boltzmann hierarchy. Depending upon the frequency channel and upon the magnetic field intensity, the VT power spectra of the circular polarization can even be seven orders of magnitude larger than a putative B-mode polarization stemming from the lensing of the primary anisotropies. Specific programs aimed at the direct detection of the V-mode polarization of the Cosmic Microwave Background could provide a new observational tool for the scrutiny of pre-decoupling physics.

The large-scale curvature perturbations induced by spectator anisotropic stresses are analyzed across the matter-radiation transition. It is assumed that the anisotropic stress is associated with a plasma component whose energy density is subdominant both today and prior to photon decoupling. The enforcement of the momentum constraint and the interplay with the neutrino anisotropic stress determine the regular initial conditions of the Einstein-Boltzmann hierarchy. The Cosmic Microwave Background observables have shapes and phases which differ both from the ones of the conventional adiabatic mode as well as from their non-adiabatic counterparts.

The magnetogenesis scenarios triggered by the early variation of the gauge coupling are critically analyzed. In the absence of sources, it is shown that the electric and magnetic power spectra can be explicitly computed by means of electric-magnetic duality transformations. The remnants of a pre-inflationary expansion and the reheating process break explicitly electric-magnetic duality by inducing Ohmic currents. The generation of large-scale magnetic field and the physical distinction between electric and magnetic observables stems, in this class of models, from the final value reached by the conductivity of the plasma right after inflation. Specific numerical examples are given. The physical requirements of viable magnetogenesis scenarios are spelled out.

The impact of large-scale magnetic fields on the interplay between the ordinary and integrated Sachs-Wolfe effects is investigated in the presence of a fluctuating dark energy component. The modified initial conditions of the Einstein-Boltzmann hierarchy allow for the simultaneous inclusion of dark energy perturbations and of large-scale magnetic fields. The temperature and polarization angular power spectra are compared with the results obtained in the magnetized version of the (minimal) concordance model. Purported compensation effects arising at large scales are specifically investigated. The fluctuating dark energy component modifies, in a computable manner, the shapes of the 1- and 2-$\sigma$ contours in the parameter space of the magnetized background. The allowed spectral indices and magnetic field intensities turn out to be slightly larger than those determined in the framework of the magnetized concordance model where the dark energy fluctuations are absent.

Ohmic currents induced prior to decoupling are investigated in a standard transport model accounting both for the expansion of the background geometry as well as of its relativistic inhomogeneities. The relative balance of the Ohmic electric fields in comparison with the Hall and thermoelectric contributions is specifically addressed. The impact of the Ohmic currents on the evolution of curvature perturbations is discussed numerically and it is shown to depend explicitly upon the evolution of the conductivity.

Pre-decoupling magnetic fields affect the scalar modes of the geometry and produce observable effects which can be constrained also through the use of current (as opposed to forthcoming) data stemming from the Cosmic Microwave Background observations. The dependence of the temperature and polarization angular power spectra upon the parameters of an ambient magnetic field is encoded in the scaling properties of a set of basic integrals whose derivation is simplified in the limit of small angular scales. The magnetically-induced distortions patterns of the relevant observables can be computed analytically by employing scaling considerations which are corroborated by numerical results.

The temperature and polarization inhomogeneities of the Cosmic Microwave Background might bear the mark of pre-decoupling magnetism. The parameters of a putative magnetized background are hereby estimated from the observed temperature autocorrelation as well as from the measured temperature-polarization cross-correlation.

Stochastic backgrounds or relic gravitons, if ever detected, will constitute a prima facie evidence of physical processes taking place during the earliest stages of the evolution of the plasma. The essentials of the stochastic backgrounds of relic gravitons are hereby introduced and reviewed. The pivotal observables customarily employed to infer the properties of the relic gravitons are discussed both in the framework of the $\Lambda$CDM paradigm as well as in neighboring contexts. The complementarity between experiments measuring the polarization of the Cosmic Microwave Background (such as, for instance, WMAP, Capmap, Quad, Cbi, just to mention a few) and wide band interferometers (e.g. Virgo, Ligo, Geo, Tama) is emphasized. While the analysis of the microwave sky strongly constrains the low-frequency tail of the relic graviton spectrum, wide-band detectors are sensitive to much higher frequencies where the spectral energy density depends chiefly upon the (poorly known) rate of post-inflationary expansion.

Possible deviations from a radiation-dominated evolution, occurring prior the synthesis of light nuclei, impacted on the spectral energy density of high-frequency gravitons. For a systematic scrutiny of this situation, the $\Lambda$CDM paradigm must be complemented by (at least two) physical parameters describing, respectively, a threshold frequency and a slope. The supplementary frequency scale sets the lower border of a high-frequency domain where the spectral energy grows with a slope which depends, predominantly, upon the total sound speed of the plasma right after inflation. While the infra-red region of the graviton energy spectrum is nearly scale-invariant, the expected signals for typical frequencies larger than 0.01 nHz are hereby analyzed in a model-independent framework by requiring that the total sound speed of the post-inflationary plasma be smaller than the speed of light. Current (e.g. low-frequency) upper limits on the tensor power spectra (determined from the combined analysis of the three large-scale data sets) are shown to be compatible with a detectable signal in the frequency range of wide-band interferometers. In the present context, the scrutiny of the early evolution of the sound speed of the plasma can then be mapped onto a reliable strategy of parameter extraction including not only the well established cosmological observables but also the forthcoming data from wide band interferometers.

Even in the absence of a sizable tensor contribution, a B-mode polarization can be generated because of the competition between a pseudo-scalar background and pre-decoupling magnetic fields. By investigating the dispersion relations of a magnetoactive plasma supplemented by a pseudo-scalar interaction, the total B-mode polarization is shown to depend not only upon the plasma and Larmor frequencies but also on the pseudo-scalar rotation rate. If the (angular) frequency channels of a given experiment are larger than the pseudo-scalar rotation rate, the only possible source of (frequency dependent) B-mode autocorrelations must be attributed to Faraday rotation. In the opposite case the pseudo-scalar contribution dominates and the total rate becomes, in practice, frequency-independent. The B-mode cross-correlations can be used, under certain conditions, to break the degeneracy by disentangling the two birefringent contributions.

Polarimetry of the Cosmic Microwave Background (CMB) represents one of the possible diagnostics aimed at testing large-scale magnetism at the epoch of the photon decoupling. The propagation of electromagnetic disturbances in a magnetized plasma leads naturally to a B-mode polarization whose angular power spectrum is hereby computed both analytically and numerically. Combined analyses of all the publicly available data on the B-mode polarization are presented, for the first time, in the light of the magnetized $\Lambda$CDM scenario. Novel constraints on pre-equality magnetism are also derived in view of the current and expected sensitivities to the B-mode polarization.

Recent analyses of the WMAP 5-year data constrain possible non-adiabatic contributions to the initial conditions of CMB anisotropies. Depending upon the early dynamics of the plasma, the amplitude of the entropic modes can experience a different suppression by the time of photon decoupling. Explicit examples of the latter observation are presented both analytically and numerically when the post-inflationary dynamics is dominated by a stiff contribution.

Absent any indirect tests on the thermal history of the Universe prior to the formation of light nuclear elements, it is legitimate to investigate situations where, before nucleosyntheis, the sound speed of the plasma was larger than $c/\sqrt{3}$, at most equalling the speed of light $c$. In this plausible extension of the current cosmological paradigm, hereby dubbed Tensor-$\Lambda$CDM (i.e. T$\Lambda$CDM) scenario, high-frequency gravitons are copiously produced. Without conflicting with the bounds on the tensor to scalar ratio stemming from the combined analysis of the three standard cosmological data sets (i.e. cosmic microwave background anisotropies, large-scale structure and supenovae), the spectral energy density of the relic gravitons in the T$\Lambda$CDM scenario can be potentially observable by wide-band interferometers (in their advanced version) operating in a frequency window which ranges between few Hz and few kHz.

The high- and low-frequency descriptions of the pre-decoupling plasma are deduced from the Vlasov-Landau treatment generalized to curved space-times and in the presence of the relativistic fluctuations of the geometry. It is demonstrated that the interplay between one-fluid and two-fluid treatments is mandatory for a complete and reliable calculation of the polarization observables. The Einstein-Boltzmann hierarchy is generalized to handle the dispersive propagation of the electromagnetic disturbances in the pre-decoupling plasma. Given the improved physical and numerical framework, the polarization observables are computed within the magnetized $\Lambda$CDM paradigm (m$\Lambda$CDM). In particular, the Faraday-induced B-mode is consistently estimated by taking into account the effects of the magnetic fields on the initial conditions of the Boltzmann hierarchy, on the dynamical equations and on the dispersion relations. The complete calculations of the angular power spectra constitutes the first step for the derivation of magnetized maps of the CMB temperature and polarization which are here obtained for the first time and within the minimal m$\Lambda$CDM model. The obtained results set the ground for direct experimental scrutiny of large-scale magnetism via the low and high frequency instruments of the Planck explorer satellite.

The complete calculation of the CMB polarization observables (i.e. E- and B-modes) is reported within the conventional $\Lambda$CDM paradigm supplemented by a stochastic magnetic field. Intriguing perspectives for present and forthcoming CMB polarization experiments are outlined.

The most general initial conditions of CMB anisotropies, compatible with the presence of pre-equality magnetic fields, are derived. When the plasma is composed by photons, baryons, electrons, CDM particles and neutrinos, the initial data of the truncated Einstein-Boltzmann hierarchy contemplate one magnetized adiabatic mode and four (magnetized) non-adiabatic modes. After obtaining the analytical form of the various solutions, the Einstein-Boltzmann hierarchy is numerically integrated for the corresponding sets of initial data. The TT, TE and EE angular power spectra are illustrated and discussed for the magnetized generalization of the CDM-radiation mode, of the baryon-radiation mode and of the non-adiabatic mode of the neutrino sector. Mixtures of initial conditions are examined by requiring that the magnetized adiabatic mode dominates over the remaining non-adiabatic contributions. In the latter case, possible degeneracies between complementary sets of initial data might be avoided through the combined analysis of the TT, TE and EE angular power spectra at high multipoles (i.e. $\ell >1000$).

Large-scale magnetic fields affect the scalar modes of the geometry whose ultimate effect is to determine the anisotropies of the Cosmic Microwave Background (CMB in what follows). For the first time, a consistent numerical approach to the magnetized CMB anisotropies is pursued with the aim of assessing the angular power spectra of temperature and polarization when the scalar modes of the geometry and a stochastic background of inhomogeneous magnetic fields are simultaneously present in the plasma. The effects related to the magnetized nature of the plasma are taken into account both at the level of the dynamical equations and at the level of the initial conditions of the Einstein-Boltzmann hierarchy. The temperature and polarization observables are exploited to infer the peculiar signatures of a pre-equality magnetic field. Using the extrapolated best fit to the three year WMAP data the increase and distortions of the first seven peaks in the TT autocorrelations are monitored for different values of the regularized magnetic field intensity and for the physical range of spectral indices. Similar analyses are also conducted for the first few anticorrelation (and corrrelation) peaks of the TE power spectra. Possible interesting degeneracies and stimulating perspectives are pointed out and explored.

The standard $\Lambda$CDM paradigm is complemented with a magnetized contribution whose effects on the anisotropies of the Cosmic Microwave Background (CMB) are assessed by means of a dedicated numerical approach. The accuracy on the temperature and polarization correlations stems from the inclusion of the large-scale magnetic fields both at the level of the initial conditions and at the level of the Einstein-Boltzmann hierarchy which is consistently embedded in a generalized magnetohydrodynamical framework. Examples of the calculations of the temperature and polarization angular power spectra are illustrated and discussed. The reported results and the described numerical tools set the ground for a consistent inclusion of a magnetized contribution in current strategies of cosmological parameter estimation.

A viable class of magnetogenesis models can be constructed by coupling the kinetic term of the hypercharge to a spectator field whose dynamics does not affect the inflationary evolution. The magnetic power spectrum is explicitly related to the power spectrum of (adiabatic) curvature inhomogeneities when the quasi-de Sitter stage of expansion is driven by a single scalar degree of freedom. Depending upon the value of the slow-roll parameters, the amplitude of smoothed magnetic fields over a (comoving) Mpc scale can be as large as 0.01--0.1 nG at the epoch of the gravitational collapse of the protogalaxy. The contributions of the magnetic fields to the Sachs-Wolfe plateau and to the temperature autocorrelations in the Doppler region compare favourably with the constraints imposed by galactic magnetogenesis. Stimulating lessons are drawn on the interplay between magnetogenesis models and their possible CMB signatures.

Time-dependent soliton solutions are explicitly derived in a five-dimensional theory endowed with one (warped) extra-dimension. Some of the obtained geometries, everywhere well defined and technically regular, smoothly interpolate between two five-dimensional anti-de Sitter space-times for fixed value of the conformal time coordinate. Time dependent solutions containing both topological and non-topological sectors are also obtained. Supplementary degrees of freedom can be also included and, in this case, the resulting multi-soliton solutions may describe time-dependent kink-antikink systems.

It is shown that gravitating magnetic fields affect the evolution of curvature perturbations in a way that is reminiscent of a pristine non-adiabatic pressure fluctuation. The gauge-invariant evolution of curvature perturbations is used to constrain the magnetic power spectrum. Depending on the essential features of the thermodynamic history of the Universe, the explicit derivation of the bound is modified. The theoretical uncertainty in the constraints on the magnetic energy spectrum is assessed by comparing the results obtained in the case of the conventional thermal history with the estimates stemming from less conventional (but phenomenologically allowed) post-inflationary evolutions.

The cosmic microwave background (CMB) temperature autocorrelations, induced by a magnetized adiabatic mode of curvature inhomogeneities, are computed with semi-analytical methods. As suggested by the latest CMB data, a nearly scale-invariant spectrum for the adiabatic mode is consistently assumed. In this situation, the effects of a fully inhomogeneous magnetic field are scrutinized and constrained with particular attention to harmonics which are relevant for the region of Doppler oscillations. Depending on the parameters of the stochastic magnetic field a hump may replace the second peak of the angular power spectrum. Detectable effects on the Doppler region are then expected only if the magnetic power spectra have quasi-flat slopes and typical amplitude (smoothed over a comoving scale of Mpc size and redshifted to the epoch of gravitational collapse of the protogalaxy) exceeding 0.1 nG. If the magnetic energy spectra are bluer (i.e. steeper in frequency) the allowed value of the smoothed amplitude becomes, comparatively, larger (in the range of 20 nG). The implications of this investigation for the origin of large-scale magnetic fields in the Universe are discussed. Connections with forthcoming experimental observations of CMB temperature fluctuations are also suggested and partially explored.

The aim of these lectures is to introduce some basic problems arising in gravitation and modern cosmology. All along the discussion the guiding theme is provided by the phenomenological and theoretical properties of the Cosmic Microwave Background (CMB). These lectures have been prepared for a regular Phd course of the University of Milan-Bicocca.

The photon splitting gamma -> gamma gamma in a time-independent and inhomogeneous magnetized background is considered when neutral and ultralight spin-0 particles are coupled to two-photons. Depending on the inhomogeneity scale of the external field, resonant photon splitting can occur. If an optical laser crosses a magnetic field of few Tesla with typical inhomogeneity scale of the order of the meter, a potentially observable rate of photon splittings is expected for the PVLAS range of couplings and masses.

The main motivations and challenges related with the physics of large-scale magnetic fields are briefly analyzed. The interplay between large-scale magnetic fields and scalar CMB anisotropies is addressed with specific attention on recent progresses.

A class of five-dimensional warped solutions is presented. The geometry is everywhere regular and tends to five-dimensional anti-de Sitter space for large absolute values of the bulk coordinate. The physical features of the solutions change depending on the value of an integer parameter. In particular, a set of solutions describes generalized gravitating kinks where the scalar field interpolates between two different minima of the potential. The other category of solutions describes instead gravitating defects where the scalar profile is always finite and reaches the same constant asymptote both for positive and negative values of the bulk coordinate. In this sense the profiles are non-topological. The physical features of the zero modes are discussed.

The tight coupling expansion, appropriately generalized to include large-scale magnetic fields, allows the estimate of the brightness perturbations of CMB anisotropies for typical wavelengths that are larger than the Hubble radius after matter-radiation equality. After discussing the basic features of the the pre-decoupling initial conditions in the presence of fully inhomogeneous magnetic fields, the tight coupling expansion is studied both analytically and numerically. From the requirement that the amplitudes and phases of Sakharov oscillations are (predominantly) adiabatic and from the inferred value of the plateau in the temperature autocorrelation, the effects of the magnetized contribution can be systematically investigated and constrained.

An appropriate gauge-invariant framework for the treatment of magnetized curvature and entropy modes is developed. It is shown that large-scale magnetic fields, present after neutrino decoupling, affect curvature and entropy perturbations. The evolution of different magnetized modes is then studied across the matter-radiation transition both analytically and numerically. From the observation that, after equality (but before decoupling) the (scalar) Sachs-Wolfe contribution must be (predominantly) adiabatic, constraints on the magnetic power spectra are deduced. The present results motivate the experimental analysis of more general initial conditions of CMB anisotropies (i.e. mixtures of magnetized adiabatic and isocurvature modes during the pre-decoupling phase). The role of the possible correlations between the different components of the fluctuations is partially discussed.

Large-scale magnetic fields can affect scalar cosmological perturbations whose evolution is described in the conformally Newtonian gauge and within the tight coupling approximation. The magnetized curvature perturbations present after matter radiation equality (and prior to decoupling) are computed in terms of an appropriate transfer matrix allowing a general estimate of the Sachs-Wolfe plateau. From the observation that CMB initial conditions should be (predominantly) adiabatic, the contribution of the magnetic field intensity can be constrained.

The dynamical effects of the tensor modes of the geometry are investigated in the context of curvature bounces. Since the bouncing behaviour implies sharp deviations from a radiation-dominated evolution, significant back-reaction effects of relic gravitons may be expected at short wavelengths. After developing a general iterative framework for the calculation of dynamical back-reaction effects, explicit analytical and numerical examples are investigated for different parametrizations of the energy-momentum pseudo-tensor(s) of the produced gravitons. The reported results suggest that dynamical back-reaction effects are a necessary ingredient for a consistent description of bouncing models at late times.

Possible effects of large-scale magnetic fields on the Cosmic Microwave Background (CMB) are reviewed. Depending on the specific branch of the spectrum of plasma excitations, magnetic fields are treated either within a two-fluid plasma description or within an effective (one-fluid) approach. The uniform field approximation is contrasted with the fully inhomogeneous field approximation. It is argued that the interplay between CMB physics and large-scale magnetic fields will represent a rather interesting cross-disciplinary arena along the next few years.

Interacting fluids, endowed with bulk viscous stresses, are discussed in a unified perspective with the aim of generalizing the treatment of cosmological perturbation theory to the case where both fluctuating decay rates and fluctuating bulk viscosity coefficients are simultaneously present in the relativistic plasma. A gauge-invariant treatment of the qualitatively new phenomena arising in this context is provided. In a complementary approach, faithful gauge-fixed descriptions of the gravitational and hydrodynamical fluctuations are developed and exploited. To deepen the interplay between bulk viscous stresses and fluctuating decay rates, illustrative examples are proposed and discussed both analytically and numerically. Particular attention is paid to the coupled evolution of curvature and entropy fluctuations when, in the relativistic plasma, at least one of the interacting fluids possesses a fluctuating bulk viscosity coefficient. It is argued that this class of models may be usefully employed as an effective description of the decay of the inflaton as well as of other phenomena involving imperfect relativistic fluids.

Motivated by recent claims stating that the acceleration of the present Universe is due to fluctuations with wavelength larger than the Hubble radius, we present a general analysis of various perturbative solutions of fully inhomogeneous Einstein equations supplemented by a perfect fluid. The equivalence of formally different gradient expansions is demonstrated. If the barotropic index vanishes, the deceleration parameter is always positive semi-definite.

We investigate the way big rips are approached in a fully inhomogeneous description of the space-time geometry. If the pressure and energy densities are connected by a (supernegative) barotropic index, the spatial gradients and the anisotropic expansion decay as the big rip is approached. This behaviour is contrasted with the usual big-bang singularities. A similar analysis is performed in the case of sudden (quiescent) singularities and it is argued that the spatial gradients may well be non-negligible in the vicinity of pressure singularities.

Exact results stemming directly from Einstein equations imply that inhomogeneous Universes endowed with vanishing pressure density can only decelerate, unless the energy density of the Universe becomes negative. Recent proposals seem to argue that inhomogeneous (but isotropic) space-times, filled only with incoherent matter,may turn into accelerated Universes for sufficiently late times. To scrutinize these scenarios, fully inhomogeneous Einstein equations are discussed in the synchronous system. In a dust-dominated Universe, the inhomogeneous generalization of the deceleration parameter is always positive semi-definite implying that no acceleration takes place.

This review presents, in a self-consistent manner, those analytical tools that are relevant for the analysis of the physics of CMB anisotropies generated in different theoretical models of the early Universe. After introducing the physical foundations of the Sachs-Wolfe effect, the origin and evolution of the scalar, tensor and vector modes of the geometry is treated in both gauge-invariant and gauge-dependent descriptions. Some of the recent progresses in the theory of cosmological perturbations are scrutinized with particular attention to their implications for the adiabatic and isocurvature paradigms, whose description is reviewed both within conventional fluid approaches and within the the Einstein-Boltzmann treatment. Open problems and theoretical challenges for a unified theory of the early Universe are outlined in light of their implications for the generation of large-scale anisotropies in the CMB sky and in light of the generation of stochastic backgrounds of relic gravitons between few Hz and the GHz.

The evolution of the rotational inhomogeneities is investigated in the specific framework of four-dimensional pre-big bang models. While minimal (dilaton-driven) scenarios do not lead to rotational fluctuations, in the case of non-minimal (string-driven) models, fluid sources are present in the pre-big bang phase. The rotational modes of the geometry, coupled to the divergenceless part of the velocity field, can then be amplified depending upon the value of the barotropic index of the perfect fluids. In the light of a possible production of rotational inhomogeneities, solutions describing the coupled evolution of the dilaton field and of the fluid sources are scrutinized in both the string and Einstein frames. In semi-realistic scenarios, where the curvature divergences are regularized by means of a non-local dilaton potential, the rotational inhomogeneities are amplified during the pre-big bang phase but they decay later on. Similar analyses can also be performed when a contraction occurs directly in the string frame metric.

This paper introduces a systematic treatment of the linear theory of scalar gravitational perturbations in the presence of a fully inhomogeneous magnetic field. The analysis is conducted both in the synchronous and in the conformally Newtonian gauges. The cosmological plasma is assumed to be composed of cold dark mattter, baryons, photons, neutrinos. The problem of super-horizon initial conditions for the fluid variables of the various species and for the coupled system of Boltzmann-Einstein equations is discussed in the presence of an inhomogeneous magnetic field. The tight coupling approximation for the Boltzmann hierarchy is extended to the case where gravitating magnetic fields are included.

The polarization plane of the cosmic microwave background radiation can be rotated either in a magnetized plasma or in the presence of a quintessential background with pseudoscalar coupling to electromagnetism. A unified treatment of these two phenomena is presented for cold and warm electron-ion plasmas at the pre-recombination epoch. The electron temperature is only relevant to the relativistic correction of the cold plasma results. The spectrum of plasma excitations is obtained from a generalized Appleton--Hartree equation, describing simultaneously the high-frequency propagation of electromagnetic waves in a magnetized plasma with a dynamical quintessence field. It is shown that these two effects are comparable for the plausible range of parameters allowed by present constraints. It is then argued that the generalized expressions derived in the present study may be relevant for direct searches of a possible rotation of the cosmic microwave background polarization.

It is argued that in the case of a smooth transition across a (dilaton-driven) curvature bounce the growing mode of the vector fluctuations matches continuously with a decaying mode at later times. Analytical examples of this observation are given both in the presence and in the absence of fluid sources. In the case of multidimensional bouncing models the situation is different, since the system of differential equations describing the vector modes of the geometry has a richer structure. The amplification of the vector modes of the geometry is specifically investigated in a regular five-dimensional bouncing curvature model where scale factors of the external and internal manifolds evolve at a dual rate. Vector fluctuations, in this case, can be copiously produced and are continuous across the bounce. The relevance of these results is critically illustrated.

Cosmology, high-energy physics and astrophysics are converging on the study of large-scale magnetic fields. While the experimental evidence for the existence of large-scale magnetization in galaxies, clusters and superclusters is rather compelling, the origin of the phenomenon remains puzzling especially in light of the most recent observations. The purpose of the present review is to describe the physical motivations and some of the open theoretical problems related to the existence of large-scale magnetic fields.

String-inspired cosmologies, whereby a non-singular curvature bounce is induced by a general-covariant, $T$-duality-invariant, non-local dilaton potential, are used to study numerically how inhomogeneities evolve and to compare the outcome with analytic expressions obtained through different matching conditions across the bounce. Good agreement is found if continuity across the bounce is assumed to hold for $\cal{R}$, the curvature perturbation on comoving hypersurfaces, rather than for the Bardeen potential.

The ratio of the curvaton energy density to that of the dominant component of the background sources may be constant during a significant period in the evolution of the Universe. The possibility of having tracking curvatons, whose decay occurs prior to the nucleosynthesis epoch, is studied. It is argued that the tracking curvaton dynamics is disfavoured since the value of the curvature perturbations prior to curvaton decay is smaller than the value required by observations. It is also argued, in a related context, that the minimal inflationary curvature scale compatible with the curvaton paradigm may be lowered in the case of low-scale quintessential inflation.

Quantum-mechanical initial conditions for the fluctuations of the geometry can be assigned in excess of a given physical wavelength. The two-point functions of the scalar and tensor modes of the geometry will then inherit corrections depending on which Hamiltonian is minimized at the initial stage of the evolution. The energy density of the background geometry is compared with the energy-momentum pseudo-tensor of the fluctuations averaged over the initial states, minimizing each different Hamiltonian. The minimization of adiabatic Hamiltonians leads to initial states whose back-reaction on the geometry is negligible. The minimization of non-adiabatic Hamiltonians, ultimately responsible for large corrections in the two-point functions, is associated with initial states whose energetic content is of the same order as the energy density of the background.

We present a comprehensive study of the cosmological solutions of 6D braneworld models with azimuthal symmetry in the extra dimensions, moduli stabilization by flux or a bulk scalar field, and which contain at least one 3-brane that could be identified with our world. We emphasize an unusual property of these models: their expansion rate depends on the 3-brane tension either not at all, or in a nonstandard way, at odds with the naive expected dimensional reduction of these systems to 4D general relativity at low energies. Unlike other braneworld attempts to find a self-tuning solution to the cosmological constant problem, the apparent failure of decoupling in these models is not associated with the presence of unstabilized moduli; rather it is due to automatic cancellation of the brane tension by the curvature induced by the brane. This provides some corroboration for the hope that these models provide a distinctive step toward understanding the smallness of the observed cosmological constant. However, we point out some challenges for obtaining realistic cosmology within this framework.

We complement the low-energy gravi-dilaton effective action of string theory with a non-local, general-covariant dilaton potential, and obtain homogeneous solutions describing a non-singular (bouncing-curvature) cosmology. We then compute, both analytically and numerically, the spectrum of amplified scalar and tensor perturbations, and draw some general lessons on how to extract observable consequences from pre-big bang and ekpyrotic scenarios.

The interplay among the possible variation of gauge coupling and the inflationary dynamics is investigated in a simplified toy model. Depending upon various parameters (scalar mass, curvature scale at the end of inflation and at the onset of the radiation epoch), the two-point function of the magnetic inhomogeneities grows during the de Sitter stage of expansion and, consequently, large scale magnetic fields are generated. The requirements coming from inflationary magnetogenesis are examined together with the theoretical constraints stemming from the explicit model of the evolution of the gauge coupling.

The so-called curvaton mechanism --a way to convert isocurvature perturbations into adiabatic ones-- is investigated both analytically and numerically in a pre-big bang scenario where the role of the curvaton is played by a sufficiently massive Kalb--Ramond axion of superstring theory. When combined with observations of CMBR anisotropies at large and moderate angular scales, the present analysis allows us to constrain quite considerably the parameter space of the model: in particular, the initial displacement of the axion from the minimum of its potential and the rate of evolution of the compactification volume during pre-big bang inflation. The combination of theoretical and experimental constraints favours a slightly blue spectrum of scalar perturbations, and/or a value of the string scale in the vicinity of the SUSY-GUT scale.

Large scale magnetic fields represent a triple point where cosmology, high-energy physics and astrophysics meet for different but related purposes. After reviewing the implications of large scale magnetic fields in these different areas, the role of primordial magnetic fields is discussed in various physical processes occurring prior to the decoupling epoch with particular attention to the big bang nucleosynthesis (BBN) epoch and to the electroweak (EW) epoch. The generation of matter--antimatter isocurvature fluctuations, induced by hypermagnetic fields, is analyzed in light of a possible increase of extra-relativistic species at BBN. It is argued that stochastic GW backgrounds can be generated by hypermagnetic fields at the LISA frequency. The problem of the origin of large scale magnetic fields is also scrutinized.

We present the results of a detailed study of how isocurvature axion fluctuations are converted into adiabatic metric perturbations through axion decay, and discuss the constraints on the parameters of pre-big bang cosmology needed for consistency with present CMB-anisotropy data. The large-scale normalization of temperature fluctuations has a non-trivial dependence both on the mass and on the initial value of the axion. In the simplest, minimal models of pre-big bang inflation, consistency with the COBE normalization requires a slightly tilted (blue) spectrum, while a strictly scale-invariant spectrum requires mild modifications of the minimal backgrounds at large curvature and/or string coupling.

Provided that matter-antimatter domains are present at the onset of big bang nucleosynthesis (BBN), the number of allowed additional relativistic species increases, compared to the standard scenario when matter-antimatter domains are absent. The extra relativistic species may take the form of massless fermions or even massless bosons, like relic gravitons. The number of additional degrees of freedom compatible with BBN depends, in this framework, upon the typical scale of the domains and the antimatter fraction. Since the presence of matter-antimatter domains allows a reduction of the neutron to proton ratio prior to the formation of 4He, large amounts of radiation-like energy density are allowed. Our results are compared with other constraints on the number of supplementary relativistic degrees of freedom and with different nonstandard BBN scenarios where the reduction of the neutron to proton ratio occurs via a different physical mechanism. The implications of our considerations for the upper limits on stochastic gravitational waves backgrounds of cosmological origin are outlined.

Problems and perspectives concerning the existence of large-scale magnetic fields are described. Heeding observations, possible origins and implications of magnetic fields in spiral galaxies and in regular clusters are scrutinized in different cosmological frameworks including scenarios where conformal invariance is broken because of the evolution of the gauge couplings.The generation of the BAU via hypermagnetic knots is reviewed. It is also argued that stochastic GW backgrounds can be generated in the LISA frequency range thanks to the presence of hypermagnetic fields.

The sensitivity achievable by a pair of VIRGO detectors to stochastic and isotropic gravitational wave backgrounds produced in pre-big-bang models is discussed in view of the development of a second VIRGO interferometer. We describe a semi-analytical technique allowing to compute the signal-to-noise ratio for (monotonic or non-monotonic) logarithmic energy spectra of relic gravitons of arbitrary slope. We apply our results to the case of two correlated and coaligned VIRGO detectors and we compute their achievable sensitivities. We perform our calculations both for the usual case of minimal string cosmological scenario and in the case of a non-minimal scenario where a long dilaton dominated phase is present prior to the onset of the ordinary radiation dominated phase. In this framework, we investigate possible improvements of the achievable sensitivities by selective reduction of the thermal contributions (pendulum and pendulum's internal modes) to the noise power spectra of the detectors. Since a reduction of the shot noise does not increase significantly the expected sensitivity of a VIRGO pair (in spite of the relative spatial location of the two detectors) our findings support the experimental efforts directed towards a substantial reduction of thermal noise.

It is shown that the evolution of the (Abelian) gauge coupling during an inflationary phase of de Sitter type drives the growth of the two-point function of the magnetic inhomogeneities. After examining the constraints on the variation of the gauge coupling arising in a standard model of inflationary and post-inflationary evolution, magnetohydrodynamical equations are generalized to the case of time evolving gauge coupling. It is argued that large scale magnetic fields can be copiously generated. Other possible implications of the model are outlined.

In this note we argue that the breaking of conformal invariance because of the coupling of a charged scalar field to gravity is not sufficient for the production of seed galactic magnetic fields during inflation.

The interplay between different ground based detectors and stochastic backgrounds of relic GW is described. A simultaneous detection of GW in the kHz and in the MHz--GHz region can point towards a cosmological nature of the signal. The sensitivity of a pair of VIRGO detectors to string cosmological models is presented. The implications of microwave cavities for stochastic GW backgrounds are discussed.

The dynamical evolution of internal space-like dimensions breaks the invariance of the Maxwell's equations under Weyl rescaling of the (conformally flat) four-dimensional metric. Depending upon the number and upon the dynamics of internal dimensions large scale magnetic fields can be created. The requirements coming from magnetogenesis together with the other cosmological constraints are examined under the assumption that the internal dimensions either grow or shrink (in conformal time) prior to a radiation dominated epoch. If the internal dimensions are growing the magnitude of the generated magnetic fields can seed the galactic dynamo mechanism.

The hot plasma above the electroweak scale contains (hyper) charged scalar particles which are coupled to Abelian gauge fields. Scalars may interact with gravity in a non-conformally invariant way and thus their fluctuations can be amplified during inflation. These fluctuations lead to creation of electric currents and produce inhomogeneous distribution of charge density, resulting in the generation of cosmological magnetic fields. We address the question whether these fields can be coherent at large scales so that they may seed the galactic magnetic fields. Depending upon the mass of the charged scalar and upon various cosmological (critical fraction of energy density in matter, Hubble constant) and particle physics parameters we found that the magnetic fields generated in this way are much larger than vacuum fluctuations. However, their amplitude on cosmological distances is found to be too small for seeding the galactic magnetic fields.

Physical scenarios, leading to highly energetic stochastic gravitational waves backgrounds (for frequencies ranging from the $\mu$Hz up to the GHz) are examined. In some cases the typical amplitude of the logarithmic energy spectrum can be even eight orders of magnitude larger than the ordinary inflationary prediction. Scaling violations in the frequency dependence of the energy density of the produced gravitons are discussed.

The sensitivity of a pair of VIRGO interferometers to gravitational waves backgrounds of cosmological origin is analyzed for the cases of maximal and minimal overlap of the two detectors. The improvements in the detectability prospects of scale-invariant and non-scale-invariant logarithmic energy spectra of relic gravitons are discussed.

The sensitivity achievable by a pair of VIRGO detectors to stochastic and isotropic gravitational wave backgrounds of cosmological origin is discussed in view of the development of a second VIRGO interferometer. We describe a semi-analytical technique allowing to compute the signal-to-noise ratio for (monotonic or non-monotonic) logarithmic energy spectra of relic gravitons of arbitrary slope. We apply our results to the case of two correlated and coaligned VIRGO detectors and we compute their achievable sensitivities. The maximization of the overlap reduction function is discussed. We focus our attention on a class of models whose expected sensitivity is more promising, namely the case of string cosmological gravitons. We perform our calculations both for the case of minimal string cosmological scenario and in the case of a non-minimal scenario where a long dilaton dominated phase is present prior to the onset of the ordinary radiation dominated phase. In this framework, we study possible improvements of the achievable sensitivities by selective reduction of the thermal contributions (pendulum and pendulum's internal modes) to the noise power spectra of the detectors. Since a reduction of the shot noise does not increase significantly the expected sensitivity of a VIRGO pair (in spite of the relative spatial location of the two detectors) our findings support the experimental efforts directed towards a substantial reduction of thermal noise.

At finite hyperconductivity and finite fermionic density the flux lines of long range hypermagnetic fields may not have a topologically trivial structure. The combined evolution of the chemical potentials and of pseudoscalar fields (like the axial Higgs), possibly present for temperatures in the TeV range, can twist the hypercharge flux lines, producing, ultimately, hypermagnetic knots (HK). The dynamical features of the HK depend upon the various particle physics parameters of the model (pseudoscalar masses and couplings, strength of the electroweak phase transition, hyperconductivity of the plasma) and upon the magnitude of the primordial flux sitting in topologically trivial configurations of the hypermagnetic field. We study different cosmological scenarios where HK can be generated. We argue that the fermionic number sitting in HK can be released producing a seed for the Baryon Asymmetry of the Universe (BAU) provided the typical scale of the knot is larger than the diffusivity length scale. We derive constraints on the primordial hypermagnetic flux required by our mechanism and we provide a measure of the parity breaking by connecting the degree of knottedness of the flux lines to the BAU. We rule out the ordinary axion as a possible candidate for production (around temperatures of the order of the GeV) of {\em magnetic} knots since the produced {\em electromagnetic} helicity is negligible (for cosmological standard) if the initial amplitude of the axion oscillations is of the order of the Peccei-Quinn breaking scale.

There are no reasons why the energy spectra of the relic gravitons, amplified by the pumping action of the background geometry, should not increase at high frequencies. A typical example of this behavior are quintessential inflationary models where the slopes of the energy spectra can be either blue or mildly violet. In comparing the predictions of scenarios leading to blue and violet graviton spectra we face the problem of correctly deriving the sensitivities of the interferometric detectors. Indeed, the expression of the signal-to-noise ratio not only depends upon the noise power spectra of the detectors but also upon the spectral form of the signal and, therefore, one can reasonably expect that models with different spectral behaviors will produce different signal-to-noise ratios. By assuming monotonic (blue) spectra of relic gravitons we will give general expressions for the signal-to-noise ratio in this class of models. As an example we studied the case of quintessential gravitons. The minimum achievable sensitivity to $h^2_{0} \Omega_{GW}$ of different pairs of detectors is computed, and compared with the theoretical expectations.

Topologically non-trivial configurations of the hypermagnetic flux lines lead to the formation of hypermagnetic knots (HK) whose decay might seed the Baryon Asymmetry of the Universe (BAU).HK can be dynamically generated provided a topologically trivial (i.e. stochastic) distribution of flux lines is already present in the symmetric phase of the electroweak (EW) theory. In spite of the mechanism generating the HK, their typical size must exceed the diffusivity length scale. In the minimal standard model (MSM) (but not necessarily in its supersymmetric extension) HK are washed out. A classical hypermagnetic background in the symmetric phase of the EW theory can produce interesting amounts of gravitational radiation.

If the energy momentum tensor contains bulk viscous stresses violating the dominant energy condition (DOC) the energy spectra of the relic gravitons (produced at the time of the DOC's violation) increase in frequency in a calculable way. In a general relativistic context we give examples where the DOC is only violated for a limited amount of time after which the ordinary (radiation dominated) evolution takes place. We connect our discussion to some recent remarks of Grishchuk concerning the detectability of the stochastic gravitational wave background by the forthcoming interferometric detectors.

The energy spectra of gravitational waves (GW) produced in quintessential inflationary models increase in frequency and exhibit a sharp spike around 170 GHz where the associated fraction of critical energy density today stored in relic gravitons is of the order of $10^{-6}$. We contrast our findings with the spectra of ordinary inflationary models and we comment about possible detetction strategies of the spike.

A large class of quintessential inflationary models, recently proposed by Peebles and Vilenkin, leads to post-inflationary phases whose effective equation of state is stiffer than radiation. The expected gravitational waves logarithmic energy spectra are tilted towards high frequencies and characterized by two parameters: the inflationary curvature scale at which the transition to the stiff phase occurs and the number of (non conformally coupled) scalar degrees of freedom whose decay into fermions triggers the onset of a gravitational reheating of the Universe. Depending upon the parameters of the model and upon the different inflationary dynamics (prior to the onset of the stiff evolution) the relic gravitons energy density can be much more sizeable than in standard inflationary models, for frequencies larger than 1 Hz. We estimate the required sensitivity for detection of the predicted spectral amplitude and show that the allowed region of our parameter space leads to a signal smaller (by one 1.5 orders of magnitude) than the advanced LIGO sensitivity at a frequency of 0.1 KHz. The maximal signal, in our context, is expected in the GHz region where the energy density of relic gravitons in critical units (i.e. $h_0^2 \Omega_{GW}$) is of the order of $10^{-6}$, roughly eight orders of magnitude larger than in ordinary inflationary models. Smaller detectors (not necessarily interferometers) can be relevant for detection purposes in the GHz frequency window. We suggest/speculate that future measurements through microwave cavities can offer interesting perspectives.

The dynamical symmetries of hot and electrically neutral plasmas in a highly conducting medium suggest that, after the epoch of the electron-positron annihilation, magnetohydrodynamical configurations carrying a net magnetic helicity can be present. The simultaneous conservation of the magnetic flux and helicity implies that the (divergence free) field lines will possess inhomogeneous knot structures acting as source seeds in the evolution equations of the scalar, vector and tensor fluctuations of the background geometry. We give explicit examples of magnetic knot configurations with finite energy and we compute the induced metric fluctuations. Since magnetic knots are (conformally) coupled to gravity via the vertex dictated by the equivalence principle, they can imprint spikes in the gravitational wave spectrum for frequencies compatible with the typical scale of the knot corresponding, in our examples, to a present frequency range of $10^{-11}$--$10^{-12}$ Hertz. At lower frequencies the spectrum is power-suppressed and well below the COBE limit. For smaller length scales (i.e. for larger frequencies) the spectrum is exponentially suppressed and then irrelevant for the pulsar bounds. Depending upon the number of knots of the configuration, the typical amplitude of the gravitational wave logarithmic energy spectrum (in critical units) can be even four orders of magnitude larger than the usual flat (inflationary) energy spectrum generated thanks to the parametric amplification of the vacuum fluctuations.

We point out that the existence of post-inflationary phases stiffer than radiation leads to the production of stochastic gravitational waves (GW) backgrounds whose logarithmic energy spectra (in critical units) are typically "blue" at high frequencies. The maximal spectral slope (for present frequencies larger than $10^{-16}$ Hz) is of order one and it is related to the maximal sound velocity of the stiff plasma governing the evolution of the geometry. The duration of the stiff phase is crucially determined by the backreaction of the GW leaving the horizon during the de Sitter phase and re-entering during the stiff phase. Therefore, the maximal (inflationary) curvature scale has to be fine-tuned to a value smaller than the limits set by the large scale measurements ($H_{dS}\laq 10^{-6}~M_{P}$) in order to a have a sufficiently long stiff phase reaching an energy scale of the order of the $1~{\rm TeV}$ and even lower if we want the stiff phase to touch the hadronic era (corresponding to $T_{had}\sim 140~{\rm MeV}$). By looking more positively at our exercise we see that, if an inflationary phase is followed by a stiff phase, there exist an appealing possibility of "graviton reheating" whose effective temperature can be generally quite low.

The high-temperature plasma above the electroweak scale $\sim 100$ GeV may have contained a primordial hypercharge magnetic field whose anomalous coupling to the fermions induces a transformation of the hypermagnetic energy density into fermionic number. In order to describe this process, we generalize the ordinary magnetohydrodynamical equations to the anomalous case. We show that a not completely homogeneous hypermagnetic background induces fermion number fluctuations, which can be expressed in terms of a generic hypermagnetic field configuration. We argue that, depending upon the various particle physics parameters involved in our estimate (electron Yukawa coupling, strength of the electroweak phase transition) and upon the hypermagnetic energy spectrum, sizeable matter-antimatter fluctuations can be generated in the plasma. These fluctuations may modify the predictions of the standard Big Bang nucleosynthesis (BBN). We derive constraints on the magnetic fields from the requirement that the homogeneous BBN is not changed. We analyse the influence of primordial magnetic fields on the electroweak phase transition and show that some specific configurations of the magnetic field may be converted into net baryon number at the electroweak scale.

We show that the presence of primordial stochastic (hypercharge) magnetic fields before the electroweak (EW) phase transition induces isocurvature fluctuations (baryon number inhomogeneities). Depending on the details of the magnetic field spectrum and on the particle physics parameters (such as the strength of the EW phase transition and electron Yukawa couplings) these fluctuations may survive until the Big Bang nucleosynthesis (BBN). Their lenghtscale may exceed the neutron diffusion length at that time, while their magnitude can be so large that sizable antimatter domains are present. This provides the possibility of a new type of initial conditions for non-homogeneous BBN or, from a more conservative point of view, stringent bounds on primordial magnetic fields.

A magnetic field, coherent over the horizon size at the decoupling and strong enough to rotate the polarization plane of the CMBR, can be generated from the electromagnetic vacuum fluctuations amplified by the space-time evolution of the dilaton coupling. The possible relevance of this result for superstring inspired cosmological models is discussed. Particular attention will be paid to the connection between Faraday rotation signals and stochastic gravity-wave backgrounds.

We discuss a gauge invariant approach to the theory of cosmological perturbations in a higher-dimensonal background. We find the normal modes which diagonalize the perturbed action, for a scalar field minimally coupled to gravity, in a higher-dimensional manifold M of the Bianchi-type I, under the assumption that the translations along an isotropic spatial subsection of M are isometries of the full, perturbed background. We show that, in the absence of scalar field potential, the canonical variables for scalar and tensor metric perturbations satisfy exactly the same evolution equation, and we discuss the possible dependence of the spectrum on the number of internal dimensions.

We discuss some perturbative techniques suitable for the gauge-invariant treatment of the scalar and tensor inhomogeneities of an anisotropic and homogeneous background geometry whose spatial section naturally decomposes into the direct product of two maximally symmetric Eucledian manifolds, describing a general situation of dimensional decoupling in which $d$ external dimensions evolve (in conformal time) with scale factor $a(\eta)$ and $n$ internal dimensions evolve with scale factor $b(\eta)$. We analyze the growing mode problem which typically arises in contracting backgrounds and we focus our attention on the situation where the amplitude of the fluctuations not only depends on the external space-time but also on the internal spatial coordinates. In order to illustrate the possible relevance of this analysis we compute the gravity waves spectrum produced in some highly simplified model of cosmological evolution and we find that the spectral amplitude, whose magnitude can be constrained by the usual bounds applied to the stochastic gravity waves backgrounds, depends on the curvature scale at which the compactification occurs and also on the typical frequency of the internal excitations.

A spectrum of relic stochastic gravitational radiation, strongly tilted towards high frequencies, and characterized by two basic parameters is shown to emerge in a class of string theory models. We estimate the required sensitivity for detection of the predicted gravitational radiation and show that a region of our parameter space is within reach for some of the plannedgravitational-wave detectors.

A large class of string-cosmology backgrounds leads to a spectrum of relic stochastic gravitational waves, strongly tilted towards high frequencies, and characterized by two basic parameters of the cosmological model. We estimate the required sensitivity for detection of the predicted gravitational radiation and show that a region of our parameter space is within reach for some of the planned gravitational-wave detectors.

In the inflationary scenarios suggested by string theory, the vacuum fluctuations of the electromagnetic field can be amplified by the time-evolution of the dilaton background, and can grow large enough to explain both the origin of the cosmic magnetic fields and of the observed CMB anisotropy. The normalization of the perturbation spectrum is fixed, and implies a relation between the perturbation amplitude at the COBE scale and the spectral index $n$. Working within a generic two-parameter family of backgrounds, a large scale anisotropy $\Delta T/T\simeq 10^{-5}$ is found to correspond to a spectral index in the range $n\simeq 1.11 - 1.17$.

Sufficiently large seeds for generating the observed (inter)galactic magnetic fields emerge naturally in string cosmology from the amplification of electromagnetic vacuum fluctuations due to a dynamical dilaton background. The success of the mechanism depends crucially on two features of the so-called pre-big-bang scenario, an early epoch of dilaton-driven inflation at very small coupling, and a sufficiently long intermediate stringy era preceding the standard radiation-dominated evolution.

We compute the spectrum of scalar and tensor metric perturbations generated, as amplified vacuum fluctuations, during an epoch of dilaton-driven inflation of the type occurring naturally in string cosmology. In the tensor case the computation is straightforward while, in the scalar case, it is made delicate by the appearance of a growing mode in the familiar longitudinal gauge. In spite of this, a reliable perturbative calculation of perturbations far outside the horizon can be performed by resorting either to appropriate gauge invariant variables, or to a new coordinate system in which the growing mode can be "gauged down". The simple outcome of this complicated analysis is that both scalar and tensor perturbations exhibit nearly Planckian spectra, whose common "temperature" is related to some very basic parameters of the string-cosmology background.

The entropy growth in a cosmological process of pair production is completely determined by the associated squeezing parameter, and is insensitive to the number of particles in the initial state. The total produced entropy may represent a significant fraction of the entropy stored today in the cosmic black-body radiation, provided pair production originates from a change in the background metric at a curvature scale of the Planck order.

We consider the stimulated emission of gravitons from an initial state of thermal equilibrium, under the action of the cosmic gravitational background field. We find that the low-energy graviton spectrum is enhanced if compared with spontaneous creation from the vacuum; as a consequence, the scale of inflation must be lowered, in order not to exceed the observed CMB quadrupole anisotropy. This effect is particularly important for models based on a symmetry-breaking transition which require, as initial condition, a state of thermal equilibrium at temperatures of the order of the inflation scale.