A covariant modified gravity (MOG) is formulated by adding to general relativity two new degrees of freedom, a scalar field gravitational coupling strength $G= 1/\chi$ and a gravitational spin 1 vector field $\phi_\mu$. The $G$ is written as $G=G_N(1+\alpha)$ where $G_N$ is Newton's constant, and the gravitational source charge for the vector field is $Q_g=\sqrt{\alpha G_N}M$, where $M$ is the mass of a body. Cosmological solutions of the theory are derived in a homogeneous and isotropic cosmology. Black holes in MOG are stationary as the end product of gravitational collapse and are axisymmetric solutions with spherical topology. It is shown that the scalar field $\chi$ is constant everywhere for an isolated black hole with asymptotic flat boundary condition. A consequence of this is that the scalar field loses its monopole moment radiation.

We present calculations for the anticipated shadow sizes of Sgr A* and the supermassive black hole in the galaxy M87 in the context of the MOG modified theory of gravity (also known as Scalar-Tensor-Vector-Gravity, or STVG). We demonstrate that mass estimates derived from stellar and gas dynamics in the vicinity of these black holes are the Newtonian masses of the black holes even in the MOG theory. Consequently, shadow sizes increase as a function of the key dimensionless MOG parameter alpha that characterizes the variable gravitational coupling coefficient G, and may offer an observational means to distinguish the MOG theory from standard general relativity.

A modified gravitational theory explains early universe and late time cosmology, galaxy and galaxy cluster dynamics. The modified gravity (MOG) theory extends general relativity (GR) by three extra degrees of freedom: a scalar field $G$, enhancing the strength of the Newtonian gravitational constant $G_N$, a gravitational, spin 1 vector graviton field $\phi_\mu$, and the effective mass $\mu$ of the ultralight spin 1 graviton. For $t < t_{\rm rec}$, where $t_{\rm rec}$ denotes the time of recombination and re-ionization, the density of the vector graviton $\rho_\phi > \rho_b$, where $\rho_b$ is the density of baryons, while for $t > t_{\rm rec}$ we have $\rho_b > \rho_\phi$. The matter density is parameterized by $\Omega_M=\Omega_b+\Omega_\phi+\Omega_r$ where $\Omega_r=\Omega_\gamma+\Omega_\nu$. For the cosmological parameter values obtained by the Planck Collaboration, the CMB acoustical oscillation power spectrum, polarization and lensing data can be fitted as in the $\Lambda$CDM model. When the baryon density $\rho_b$ dominates the late time universe, MOG explains galaxy rotation curves, the dynamics of galaxy clusters, galaxy lensing and the galaxy clusters matter power spectrum without dominant dark matter.

The equation of motion in the generally covariant modified gravity (MOG) theory leads, for weak gravitational fields and non-relativistic motion, to a modification of Newton's gravitational acceleration law. In addition to the metric $g_{\mu\nu}$, MOG has a vector field $\phi_\mu$ that couples with gravitational strength to all baryonic matter. The gravitational coupling strength is determined by the MOG parameter $\alpha$, while parameter $\mu$ is the small effective mass of $\phi_\mu$. The MOG acceleration law has been demonstrated to fit a wide range of galaxies, galaxy clusters and the Bullet Cluster and Train Wreck Cluster mergers. For the SPARC sample of rotationally supported spiral and irregular galaxies, McGaugh et al. [24] (MLS) have found a radial acceleration relation (RAR) that relates accelerations derived from galaxy rotation curves to Newtonian accelerations derived from galaxy mass models. Using the same SPARC galaxy data, mass models independently derived from that data, and MOG parameters $\alpha$ and $\mu$ that run with galaxy mass, we demonstrate that adjusting galaxy parameters within $\pm 1$-sigma bounds can yield MOG predictions consistent with the given rotational velocity data. Moreover, the same adjusted parameters yield a good fit to the RAR of MLS, with the RAR parameter $a_0=(5.4\pm .3)\times 10^{-11}\,{\rm m/s^2}$.

We consider the case of rotating black holes in a dark-matter-emulating theory of gravity called MOG. The latter introduces a gravitational vector field with an associated gravitational charge proportional to the black hole mass and a scalar field in place of the gravitational constant. The resulting black hole metrics resemble the Kerr-Newman geometry and enjoy superradiant scattering. MOG, however, presents important new features. By studying the scattering of a scalar field, we show that there is a marked reduction of the critical frequency of mode amplification. This corresponds to saying that the superradiance peak frequency is red shifted. Analyses of the reflected energy flux also show that MOG black holes are fainter with respect to the standard ones. The proposed results pave the way for testing MOG against astronomical observations.

We investigate the solution of Maxwell's equations in curved spacetime within the framework of Modified Gravity (MOG). We show that besides the null-geodesic treatment of photons in MOG, using Maxwell's equations and covariant coupling with the extra vector sector of gravitation in MOG, we can extract the equation for the propagation of light. We obtain Fermat's potential and calculate the deflection angle of light during lensing from a point-like star. Our results show that the deflection angle obtained from the solution of the wave equation in MOG for the large-scale structures with larger impact parameter of light rays is proportional to that of General Relativity (GR). For solar mass stars the deflection angle agrees with the prediction of GR. However, for the compact structures like the supermassive black hole Sagittarius A* at the centre of the Milky Way, the prediction for the deflection angle is larger than GR, which can be tested in future observations.

We model the velocity dispersion of the ultra-diffuse galaxy NGC 1052-DF2 using Newtonian gravity and modified gravity (MOG). The velocity dispersion predicted by MOG is higher than the Newtonian gravity prediction, but it is fully consistent with the observed velocity dispersion that is obtained from the motion of 10 globular clusters.

We describe an independent analysis of LIGO data for black hole coalescence events. Gravitational wave strain waveforms are extracted directly from the data using a filtering method that exploits the observed or expected time-dependent frequency content. Statistical analysis of residual noise, after filtering out spectral peaks (and considering finite bandwidth), shows no evidence of non-Gaussian behaviour. There is also no evidence of anomalous causal correlation between noise signals at the Hanford and Livingston sites. The extracted waveforms are consistent with black hole coalescence template waveforms provided by LIGO. Simulated events, with known signals injected into real noise, are used to determine uncertainties due to residual noise and demonstrate that our results are unbiased. Conceptual and numerical differences between our RMS signal-to-noise ratios (SNRs) and the published matched-filter detection SNRs are discussed.

A regular rotating MOG compact object is derived that reduces to the Kerr black hole when the parameter $\alpha=0$. Physical consequences of the dark compact object, which is regular everywhere in spacetime for $\alpha > \alpha_{\rm crit}=0.674$ and is a rotating Kerr-MOG black hole for $\alpha < \alpha_{\rm crit}$ are investigated.

We investigate gravitational lensing in the context of the MOG modified theory of gravity. Using a formulation of the theory with no adjustable or fitted parameters, we present the MOG equations of motion for slow, nonrelativistic test particles and for ultrarelativistic test particles, such as rays of light. We demonstrate how the MOG prediction for the bending of light can be applied to astronomical observations. Our investigation first focuses on a small set of strong lensing observations where the properties of the lensing objects are found to be consistent with the predictions of the theory. We also present an analysis of the colliding clusters 1E0657-558 (known also as the Bullet Cluster) and Abell 520; in both cases, the predictions of the MOG theory are in good agreement with observation.

A major hurdle for modified gravity theories is to explain the dynamics of galaxy clusters. A case is made for a generalized gravitational theory called Scalar-Tensor-Vector-Gravity (STVG) or MOG (modified gravity) to explain merging cluster dynamics. The paper presents the results of a re-analysis of the Bullet Cluster, as well as an analysis of the Train Wreck Cluster in the weak gravitational field limit without dark matter. The King-$\beta$ model is used to fit the X-ray data of both clusters, and the $\kappa$-maps are computed using the parameters of this fit. The amount of galaxies in the clusters is estimated by subtracting the predicted $\kappa$-map from the $\kappa$-map data. The estimate for the Bullet Cluster is that $14.1\%$ of the cluster is composed of galaxies. For the Train Wreck Cluster, if the Jee et al. data are used, $25.7\%$ of the cluster is composed of galaxies. The baryon matter in the galaxies and the enhanced strength of gravitation in MOG, shift the lensing peaks making them offset from the gas. The work demonstrates that this generalized gravitational theory can explain merging cluster dynamics without dark matter.

Modified gravity (MOG) is a covariant, relativistic, alternative gravitational theory whose field equations are derived from an action that supplements the spacetime metric tensor with vector and scalar fields. Both gravitational (spin 2) and electromagnetic waves travel on null geodesics of the theory's one metric. Despite a recent claim to the contrary, MOG satisfies the weak equivalence principle and is consistent with observations of the neutron star merger and gamma ray burster event GW170817/GRB170817A.

