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.

In this article we construct an extended relativistic f(R) theory of gravity with matter-curvature couplings F(R,L_matt) for which its weak field limit of approximation recovers the simplest version of MOND. We do this by (a) performing an order of magnitude approach and (b) by perturbing the resulting field equations of the theory to the weakest field limit of approximation. We also compute the geodesic equation of the resulting theory and show that it has an extra force, a fact that commonly appears in general matter-curvature couplings.

In this article we develop a Primitive Variable Recovery Scheme (PVRS) to solve any system of coupled differential conservative equations. This method obtains directly the primitive variables applying the chain rule to the time term of the conservative equations. With this, a traditional finite volume method for the flux is applied in order avoid violation of both, the entropy and "Rankine-Hugoniot" jump conditions. The time evolution is then computed using a forward finite difference scheme. This numerical technique evades the recovery of the primitive vector by solving an algebraic system of equations as it is often used and so, it generalises standard techniques to solve these kind of coupled systems. The article is presented bearing in mind special relativistic hydrodynamic numerical schemes with an added pedagogical view in the appendix section in order to easily comprehend the PVRS. We present the convergence of the method for standard shock-tube problems of special relativistic hydrodynamics and a graphical visualisation of the errors using the fluctuations of the numerical values with respect to exact analytic solutions. The PVRS circumvents the sometimes arduous computation that arises from standard numerical methods techniques, which obtain the desired primitive vector solution through an algebraic polynomial of the charges.

Removing outbursts from multiwavelength light curves of the blazar Mrk~421, we construct outburstless time series for this system. A model-independent power spectrum light curve analysis in the optical, hard X-ray and gamma-rays for this outburstless state and also the full light-curves, show clear evidence for a periodicity of ~ 310 days across all wavelengths studied. A subsequent full maximum likelihood analysis fitting an eclipse model confirms a periodicity of 310 \pm 1 days. The power spectrum of the signal in the outburstless state of the source does not follow a flicker noise behaviour and so, the system producing it is not self-organised. This and the fact that the periodicity is better defined in the outburstless state, strongly suggests that it is not produced by any internal physical processes associated to the central engine. The simplest physical mechanism to which this periodicity could be ascribed is a dynamical effect produced by an orbiting supermassive black hole companion eclipsing the central engine. Interestingly, the optimal eclipse model infers a brightness enhancement of (136.4 \pm 20 )%, suggesting an eclipse resulting in a gravitational lensing brightening. Consisting with this interpretation, the eclipse occurs for only ( 9.7 \pm 0.2)% of the orbital period.

Under natural assumptions on the thermodynamic properties of space and time with the holographic principle we reproduce a MOND-like behaviour of gravity on particular scales of mass and length, where Newtonian gravity requires a modification or extension if no dark matter component is introduced in the description of gravitational phenomena. The result is directly obtained with the assumption that a fundamental constant of nature with dimensions of acceleration needs to be introduced into gravitational interactions. This in turn allows for modifications or extensions of the equipartion law and/or the holographic principle. In other words, MOND-like phenomenology can be reproduced when appropriate generalised concepts at the thermodynamical level of space and/or at the holographic principle are introduced. Thermodynamical modifications are reflected in extensions to the equipartition law which occur when the temperature of the system drops below a critical value, equals to Unruh's temperature evaluated at the acceleration constant scale introduced for the description of the gravitational phenomena. Our calculations extend the ones by Verlinde (2011) in which Newtonian gravity is shown to be an emergent phenomenon, and together with it reinforces the idea that gravity at all scales is emergent.

Knots or blobs observed in astrophysical jets are commonly interpreted as shock waves moving along them. Long time observations of the HST-1 knot inside the jet of the galaxy M87 have produced detailed multi-wavelength light curves. In this article, we model these light curves using the semi-analytical approach developed by Mendoza et al. (2009). This model was developed to account for the light curves of working surfaces moving along relativistic jets. These working surfaces are generated by periodic oscillations of the injected flow velocity and mass ejection rates at the base of the jet. Using genetic algorithms to fit the parameters of the model, we are able to explain the outbursts observed in the light curves of the HST-1 knot with an accuracy greater than a 2 sigma statistical confidence level.

The following is a comment on the recent letter by Iocco et al. (2015, arXiv:1502.03821) where the authors claim to have found "...convincing proof of the existence of dark matter...". The letter in question presents a compilation of recent rotation curve observations for the Milky Way, together with Newtonian rotation curve estimates based on recent baryonic matter distribution measurements. A mismatch between the former and the latter is then presented as "evidence for dark matter". Here we show that the reported discrepancy is the well known gravitational anomaly which consistently appears when dynamical accelerations approach the critical Milgrom acceleration a_0 = 1.2 \times 10^{-10} m / s^2. Further, using a simple modified gravity force law, the baryonic models presented in Iocco et al. (2015), yield dynamics consistent with the observed rotation values.

