In this article, we present the results of a fourth order perturbation analysis of the metric theory of gravity $f(\chi) = \chi^{3/2}$, with $\chi$ a suitable dimensionless Ricci scalar. Such model corresponds to a specific $f(R)$ metric theory of gravity, where the mass of the system is included into the gravitational field's action. In previous works we have shown that, up to the second order in perturbations, this theory reproduces flat rotation curves of galaxies and the details of the gravitational lensing in individual, groups and clusters of galaxies. Here, leaving fixed the results from our previous works, we show that the theory reproduces the dynamical masses of 12 Chandra X-ray galaxy clusters, without the need of dark matter, through the metric coefficients up to the fourth order of approximation. In this sense, we calculate the first relativistic correction of the $f(\chi)$ metric theory and apply it to fit the dynamical masses of the clusters of galaxies.

Cold dark matter (CDM) has shown to be an excellent candidate for the dark matter (DM) of the Universe at large scales, however it presents some challenges at the galactic level. The scalar field dark matter (SFDM), also called fuzzy, wave, Bose-Einstein condensate or ultra-light axion DM, is identical to CDM at cosmological scales but different at the galactic ones. SFDM forms core halos, it has a natural cut-off in its matter power spectrum and it predicts well-formed galaxies at high redshifts. In this work we reproduce the rotation curves of high-resolution low surface brightness (LSB) and SPARC galaxies with two SFDM profiles: (1)~The soliton+NFW profile in the fuzzy DM (FDM) model, arising empirically from cosmological simulations of real, non-interacting scalar field (SF) at zero temperature, and (2)~the multistate SFDM (mSFDM) profile, an exact solution to the Einstein-Klein-Gordon equations for a real, self-interacting SF, with finite temperature into the SF potential, introducing several quantum states as a realistic model for a SFDM halo. From the fits with the soliton+NFW profile, we obtained for the boson mass $0.212< m_\psi/(10^{-23} \mathrm{eV}/c^2)<27.0$ and for the core radius $0.326< r_c/\mathrm{kpc}<8.96$. From the combined analysis with the LSB galaxies, we obtained $m_\psi = 0.554 \times10^{-23}\mathrm{eV}$, a result in tension with the severe cosmological constraints. Also, we show the analytical mSFDM model fits the observations as well as or better than the empirical soliton+NFW profile, and it reproduces naturally the wiggles present in some galaxies, being a theoretically motivated framework additional or alternative to the FDM profile.

One alternative to the cold dark matter (CDM) paradigm is the scalar field dark matter (SFDM) model, which assumes dark matter is a spin-0 ultra-light scalar field (SF) with a typical mass $m\sim10^{-22}\mathrm{eV}/c^2$ and positive self-interactions. Due to the ultra-light boson mass, the SFDM could form Bose-Einstein condensates (BEC) in the very early Universe, which are interpreted as the dark matter haloes. Although cosmologically the model behaves as CDM, they differ at small scales: SFDM naturally predicts fewer satellite haloes, cores in dwarf galaxies and the formation of massive galaxies at high redshifts. The ground state (or BEC) solution at zero temperature suffices to describe low-mass galaxies but fails for larger systems. A possible solution is adding finite-temperature corrections to the SF potential which allows combinations of excited states. In this work, we test the finite-temperature multistate SFDM solution at galaxy cluster scales and compare our results with the Navarro-Frenk-White (NFW) and BEC profiles. We achieve this by fitting the mass distribution of 13 \textit{Chandra} X-ray clusters of galaxies, excluding the region of the brightest cluster galaxy. We show that the SFDM model accurately describes the clusters' DM mass distributions offering an equivalent or better agreement than the NFW profile. The complete disagreement of the BEC model with the data is also shown. We conclude that the theoretically motivated multistate SFDM profile is an interesting alternative to empirical profiles and ad hoc fitting-functions that attempt to couple the asymptotic NFW decline with the inner core in SFDM.

In this article we perform a second order perturbation analysis of the gravitational metric theory of gravity $ f(\chi) = \chi^{3/2} $ developed by Bernal et al. (2011). We show that the theory accounts in detail for two observational facts: (1) the phenomenology of flattened rotation curves associated to the Tully-Fisher relation observed in spiral galaxies, and (2) the details of observations of gravitational lensing in galaxies and groups of galaxies, without the need of any dark matter. We show how all dynamical observations on flat rotation curves and gravitational lensing can be synthesised in terms of the empirically required metric coefficients of any metric theory of gravity. We construct the corresponding metric components for the theory presented at second order in perturbation, which are shown to be perfectly compatible with the empirically derived ones. It is also shown that under the theory being presented, in order to obtain a complete full agreement with the observational results, a specific signature of Riemann's tensor has to be chosen. This signature corresponds to the one most widely used nowadays in relativity theory. Also, a computational program, the MEXICAS (Metric EXtended-gravity Incorporated through a Computer Algebraic System) code, developed for its usage in the Computer Algebraic System (CAS) Maxima for working out perturbations on a metric theory of gravity, is presented and made publicly available.