A promising signature of coalescing black holes is their spin angular distribution. We consider the aLIGO collaboration gravitational wave measurements of the binary black hole spins and the predicted modified gravity (MOG) preference for misaligned spins of the coalescing black holes. In MOG, during the merger of two black holes, the enhanced strength of gravitation reduces the effective spin parameter $\chi_{\rm eff}\sim 0$ in agreement with the measured spin misalignment of the merging black holes observed in the gravitational wave detections by the aLIGO collaboration.

The weak lensing of cosmic web filaments is investigated in modified gravity (MOG) and it is demonstrated that the detected galaxies and baryonic gas of order $5-15 \%$ in filaments can with the enhanced value of the strength of gravity agree with the lensing data for filaments without dark matter.

The galaxy cluster system Abell 1689 has been well studied and yields good lensing and X-ray gas data. Modified gravity (MOG) is applied to the cluster Abell 1689 and the acceleration data is well fitted without assuming dark matter. Newtonian dynamics and Modified Newtonian dynamics (MOND) are shown not to fit the acceleration data, while a dark matter model based on the Navarro-Frenk-White (NFW) mass profile is shown to fit the acceleration data for the radial range $r > 200$ kpc.

The nature of gravitational waves in a generalized gravitation theory is investigated. The linearized field equations and the metric tensor quadrupole moment power and the decrease in radius of an inspiralling binary system of two compact objects are derived. The generalized Kerr metric describing a spinning black hole is determined by its mass $M$ and the spin parameter $a=cS/GM^2$. The LIGO-Virgo collaboration data is fitted with smaller binary black hole masses in agreement with the current electromagnetic, observed X-ray binary upper bound for a black hole mass, $M\lesssim 10M_\odot$.

The equation of motion in the generally covariant modified gravity (MOG) theory leads for weak gravitational fields and the non-relativistic limit to a modification of the Newtonian gravitational acceleration law, expressed in terms of two parameters $\alpha$ and $\mu$. The parameter $\alpha$ determines the strength of the gravitational field and $\mu$ is the effective mass of the vector field $\phi_\mu$, coupled with gravitational strength to baryonic matter. The MOG acceleration law for weak field gravitation and non-relativistic particles has been demonstrated to fit a wide range of galaxies, galaxy clusters and the Bullet Cluster and Train Wreck Cluster mergers. We demonstrate that the MOG acceleration law for a point mass source is in agreement with the McGaugh et al., correlation between the radial acceleration traced by galaxy rotation curves and the distribution of baryonic matter for the SPARC sample of 153 rotationally supported spiral and irregular galaxies.

An alternative to the postulate of dark energy required to explain the accelerated expansion of the universe is to adopt an inhomogeneous cosmological model to explain the supernovae data without dark energy. We adopt a void cosmology model, based on the inhomogeneous Lema\^{i}tre-Tolman-Bondi solution of Einstein's field equations. The model can resolve observational anomalies in the $\Lambda CDM$ model, such as the discrepancy between the locally measured value of the Hubble constant, $H_0=73.24\pm 1.74\,{\rm km}\,{\rm s}^{-1}\,{\rm Mpc}^{-1}$, and the $H_0=66.93\pm 0.62\,{\rm km}\,{\rm s}^{-1}\,{\rm Mpc}^{-1}$ determined by the Planck satellite data and the $\Lambda CDM$ model, and the lithium $^{7}{\rm Li}$ problem, which is a $5\sigma$ mismatch between the theoretical prediction for the $^{7}{\rm Li}$ from big bang nucleosynthesis and the value that we observe locally today at $z=0$. The void model can also resolve the tension between the number of massive clusters derived from the Sunyaev-Zel'dovich effect by the Planck satellite and the number expected from the CMB anisotropies, and the CMB weak lensing anomaly observed in the Planck data. The cosmological Copernican principle and the time and position today coincidence conundrums in the $\Lambda CDM$ and void cosmological models are discussed.

A variable speed of light (VSL) cosmology is described in which the causal mechanism of generating primordial perturbations is achieved by varying the speed of light in a primordial epoch. This yields an alternative to inflation for explaining the formation of the cosmic microwave background (CMB) and the large scale structure (LSS) of the universe. The initial value horizon and flatness problems in cosmology are solved. The model predicts primordial scalar and tensor fluctuation spectral indices $n_s=0.96$ and $n_t=- 0.04$, respectively. We make use of the $\delta{\cal N}$ formalism to identify signatures of primordial nonlinear fluctuations, and this allows the VSL model to be distinguished from inflationary models. In particular, we find that the parameter $f_{\rm NL}=5$ in the variable speed of light cosmology. The value of the parameter $g_{\rm NL}$ evolves during the primordial era and shows a running behavior.

Deviations from the standard $\Lambda$CDM model motivate an interpretation of early universe cosmology using the Scalar-Tensor-Vector-Gravity (STVG) theory. A constraint analysis carried out by Valentino, Melchiorri and Silk, revealed deviations from the growth of structure predicted by General Relativity, and a lensing anomaly in the angular CMB power spectrum data. The modified gravity (MOG) theory resolves the lensing deviation from the standard model and provides an explanation of the CMB and structure growth data.

A variable speed of light (VSL) cosmology is described in which the causal mechanism of generating primordial perturbations is achieved by varying the speed of light in a primordial epoch. This yields an alternative to inflation for explaining the formation of the cosmic microwave background (CMB) and the large scale structure (LSS) of the universe. We make use of the $\delta{\cal N}$ formalism to identify signatures of primordial nonlinear fluctuations, and this allows the VSL model to be distinguished from inflationary models. In particular, we find that the parameter $f_{\rm NL}=5$ in the variable speed of light cosmology. The value of the parameter $g_{\rm NL}$ evolves during the primordial era and shows a running behavior.

A modified gravitational theory is developed in which the gravitational coupling constants $G$ and $Q$ and the effective mass $m_\phi$ of a repulsive vector field run with momentum scale $k$ or length scale $\ell =1/k$, according to a renormalization group flow. The theory can explain cosmological early universe data with a dark hidden photon and late time galaxy and cluster dynamics without dark matter. The theory agrees with solar system and binary pulsar observations.

The field equations for Scalar-Tensor-Vector-Gravity (STVG) or modified gravity (MOG) have a static, spherically symmetric black hole solution determined by the mass $M$ with two horizons. The strength of the gravitational constant is $G=G_N(1+\alpha)$ where $\alpha$ is a parameter. A regular singularity-free MOG solution is derived using a nonlinear field dynamics for the repulsive gravitational field component and a reasonable physical energy-momentum tensor. The Kruskal-Szekeres completion of the MOG black hole solution is obtained. The Kerr-MOG black hole solution is determined by the mass $M$, the parameter $\alpha$ and the spin angular momentum $J=Ma$. The equations of motion and the stability condition of a test particle orbiting the MOG black hole are derived, and the radius of the black hole photosphere and the shadows cast by the Schwarzschild-MOG and Kerr-MOG black holes are calculated. A traversable wormhole solution is constructed with a throat stabilized by the repulsive component of the gravitational field.

The shadows cast by non-rotating and rotating modified gravity (MOG) black holes are determined by the two parameters mass $M$ and angular momentum $J=Ma$. The sizes of the shadows cast by the spherically symmetric static Schwarzschild-MOG and Kerr-MOG rotating black holes increase significantly as the free parameter $\alpha$ is increased from zero. The Event Horizon Telescope (EHT) shadow image measurements can determine whether Einstein's general relativity is correct or whether it should be modified in the presence of strong gravitational fields.

Galaxy rotation curves determined observationally out to a radius well beyond the galaxy cores can provide a critical test of modified gravity models without dark matter. The predicted rotational velocity curve obtained from Scalar-Vector-Tensor Gravity (STVG or MOG) is in excellent agreement with data for the Milky Way without a dark matter halo, with a mass of $5\times 10^{10}\,M_{\odot}$. The velocity rotation curve predicted by modified Newtonian Dynamics (MOND) does not agree with the data.

A conformal coupling of the metric in the Jordan frame to the energy-momentum tensor, screens the scalar field gravitational coupling strength $G$ in modified gravity (MOG). The scalar field acquires a mass which depends on the local matter density: the scalar field particle is massive for the Sun and earth, where the density is high compared to low density environments in cosmology and astrophysics. Together with the screening of the vector field $\phi_\mu$, this guarantees that solar system tests of gravity are satisfied. The conformal metric is coupled to the electromagnetic matter field and energy-momentum tensor, screening $G$ for the Sun and the deflection of light by the Sun and the Shapiro time delay in MOG are in agreement with general relativity. For galaxies and galactic clusters the enhanced gravitational coupling constant $G$ leads to agreement with gravitational lensing without dark matter. For compact binary pulsars the screening of $G$ removes the monopole and dipole gravitational radiation modes in agreement with the binary pulsar timing data.