From 1975 to 1976, an outburst was detected in the light curve of the X-ray transient A0620-00 using the Ariel V and SAS-3 experiments. In this letter we model the outburst with the hydrodynamical model proposed by Mendoza et al. (2009). The physical model is constructed assumming basic mass and momentum conservation laws associated to the motion of the shock waves developed inside the expanding relativistic jet of the source. These internal shock waves are produced as a result of periodic variations of the inyected mass and velocity of the flow at the base of the jet. The observations of this X-ray light curve present two clear bumps. The first one is modelled assuming periodic variations of the inyected velocity at the base of the jet, while the second one can either be modelled by a further velocity oscillations, or by a periodic variation of the mass injection rate at the base of the jet at a latter time. The fitting of the data fixes different parameters of the model, such as the mean mass injection rate at the base of the jet and the oscillation frequency of the flow as measured on the rest frame of the central source.

We give precise details to support that observations of gravitational lensing at scales of individual, groups and clusters of galaxies can be understood in terms of non-Newtonian gravitational interactions with a relativistic structure compatible with the Einstein Equivalence Principle. This result is derived on very general grounds without knowing the underlying structure of the gravitational field equations. As such, any developed gravitational theory built to deal with these astrophysical scales needs to reproduce the obtained results of this article.

A physical description of the formation and propagation of the working surface inside the relativistic jet of the Blazar PKS 1510-089 are used to model its gamma-ray variability light curve using FERMI-LAT data from 2008 to 2012. The physical model is based on conservation laws of mass and momentum at the working surface as explained by Mendoza et al. (2009). The hydrodynamical description of the working surface is parametrised by the initial velocity and mass injection rate at the base of the jet. We show that periodic variations on the injected velocity profiles are able to account for the observed luminosity. With this, we are able to obtain mass ejection rates of the central engine which are injected at the base of the jet, and oscillation frequencies of the flow, amongst other physical parameters.

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.

In this chapter it is shown how the introduction of a fundamental constant of nature with dimensions of acceleration into the theory of gravity makes it possible to extend gravity in a very consistent manner.

Introducing a fundamental constant of nature with dimensions of acceleration into the theory of gravity makes it possible to extend gravity in a very consistent manner. At the non-relativistic level a MOND-like theory with a modification in the force sector is obtained, which is the limit of a very general metric relativistic theory of gravity. Since the mass and length scales involved in the dynamics of the whole universe require small accelerations of the order of Milgrom's acceleration constant a_0, it turns out that the relativistic theory of gravity can be used to explain the expansion of the universe. In this work it is explained how to use that relativistic theory of gravity in such a way that the overall large-scale dynamics of the universe can be treated in a pure metric approach without the need to introduce dark matter and/or dark energy components.

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.

We construct a general relativistic model for the accretion flow of a rotating finite cloud of non-interacting particles infalling onto a Schwarzschild black hole. The streamlines start at a spherical shell, where boundary conditions are fixed, and are followed down to the point at which they either cross the black hole horizon or become incorporated into an equatorial thin disc. Analytic expressions for the streamlines and the velocity field are given, in terms of Jacobi elliptic functions, under the assumptions of stationarity and ballistic motion. A novel approach allows us to describe all of the possible types of orbit with a single formula. A simple numerical scheme is presented for calculating the density field. This model is the relativistic generalisation of the Newtonian one developed by Mendoza, Tejeda, Nagel, 2009 and, due to its analytic nature, it can be useful in providing a benchmark for general relativistic hydrodynamical codes and for exploring the parameter space in applications involving accretion onto black holes when the approximations of steady state and ballistic motion are reasonable ones.

The motion of knots in astrophysical jets is commonly interpreted as shock waves moving along it. Observations of the HST-1 knot during extended periods of time have produced complicated light curves on many wavelengths which are difficult to account using standard hydrodynamical models. Here we reproduce these light curves using the semi-analytical approach developed by Mendoza et al. (2009), developed to reproduce light curves of working surfaces moving along relativistic jets. These working surfaces are generated by periodic oscillations of the injected flow velocity and discharge at the base of the jet. In particular, we use this approach to reproduce the exotic observed features of the light curves of the HST-1 knot in M87. We show that the complicated fits to the light curves are reproduced with high accuracy in all wavelengths. As a bonus, we show that this model is also able to reproduce the light curve of the micro-qsr A06200-00 with high accuracy.

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.

