Effective horizons, junction conditions and large-scale magnetism

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 qualitative description of large-scale cosmological perturbations [1,2,3] stipulates that a given wavelength exits the Hubble radius at some typical conformal time τ ex during an inflationary stage of expansion and approximately reenters at τ re , when the Universe still expands but in a decelerated manner. By a mode being beyond the horizon we only mean that the physical wavenumber is much less than the expansion rate: this does not necessarily have anything to do with causality [2]. Indeed, the initial conditions of the Einstein-Boltzmann hierarchy (mandatory for the calculation of the temperature and polarization anisotropies) are set when the relevant modes are larger than the Hubble radius prior to matter-radiation equality [3]. Similarly the physical wavenumbers of the hyperelectric and hypermagnetic fields can be much smaller than the rate of variation of the susceptibility (χ in what follows) which now plays the role of the effective horizon. The junction conditions of the gauge power spectra will be derived in general terms and then corroborated by specific examples with particular attention to the the case of sudden (but continuous) postinflationary transitions. Using the obtained results the gauge power spectra will be computed in the case of generalized quantum mechanical initial conditions of the hypercharge field.
The four-dimensional action discussed in [4] concisely summarizes a large class of magnetogenesis scenarios and it can be written, for the present ends, as 2 : where g denotes the determinant of the four-dimensional metric 3 ; Y αβ and Y αβ are, respectively, the gauge field strength and its dual. While the two symmetric tensors M ρ σ and N ρ σ parametrize, in full generality, the dependence upon the electric and magnetic susceptibilities, Eq. (1) includes, as a special case, the derivative couplings typical of the relativistic theory of Casimir-Polder and Van der Waals interactions [5]. Even though the whole discussion could be carried on in the case of unequal magnetic and electric susceptibilities by using the results reported in [4], for the sake of simplicity the attention will now be focussed on the case M ρ σ = N ρ σ = (λ/2)δ ρ σ . In this instance the evolution equations derived from Eq. (1) are: where, as already mentioned, χ = √ λ represents the susceptibility and the prime denotes a derivation with respect to the conformal time coordinate. As implied by the duality 2 We shall be working in a conformally flat background geometry g µν = a 2 (τ )η µν where η µν denotes the Minkowski metric, a(τ ) is the scale factor and τ parametrizes the conformal time coordinate. The components of the Abelian field strength appearing in Eq. (1) are Y 0i = e i /a 2 and Y ij = − ijk b k /a 2 . The canonical electric and magnetic fields of Eq. (2) are defined as B = a 2 √ λ b and E = a 2 √ λ e. 3 We shall be working in a conformally flat background geometry g µν = a 2 (τ )η µν where η µν denotes the Minkowski metric, a(τ ) is the scale factor and τ parametrizes the conformal time coordinate. The components of the Abelian field strength appearing in Eq. (1) are Y 0i = e i /a 2 and Y ij = − ijk b k /a 2 . The canonical electric and magnetic fields of Eq. (2) are defined as B = a 2 √ λ b and E = a 2 √ λ e. symmetry [6], when χ → 1/χ (i.e. F → −F) the two equations appearing in Eq. (2) are interchanged provided E → − B and B → E. Equations (1) and (2) contain, as a particular case, a class of magnetogenesis models based on the evolution of the inflaton or of some other spectator field (see, e. g. [7,8,9] for an incomplete list of references 4 ). Various scenarios aim at producing magnetic fields with approximate intensities of a few hundredths of a nG (1 nG = 10 −9 G) and over typical comoving scales between few Mpc and 100 Mpc. When the intensities are much lower than O(10 −3 ) nG a dynamo action (of some sort) seems mandatory (see, for instance, Ref. [11] for a time ordered but still incomplete list of review articles).
In conformally flat backgrounds geometries of Friedmann-Robertson-Walker type the Coulomb gauge condition (i.e. Y 0 = 0 and ∇ · Y = 0) is preferable since it is preserved (unlike the Lorentz gauge condition) by a conformal rescaling of the metric; with this choice, when the electric and magnetic susceptibilities coincide, Eq. (1) reduces to: where A = λ/(4π) Y . In terms of A and of its conjugate momentum Π the canonical Hamiltonian derived from the action (3) is: In terms of the normal modes A he hyperelectric and hypermagnetic fields of Eq. (2) are defined, respectively, as E = − Π(τ, x) and B = ∇ × A; the corresponding field operators in the Heisenberg description are: where the sum is performed over the physical polarizations e (α) i while the mode functions f k and g k obey, in the absence of conductivity, the following pair of dual equations: From Eq. (7) two (decoupled) second-order differential equations can be derived for f k (i.e. f k + [k 2 − χ /χ]f k = 0) and for g k (i.e. g k + [k 2 − χ(1/χ) ]g k = 0). To guarantee the correct formulation of the Cauchy problem the initial conditions must be assigned in agreement with g k = f k − Ff k but without the continuity of χ (and of its first derivative) the pump fields χ /χ and χ(χ −1 ) will be singular at the transition points. Moreover, to enforce the canonical form of the commutation relations 5 the Wronskian W(τ ) = f k (τ )g * k (τ )−f * k (τ )g k (τ ) (conserved and invariant under duality) must be normalized as W(τ ) = i.
Using Eqs. (5) and (6) the magnetic and electric field operators are expressible in Fourier space and their expectation values at coincident conformal times are: where P B (k, τ ) and P E (k, τ ) denote respectively the hypermagnetic and the hyperelectric power spectra 6 that shall now be derived without relying on the nature of the transition but only on the overall continuity and differentiability of the evolution. For this purpose Eq. (7) can be transformed into a pair of integral equations with initial conditions assigned at τ ex : Depending on the evolution of the background (either before or after τ i ) the condition k 2 = |χ /χ| defines either τ ex or τ re ; the latter condition can also be dubbed, after simple algebra, as: where the overdot denotes a derivation with respect to the cosmic time coordinate; moreover H =ȧ/a = H/a is the Hubble rate while F =χ/χ = F/a is the rate of variation of χ. In Eq. (12) F is the analog of the conventional slow-roll parameter (i.e. H = −Ḣ/H 2 ). Inside the effective horizon (i.e. k/F 1) the initial conditions for τ ≤ −τ ex appearing in Eqs. (10) and (11) Defining, as usual, 3 is the transverse generalization of the Dirac delta function. 6 As in the case of Eq. (2), when F → −F (i.e. χ → 1/χ), the equations of Eq. (7) are interchanged provided f k → g k /k and g k → −kf k (see Eq. (7)). Under the same duality transformation Eqs. (8)- (9) imply that P B → P E and vice versa. Again this symmetry [6] is verified provided χ and χ are simultaneously continuous everywhere and, in particular, across the inflationary boundary.
Eq. (7) for k F. The vacuum Cauchy data correspond to b + (k) → 1 and b − (k) → 0. Conversely when b + (k) = 1 and b − (k) = 0 the mode functions for τ < τ ex correspond to an initial state whose average multiplicity is |b − (k)| 2 . The iterative solution of Eqs. (10) and (11) can then be obtained to the wanted order in k 2 τ 2 but the lowest order solution reduces to the evaluation of the following pair of (dual) integrals: which are both defined provided the integrand is (at least) continuous. Because of this property the evolution across the inflationary boundary can be globally described by introducing the following averages of F and of H/F namely Recalling Eq. (14) and that χ and F are continuous everywhere (and in particular across the inflationary boundary), after two integrations by parts the integral I f (τ ex , τ re ) becomes: According to Eq. (15) χ cannot freeze instantaneously to a constant value after the end of inflation, as observed in explicit numerical integrations (see, for instance, the last paper of Ref. [8]). We now observe that I f (τ ex , τ re ) multiplies the second term at the right hand side of Eq. (15). As a consequence the wanted integral (appearing both at the right and at the left hand side of Eq. (15)) can be solely expressed in terms of F and (H/F ) and f k (τ ) is explicitly given by: We shall now parametrize the inflationary evolution of the susceptibility as 7 χ inf (τ ) = marks the end of the inflationary phase). Conversely for τ ≥ −τ i we shall just assume the continuity of χ and F; with these premises we obtain, quite generically 8 that 2 − F + (H/F ) = 0. All in all we can then say that the continuity properties of the transition imply that the hypermagnetic power spectra of Eqs. (8)-(9) are given by: where, for the sake of simplicity, the vacuum initial conditions have been imposed by setting |b + (k)| = 1 and |b − (k)| = 0. Except for specific values of δ (possibly leading to conspiratorial cancellations) Eq. (17) determines the hypermagnetic power spectra up to overall factors of order 1. With similar considerations the electric power spectra can also be derived; furthermore, since under duality χ → 1/χ and F → −F, we will also have that ν → 1 − ν implying, from Eq. (17) and only using the symmetries of the problem, that P E (k, τ re ) ∝ |kτ i | 7−2ν . Even if they will not be directly relevant for the present discussion we explicitly verified that the method described here correctly leads to the electric power spectra implied by the duality symmetry [6]. Indeed, as soon as the gauge modes reenter the effective horizon duality is explicitly broken since the evolution equations will only contain electric (and not magnetic) sources [12]. When the extrinsic curvature, the susceptibility and the effective horizon are (simultaneously and explicitly) continuous across the inflationary boundary, the general derivation leading to Eq. (17) One of the simplest situations compatible with a sudden transition stipulates that the susceptibility approaches exponentially its (constant) asymptotic value; the explicit expressions of χ inf (τ ) and χ rad (τ ) are given, respectively, by: where we defined, for the sake of conciseness, C = [1 − (1 − 2ν)/(2β)] and D = (1 − C).
Equations (18) and (19) imply the continuity of where the shorthand notations z re = β(τ re /τ i + 1) and kτ ex = q(ν) = ν 2 − 1/4 have been adopted. The second turning point is determined by the analog of Eq. (12) (i.e. χ rad /χ rad = k 2 ) implying 11 e −zre = (C/D)k 2 τ 2 i /β 2 . Therefore we will have that τ re /τ i will be given by: From Eqs. (20) and (21) the final expressions of the mode functions are given by: As long as τ ≥ τ re the presence of the conductivity σ breaks the explicit duality symmetry so that the second equation of Eq. (7) will be replaced by g k = −k 2 f k − Fg k − 4πσg k . In explicit numerical integrations σ smoothly increases (see e.g. the last papers of Refs. [6] and [8]) and the equation obeyed by f k will then be: where f k = exp [2π τ σ(τ )dτ ] f k (τ ). The structure of the turning points implied by Eq. (23) is different, namely, (k 2 − 4πσ 2 )τ 2 i = e −zre [β 2 − 4πστ i β]/(C + e −zre ). According to Eq. (23) the turning point is predominantly fixed by the largeness of στ i rather than by the smallness of kτ i : while kτ i is (at most) of order 1 (and it is much smaller than 1 for the galactic scale) we have instead that στ i 1, as already stressed in explicit numerical integrations of the of the power spectra 12 (see, in particular, the last paper of [8]).
The hypermagnetic power spectra for generalized quantum mechanical Cauchy data can therefore be expressed as: where N B depends on the suddenness of the transition (parametrized, in the above example, by the value of β); when β 1, we will have that N B = (4π 2 ) −1 + O(β −1 ). Equations (17) and (24) are clearly compatible in the case b − (k) → 0 and b + (k) → 1. The relation between comoving and physical power spectra (i.e. P B (k, τ ) = a 4 (τ )χ 2 (τ )P b (k, τ )) follows from the relation between the physical and comoving field operators 13 , i.e. respectively B and b (see Eq. (2) and discussion therein). In the spirit of the present discussion it is therefore interesting to compute the explicit relation between P b (k, τ re ) and P b (k, τ i ) that is given by: where K(ν, β) = 2β/(2β − 1 + 2ν) is a numerical factor of order 1; in Eq. (25) we also used that (τ re /τ i ) = −(2/β) ln kτ i + (3/β) ln β as it follows from Eq. (21) for β ≥ 1. Since the value of where A R is the amplitude of the scalar power spectrum at the pivot scale k p = 0.002 Mpc −1 ; the latter scale also defines the tensor to scalar ratio Depending on the total number of efolds N t , the initial state can influence the late-time spectrum. In this respect the critical number of efolds is given by N c and it is defined as e Nc = (2π Ω R0 A R r T ) 1/4 (M P /H 0 )/4 where Ω R0 is the present energy density of radiation in critical units, H −1 0 is the Hubble radius today and r T is the conventional tensor to scalar ratio 15 . If N t = N c the inflationary event horizon (redshifted at the present epoch) coincides with the Hubble radius today: 13 This relation has nothing to do with the approximate flux conservation during a radiation-dominated stage of expansion, as it is erroneously stated by some. This relation stems directly from the properties of the canonical normal modes of the action (see also Eqs. (2) and (3)).
14 The natural logarithm of x will be denoted hereunder by ln x; common logarithms will be instead denoted by log x. 15 The tensor to scalar ratio itself is affected by the evolution of the large-scale gauge fields (see last paper of [9]). The latter observation implies that r T is bounded from below by the dominance of the adiabatic contribution and it cannot be smaller than 10 −3 , at least in the case of single-field inflationary models. We should therefore assume that 0.001 < r T < 0.1.
When N t > N c the redshifted value of the inflationary event horizon is larger than the present value of the Hubble radius and for N t N c we plausibly expect (at least in conventional inflationary models, that any finite portion of the Universe gradually loses the memory of an initially imposed anisotropy or inhomogeneity so that the Universe attains the observed regularity regardless of the initial boundary conditions). To investigate the role of the initial  Figure 1: The parameter space is described in the (ν, ζ) plane for two different values of the number of efolds and for two complementary comoving scales. The labels appearing on the various isospectral curves denote te common logarithm of P B (k, τ c ) expressed in Gauss. data we shall now define the protoinflationary boundary τ * as the approximate moment at which the background starts inflating. We can then assume, without loss of generality, that |b − (k)| 2 = (k/k * ) ζ where k * τ −1 * . The scale k * is, by definition, the maximal wavenumber of the initial spectrum. The energy density of the initial state at τ * is approximately dρ * 2|b − (k)| 2 d 3 k/(2π) 3 . If ζ > −4 the total energy density is is dominated by the largest scale (i.e. ρ * = O(k 4 * )). To avoid an excessive contribution of the initial energy density to the dynamics of the background we must require 8πρ * where Q i = H * a * / √ H i M P < 1 if the energetic content of the initial state. In practice we can choose, for instance, Q i = O(10 −4 ) as a reasonable fiducial value. Note that when the initial state is thermal |b − (k)| 2 → 1/(e k/k T − 1) (as implied by the Bose-Einstein occupation number) where k T coincides, in the present units, with the putative comoving temperature of the initial state. This case will not be explicitly discussed here 16 but it has been carefully scrutinized in Ref. [12].
The two-point function if the magnetic field intensity at coincident times is derived from Eq. (8) and it is implying 17 that P B (k, τ c ) and P B (k, τ c ) measure, respectively, the energy density the field intensity over the typical wavenumber 1/r. For the applications to magnetogenesis problems [10,11] it is useful to compute the physical power spectrum expressed in units of Gauss (G, in what follows) and at the time of the gravitational collapse of the protogalaxy. The result of this calculation is where, as already mentioned, P B (k, τ ) = a 2 (τ )χ(τ ) P b (k, τ ) for τ < τ re and P B (k, τ ) = a 2 (τ ) P b (k, τ ) for τ > τ re ; the function M(ν, ζ, k, k * ) appearing in Eq. (30) is given by In Eq. (30) we used that the non-screened vector modes of the hypercharge field the project on the electromagnetic fields through the cosine of the Weinberg angle cos θ W .
Equations (30) and (31) shall now be analyzed in the (ν, ζ) plane illustrated in Fig. 1 where we plot the isospectral lines (i.e. the lines of the parameter space over which the magnetic power spectrum is approximately constant) for few physically meaningful choices of the parameters. When ζ → 0 the scale-invariant limit of the magnetic power spectrum corresponds to ν → 5/2 = 2.5 with typical amplitude O(10 −10 )G: this limit can be verified in both plots of Fig. 1. When ζ = 0, quasi-flat power spectra can be obtained as long as ν remains in the region where log [ P B (k, τ c )] = O(−10). In the left plot of Fig. 1 we illustrated the benchmark scale of the gravitational collapse of the protogalaxy (i.e. k = O(1) Mpc −1 ) when the number of efolds is just critical (i.e. N t = N c ); in the right plot we instead considered a larger typical length-scale (i.e. smaller wavenumber k = O(10 2 ) Mpc −1 ) and a larger total number of efolds (i.e. N t = 70 > N c ). What matters is the inclination of the (almost straight) isospectral lines for ζ > 0. Figure 1 shows that for N t > N c the inclination diminishes and the lines become more and more vertical when N t N c . The latter observation shows that as long as N t = O(N c ) the parameter space of inflationary magnetogenesis is comparatively larger than in the case ζ = 0: we pass from a point (i.e. ζ = 0 and ν 2.5) to a whole isospectral line in the (ν, ζ) plane.
17 Equations (5) and (8)  If P b (k c , τ c ) is approximately larger than about 10 −24 G (but still much smaller than O(10 −10 ) G) the observed galactic field intensity can only be reached if the fields are amplified (for τ > τ re ) by the combined action of the gravitational collapse and of the galactic rotation. The latter effect may hopefully transform, under various conditions, the kinetic energy of the plasma into magnetic energy [11]. The most optimistic estimates for the required initial conditions are derived by assuming that every rotation of the galaxy would increase the magnetic field of one efold. The number of galactic rotations since the collapse of the protogalaxy can be between 30 and 35, leading approximately to a purported growth of 13 orders of magnitude. If the dynamo action is totally absent, the required field should be O(10 −11 ) G. In this case, during the collapse of the protogalaxy, the magnetic field will increase by about 5 orders of magnitude. In the literature it is sometimes practical to refer to some hypothetical seed field supposedly present at the time of the collapse of the protogalaxy. By definition B seed = P b (k c , τ c ) where following the standard conventions [11] we took k c = O(1) Mpc −1 . From the above considerations we have therefore that O(10 −24 ) G ≤ B seed (τ c ) ≤ O(10 −12 )G.
After the gauge modes reenter the effective horizon the approximate flux conservation implies B seed (τ c ) = B seed (τ re )(H 0 /M P ) Ω R0 /[π H A R ] and Eq. (32), at τ re can be written, up to the insignificant logarithmic corrections discussed above, as: O(10 33 ) G ≤ B seed (τ re ) ≤ O(10 45 )G.
Equations (32) and (33) are insensitive to the properties of the initial state but they depend on the postinflationary thermal history, as already discussed in the past [12]. In summary the junction conditions for the gauge fields are compatible with sudden and delayed transitions of the effective horizon. The dynamical evolution of the gauge modes has been rephrased in terms of a pair of integral equations related by duality transformations. After showing how the continuity of the susceptibility and of its first derivative determines the hypermagnetic and hyperelectric power spectra, explicit examples of smooth transitions have been proposed to corroborate the analytic discussion. The general arguments based on the continuity of the effective horizon are valid up to logarithmic corrections which are numerically not significant when the gauge modes reenter the effective horizon. Moreover, after reentry these corrections are anyway overwhelmed by the dominance of the conductivity. As long as the total duration of the inflationary phase is nearly minimal the spectral slopes may be directly affected by the properties of the initial state. In the latter case case the parameter space of quasi-flat spectra gets larger. Conversely, when the number of efolds increases beyond a certain critical value, the present findings reproduce the previous results since the effects of the initial state are exponentially suppressed.