An important piece of evidence for dark matter is the need to explain the growth of structure from the time of horizon entry and radiation-matter equality to the formation of stars and galaxies. This cannot be explained by using general relativity without dark matter. So far, dark matter particles have not been detected in laboratory measurements or at the LHC. We demonstrate that enhanced structure growth can happen in a modified gravity theory (MOG). The vector field and particle introduced in the theory to explain galaxy and cluster dynamics plays an important role in generating the required structure growth. The particle called the phion (a light hidden photon) is neutral and is a dominant, pressureless component in the MOG Friedmann equations, before the time of decoupling. The dominant energy density of the phion particle in the early universe, generates an explanation for the growth of density perturbations. The angular acoustical power spectrum due to baryon-photon pressure waves is in agreement with the Planck 2013 data. As the universe expands and large scale structures are formed, the density of baryons dominates and the rotation curves of galaxies and the dynamics of clusters are explained in MOG, when the phion particle in the present universe is ultra-light. The matter power spectrum determined by the theory is in agreement with current galaxy redshift surveys.

The quantum gravity effects of vacuum polarization of gravitons propagating in a curved spacetime cause the quantum vacuum to act as a dispersive medium with a refractive index. Due to this dispersive medium gravitons acquire superluminal velocities. The dispersive medium is produced by higher derivative curvature contributions to the effective gravitational action. It is shown that in a Friedmann-Lema\^{i}tre-Robertson-Walker spacetime in the early universe near the Planck time $t_{\rm PL}\gtrsim 10^{-43}\,{\rm sec}$, the speed of gravitational waves $c_g\gg c_{g0}=c_0$, where $c_{g0}$ and $c_0$ are the speeds of gravitational waves and light today. The large speed of gravitational waves stretches their wavelengths to super-horizon sizes, allowing them to be observed in B-polarization experiments.

We apply the weak field approximation limit of the covariant Scalar-Tensor-Vector Gravity (STVG) theory, so-called MOdified gravity (MOG), to the dynamics of clusters of galaxies by using only baryonic matter. The MOG effective gravitational potential in the weak field approximation is composed of an attractive Newtonian term and a repulsive Yukawa term with two parameters $\alpha$ and $\mu$. The numerical values of these parameters have been obtained by fitting the predicted rotation curves of galaxies to observational data, yielding the best fit result: $\alpha = 8.89 \pm 0.34$ and $\mu= 0.042\pm 0.004$ kpc$^{-1}$~\cite{rah13}. We extend the observational test of this theory to clusters of galaxies, using data for the ionized gas and the temperature profile of nearby clusters obtained by the Chandra X-ray telescope. Using the MOG virial theorem for clusters, we compare the mass profiles of clusters from observation and theory for eleven clusters. The theoretical mass profiles for the inner parts of clusters exceed the observational data. However, the observational data for the inner parts of clusters (i.e., $r<0.1 r_{500}$) is scattered, but at distances larger than $\sim 300$ kpc, the observed and predicted mass profiles converge. Our results indicate that MOG as a theory of modified gravity is compatible with the observational data from the the solar system to Mega parsec scales without invoking dark matter.

As an alternative to dark matter models, MOdified Gravity (MOG) theory can compensate for dark matter by a covariant modification of Einstein gravity. The theory introduces two additional scalar fields and one vector field. The aim is to explain the dynamics of astronomical systems based only on their baryonic matter. The effect of the vector field in the theory resembles a Lorentz force where each mass has a charge proportional to the inertial mass. In this work, we obtain the weak field approximation of MOG by perturbing the metric and the fields around Minkowski space-time. We derive an effective gravitational potential which yields the Newtonian attractive force plus a repulsive Yukawa force. This potential, in addition to the Newtonian gravitational constant, $G_N$, has two additional constant parameters $\alpha$ and $\mu$. We use the THINGS catalog of galaxies and fix the two parameters $\alpha$ and $\mu$ of the theory to be $\alpha =8.89 \pm 0.34$ and $\mu =0.04 \pm 0.004 {\rm kpc}^{-1}$. We then apply the effective potential with the fixed universal parameters to the Ursa-Major catalog of galaxies and obtain good fits to galaxy rotation curve data with an average value of $\bar{\chi^2} = 1.07 $. In the fitting process, only the stellar mass-to-light ratio $(M/L)$ of the galaxies is a free parameter. As predictions of MOG, our derived $M/L$ is shown to be correlated with the color of galaxies, and we fit the Tully-Fisher relation for galaxies. As an alternative to dark matter, introducing an effective weak field potential for MOG opens a new window to the astrophysical applications of the theory.

A bimetric gravity model with a variable speed of light is shown to be in agreement with the results reported from the Planck satellite in 2013. The predicted scalar mode spectral index is $n_s\approx 0.96$ and its running is $\alpha_s\approx 8\times 10^{-4}$ when the fundamental length scale $\ell_0$ in the model is fixed to be $\ell_0\approx 10^5\ell_P$, where $\ell_P$ is the Planck length $\ell_P=1.62\times 10^{-33}\,{\rm cm}$, giving the observed CMB fluctuations: $\delta_H\approx 10^{-5}$. The enlarged lightcone ensures that horizon and flatness problems are solved. The model is free from many of the fine-tuning problems of the inflationary models and the fluctuations that form the seeds of structure formation do not lead to a chaotic inhomogeneous universe and the need for a multiverse.

The Jordan-Brans-Dicke theory of gravitation, which promotes the gravitational constant to a dynamical scalar field, predicts a value for the Eddington-Robertson post-Newtonian parameter gamma that is significantly different from the general relativistic value of unity. This contradicts precision solar system measurements that tightly constrain gamma around 1. We consider a modification of the theory, in which the scalar field is sourced explicitly by matter. We find that this leads to a modified expression for the gamma-parameter. In particular, a specific choice of the scalar current yields gamma=1, just as in general relativity, while the weak equivalence principle is also satisfied. This result has important implications for theories that mimic Jordan-Brans-Dicke theory in the post-Newtonian limit in the solar system, including our scalar-tensor-vector modified gravity theory (MOG).

We explore the cosmological consequences of Modified Gravity (MOG), and find that it provides, using a minimal number of parameters, good fits to data, including CMB temperature anisotropy, galaxy power spectrum, and supernova luminosity-distance observations without exotic dark matter. MOG predicts a bouncing cosmology with a vacuum energy term that yields accelerating expansion and an age of ~13 billion years.

Our modified gravity theory (MOG) is a gravitational theory without exotic dark matter, based on an action principle. MOG has been used successfully to model astrophysical phenomena such as galaxy rotation curves, galaxy cluster masses, and lensing. MOG may also be able to account for cosmological observations. We assume that the MOG point source solution can be used to describe extended distributions of matter via an appropriately modified Poisson equation. We use this result to model perturbation growth in MOG and find that it agrees well with the observed matter power spectrum at present. As the resolution of the power spectrum improves with increasing survey size, however, significant differences emerge between the predictions of MOG and the standard LCDM model, as in the absence of exotic dark matter, oscillations of the power spectrum in MOG are not suppressed. We can also use MOG to model the acoustic power spectrum of the cosmic microwave background. A suitably adapted semi-analytical model offers a first indication that MOG may pass this test, and correctly model the peak of the acoustic spectrum.

We comment on arXiv:1112.1320 and point out that baryonic oscillations of the matter power spectrum, while predicted by theories that do not incorporate collisionless cold dark matter, are strongly suppressed by the statistical window function that is used to process finite-sized galaxy samples. We assert that with present-day data sets, the slope of the matter power spectrum is a much stronger indicator of a theory's validity. We also argue that MOND should not be used as a strawman theory as it is not in general representative of modified gravity theories; some theories, notably our scalar-vector-tensor MOdified Gravity (MOG), offer much more successful predictions of cosmological observations.

Our modified gravity theory (MOG) was used successfully in the past to explain a range of astronomical and cosmological observations, including galaxy rotation curves, the CMB acoustic peaks, and the galaxy mass power spectrum. MOG was also used successfully to explain the unusual features of the Bullet Cluster 1E0657-558 without exotic dark matter. In the present work, we derive the relativistic equations of motion in the spherically symmetric field of a point source in MOG and, in particular, we derive equations for light bending and lensing. Our results also have broader applications in the case of extended distributions of matter, and they can be used to validate the Bullet Cluster results and provide a possible explanation for the merging clusters in Abell 520.

We compare the MOG and MOND predictions of the Tully-Fisher relationship to observations using gas-rich galaxies for which the baryonic mass can be estimated accurately. We find that both theories are in good agreement with the data. Future observations of large, gas-rich galaxies may offer a means to distinguish between these two theories.