We present heuristic arguments suggesting that if EM waves with wavelengths somewhat larger than the Schwarzschild radius of a black hole were fully absorbed by it, the second law of thermodynamics would be violated, under the Bekenstein interpretation of the area of a black hole as a measure of its entropy. Thus, entropy considerations make the well known fact that large wavelengths are only marginally absorbed by black holes, a natural consequence of thermodynamics. We also study numerically the ingoing radial propagation of a scalar field wave in a Schwarzschild metric, relaxing the standard assumption which leads to the eikonal equation, that the wave has zero spatial extent. We find that if these waves have wavelengths larger that the Schwarzschild radius, they are very substantially reflected, fully to numerical accuracy. Interestingly, this critical wavelength approximately coincides with the one derived from entropy considerations of the EM field, and is consistent with well known limit results of scattering in the Schwarzschild metric. The propagation speed is also calculated and seen to differ from the value $c$, for wavelengths larger than $R_{s}$, in the vicinity of $R_{s}$. As in all classical wave phenomena, whenever the wavelength is larger or comparable to the physical size of elements in the system, in this case changes in the metric, the zero extent 'particle' description fails, and the wave nature becomes apparent.

We construct a steady analytic accretion flow model for a finite rotating gas cloud that accretes material to a central gravitational object. The pressure gradients of the flow are considered to be negligible and so, the flow is ballistic. We also assume a steady flow and consider the particles at the boundary of the spherical cloud to be rotating as a rigid body, with a fixed amount of inwards radial velocity. This represents a generalisation to the traditional infinite gas cloud model described by Ulrich (1976). We show that the streamlines and density profiles obtained deviate largely from the ones calculated by Ulrich. The extra freedom in the choice of the parameters on the model can naturally account for the study of protostars formed in dense clusters by triggered mechanisms, where a wide variety of external physical mechanisms determine the boundary conditions. Also, as expected, the model predicts the formation of an equatorial accretion disc about the central object with a radius different from the one calculated by Ulrich (1976).

We present a ballistic description of the formation and propagation of the working surface of a relativistic jet. Using simple laws of conservation of mass and linear momentum at the working surface, we obtain a full description of the working surface flow parametrised by the initial velocity and mass injection rate. This allows us to compute analytically the energy release at any time in the working surface. We compare this model with the results obtained numerically through a new hydrodynamical code applied to the propagation of a relativistic fluid in one dimension in order to test the limits of our study. Finally, we compare our analytical results with observed light curves of five long gamma ray bursts and show that our model is in very good agreement with observations using simple periodic variations of the injected velocity profiles. This simple method allows us to recover initial mass discharge and energy output ejected during the burst.

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)).

We investigate in detail two physical properties of the metric f(R) theory developed by Sobouti (2007). We first look for the possibility of producing gravitational waves that travel at the speed of light. We then check the possibility of producing extra bending in the lenses produced by the theory. We do this by using standard weak field approximations to the gravitational field equations that appear in Sobouti's theory. We show in this article that the metric theory of gravitation proposed by Sobouti (2007) predicts the existence of gravitational waves travelling at the speed of light in vacuum. In fact, this is proved in general terms for all metric theories of gravity which can be expressed as powers of Ricci's scalar. We also show that an extra additional lensing as compared to the one predicted by standard general relativity is produced. These two points are generally considered to be of crucial importance in the development of relativistic theories of gravity that could provide an alternative description to the dark matter paradigm.

We construct a simple accretion model of a rotating pressureless gas sphere onto a Schwarzschild black hole. We show how to build analytic solutions in terms of Jacobi elliptic functions. This construction represents a general relativistic generalisation of the Newtonian accretion model first proposed by Ulrich (1976). In exactly the same form as it occurs for the Newtonian case, the flow naturally predicts the existence of an equatorial rotating accretion disc about the hole. However, the radius of the disc increases monotonically without limit as the flow reaches its maximum allowed angular momentum.

We show that if the total internal energy of a black hole is constructed as the sum of $N$ photons all having a fixed wavelength chosen to scale with the Schwarzschild radius as $\lambda=\alpha R_{s}$, then $N$ will scale with $R_{s}^{2}$. A statistical mechanical calculation of the configuration proposed yields (\alpha = 4 \pi^2 / \ln(2)) and a total entropy of the system $S=k_{B}N \ln(2)$, in agreement with the Bekenstein entropy of a black hole . It is shown that the critical temperature for Bose-Einstein condensation for relativistic particles of $\lambda=\alpha R_{s}$ is always well below the Hawking temperature of a black hole, in support of the proposed internal configuration. We then examine our results from the point of view of recent loop quantum gravity ideas and find that a natural consistency of both approaches appears. We show that the Jeans criterion for gravitational instability can be generalised to the special and general relativistic regimes and holds for any type of mass--energy distribution.

