Scalar-Tensor-Vector modified gravity in light of the Planck 2018 data

The recent data release by the Planck satellite collaboration presents a renewed challenge for modified theories of gravitation. Such theories must be capable of reproducing the observed angular power spectrum of the cosmic microwave background radiation. For modified theories of gravity, an added challenge lies with the fact that standard computational tools do not readily accommodate the features of a theory with a variable gravitational coupling coefficient. An alternative is to use less accurate but more easily modifiable semianalytical approximations to reproduce at least the qualitative features of the angular power spectrum. We extend a calculation that was used previously to demonstrate compatibility between the Scalar-Tensor-Vector-Gravity (STVG) theory, also known by the acronym MOG, and data from the Wilkinson Microwave Anisotropy Probe (WMAP) to show consistency between the theory and the newly released Planck 2018 data. We find that within the limits of this approximation, the theory accurately reproduces the features of the angular power spectrum.


I. INTRODUCTION
Though highly isotropic, the cosmic microwave background (CMB) shows small temperature fluctuations as a function of sky direction. The magnitude of these fluctuations depends on angular size. This location and size of these peaks is an important prediction of the standard model of cosmology, which has been confirmed by increasingly accurate experiments, such as the Boomerang experiment [1], the Wilkinson Microwave Anisotropy Probe (WMAP, [2]), and the Planck satellite [3].
The angular power spectrum of the CMB can be calculated in a variety of ways. The preferred method is to use numerical software, such as CMBFAST [4]. Unfortunately, such software packages cannot easily be adapted for use with a variable-G theory of gravitation, such as Scalar-Tensor-Vector-Gravity (STVG [5]) theory, also known as MOdified Gravity (MOG).
There are alternative methods of calculation, which, though somewhat less accurate, nonetheless capture the essential qualitative features of the CMB angular power spectrum. The advantage of such calculations is that the physics is transparent, not obscured by "black box" computer code, and that the calculations can be adapted with relative ease to accommodate a different theory. One such method is the semianalytical approximation presented in Ref. [6].
We have, in fact, used this approximation in the past showing agreement between the predictions of the MOG theory and the WMAP results [7,8]. In light of the recently published Planck 2018 results, it is important to revisit and refine this computation, and also extend it to high values of the multipole index ℓ.
We begin in Sec. II with a brief introduction to the MOG theory and its acceleration law, which gives rise to the theory's effective gravitational coupling parameter. In Sec. III we introduce the angular power spectrum and its semianalytical approximation, as presented in Ref. [6] (which the interested reader is advised to consult for details.) We adapt the calculation to the MOG theory and show results. We conclude by presenting our discussion and conclusions in IV.

II. THE MOG THEORY
Our MOG modified gravity theory, also known as Scalar-Tensor-Vector-Gravity (STVG [5]), is a relativistic theory of gravitation based on an action principle. In addition to the metrical field of gravitation, the theory introduces a repulsive vector field of finite range. The gravitational constant and the vector field range (mass) parameter are promoted to dynamical (massless) scalar fields. Within the range of the vector field, the theory replicates Newtonian gravitation; outside this range, in the absence of the repulsive force, gravitation is stronger. By this feature, the theory successfully accounts for galaxy rotation curves [9][10][11][12][13], the matter power spectrum [8], while also remaining consistent with recent gravitational wave data [14].
In the weak field, low-velocity regime, the MOG theory yields a simple gravitational acceleration law [15]. For a point source of gravitation characterized by mass M at the origin, the gravitational acceleration at position r is given byr with where G N is Newton's constant of gravitation, the dimensionless quantity α = (G−G N )/G N (i.e., G = (1+α)G N ) characterizes the difference between the theory's variable gravitational coupling coefficient G and G N , and µ is the mass of the vector field.
In an approximately homogeneous and isotropic universe, α and µ can be taken as constants. Consequently, at distance scales characterized by µr ≫ 1, G eff ∼ G N (1 + α) can be treated as constant as well.
The Friedmann equations that describe a homogeneous and isotropic universe remain valid in the MOG theory [7,8,16,17], with only trivial modifications. Not surprising, given that these equations can also be heuristically derived from the Newtonian theory [6,18]. The equations read (using c = 1): The critical density, characterized by k = 0, Λ = 0, is given by where H =ȧ/a. Note the presence of a factor of 1/(1 + α) in this definition of ρ crit . Consequently, for a given baryon density ρ b , the corresponding density parameter Ω b is inflated by this same factor: Often, in cosmological calculations, Ω b is used to represent the baryon density in equations describing both gravitational and nongravitational interactions. Clearly, this convenience is lost in the case of modified gravity, in the presence of the (1 + α) factor, which applies only to gravitational interactions.

III. MODELING THE COSMIC MICROWAVE BACKGROUND (CMB)
It is surprisingly difficult to analyze high quality Cosmic Microwave Background (CMB) data sets from the perspective of a modified gravity theory, such as MOG. The main reason for this difficulty lies with the fact that, as we alluded to above, the dimensionless density parameter Ω b is used to represent baryonic matter in calculations that involve gravity as well as calculations that represent nongravitational physics.
Why is this a problem? Consider the definition: In the standard theory, this expression will suffice. But what about a theory, such as MOG, with a variable gravitational coefficient G = G eff = (1+α)G N ? Clearly, when the context is gravitational, the product Gρ b accurately reflects the gravitational contribution of baryonic matter.
But when, e.g., the pressure of the medium is considered, Ω b is not supposed to be scaled in this manner (pressure does not increase just because gravitation is stronger). Disentangling these issues in computer codes that have been in use for years or decades, written or rewritten by multiple authors, perhaps even machine translated from one programming language to another (e.g., from FOR-TRAN to C) is a daunting task.
Without access to a standard suite of computer programs that can reliably and provably deal with a variable-G modified theory of gravity, we opted for another approach: use a semianalytical approximation that is sufficiently accurate to reproduce the key qualitative features of the CMB angular power spectrum, and perhaps even allow us to make some cautious predictions.
Such an approximation method was published by Mukhanov [6]. We previously used this approximation method in the context of WMAP results, showing that MOG indeed fits the angular power spectrum well. In light of the recent release of Planck 2018 data, we found it imperative to revisit and, if necessary, refine this calculation and compare Planck results against MOG predictions.