The Sloan Digital Sky Survey (SDSS) provides data on several hundred thousand galaxies. Precise location of these galaxies in the sky, along with information about their luminosities and line-of-sight (Doppler) velocities allows one to construct a three-dimensional map of their location and estimate their line-of-sight velocity dispersion. This information, in principle, allows one to test dynamical gravity models, specifically models of satellite galaxy velocity dispersions near massive hosts. A key difficulty is the separation of true satellites from interlopers. We sidestep this problem by not attempting to derive satellite galaxy velocity dispersions from the data, but instead incorporate an interloper background into the mathematical models and compare the result to the actual data. We find that due to the presence of interlopers, it is not possible to exclude several gravitational theories on the basis of the SDSS data.

Modified Gravity (MOG) has been used successfully to explain the rotation curves of galaxies, the motion of galaxy clusters, the Bullet Cluster, and cosmological observations without the use of dark matter or Einstein's cosmological constant. We review the main theoretical ideas and applications of the theory to astrophysical and cosmological data.

In this Letter we derive the gravity field equations by varying the action for an ultraviolet complete quantum gravity. Then we consider the case of a static source term and we determine an exact black hole solution. As a result we find a regular spacetime geometry: in place of the conventional curvature singularity extreme energy fluctuations of the gravitational field at small length scales provide an effective cosmological constant in a region locally described in terms of a de Sitter space. We show that the new metric coincides with the noncommutative geometry inspired Schwarzschild black hole. Indeed, we show that the ultraviolet complete quantum gravity, generated by ordinary matter is the dual theory of ordinary Einstein gravity coupled to a noncommutative smeared matter. In other words we obtain further insights about that quantum gravity mechanism which improves Einstein gravity in the vicinity of curvature singularities. This corroborates all the existing literature in the physics and phenomenology of noncommutative black holes.

We apply our scalar-tensor-vector (STVG) modified gravity theory (MOG) to calculate the infall velocities of the two clusters constituting the Bullet Cluster 1E0657-06. In the absence of an applicable two-body solution to the MOG field equations, we adopt an approximate acceleration formula based on the spherically symmetric, static, vacuum solution of the theory in the presence of a point source. We find that this formula predicts an infall velocity of the two clusters that is consistent with estimates based on hydrodynamic simulations.

Two possible explanations for the type SNe Ia supernovae observations are a nonlinear, underdense void embedded in a matter-dominated Einstein-de Sitter spacetime or dark energy in the $\Lambda$CDM model. Both of these alternatives are faced with Copernican fine-tuning problems. A case is made for the void scenario that avoids introducing undetected dark energy.

MOG is a fully relativistic modified theory of gravity based on an action principle. The MOG field equations are exactly solvable numerically in two important cases. In the spherically symmetric, static case of a gravitating mass, the equations also admit an approximate solution that closely resembles the Reissner-Nordstrom metric. Furthermore, for weak gravitational fields, a Yukawa-type modification to the Newtonian acceleration law can be obtained, which can be used to model a range of astronomical observations. Without nonbaryonic dark matter, MOG provides good agreement with the data for galaxy rotation curves, galaxy cluster masses, and gravitational lensing, while predicting no appreciable deviation from Einstein's predictions on the scale of the solar system. Another solution of the field equations is obtained for the case of a a spatially homogeneous, isotropic cosmology. MOG predicts an accelerating universe without introducing Einstein's cosmological constant; it also predicts a CMB acoustic power spectrum and a mass power spectrum that are consistent with observations without relying on non-baryonic dark matter. Increased sensitivity in future observations or space-based experiments may be sufficient to distinguish MOG from other theories, notably the LCDM "standard model" of cosmology.

Modified gravity theory is known to violate Birkhoff's theorem. We explore a key consequence of this violation, the effect of distant matter in the Universe on the motion of test particles. We find that when a particle is accelerated, a force is experienced that is proportional to the particle's mass and acceleration and acts in the direction opposite to that of the acceleration. We identify this force with inertia. At very low accelerations, our inertial law deviates slightly from that of Newton, yielding a testable prediction that may be verified with relatively simple experiments. Our conclusions apply to all gravity theories that reduce to a Yukawa-like force in the weak field approximation.

Modified Gravity (MOG) has been used successfully to explain the rotation curves of galaxies, the motion of galaxy clusters, the Bullet Cluster, and cosmological observations without the use of dark matter or Einstein's cosmological constant. We now have the ability to demonstrate how these solutions can be obtained directly from the action principle, without resorting to the use of fitted parameters or empirical formulae. We obtain numerical solutions to the theory's field equations that are exact in the sense that no terms are omitted, in two important cases: the spherically symmetric, static vacuum solution and the cosmological case of an homogeneous, isotropic universe. We compare these results to selected astrophysical and cosmological observations.

A modified gravity (MOG) theory, that has been successfully fitted to galaxy rotational velocity data, cluster data, the Bullet Cluster 1E0657-56 and cosmological observations, is shown to be in good agreement with the motion of satellite galaxies around host galaxies at distances 50-400 kpc.

Globular clusters (GCs) in the Milky Way have characteristic velocity dispersions that are consistent with the predictions of Newtonian gravity, and may be at odds with Modified Newtonian Dynamics (MOND). We discuss a modified gravity (MOG) theory that successfully predicts galaxy rotation curves, galaxy cluster masses and velocity dispersions, lensing, and cosmological observations, yet produces predictions consistent with Newtonian theory for smaller systems, such as GCs. MOG produces velocity dispersion predictions for GCs that are independent of the distance from the galactic center, which may not be the case for MOND. New observations of distant GCs may produce strong criteria that can be used to distinguish between competing gravitational theories.

A detailed analysis of the November 15, 2006 data release (Clowe et al., 2006) X-ray surface density Sigma-map and the strong and weak gravitational lensing convergence kappa-map for the Bullet Cluster 1E0657-558 is performed and the results are compared with the predictions of a modified gravity (MOG) and dark matter. Our surface density Sigma-model is computed using a King beta-model density, and a mass profile of the main cluster and an isothermal temperature profile are determined by the MOG. We find that the main cluster thermal profile is nearly isothermal. The MOG prediction of the isothermal temperature of the main cluster is T = 15.5 +- 3.9 keV, in good agreement with the experimental value T = 14.8{+2.0}{-1.7} keV. Excellent fits to the two-dimensional convergence kappa-map data are obtained without non-baryonic dark matter, accounting for the 8-sigma spatial offset between the Sigma-map and the kappa-map reported in Clowe et al. (2006). The MOG prediction for the kappa-map results in two baryonic components distributed across the Bullet Cluster 1E0657-558 with averaged mass-fraction of 83% intracluster medium (ICM) gas and 17% galaxies. Conversely, the Newtonian dark matter kappa-model has on average 76% dark matter (neglecting the indeterminant contribution due to the galaxies) and 24% ICM gas for a baryon to dark matter mass-fraction of 0.32, a statistically significant result when compared to the predicted Lambda-CDM cosmological baryon mass-fraction of 0.176{+0.019}{-0.012} (Spergel et al., 2006).

An exact inhomogeneous solution of Einstein's field equations is shown to be able to inflate in a non-uniform way in the early universe and explain anomalies in the WMAP power spectrum data. It is also possible for the model to explain the accelerated expansion of the universe by late-time inhomogeneous structure.

In a recent note Ellis criticizes varying speed of light theories on the grounds of a number of foundational issues. His reflections provide us with an opportunity to clarify some fundamental matters pertaining to these theories.

A fully relativistic modified gravitational theory including a fifth force skew symmetric field is fitted to the Pioneer 10/11 anomalous acceleration. The theory allows for a variation with distance scales of the gravitational constant G, the fifth force skew symmetric field coupling strength omega and the mass of the skew symmetric field mu=1/lambda. A fit to the available anomalous acceleration data for the Pioneer 10/11 spacecraft is obtained for a phenomenological representation of the "running" constants and values of the associated parameters are shown to exist that are consistent with fifth force experimental bounds. The fit to the acceleration data is consistent with all current satellite, laser ranging and observations for the inner planets.

Gravitational lensing in a modified gravity (MOG) is derived and shown to describe lensing without postulating dark matter. The recent data for merging clusters identified with the interacting cluster 1E0657-56 is shown to be consistent with a weak lensing construction based on MOG without exotic dark matter.

The inhomogeneous distribution of matter in the non-linear regime of galaxies, clusters of galaxies and voids is described by an exact, spherically symmetric inhomogeneous solution of Einstein's gravitational field equations, corresponding to an under-dense void. The solution becomes the homogeneous and isotropic Einstein-de Sitter solution for a red shift $z > 10-20$, which describes the matter dominated CMB data with small inhomogeneities $\delta\rho/\rho\sim 10^{-5}$. A spatial volume averaging of physical quantities is introduced and the averaged time evolution expansion parameter $\theta$ in the Raychoudhuri equation can give rise in the late-time universe to a volume averaged deceleration parameter $<q>$ that is negative for a positive matter density. This allows for a region of accelerated expansion which does not require a negative pressure dark energy or a cosmological constant. A negative deceleration parameter can be derived by this volume averaging procedure from the Lema\^{i}tre-Tolman-Bondi open void solution, which describes the late-time non-linear regime associated with galaxies and under-dense voids and solves the ``coincidence'' problem.