The idea of a unified model for all astrophysical jets has been considered for quite some time. We present here a few scaling laws relevant to all type of astrophysical jets, analogous to those of \citet{sams96} which are widely used for astrophysical black holes. We use Buckingham's $ \Pi $ theorem of dimensional analysis to obtain a family of dimensional relations among the physical quantities associated to astrophysical jets.

Self-similar solutions to the problem of a strong imploding relativistic shock wave are calculated. These solutions represent the relativistic generalisation of the Newtonian Gouderley-Landau-Stanyukovich problem of a strong imploding spherical shock wave converging to a centre. The solutions are found assuming that the pre-shocked flow has a uniform density, and are accurate for sufficiently large times after the formation of the shock wave.

The idea of a unified model for quasars and micro-quasars has been considered for a long time, despite the different environments and physical conditions where both classes of objects reside. Here we show the existence of a simple scaling law, relating the maximum size of a jet to the properties of the gas medium into which it expands. This appears to be valid for all types of hydrodynamical jets, and can be thought of as a broad unified model. The expansion velocity of the jet and the physical properties of the surrounding gas combine in such a way that a limit to the maximum extent of jets at different scales can be obtained.

The idea of a unified model for all astrophysical jets has been considered for some time now. We present here some hydrodynamical scaling laws relevant for all type of astrophysical jets, analogous to those of Sams et al. (1996). We use Buckingham's Pi theorem of dimensional analysis to obtain a family of dimensional relations among the physical quantities associated with the jets.

We present a ballistic description of the propagation of the working surface of a relativistic jet. Using simple laws of conservation of mass and linear momentum at the working surface, we obtain a full description of the jet flow parametrised by the initial velocity and mass injection. This analysis will soon be applied to particular cases of time-dependent injection of mass and velocity into the jet.

Astronomical observations are suggesting that the fine structure constant varies cosmologically. We present an analysis on the consequences that these variations might induce on the electromagnetic field as a whole. We show that under these circumstances the electrodynamics in vacuum could be described by two fields, the ``standard'' Maxwell's field and a new scalar field. We provide a generalised Lorentz force which can be used to test our results experimentally.

The radio/optical alignment effect for small powerful radio galaxies has been shown to be produced by shock waves formed by the interaction of the head of the jet and/or cocoon with clouds embedded in the interstellar/intergalactic medium. We present here preliminary results of analytical and numerical solutions that have been made to account for the production of implosive shock waves induced by embedding cold clouds in the radio lobe of expanding powerful radio sources.

Molecular clouds in the interstellar medium suffer gravitational instabilities that lead to the formation of one or multiple stars. A recently formed star inside a cold cloud communicates its gravitational force to the surrounding environment and soon an accretion flow falling into the star develops. After their formation, all stars soon eject a wind of gas that interacts with the external accretion flow. This interaction produces a shock wave that evolves with time. The work presented in this article formulates a simple prescription for the evolution of this interaction. With the aid of this model we construct a few radio continuum maps of the source.

We discuss the circumstances under which the bending of a jet can generate an internal shock wave. The analysis is carried out for relativistic and non-relativistic astrophysical jets. The calculations are done by the method of characteristics for the case of steady simple waves. This generalises the non-relativistic treatment first used by Icke (1991). We show that it is possible to obtain an upper limit to the bending angle of a jet in order not to create a shock wave at the end of the curvature. This limiting angle has a value of 75 degrees for non-relativistic jets with a polytropic index k = 4/3, 135 degrees for non-relativistic jets with k = 5/3) and 50 degrees for relativistic jets with k = 5/3. We also discuss under which circumstances jets will form internal shock waves for smaller deflection angles.

A non--relativistic and relativistic model in which astrophysical jets are deflected on passing through an isothermal high density region is analysed. The criteria for the stability of jets due to the formation of internal shocks are discussed.

The model first introduced by Raga & Canto (1996) in which astrophysical jets are deflected on passing through an isothermal high density region is generalised by taking into account gravitational effects on the motion of the jet as it crosses the high density cloud. The problem is also generalised for relativistic jets in which gravitational effects induced by the cloud are neglected. Two further cases, classical and relativistic, are discussed for the cases in which the jet is deflected on passing through the interstellar gas of a galaxy in which a dark matter halo dominates the gravitational potential. The criteria for the stability of jets due to the formation of internal shocks are discussed.