Using dimensional analysis techniques we present an extension of Newton's gravitational theory built under the assumption that Milgrom's acceleration constant is a fundamental quantity of nature. The gravitational force converges to Newton's gravity and to a MOND-like description in two different mass and length regimes. It is shown that a modification on the force sector (and not in the dynamical one as MOND does) is more convenient and can reproduce and predict different phenomena usually ascribed to dark matter at the non-relativistic level.

We show that the Modified Newtonian Dynamics (MOND) regime can be fully recovered as the weak-field limit of a particular theory of gravity formulated in the metric approach. This is possible when Milgrom's acceleration constant is taken as a fundamental quantity which couples to the theory in a very consistent manner. As a consequence, the scale invariance of the gravitational interaction is naturally broken. In this sense, Newtonian gravity is the weak-field limit of general relativity and MOND is the weak-field limit of that particular extended theory of gravity.

We show that the Modified Newtonian Dynamics (MOND) regime can be fully recovered as the weak-field limit of a particular theory of gravity formulated in the metric approach. This is possible when Milgrom's acceleration constant is taken as a fundamental quantity which couples to the theory in a very consistent manner. As a consequence, the scale invariance of the gravitational interaction is naturally broken. In this sense, Newtonian gravity is the weak-field limit of general relativity and MOND is the weak-field limit of that particular extended theory of gravity. We also prove that a Noether's symmetry approach to the problem yields a conserved quantity coherent with this relativistic MONDian extension.

Some quantum-cosmic scaling relations indicate that the MOND acceleration parameter a_0 could be a fundamental quantity ruling the self-gravitating structures, ranging from stars and globular clusters up to superclusters of galaxies and the whole observed universe. We discuss such coincidence relations starting from the Dirac quantization condition ruling the masses of primordial black holes.

In the pursuit of a general formulation for a modified gravitational theory at the non-relativistic level and as an alternative to the dark matter hypothesis, we construct a model valid over a wide variety of astrophysical scales. Through the inclusion of Milgrom's acceleration constant into a gravitational theory, we show that very general formulas can be constructed for the acceleration felt by a particle. Dimensional analysis shows that this inclusion naturally leads to the appearance of a mass-length scale in gravity, breaking its scale invariance. A particular form of the modified gravitational force is constructed and tested for consistency with observations over a wide range of astrophysical environments, from solar system to extragalactic scales. We show that over any limited range of physical parameters, which define a specific class of astrophysical objects, the dispersion velocity of a system must be a power law of its mass and size. These powers appear linked together through a natural constraint relation of the theory. This yields a generalised gravitational equilibrium relation valid for all astrophysical systems. A general scheme for treating spherical symmetrical density distributions is presented, which in particular shows that the fundamental plane of elliptical galaxies, the Newtonian virial equilibrium, the Tully-Fisher and the Faber-Jackson relations, as well as the scalings observed in local dwarf spheroidal galaxies, are nothing but particular cases of that relation when applied to the appropriate mass-length scales. We discuss the implications of this approach for a modified theory of gravity and emphasise the advantages of working with the force, instead of altering Newton's second law of motion, in the formulation of a gravitational theory.

We use a specific form of the interpolation function in the MOND formalism, which optimally accounts for the internal structure of dwarf spheroidal (dSph) galaxies, to explore the consequences it has on the scaling relations seen in these systems.}} {{The particular form of the interpolation function we used leads to a law of gravity that does not degrade the good fit of the MOND proposal on galactic scales, and in fact, slightly improves the accordance with observations at dSph scales. This formalism yields a good description of gravitational phenomena without the need of invoking any still undetected and hypothetically dominant dark matter, in the weak field regime probed by local dSph galaxies.}} {{Isothermal equilibrium density profiles then yield projected surface density profiles for the local dSph galaxies in very good agreement with observational determinations, for values of the relevant parameters as inferred from recent observations of these Galactic satellites. The observed scaling relations for these systems are also naturally accounted for within the proposed scheme, including a previously unrecognised correlation of the inferred mass-to-light ratios of local dSph's with the ages of their stellar populations, which is natural in modified gravity schemes in the absence of dark matter.}} {The results shed some light on the form that the MOND interpolating function may have in the most challenging regime, which occurs at moderate accelerations and intermediate mass--weighted lengths.

Recently, Y. Sobouti (2007) has provided a metric theory f(R) that can account for certain dynamical anomalies observed in spiral galaxies. Mendoza & Rosas-Guevara (2007) have shown that in this theory there is an extra-bending as compared to standard general relativity. In the present work we have developed in more specific detail this additional lensing effect and we have made evaluations of the alpha parameter used in the model adjusting the theory to observations in X-rays of 13 clusters of galaxies with gravitational lensing (Hoekstra (2007)).