The inhomogeneous distribution of matter in the non-linear regime of galaxies, clusters of galaxies and voids is described by an exact, spherically symmetric inhomogeneous solution of Einstein's gravitational field equations, corresponding to an under-dense void. The solution becomes the homogeneous and isotropic Einstein-de Sitter solution for a red shift $z > 10-20$, which describes the matter dominated CMB data with small inhomogeneities. A spatial volume averaging of physical quantities is introduced and the averaged time evolution expansion parameter $\theta$ can give rise in the late-time universe to a volume averaged deceleration parameter $<q>$ that is negative for a positive matter density. This allows for a region of accelerated expansion which does not require a cosmological constant. A negative deceleration parameter can be derived by this volume averaging procedure from the Lema\^{i}tre-Tolman-Bondi open void solution, which describes the late-time non-linear regime associated with galaxies and under-dense voids and solves the ``coincidence'' problem.

A modified gravity (MOG) possesses a light, neutral vector particle called a ``phion'' associated with a vector field $\phi^\mu$, which forms a cold fluid of Bose-Einstein condensates before recombination with zero pressure and zero shear viscosity. The energy density associated with this Bose-Einstein condensate fluid dominates the energy density before recombination and produces a density parameter, $\Omega_\phi\sim 0.3$, that together with the fractional baryon density $\Omega_b\sim 0.04$, and a cosmological constant parameter $\Omega_\Lambda\sim 0.7$ yields an approximate fit to the data for the acoustical oscillations in the CMB power spectrum. The quantum phion condensate fluid is abundant well before recombination and can clump and form the primordial structure for galaxies. At late times in the expanding universe, in local bound systems such as galaxies ordinary baryonic matter dominates the matter density. For galactic systems in the present epoch, the modified Newtonian acceleration law determined by MOG describes well galaxy rotation curve data and X-ray cluster mass profile data.

A covariant scalar-tensor-vector gravity theory is developed which allows the gravitational constant $G$, a vector field coupling $\omega$ and the vector field mass $\mu$ to vary with space and time. The equations of motion for a test particle lead to a modified gravitational acceleration law that can fit galaxy rotation curves and cluster data without non-baryonic dark matter. The theory is consistent with solar system observational tests. The linear evolutions of the metric, vector field and scalar field perturbations and their consequences for the observations of the cosmic microwave background are investigated.

We consider a cosmology in which a spherically symmetric large scale inhomogeneous enhancement or a void are described by an inhomogeneous metric and Einstein's gravitational equations. For a flat matter dominated universe the inhomogeneous equations lead to luminosity distance and Hubble constant formulas that depend on the location of the observer. For a general inhomogeneous solution, it is possible for the deceleration parameter to differ significantly from the FLRW result. The deceleration parameter $q_0$ can be interpreted as $q_0 > 0$ ($q_0=1/2$ for a flat matter dominated universe) in a FLRW universe and be $q_0 < 0$ as inferred from the inhomogeneous enhancement that is embedded in a FLRW universe. A spatial volume averaging of local regions in the backward light cone has to be performed for the inhomogeneous solution at late times to decide whether the decelerating parameter $q$ can be negative for a positive energy condition. The CMB temperature fluctuations across the sky can be unevenly distributed in the northern and southern hemispheres in the inhomogeneous matter dominated solution, in agreement with the analysis of the WMAP power spectrum data by several authors. The model can possibly explain the anomalous alignment of the quadrupole and octopole moments observed in the WMAP data.

We apply the modified acceleration law obtained from Einstein gravity coupled explaining galaxy rotation curves without exotic dark matter. Our sample of galaxies includes low surface brightness (LSB) and high surface brightness (HSB) galaxies, and an elliptical galaxy. In those cases where photometric data are available, a best fit via the single parameter (M/L)_{stars} to the luminosity of the gaseous (HI plus He) and luminous stellar disks is obtained. Additionally, a best fit to the rotation curves of galaxies is obtained in terms of a parametric mass distribution (independent of luminosity observations) -- a two parameter fit to the total galactic mass, (or mass-to-light ratio M/L), and a core radius associated with a model of the galaxy cores using a nonlinear least-squares fitting routine including estimated errors. The fits are compared to those obtained using Milgrom's phenomenological MOND model and to the predictions of the Newtonian-Kepler acceleration law.

We apply the modified acceleration law obtained from Einstein gravity coupled to a massive skew symmetric field, F_{\mu\nu\lambda}, to the problem of explaining X-ray galaxy cluster masses without exotic dark matter. Utilizing X-ray observations to fit the gas mass profile and temperature profile of the hot intracluster medium (ICM) with King beta-models, we show that the dynamical masses of the galaxy clusters resulting from our modified acceleration law fit the cluster gas masses for our sample of 106 clusters without the need of introducing a non-baryonic dark matter component. We are further able to show for our sample of 106 clusters that the distribution of gas in the ICM as a function of radial distance is well fit by the dynamical mass distribution arising from our modified acceleration law without any additional dark matter component. In previous work, we applied this theory to galaxy rotation curves and demonstrated good fits to our sample of 101 LSB, HSB and dwarf galaxies including 58 galaxies that were fit photometrically with the single parameter (M/L)_{stars}. The results there were qualitatively similar to those obtained using Milgrom's phenomenological MOND model, although the determined galaxy masses were quantitatively different and MOND does not show a return to Keplerian behavior at extragalactic distances. The results here are compared to those obtained using Milgrom's phenomenological MOND model which does not fit the X-ray galaxy cluster masses unless an auxiliary dark matter component is included.

Einstein gravity coupled to a massive skew symmetric field F_{\mu\nu\lambda} leads to an acceleration law that modifies the Newtonian law of attraction between particles. We use a framework of non-perturbative renormalization group equations as well as observational input to characterize special renormalization group trajectories to allow for the running of the effective gravitational coupling G and the coupling of the skew field to matter. The latter lead to an increase of Newton's constant at large galactic and cosmological distances. For weak fields a fit to the flat rotation curves of galaxies is obtained in terms of the mass (mass-to-light ratio M/L) of galaxies. The fits assume that the galaxies are not dominated by exotic dark matter and that the effective gravitational constant G runs with distance scale. The equations of motion for test particles yield predictions for the solar system and the binary pulsar PSR 1913+16 that agree with the observations. The gravitational lensing of clusters of galaxies can be explained without exotic dark matter. An FLRW cosmological model with an effective G=G(t) running with time can lead to consistent fits to cosmological data without assuming the existence of exotic cold dark matter.

We describe the universe as a local, inhomogeneous spherical bubble embedded in a flat matter dominated FLRW universe. Generalized exact Friedmann equations describe the expansion of the universe and an early universe inflationary de Sitter solution is obtained. A non-perturbative expression for the deceleration parameter q is derived that can possibly describe the acceleration of the universe without dark energy, due to the effects associated with very long wave length super-horizon inflationary perturbations. The suggestion by Kolbe et al. [9] that long wave length super-horizon inflationary modes can affect a local observable through inhomogeneities is considered in the light of our exact inhomogeneous model.

A scalar-tensor bimetric gravity model of early universe cosmology is reviewed. The metric frame with a variable speed of light (VSL) and a constant speed of gravitational waves is used to describe a Friedmann-Robertson-Walker universe. The Friedmann equations are solved for a radiation dominated equation of state and the power spectrum is predicted to be scale invariant with a scalar mode spectral index $n_s=0.97$. The scalar modes are born in a ground state superhorizon and the fluctuation modes are causally connected by the VSL mechanism. The cosmological constant is equated to zero and there is no significant dependence on the scalar field potential energy. A possible way of distinguishing the metric gravity model from standard inflationary models is discussed.

The problem of explaining the acceleration of the expansion of the universe and the observational and theoretical difficulties associated with dark matter and dark energy are discussed. The possibility that Einstein gravity does not correctly describe the large-scale structure of the universe is considered and an alternative gravity theory is proposed as a possible resolution to the problems.

The nonsymmetric gravitational theory predicts an acceleration law that modifies the Newtonian law of attraction between particles. For weak fields a fit to the flat rotation curves of galaxies is obtained in terms of the mass (mass-to-light ratio M/L) of galaxies. The fits assume that the galaxies are not dominated by exotic dark matter. The equations of motion for test particles reduce for weak gravitational fields to the GR equations of motion and the predictions for the solar system and the binary pulsar PSR 1913+16 agree with the observations. The gravitational lensing of clusters of galaxies can be explained without exotic dark matter.