A. Semi-analytical estimation of CMB anisotropies
The general expression for the cosmic mean of the CMB temperature autocorrelation function, expressed in terms of multipoles C ℓ (with the monopole and dipole components, ℓ = 0, 1, excluded), can be written as (see Eq. (9.38) in [6] and the discussion therein for details): where Φ k is the Fourier-decomposition of the gravitational potential Φ with respect to wavenumber k, η r is the conformal time at recombination and η 0 corresponds to the present time. The quantity δ is the fractional energy fluctuation of radiation, defined using the 00-component of the radiation energy-momentum tensor before recombination as T 0 0 = ǫ(1 + δ) where ǫ is the radiation energy density; δ ′ is the derivative with respect to conformal time and and j ℓ are the spherical Bessel functions.
The non-oscillatory part, in turn, is split into a sum: where The parameters that occur in these expressions are as follows. First, the baryon density parameter: where where Ω m ≃ 0.3 is the total matter content (baryonic matter, neutrinos, and cold dark matter). The freestreaming and Silk damping scales are determined, respectively, by Lastly, the location of the acoustic peaks is determined by the parameter 1 .
Finally, we note that the calculated result for C ℓ assumes scale invariance. For small deviations from scale invariance characterized as usual by the parameter n s (with |n s − 1| ≪ 1), the result is scaled: The quality of this approximation is demonstrated in Fig. 1 (top left), which shows the estimated angular power spectrum using nominal parameters (H 0 ∼ 67.4 km/s/Mpc, h 2 Ω b = 0.0224, Ω m = 0.315, n s = 0.965 with spatially flat cosmology, Ω Λ = 1 − Ω m ), against Planck 2018 data 2 from http://pla.esac.esa.int/pla/#cosmology.
The quality of this fit improves significantly if we allow some of the parameters to vary. For instance, using a simple least squares fit, we obtain h 2 Ω b = 0.0187, n s = 0.965 (see Fig. 1 top right).

B. The MOG CMB spectrum
What are the key differences between the MOG theory and standard cosmology?
1 Note the slight changes in the coefficients in (22) and (24) compared to the value published in [6]. We used best fit values for these coefficients from Mukhanov's approximation using the Planck collaboration's best estimates for the parameters of the standard ΛCDM cosmology. In the standard ΛCDM model in the early universe, there are two main sources of gravitation: baryonic matter and collisionless cold dark matter (CDM). The distribution of matter in the universe is still largely homogeneous, and the gravitational field is determined by the sum Ω m = Ω b + Ω DM .
In the MOG theory, Ω DM is, of course, absent. However, the gravitational coupling parameter is no longer Newton's constant. In the late time universe, we expect the gravitational coupling parameter to vary from region to region (an essential feature of the MOG theory that accounts for its ability to model phenomena such as galaxy rotation curves successfully.) In the early, mostly homogeneous universe, we expect little variation in the value of G; however, G = G N .
This means that gravitational interactions are "enhanced" by the factor 1 + α, defined by the relationship G = G eff = (1 + α)G N . When computing results such as the angular power spectrum, this must be taken into account.
This actually leads to a fairly simple prescription. In the formulation presented in the previous subsection, the density parameter for matter, Ω m , must be replaced by These changes are, of course, trivial. Ω b appears only in Eq. (20). For otherwise identical parameterization, we expect identical results.
Instead, we opted to relax the parameter space further as we investigate the MOG solution. Fig. 1, bottom left, was obtained by fitting the values of h 2 Ω b = 0.0197, α = 5.27 and n s = 0.951.
As we explored the parameter space, it became evident that there is significant degeneracy with respect to the value of H 0 . Fig. 1 (lower right) shows another fit, after setting H 0 = 73 km/s/Mpc, resulting in h 2 Ω b = 0.0199, α = 4.75 and n s = 0.949. We believe that this degeneracy demonstrates the limit of the Mukhanov approximation.

IV. CONCLUSIONS
Recently, the Planck collaboration released a data set characterizing the cosmic microwave background's angu-lar power spectrum in more detail than anything previously published. This data release raises the bar for modified theories of gravitation that compete with the standard ΛCDM model as potentially viable representations of the evolution of, and structure formation in the universe.
We investigated in particular the behavior of Scalar-Tensor-Vector-Gravity, also known by the acronym MOG, in the light of these new data. A key feature of the MOG theory is the presence of a variable gravitational coupling coefficient, which makes the task of adapting existing numerical models of the CMB or structure formation difficult. Large numerical code bases that are opaque and often use the dimensionless density parameters Ω b , Ω m , etc., to model both gravitational and nongravitational interactions, cannot be easily modified.
Instead, in this paper, we revived a model that we first employed in the wake of the WMAP data release. Extending the calculations to higher multipoles (up to ℓ = 2500) we were able to demonstrate that the MOG theory correctly reproduces the qualitative features of the CMB, and that within the limits of the approximation, also produces good quantitative fits. At the same time, we also saw the limitations of the method, notable among them a degeneracy with respect to H 0 . This leads us to conclude that, for instance, to decide whether or not the MOG theory can offer a better resolution to the Hubble tension (see [19] for an up-to-date review), more sophisticated methods will be required.