We study the gravitomagnetism in the TeVeS theory. We compute the gravitomagnetic field that a slow-moving mass distribution produces in its Newtonian regime. We report that the consistency between the TeVeS gravitomagnetic field and that predicted by the Einstein-Hilbert theory leads to a relation between the vector and scalar coupling constants of the theory. We observe that requiring consistency between the near horizon geometry of a black hole in TeVeS and the image of the black hole taken Event Horizon Telescope leads to another relation between the coupling constants of the TeVeS theory and enable us to identify the coupling constants of the theory.

We introduce the Gauge Vector-Tensor (GVT) theory by extending the AQUAL's approach to the GravitoElectroMagnetism (GEM) approximation of gravity. GVT is a generally covariant theory of gravity composed of a pseudo Riemannian metric and two U(1) gauge connections that reproduces MOND in the limit of very weak gravitational fields while remains consistent with the Einstein-Hilbert gravity in the limit of strong and Newtonian gravitational fields. GVT also provides a simple framework to study the GEM approximation to gravity. We illustrate that the gravitomagnetic force at the edge of a galaxy can be in accord with either GVT or \Lambda CDM but not both. We also study the physics of the GVT theory around the gravitational saddle point of the Sun and Jupiter system. We notice that the conclusive refusal of the GVT theory demands measuring either both of the gravitoelectric and gravitomagnetic fields inside the Sun-Jupiter MOND window, or the gravitoelectric field inside two different solar GVT MOND windows. The GVT theory, however, will be favored by observing an anomaly in the gravitoelectric field inside a single MOND window.

The MOND paradigm to the missing mass problem requires introducing a functional that is to be identified through observations and experiments. We consider AQUAL theory as a realization of the MOND. We show that the accurate value of the Earth GM measured by the Lunar Laser Ranging and that by various artificial Earth satellites, including the accurate tracking of the LAGEOS satellites, constrain this functional such that some of the chosen/proposed functional are refuted.

We consider the AQUAL theory - a theory of modified gravity capable of resolving the missing mass problem - and study its predictions for micro gravity tests at the gravitational saddle points of the Solar system. We report that the AQUAL model enhances the gravity at the sub-micrometer ranges around the gravitational saddle points in a way that so far has been unnoticed. This enhancement can be measured. We, therefore, call for moving toward implementing micrometer gravity tests within the Solar gravitational saddle points.

We calculate the birefringent in the vacuum for light at the leading and sub-leading orders for the CPT-even part of the SME. We report that all the LIV coefficients absent in the leading order, but the isotropic one, contributes to the sub-leading order birefringent. We consider models free of the first order birefringent. We then show that infrared, optical, and ultraviolet spectropolarimetry of cosmological sources bound the LIV coefficients to less than $ 10^{-16}$. This improves the best current bound on the parity-odd coefficients by two orders of magnitude and establishes the isotropy of the one-way light speed with the precision of $41 \frac{nm}{s}$.

We covariantly modify the Einstein-Hilbert action such that the modified action perturbatively resolves the flat rotational velocity curve of the spiral galaxies and gives rise to the Tully-Fisher relation, and dynamically generates the cosmological constant. This modification requires introducing just a single new universal parameter.

We argue that Solar system tests do not rule out 1/R gravity at least due to the reason addressed in Phys. Rev. D 74 (2006) 121501 [astro-ph/0610483] (ref. [1]) and subsequent published papers. Ref. [1] has not only modified the Einstein-Hilbert action but also has changed the boundary conditions since they altered the equations of motion. In Einstein-Hilbert action equations are second order, so the fall off of the fields suffices to single out a unique solution. In 1/R gravity the equations are fourth order, so we should impose additional boundary conditions. Perhaps the boundary condition we must impose is that the abrupt change in the second derivative of the metric near the surface of the Sun remains intact by adding `1/R' corrections to the equations of motion. The solution of 1/R gravity with this boundary condition remains consistent with the solar system tests. Ref. [1] assumes that as soon as they perturbatively modified the equations then the Ricci scalar becomes smooth on the surface of the Sun. This assumption is simply wrong because the boundary conditions and equations of motions are two different entities.

We consider corrections in the form of $\Delta L(R_{ijkl}R^{ijkl})$ to the Einstein-Hilbert Lagrangian. Then we compute the corrections to the Schwarszchild geometry due to the inclusion of this general term to the Lagrangian. We show that $\Delta L_3=\alpha_{{1/3}}(R_{ijkl}R^{ijkl})^{{1/3}}$ gives rise to a constant anomalous acceleration for objects orbiting the Sun onward the Sun. This leads to the conclusion that $\alpha_{{1/3}}=(13.91\pm 2.11) \times 10^{-26}(\frac{1}{\text{meters}})^{{2/3}}$ would have covariantly resolved the Pioneer anomaly if this value of $\alpha_{{1/3}}$ had not contradicted other observations. We notice that the experimental bounds on $\Delta L_3$ grows stronger in case we examine the deformation of the space-time geometry around objects lighter than the Sun. We therefore use the high precision measurements around the Earth (LAGEOS and LLR) and obtain a very strong constraint on the corrections in the form of $\Delta L(R_{ijkl}R^{ijkl})$ and in particular $\Delta L=\alpha_n(R_{ijkl}R^{ijkl})^n$. This bound requires $\alpha_{{1/3}}\leq6.12\times 10^{-29}(\frac{1}{\text{meters}})^{{2/3}}$. Therefore it refutes the covariant resolution of the Pioneer anomaly.

We consider extensions of the Einstein-Hilbert Lagrangian to a general functional of metric and Riemann curvature tensor. A given such Lagrangian describes two different theories depending on considering connection and metric (Palatini formulation), or only the metric (metric formulation) as independent dynamical degrees of freedom. Equivalence of the Palatini and metric formulations at the level of equations of motion, which as we will argue is a manifestation of the Equivalence Principle, is the physical criterion that restricts form of the Lagrangians of modified gravity theories. We prove that within the class of modified gravity theories we consider, only the Lovelock gravity satisfies this requirement.