A mechanism for suppressing the cosmological constant is developed, based on an analogy with a superconducting phaseshift in which free fermions coupled perturbatively to a weak gravitational field are in an unstable false vacuum state. The coupling of the fermions to the gravitational field generates fermion condensates with zero momentum and a phase transition induces a nonperturbative transition to a true vacuum state by producing a positive energy gap $\Delta$ in the vacuum energy, identified with $\sqrt{\Lambda}$, where $\Lambda$ is the cosmological constant. In the strong coupling limit a large cosmological constant induces a period of inflation in the early universe, followed by a weak coupling limit in which $\sqrt{\Lambda}$ vanishes exponentially fast as the universe expands due to the dependence of the energy gap on the density of Fermi surface fermions, $D({\epsilon})$, predicting a small cosmological constant in the present universe.

We use the idea of the symmetry between the spacetime coordinates x^\mu and the energy-momentum p^\mu in quantum theory to construct a momentum space quantum gravity geometry with a metric s_{\mu\nu} and a curvature P^\lambda_{\mu\nu\rho}. For a closed maximally symmetric momentum space with a constant 3-curvature, the volume of the p-space admits a cutoff with an invariant maximum momentum a. A Wheeler-DeWitt-type wave equation is obtained in the momentum space representation. The vacuum energy density and the self-energy of a charged particle are shown to be finite, and modifications of the electromagnetic radiation density and the entropy density of a system of particles occur for high frequencies.

The cosmological constant problem is reviewed and a possible quantum gravity resolution is proposed. A space satellite E\"otv\"os experiment for zero-point vacuum energy is proposed to see whether Casimir vacuum energy falls in a gravitational field at the same rate as ordinary matter.

The graviton is pictured as a bound state of a fermion and anti-fermion with the spacetime metric assumed to be a composite object of spinor fields, based on a globally Lorentz invariant action proposed by Hebecker and Wetterich. The additional degrees of freedom beyond those of the graviton are described by Goldstone boson gravitational degrees of freedom. If we assume that the fermion is a light neutrino with mass $m_nu\sim 10^{-3}$ eV, then we obtain the effective vacuum density ${\bar\rho}_\lambda\sim (10^{-3} eV)^4$, which agrees with the estimates for the cosmological constant from WMAP and SNIa data.

We propose that local Lorentz invariance is spontaneously violated at high energies, due to a nonvanishing vacuum expectation value of a vector field \phi^\mu, as a possible explanation of the observation of ultra-high energy cosmic rays with an energy above the GZK cutoff. Certain consequences of spontaneous breaking of Lorentz invariance in cosmology are discussed.

A possible solution to the problem of providing a spacetime description of the transmission of signals for quantum entangled states is obtained by using a bimetric spacetime structure, in which quantum entanglement measurements alter the structure of the classical relativity spacetime. A bimetric gravity theory locally has two lightcones, one which describes classical special relativity and a larger lightcone which allows light signals to communicate quantum information between entangled states, after a measurement device detects one of the entangled states. The theory would remove the tension that exists between macroscopic classical, local gravity and macroscopic nonlocal quantum mechanics.

An approximately scale invariant spectrum generating the seeds of structure formation is derived from a bimetric gravity theory. By requiring that the amplitude of the CMB fluctuations from the model matches the observed value, we determine the fundamental length scale in the model to be a factor of 10^5 times larger than the Planck length, which results in a scalar mode spectral index: n_s\approx 0.97, and its running: \alpha_s\approx -5\times 10^{-4}. This is accomplished in the variable speed of light (VSL) metric frame, in which the dynamics of perturbations of the bimetric scalar field are determined by a minimally-coupled Klein-Gordon equation, and it is assumed that modes are born in a ground state at a scale given by the fundamental length scale appearing in the bimetric structure. We show that while this is taking place for scales of interest, the background (primordial) radiation energy density is strongly suppressed as a result of the bimetric structure of the model. Nevertheless, the enlarged lightcone of matter fields ensures that the horizon and flatness problems are solved.

It is argued by Duff that only the time variation of dimensionless constants of nature is a legitimate subject of enquiry, and that dimensional constants such as c, h-bar and G... are merely human constructs whose value has no operational meaning. We refute this claim and point out that such varying dimensional "constants" can have significant physical consequences for the universe that can be directly measured in experiments. Postulating that dimensional constants vary in time can significantly change the laws of physics.

An E\"otv\"os experiment to test the weak equivalence principle (WEP) for zero-point vacuum energy is proposed using a satellite. Following the suggestion of Ross for a terrestrial experiment of this type, the acceleration of a spherical test mass of aluminum would be compared with that of a similar test mass made from another material. The estimated ratio of the zero-point vacuum energy density inside the aluminum sphere to the rest mass energy density is ~ 1.6 X 10^{-14}, which would allow a 1% resolution of a potential WEP violation observed in a satellite mission test that had a baseline sensitivity to WEP violations of ~ 10^{-16}. An observed violation of the WEP for vacuum energy density would constitute a significant clue as to the origin of the cosmological constant and the source of dark energy, and test a recently proposed resolution of the cosmological constant problem, based on a model of nonlocal quantum gravity and quantum field theory.

A nonlocal quantum gravity theory is presented which is finite and unitary to all orders of perturbation theory. Vertex form factors in Feynman diagrams involving gravitons suppress graviton and matter vacuum fluctuation loops by introducing a low-energy gravitational scale, \Lambda_{Gvac} < 2.4\times 10^{-3} eV. Gravitons coupled to non-vacuum matter loops and matter tree graphs are controlled by a vertex form factor with the energy scale, \Lambda_{GM} < 1-10 TeV.

Two variable speed of light models and their physical consequences are investigated.

It is generally believed that inflationary cosmology explains the isotropy, large scale homogeneity and flatness as well as predicting the deviations from homogeneity of our universe. We show that this is not the only cosmology which can explain successfully these features of the universe. We consider anew and modify a model in which local Lorentz invariance is spontaneously broken in the very early universe, and in this epoch the speed of light undergoes a first or second order phase transition to a value ~ 30 orders of magnitude smaller, corresponding to the presently measured speed of light. Before the phase transition at a time t ~ t_c, the entropy of the universe is reduced by many orders of magnitude, allowing for a semiclassical quantum field theory calculation of a scale invariant fluctuation spectrum. After the phase transition has occurred, the radiation density and the entropy of the universe increase hugely and the increase in the entropy follows the arrow of time determined by the spontaneously broken direction of the vev <\phi^a>_0. This solves the enigma of the arrow of time and the second law of thermodynamics. A new calculation of the primordial Gaussian and adiabatic fluctuation spectrum is carried out, leading to a scale invariant scalar component of the power spectrum. We argue that there are several attractive features of VSL theory compared to standard inflationary theory, and that it provides an alternative cosmology with potentially different predictions.

A nonlocal quantum gravity theory is presented which is finite and unitary to all orders of perturbation theory. Vertex form factors in Feynman diagrams involving gravitons suppress graviton and matter vacuum fluctuation loops by introducing a low-energy gravitational scale, Lambda_{Gvac} < 2.4 X 10^{-3} eV. Gravitons coupled to non-vacuum matter loops and matter tree graphs are controlled by a vertex form factor with the energy scale, Lambda_{GM} < 1-10 TeV. A satellite Eotvos experiment is proposed to test a violation of the equivalence principle for coupling of gravitons to pure vacuum energy compared to matter.

In the bimetric scalar-tensor gravitational theory there are two frames associated with the two metrics {\hat g}_{\mu\nu} and g_{\mu\nu}, which are linked by the gradients of a scalar field \phi. The choice of a comoving frame for the metric {\hat g}_{\mu\nu} or g_{\mu\nu} has fundamental consequences for local observers in either metric spacetimes, while maintaining diffeomorphism invariance. When the metric g_{\mu\nu} is chosen to be associated with comoving coordinates, then the speed of light varies in the frame with the metric {\hat g}_{\mu\nu}. Observers in this frame see the dimming of supernovae because of the increase of the luminosity distance versus red shift, due to an increasing speed of light in the early universe. Moreover, in this frame the scalar field \phi describes a dark energy component in the Friedmann equation for the cosmic scale without acceleration. If we choose {\hat g}_{\mu\nu} to be associated with comoving coordinates, then an observer in the g_{\mu\nu} metric frame will observe the universe to be accelerating and the supernovae will appear to be farther away. The theory predicts that the gravitational constant G can vary in spacetime, while the fine-structure constant \alpha=e^2/\hbar c does not vary. The problem of cosmological horizons as viewed in the two frames is discussed.

A finite quantum gravity theory is used to resolve the cosmological constant problem. A fundamental quantum gravity scale, \Lambda_G \leq 10^{-3} eV, is introduced above which the quantum corrections to the vacuum energy density coupled to gravity are exponentially suppressed by a graviton vertex form factor, yielding an observationally acceptable value for the particle physics contribution to the cosmological constant. Classical Einstein gravity retains its causal behavior as well as the standard agreement with observational data.

Spacetime geometry is treated as a fluctuating, stochastic quantum system allowing an effective quantum gravity solution to the cosmological constant problem. A Focker-Planck equation for the probability density of spacetime metric fluctuations is solved for an FRW universe and the most probable value for the metric scale factor R is for a very small cosmological constant.

Predictions of the CMB spectrum from a bimetric gravity theory (gr-qc/0101126) are presented. The initial inflationary period in BGT is driven by a vanishingly small speed of gravitational waves v_g in the very early universe. This initial inflationary period is insensitive to the choice of scalar field potential and initial values of the scalar field. After this initial period of inflation, v_g will increase rapidly and the effects of a potential will become important. We show that a quadratic potential introduced into BGT yields an approximately flat spectrum with inflation parameters: n_s=0.98, n_t=-0.027, alpha_s=-3.2e-4 and alpha_t=-5.0e-4, with r >= 0.014.

We construct an eleven-dimensional superspace with superspace coordinates and formulate a finite M-theory using non-anticommutative geometry. The conjectured M-theory has the correct eleven-dimensional supergravity low energy limit. We consider the problem of finding a stable finite M-theory which has de Sitter space as a natural ground state, and the problem of eliminating possible future horizons.

A conjectured finite M-theory based on eleven-dimensional supergravity formulated in a superspace with a non-anticommutative <>-product of field operators is proposed. Supermembranes are incorporated in the superspace <>-product formalism. When the deformed supersymmetry invariant action of eleven-dimensional supergravity theory is expanded about the standard supersymmetry invariant action, the spontaneously compactified M-theory can yield a four-dimensional de Sitter spacetime inflationary solution with dark energy described by the four-form F-fields. A fit to the present cosmological data for an accelerating universe is obtained from matter fields describing the dominant dark matter and the four-form F-field dark energy. Chiral fermions are obtained from the M-theory by allowing singularities in the compact internal seven-dimensional space. The possibility of obtaining a realistic M-theory containing the standard model is discussed.

The recent evidence for a cosmological evolution of the fine structure constant \alpha=e^2/\hbar c found from an analysis of absorption systems in the spectra of distant quasars, is modelled by a cosmological scenario in which it is assumed that only the speed of light varies. The model fits the spectral line data and can also lead to a solution of the initial value problems in cosmology.

An M-theory constructed in an eleven-dimensional supermanifold with a <>-product of field operators is shown to have a de Sitter space solution. Possible implications of this result for cosmology are mentioned.

The problem of explaining the acceleration of the expansion of the universe and the observational and theoretical difficulties associated with dark matter and dark energy are discussed. The possibility that GR does not correctly describe the large-scale structure of the universe is considered and an alternative gravity theory is proposed as a possible resolution to the problems.

We impose in the nonsymmetric gravitational theory, by means of Lagrange multiplier fields in the action, a set of covariant constraints on the antisymmetric tensor field. The canonical Hamiltonian constraints in the weak field approximation for the antisymmetric sector yield a Hamiltonian energy bounded from below. An analysis of the Cauchy evolution, in terms of an expansion of the antisymmetric sector about a symmetric Einstein background, shows that arbitrarily small antisymmetric Cauchy data can lead to smooth evolution.

The existence of future horizons in spacetime geometries poses serious problems for string theory and quantum field theories. The observation that the expansion of the universe is accelerating has recently been shown to lead to a crisis for the mathematical formalism of string and M-theories, since the existence of a future horizon for an eternally accelerating universe does not allow the formulation of physical S-matrix observables. Postulating that the speed of light varies in an expanding universe in the future as well as in the past can eliminate future horizons, allowing for a consistent definition of S-matrix observables.

The dynamical consequences of a bimetric scalar-tensor theory of gravity with a dynamical light speed are investigated in a cosmological setting. The model consists of a minimally-coupled self-gravitating scalar field coupled to ordinary matter fields in the standard way through the metric: $\metric_{\mu\nu}+B\partial_\mu\phi\partial_\nu\phi$. We show that in a universe with matter that has a radiation-dominated equation of state, the model allows solutions with a de Sitter phase that provides sufficient inflation to solve the horizon and flatness problems. This behaviour is achieved without the addition of a potential for the scalar field, and is shown to be largely independent of its introduction. We therefore have a model that is fundamentally different than the potential-dominated, slowly-rolling scalar field of the standard models inflationary cosmology. The speed of gravitational wave propagation is predicted to be significantly different from the speed of matter waves and photon propagation in the early universe.

In earlier published work, it was proposed that light speed was larger in the early universe by 30 orders of magnitude compared to its presently observed value. This change in the speed of light is associated with a spontaneous breaking of local Lorentz invariance in the early Universe, associated with a first order phase transition at a critical time $t=t_c$. This solves the horizon problem, leads to a mechanism of monopole suppression in cosmology and can resolve the flatness problem. it also offers the potential of solving the cosmological constant problem. After the critical time $t_c$, local Lorentz invariance is restored and light travels at its presently measured speed. We investigate the field equations in the spontaneously broken phase and study further the flatness problem and the cosmological constant problem. The entropy is shown to undergo a large increase as the light velocity goes through a phase transition. A scale invariant prediction for microwave background fluctuations originating at the time of the phase transition is discussed.

A quantum field theory formalism is reviewed that leads to a self-consistent, finite quantum gravity, Yang-Mills and Higgs theory, which is unitary and gauge invariant to all orders of perturbation theory. The gauge hierarchy problem is solved due to the exponential damping of the Higgs self-energy loop graph for energies greater than a scale $\Lambda_H\leq 1$ TeV. The cosmological constant problem is solved by introducing a fundamental quantum gravity scale, $\Lambda_G\leq 10^{-4}$ eV, above which the virtual contributions to the vacuum energy density coupled to gravity are exponentially suppressed, yielding an observationally acceptable value for the particle physics contribution to the cosmological constant. Classical Einstein gravity retains its causal behavior as well as the standard agreement with observational data. Possible experimental tests of the onset of quantum nonlocality at short distances are considered.

A brief review is given of the present observational data in cosmology. A review of a new bimetric gravity theory with multiple light cones is presented. The physical consequences of this gravity theory for the early universe are analyzed.

A vector-tensor theory of gravity that was introduced in an earlier publication is analyzed in detail and its consequences for early universe cosmology are examined. The multiple light cone structure of the theory generates different speeds of gravitational and matter wave fronts, and the contraction of these light cones produces acausal, superluminary inflation that can resolve the initial value problems of cosmology.

A gravity theory is developed with the metric ${\hat g}_{\mu\nu}= {g}_{\mu\nu}+B\partial_\mu\phi\partial_\nu\phi$. In the present universe the additional contribution from the scalar field in the metric ${\hat g}_{\mu\nu}$ can generate an acceleration in the expansion of the universe, without negative pressure and with a zero cosmological constant. In this theory, gravitational waves will propagate at a different speed from non-gravitational waves. It is suggested that gravitational wave experiments could test this observational signature.

A dynamical model for varying light velocity in cosmology is developed, based on the idea that there are two metrics in spacetime. One metric $g_{\mu\nu}$ describes the standard gravitational vacuum, and the other ${\hat g}_{\mu\nu} =g_{\mu\nu}+\beta\psi_\mu\psi_\nu$ describes the geometry through which matter fields propagate. Matter propagating causally with respect to $\hat{g}_{\mu\nu}$ can provide acausal contributions to the matter stress-energy tensor in the field equations for $g_{\mu\nu}$, which, as we explicitly demonstrate with perfect fluid and scalar field matter models, provides a mechanism for the solution of the horizon, flatness and magnetic monopole problems in an FRW universe. The field equations also provide a ``graceful exit'' to the inflationary epoch since below an energy scale (related to the mass of $\psi_\mu$) we recover exactly the standard FRW field equations.

The spontaneous breaking of local Lorentz invariance in the early Universe, associated with a first order phase transition at a critical time $t_c$, generates a large increase in the speed of light and a superluminary communication of information occurs, allowing all regions in the Universe to be causally connected. This solves the horizon problem, leads to a mechanism of monopole suppression in cosmology and can resolve the flatness problem. After the critical time $t_c$, local Lorentz (and diffeomorphism) invariance is restored and light travels at its presently measured speed. The kinematical and dynamical aspects of the generation of quantum fluctuations in the superluminary Universe are investigated. A scale invariant prediction for the fluctuation density amplitude is obtained.

Models of superdense star clusters at the center of galaxies are investigated to see whether such objects can be stable and long-lived based on evaporation and collision time-scales and stability criteria. We find that physically reasonable models of massive clusters of stellar remnants can exist with masses $\geq 10^6 M_{\odot}$, which could simulate black holes at the center of galaxies with large $M/L$ ratios and gas motions of order $\geq 10^3$ km $s^{-1}$. It follows that the evidence is not conclusive for massive dark objects at the center of galaxies being black holes.

We investigate the fate of a relativistic star cluster with a dense core which is undergoing a gravothermal catastrophe and is far from thermodynamic equilibrium. Nonlinear cooperative contributions are included in the standard transport equations for the last stage of evolution of a highly dense core of stellar remnants. We find that the core redshift does not necessarily increase without limit as the core becomes increasingly dense, preventing collapse to a black hole. In particular, the redshift can remain less than the critical value for relativistic collapse, resulting in a stable, massive dark core with a Newtonian mantle and halo.

Cosmological consequences of the nonsymmetric gravitational theory (NGT) are studied. The structure of the NGT field equations is analyzed for an inhomogeneous and anisotropic universe, based on the spherically symmetric field equations. It is assumed that the matter density and pressure are purely time dependent, and it is shown that the field equations allow open, flat and closed universes. The field equations are expanded about the Friedmann-Robertson-Walker (FRW) model for a small antisymmetric field. The observations are most in accord with the approximately spatially flat universe, which has an FRW metric to lowest order in the skew field. The recent analysis of galaxy synchroton radiation polarization data, which indicates a residual polarization after Faraday polarization is subtracted, can be explained by the birefringence of electromagnetic waves produced by the NGT spacetime.

A model of the universe as a self-organized critical system is considered. The universe evolves to a state independently of the initial conditions at the edge of chaos. The critical state is an attractor of the dynamics. Random metric fluctuations exhibit noise without any characteristic length scales, and the power spectrum for the fluctuations has a self-similar fractal behavior. In the early universe, the metric fluctuations smear out the local light cones removing the horizon problem.

A covariant Lagrangian formulation of a solution to the cosmological constant problem, based on vizualising the fluctuations of the vacuum energy as a non-equilibrium process with stochastic behaviour, is presented. The variational principle yields equations of motion for the cosmological ``constant" $\Lambda$, treated as a dynamical field, together with an equation for a Lagrange multiplier field $\phi$, and the standard Einstein field equations with a variable cosmological constant term. A stochastic model of $\Lambda$ yields a natural explanation for the smallness or zero value of the constant in the present epoch and its large value in an era of inflation in the early universe.

The fluctuations of the vacuum energy are treated as a non-equilibrium process and a stochastic model for the cosmological constant is presented, which yields a natural explanation for the smallness or zero value of the constant in the present epoch and its large value in an era of inflation in the early universe.

The gravitational wave solutions obtained from a perturbation about conformally flat backgrounds in Einstein gravity are investigated. A perturbation theory analysis of the Lesame, Ellis and Dunsby results, based on a covariant approach, shows that for gravitational waves interacting with irrotational dust, the equations are linearization unstable. The gravitational wave equations based on the Weyl curvature tensor must be solved by non-perturbative methods. The significance of this result for gravitational wave calculations and experiments is discussed.

It is shown that the new version of nonsymmetric gravitational theory (NGT) corresponds in the linear approximation to linear Einstein gravity theory and antisymmetric field equations with a non-conserved string source current. The Hamiltonian for the antisymmetric field equations is bounded from below and describes the exchange of a spin $1^+$ massive vector boson between open strings. The non-Riemannian geometrical theory is formulated in terms of a nonsymmetric fundamental tensor $g_{\mu\nu}$. The weak field limit, $g_{[\mu\nu]}\rightarrow 0$, associated with large distance scales, corresponds to the limit to a confinement region at low energies described by an effective Yukawa potential at galactic distance scales. The limit to this low-energy confinement region is expected to be singular and non-perturbative. The NGT string theory predicts that there are no black hole event horizons associated with infinite red shift null surfaces.

We establish the formalism in the nonsymmetric gravitational theory (NGT) for stellar equilibrium and gravitational collapse. We study the collapse of a pressureless, spherically symmetric dust cloud. By assuming that the interior solution is smoothly matched at the surface of the star to the quasi-static, spherically symmetric vacuum solution, we find that the star does not collapse to a black hole. It is anticipated that the final collapsed object will reach a state of equilibrium, and will emit thermal, gravitational and other forms of radiation, although the radiation may be emitted only in small amounts if the red shift from the surface of the compact object is large. No Hawking radiation is emitted and the information loss problem can be resolved at the classical level.

The analytical structure of the difference between the static vacuum solution in the nonsymmetric gravitational theory (NGT) and the Schwarzschild solution of Einstein's gravitational theory (EGT) is studied. It is proved that a smooth matching of the solutions does not exist in the range $0 < r \leq 2M$, for any non-zero values of the parameters $M$ and $s$ of the NGT solution. This means that one cannot consider the difference between the two solutions using perturbation theory in this range of $r$. Assuming that the exterior solution in gravitational collapse is a small, time dependent perturbation of the static solution for a non-zero, generic NGT source ($s\not=0$) and mass density, it is shown that the matching of the interior and exterior solutions will not lead to black hole event horizons.

In the weak field approximation, the nonsymmetric gravitational theory (NGT) has, in addition to the Newtonian gravitational potential, a Yukawa-like potential produced by the exchange of a spin $1^+$ boson between fermions. If the range $r_0=\mu^{-1}$ is $25$ kpc, then this additional potential due to the interaction with matter in the halos of galaxies can explain the flat rotation curves of galaxies and the Tully-Fisher law ($L\sim v^4$) without the dark matter hypothesis. Possible fits to clusters of galaxies and gravitational lensing observations are discussed. The results are based on a linear approximation to a new perturbatively consistent version of the NGT field equations, which does not violate the weak equivalence principle.

In the weak field approximation, the nonsymmetric gravitational theory has, in addition to the Newtonian gravitational potential, a Yukawa potential produced by the exchange of a spin $1^+$ boson between fermions. If the range $r_0$ is of order $30$ kpc, then the potential due to the interaction of known neutrinos in the halos of galaxies can explain the flat rotation curves of galaxies. The results are based on a physical linear approximation to the NGT field equations and they are consistent with equivalence principle observations, other solar system gravitational experiments and the binary pulsar data.

A dynamical model of the decaying vacuum energy is presented, which is based on Jordan-Brans-Dicke theory with a scalar field $\phi$. The solution of an evolutionary differential equation for the scalar field $\phi$ drives the vacuum energy towards a cosmological constant at the present epoch that can give for the age of the universe, $t_0\sim 13.5$ Gyr for $\Omega_0=1$, which is consistent with the age of globular clusters.

A local void in the globally Friedmann-Robertson-Walker (FRW) cosmological model is studied. The inhomogeneity is described using the Lema\^{\i}tre-Tolman-Bondi (LTB) solution with the spherically symmetric matter distribution based on the faint galaxies number counts. We investigate the effects this has on the measurement of the Hubble constant and the redshift--luminosity distance relation for moderately and very distant objects ($z \approx 0.1$ and more). The results, while fully compatible with cosmological observations, indicate that if we happened to live in such a void, but insisted on interpreting cosmological observations through the FRW model, we could get a few unexpected results. For example the Hubble constant measurement could give results depending on the separation of the source and the observer.

A massive gravitationally bound object with a radius $r \leq 2GM/c^2$, which occurs in the non-singular version of the nonsymmetric gravitational theory (NGT), replaces the black hole in Einstein gravity theory. This object is kept stable by the attractive and repulsive forces generated by NGT, as well as standard matter pressures, and is called a superdense object (SDO). The luminosity of gravitational waves emitted by a SDO with a red-shift of order unity is calculated and it is found that it could be a strong source of gravitational radiation at low frequencies. The active galactic nucleus of M87 is identified with a SDO and the recent observational results obtained by the Hubble Space Telesope are used to estimate the amplitude of gravitational radiation.

A local void in the globally Friedmann-Robertson-Walker cosmological model with the critical density ($\Omega_{0}=1$) is studied. The inhomogeneity is described using a Lema\^{\i}tre-Tolman-Bondi solution for a spherically symmetric distribution of matter. The scale of the central underdense region is $\sim 150$ Mpc. We investigate the effects this has on the cosmological time scale, the measurement of the Hubble constant and the redshift--luminosity distance for moderately and very distant objects ($z \sim 0.1$ and more). The results indicate that if we happened to live in such a void, but insisted on interpreting cosmological observations through the FRW model, we could go wrong in a few instances. For example, the Hubble constant measurement could give results depending on the separation of the source and the observer, the quasars could be younger than we think and also less distant (less energetic).

Damour, Deser and McCarthy have claimed that the nonsymmetric gravitational theory (NGT) is untenable due to curvature coupled ghost modes and bad asymptotic behavior. This claim is false for it is based on a physically inaccurate treatment of wave propagation on a curved background and an incorrect method for extracting asymptotic behavior. We show that the flux of gravitational radiation in NGT is finite in magnitude and positive in sign.