The large dark cores of common dwarf galaxies are unexplained by the standard heavy particle interpretation of dark matter. This puzzle is exacerbated by the discovery of a very large but barely visible, dark matter dominated galaxy Antlia II orbiting the Milky Way, uncovered by tracking star motions with the {\t Gaia} satellite. Although Antlia II has a low mass, its visible radius is more than double any known dwarf galaxy, with an unprecedentedly low density core. We show that Antlia II favors dark matter as a Bose-Einstein condensate, for which the ground state is a stable soliton with a core radius given by the de Broglie wavelength. The lower the galaxy mass, the larger the de Broglie wavelength, so the least massive galaxies should have the widest soliton cores of lowest density. An ultra-light boson of $m_\psi \sim 1.1 \times10^{-22}$ eV, accounts well for the large size and slowly moving stars within Antlia II, and agrees with boson mass estimates derived from the denser cores of more massive dwarf galaxies. For this very light boson, Antlia II is close to the lower limiting Jeans scale for galaxy formation permitted by the Uncertainty Principle, so other examples are expected but none significantly larger in size. This simple explanation for the puzzling dark cores of dwarf galaxies implies dark matter as an ultra-light boson, such as an axion generic in String Theory.

Dark matter as a Bose-Einstein condensate, such as the axionic scalar field particles of String Theory, can explain the coldness of dark matter on large scales. Pioneering simulations in this context predict a rich wave-like structure, with a ground state soliton core in every galaxy surrounded by a halo of excited states that interfere on the de Broglie scale. This de Broglie scale is largest for low mass galaxies as momentum is lower, providing a simple explanation for the wide cores of dwarf spheroidal galaxies. Here we extend these "wave dark matter" ($\psi$DM) predictions to the newly discovered class of "Ultra Diffuse Galaxies" (UDG) that resemble dwarf spheroidal galaxies but with more extended stellar profiles. Currently the best studied example, DF44, has a uniform velocity dispersion of $\simeq 33$km/s, extending to a least 3 kpc, that we show is reproduced by our $\psi$DM simulations with a soliton radius of $\simeq 0.5$ kpc. In the $\psi$DM context, we show the relatively flat dispersion profile of DF44 lies between massive galaxies with compact dense solitons, as may be present in the Milky Way on a scale of 100pc and lower mass galaxies where the velocity dispersion declines centrally within a wide, low density soliton, like Antlia II, of radius 3Kpc. In contrast, standard CDM requires excessive tangential stellar motions predominate to counter the inherent central cusp which would otherwise enhance the central velocity dispersion.

The lack of detection of supersymmetric particles is leading to look at alternative avenues for explaining dark matter's effects. Among them, modified theories of gravity may play an important role accounting even for both dark components needed in the standard cosmological model. Scalar-Tensor-Vector Gravity theory has been proposed to resolve the dark matter puzzle. Such a modified gravity model introduces, in its weak field limit, a Yukawa-like correction to the Newtonian potential, and is capable to explain most of the phenomenology related to dark matter at scale of galaxies and galaxy clusters. Nevertheless, some inconsistencies appears when studying systems that are supposed to be dark matter dominated such as dwarf galaxies. In this sense Antlia II, an extremely diffuse galaxy which has been recently discovered in {\em Gaia}'s second data release, may serve to probe the aforementioned theory against the need for invoking dark matter. Our analysis shows several inconsistencies and leads to argue that MOdified Gravity may not be able to shed light on the intriguing nature of dark matter.

A wavelike solution for the non-relativistic universal dark matter (wave-DM) is rapidly gaining interest, following distinctive predictions of pioneering simulations of cosmic structure as an interference pattern of coherently oscillating bosons. A prominent solitonic standing wave is predicted at the center of every galaxy, representing the ground state, that has been identified with the wide, kpc scale dark cores of common dwarf-spheroidal galaxies, providing a boson mass of, $\simeq 10^{-22}$ eV. A denser soliton is predicted for Milky Way sized galaxies where momentum is higher, so the de Broglie scale of the soliton is smaller, $\simeq 100$ pc, of mass $\simeq 10^9 M_\odot$. Here we show the central motion of bulge stars in the Milky Way implies the presence of such a dark core, where the velocity dispersion rises inversely with radius to a maximum of $\simeq 130$ km/s, corresponding to an excess central mass of $\simeq 1.5\times 10^9 M_\odot$ within $\simeq 100$ pc, favouring a boson mass of $\simeq 10^{-22}$ eV. This quantitative agreement with such a unique and distinctive prediction is therefore strong evidence for a light bosonic solution to the long standing Dark Matter puzzle, such as the axions generic in String Theory.

We analyze axion--photon mixing in the framework of quantum field theory. The condensate structure of the vacuum for mixed fields induces corrections to the oscillation formulae and leads to non-zero energy of the vacuum for the component of the photon mixed with the axion. This energy generates a new effect of the vacuum polarization and it has the state equation of the cosmological constant, $w = -1$. This result holds for any homogeneous and isotropic curved space-time, as well as for diagonal metrics. Numerical estimates of the corrections to the oscillation formulae are presented by considering the intensity of the magnetic field available in the laboratory. Moreover, we estimate the vacuum energy density induced by axion--photon mixing in the Minkowski space-time. A value compatible with that of the energy density of the universe can be obtained for axions with a mass of $(10^{-3}-10^{-2}) eV$ in the presence of the strong magnetic fields that characterize astrophysical objects such as pulsars or neutron stars. In addition, a value of the energy density less than that of the Casimir effect is obtained for magnetic fields used in experiments such as PVLAS. The vacuum polarization induced by this energy could be detected in next experiments and it might provide an indirect proof of the existence of the axion--photon mixing. The quantum field theory effects presented in this work may lead to new methods for studying axion-like particles.

The study of the dynamics of a two-body system in modified gravity constitutes a more complex problem than in Newtonian gravity. Numerical methods are typically needed to solve the equations of geodesics. Despite the complexity of the problem, the study of a two-body system in $f(R)$ gravity leads to a new exciting perspective hinting the right strategy to adopt in order to probe modified gravity. Our results point out some differences between the {\em semiclassical} (Newtonian) approach, and the {\em relativistic} (geodesic) one thus suggesting that the latter represents the best strategy for future tests of modified theories of gravity. { Finally, we have also highlighted the capability of forthcoming observations to serve as smoking gun of modified gravity revealing a departure from GR or further reducing the parameter space of $f(R)$ gravity}. \keywords{$f(R)$ gravity \and binary system \and geodesics \and precession

Decaying Dark Energy models modify the background evolution of the most common observables, such as the Hubble function, the luminosity distance and the Cosmic Microwave Background temperature-redshift scaling relation. We use the most recent observationally-determined datasets, including Supernovae Type Ia and Gamma Ray Bursts data, along with $H(z)$ and Cosmic Microwave Background temperature versus $z$ data and the reduced Cosmic Microwave Background parameters, to improve the previous constraints on these models. We perform a Monte Carlo Markov Chain analysis to constrain the parameter space, on the basis of two distinct methods. In view of the first method, the Hubble constant and the matter density are left to vary freely. In this case, our results are compatible with previous analyses associated with decaying Dark Energy models, as well as with the most recent description of the cosmological background. In view of the second method, we set the Hubble constant and the matter density to their best fit values obtained by the {\it Planck} satellite, reducing the parameter space to two dimensions, and improving the existent constraints on the model's parameters. Our results suggest that the accelerated expansion of the Universe is well described by the cosmological constant, and we argue that forthcoming observations will play a determinant role to constrain/rule out decaying Dark Energy.

We find the generalized version of the Toomre's criterion for the stability of a rotating thin disk in the context of Eddington inspired Born-Infeld (EiBI) gravity which possesses one free parameter $\chi$. To do so we use the weak field limit of the theory and find the dispersion relation for the propagation of matter density waves on the surface of a self-gravitating and differentially rotating disk. Finally we find a new version of Toomre's stability criterion for thin disks. We show that EiBI gravity with negative $\chi$ destabilizes all the rotating thin disks. On the other hand EiBI with positive $\chi$ substantially can suppress the local fragmentation, and has stabilizing effects against axi-symmetric perturbations. More specifically, we show that only an annulus remains unstable on the surface of the disk. The width of the annulus directly depends on the magnitude of $\chi$.

Alternative theories of gravity may serve to overcame several shortcomings of the standard cosmological model but, in their weak field limit, General Relativity must be recovered so as to match the tight constraints at the Solar System scale. Therefore, testing such alternative models at scales of stellar systems could give a unique opportunity to confirm or rule them out. One of the most straightforward modifications is represented by analytical $f(R)$-gravity models that introduce a Yukawa-like modification to the Newtonian potential thus modifying the dynamics of particles. Using the geodesics equations, we have illustrated the amplitude of these modifications. First, we have integrated numerically the equations of motion showing the orbital precession of a particle around a massive object. Second, we have computed an analytic expression for the periastron advance of systems having their semi-major axis much shorter than the Yukawa-scale length. Finally, we have extended our results to the case of a binary system composed of two massive objects. Our analysis provides a powerful tool to obtain constraints on the underlying theory of gravity using current and forthcoming datasets.

Light Axionic Dark Matter, motivated by string theory, is increasingly favored for the "no-WIMP era". Galaxy formation is suppressed below a Jeans scale, of $\simeq 10^8 M_\odot$ by setting the axion mass to, $m_B \sim 10^{-22}$eV, and the large dark cores of dwarf galaxies are explained as solitons on the de-Broglie scale. This is persuasive, but detection of the inherent scalar field oscillation at the Compton frequency, $\omega_B= (2.5{\rm \, months})^{-1}(m_B/10^{-22}eV)$, would be definitive. By evolving the coupled Schr\"odinger-Poisson equation for a Bose-Einstein condensate, we predict the dark matter is fully modulated by de-Broglie interference, with a dense soliton core of size $\simeq 150pc$, at the Galactic center. The oscillating field pressure induces General Relativistic time dilation in proportion to the local dark matter density and pulsars within this dense core have detectably large timing residuals, of $\simeq 400nsec/(m_B/10^{-22}eV)$. This is encouraging as many new pulsars should be discovered near the Galactic center with planned radio surveys. More generally, over the whole Galaxy, differences in dark matter density between pairs of pulsars imprints a pairwise Galactocentric signature that can be distinguished from an isotropic gravitational wave background.

We analyze the stability of self-gravitating systems which dynamics is investigated using the collisionless Boltzmann equation, and the modified Poisson equation of Eddington-inspired Born-Infield gravity. These equations provide a description of the Jeans paradigm used to determine the critical scale above which such systems collapse. At equilibrium, the systems are described using the time-independent Maxwell- Boltzmann distribution function $f_0(v)$. Considering small perturbations to this equilibrium state, we obtain a modified dispersion relation, and we find a new characteristic scale length. Our results indicate that the dynamics of the self-gravitating astrophysical systems can be fully addressed in the Eddington-inspired Born-Infield gravity. The latter modifies the Jeans instability in high densities environments while its effects become negligible in the star formation regions.

Using Planck satellite data, we construct SZ gas pressure profiles for a large, volume-complete sample of optically selected clusters. We have defined a sample of over 8,000 redMaPPer clusters from the Sloan Digital Sky Survey (SDSS), within the volume-complete redshift region 0.100 < z < 0.325, for which we construct Sunyaev-Zel'dovich (SZ) effect maps by stacking Planck data over the full range of richness. Dividing the sample into richness bins we simultaneously solve for the mean cluster mass in each bin together with the corresponding radial pressure profile parameters, employing an MCMC analysis. These profiles are well detected over a much wider range of cluster mass and radius than previous work, showing a clear trend towards larger break radius with increasing cluster mass. Our SZ-based masses fall ~24% below the mass-richness relations from weak lensing, in a similar fashion as the "hydrostatic bias" related with X-ray derived masses. We correct for this bias to derive an optimal mass-richness relation finding a slope 1.22 +/- 0.04 and a pivot mass log(M_500/M_0)= 14.432 +/- 0.041, evaluated at a richness lambda=60. Finally, we derive a tight Y_500-M_500 relation over a wide range of cluster mass, with a power law slope equal to 1.72 +/- 0.07, that agrees well with the independent slope obtained by the Planck team with an SZ-selected cluster sample, but extends to lower masses with higher precision.

In its weak field limit, Scalar-tensor-vector gravity theory introduces a Yukawa-correction to the gravitational potential. Such a correction depends on the two parameters, $\alpha$ which accounts for the modification of the gravitational constant, and $\mu^{*-1}$ which represents the scale length on which the scalar field propagates. These parameters were found to be universal when the modified gravitational potential was used to fit the galaxy rotation curves and the mass profiles of galaxy clusters, both without Dark Matter. We test the universality of these parameters using the the temperature anisotropies due to the thermal Sunyaev-Zeldovich effect. In our model the intra-cluster gas is in hydrostatic equilibrium within the modified gravitational potential well and it is described by a polytropic equation of state. We predict the thermal Sunyaev-Zeldovich temperature anisotropies produced by Coma cluster, and we compare them with those obtained using the Planck 2013 Nominal maps. In our analysis, we find $\alpha$ and the scale length, respectively, to be consistent and to depart from their universal values. Our analysis points out that the assumption of the universality of the Yukawa-correction to the gravitational potential is ruled out at more than $3.5\sigma$ at galaxy clusters scale, while demonstrating that such a theory of gravity is capable to fit the cluster profile if the scale dependence of the gravitational potential is restored.

We have constrained the spatial variation of the fine structure constant using multi-frequency measurements of the thermal Sunyaev-Zeldovich effect of 618 X-ray selected clusters. Although our results are not competitive with the ones from quasar absorption lines, we improved by a factor 10 and ~2.5 previous results from Cosmic Microwave Background power spectrum and from galaxy clusters, respectively.

We propose an improved methodology to constrain spatial variations of the fine structure constant using clusters of galaxies. We use the {\it Planck} 2013 data to measure the thermal Sunyaev-Zeldovich effect at the location of 618 X-ray selected clusters. We then use a Monte Carlo Markov Chain algorithm to obtain the temperature of the Cosmic Microwave Background at the location of each galaxy cluster. When fitting three different phenomenological parameterizations allowing for monopole and dipole amplitudes in the value of the fine structure constant we improve the results of earlier analysis involving clusters and the CMB power spectrum, and we also found that the best-fit direction of a hypothetical dipole is compatible with the direction of other known anomalies. Although the constraining power of our current datasets do not allow us to test the indications of a fine-structure constant dipole obtained though high-resolution optical/UV spectroscopy, our results do highlight that clusters of galaxies will be a very powerful tool to probe fundamental physics at low redshift.

We determine the relation between the Comptonization parameter predicted using X-ray data $Y_{C,Xray}$ and the X-ray luminosity $L_X$, both magnitudes derived from ROSAT data, with the Comptonization parameter $Y_{C,SZ}$ measured on {\it Planck} 2013 foreground cleaned Nominal maps. The 560 clusters of our sample includes clusters with masses $M\ge 10^{13}M_\odot$, one order of magnitude smaller than those used by the Planck Collaboration in a similar analysis. It also contains eight times more clusters in the redshift interval $z\le 0.3$. The prediction of the $\beta=2/3$ model convolved with the Planck antenna beam agrees with the anisotropies measured in foreground cleaned Planck Nominal maps within the X-ray emitting region, confirming the results of an earlier analysis (Atrio-Barandela et al. 2008). The universal pressure profile overestimates the signal by a 15-21\% depending on the angular aperture. We show that the discrepancy is not due to the presence of {\it cool-core} systems but it is an indication of a brake in the $L_X-M$ relation towards low mass systems. We show that relation of the Comptonization parameter averaged over the region that emits 99\% of the X-ray flux and and the X-ray luminosity is consistent with the predictions of the self-similar model. We confirm previous findings that the scaling relations studied here do not evolve with redshift within the range probed by our catalog.

In the weak field limit, analytic $f(R)$ models of gravity introduce a Yukawa-like correction to the Newtonian gravitational potential. These models have been widely tested at galactic scales and provide an alternative explanation to the dynamics of galaxies without Dark Matter. We study if the temperature anisotropies due to the thermal Sunyaev-Zeldovich effect are compatible with these Extended Theories of Gravity. We assume that the gas is in hydrostatic equilibrium within the modified Newtonian potential and it is well described by a polytropic equation of state. We particularize the model for the Coma cluster and the predicted anisotropies are compared with those measured in the foreground cleaned maps obtained using the Planck Nominal maps released in 2013. We show that the computed $f(R)$ pressure profile fits the data giving rise to competitive constraints of the Yukawa scale length $L=(2.19\pm1.02) \rm{\, Mpc}$, and of the deviation parameter $ \delta=-0.48\pm0.22$. Those are currently the tightest constraints at galaxy cluster scale, and support the idea that Extended Theories of Gravity provide an alternative explanation to the dynamics of self-gravitating systems without requiring Dark Matter.

Over the past decades, General Relativity and the concordance $\Lambda$CDM model have been successfully tested using several different astrophysical and cosmological probes based on large datasets ({\it precision cosmology}). Despite their successes, some shortcomings emerge due to the fact that General Relativity should be revised at infrared and ultraviolet limits and to the fact that the fundamental nature of Dark Matter and Dark Energy is still a puzzle to be solved. In this perspective, $f(R)$ gravity have been extensively investigated being the most straightforward way to modify General Relativity and to overcame some of the above shortcomings. In this paper, we review various aspects of $f(R)$ gravity at extragalactic and cosmological levels. In particular, we consider cluster of galaxies, cosmological perturbations, and N-Body simulations, focusing on those models that satisfy both cosmological and local gravity constraints. The perspective is that some classes of $f(R)$ models can be consistently constrained by Large Scale Structure.

We constrain the deviation of adiabatic evolution of the Universe using the data on the Cosmic Microwave Background (CMB) temperature anisotropies measured by the {\it Planck} satellite and a sample of 481 X-ray selected clusters with spectroscopically measured redshifts. To avoid antenna beam effects, we bring all the maps to the same resolution. We use a CMB template to subtract the cosmological signal while preserving the Thermal Sunyaev-Zeldovich (TSZ) anisotropies; next, we remove galactic foreground emissions around each cluster and we mask out all known point sources. If the CMB black-body temperature scales with redshift as $T(z)=T_0(1+z)^{1-\alpha}$, we constrain deviations of adiabatic evolution to be $\alpha=-0.007\pm 0.013$, consistent with the temperature-redshift relation of the standard cosmological model. This result could suffer from a potential bias $\delta\alpha$ associated with the CMB template, that we quantify it to be $|\delta\alpha|\le 0.02$ and with the same sign than the measured value of $\alpha$, but is free from those biases associated with using TSZ selected clusters; it represents the best constraint to date of the temperature-redshift relation of the Big-Bang model using only CMB data, confirming previous results.

Analytical $f(R)$-gravity models introduce Yukawa-like corrections to the Newtonian potential in the weak field limit. These models can explain the dynamics of galaxies and cluster of galaxies without requiring dark matter. To test the model, we have computed the pressure profile of 579 X-ray galaxy clusters assuming the gas is in hydrostatic equilibrium within the potential well of the modified gravitational potential. We have compared those profiles with the ones measured in the foreground cleaned SMICA released by the Planck Collaboration. Our results show that Extended Theories of Gravity explain the dynamics of self-gravitating systems at cluster scales and represent an alternative to dark matter haloes.

There are several approaches to extend General Relativity in order to explain the phenomena related to the Dark Matter and Dark Energy. These theories, generally called Extended Theories of Gravity, can be tested using observations coming from relativistic binary systems as PSR $J0348+0432$. Using a class of analytical $f(R)$-theories, one can construct the first time derivative of orbital period of the binary systems starting from a quadrupolar gravitational emission. Our aim is to set boundaries on the parameters of the theory in order to understand if they are ruled out, or not, by the observations on PSR $J0348+0432$. Finally, we have computed an upper limit on the graviton mass showing that agree with constraint coming from other observations.

Models of $f(R)$ gravity that introduce corrections to the Newtonian potential in the weak field limit are tested at the scale of galaxy clusters. These models can explain the dynamics of spiral and elliptical galaxies without resorting to dark matter. We compute the pressure profiles of 579 galaxy clusters assuming that the gas is in hydrostatic equilibrium within the potential well of the modified gravitational field. The predicted profiles are compared with the average profile obtained by stacking the data of our cluster sample in the Planck foreground clean map SMICA. We find that the resulting profiles of these systems fit the data without requiring a dominant dark matter component, with model parameters similar to those required to explain the dynamics of galaxies. Our results do not rule out that clusters are dynamically dominated by Dark Matter but support the idea that Extended Theories of Gravity could provide an explanation to the dynamics of self-gravitating systems and to the present period of accelerated expansion, alternative to the concordance cosmological model.

We study the power of PLANCK data to constrain deviations of the Cosmic Microwave Background black body temperature from adiabatic evolution using the thermal Sunyaev-Zeldovich anisotropy induced by clusters of galaxies. We consider two types of data sets: the cosmological signal is removed in the Time Ordered Information or is removed from the final maps; and two different statistical estimators, based on the ratio of temperature anisotropies at two different frequencies and on a fit to the spectral variation of the cluster signal with frequency. To test for systematics, we construct a template from clusters drawn from a hydro-simulation included in the pre-launch Planck Sky Model. We demonstrate that, using a proprietary catalog of X-ray selected clusters with measured redshifts, electron densities and X-ray temperatures, we can constrain deviations of adiabatic evolution, measured by the parameter $\alpha$ in the redshift scaling $T(z)=T_0(1+z)^{1-\alpha}$, with an accuracy of $\sigma_\alpha=0.011$ in the most optimal case and with $\sigma_\alpha=0.016$ for a less optimal case. These results represent a factor 2-3 improvement over similar measurements carried out using quasar spectral lines and a factor 6-20 with respect to earlier results using smaller cluster samples.

Dynamics and collapse of collisionless self-gravitating systems is described by the coupled collisionless Boltzmann and Poisson equations derived from $f(R)$-gravity in the weak field approximation. Specifically, we describe a system at equilibrium by a time-independent distribution function $f_0(x,v)$ and two potentials $\Phi_0(x)$ and $\Psi_0(x)$ solutions of the modified Poisson and collisionless Boltzmann equations. Considering a small perturbation from the equilibrium and linearizing the field equations, it can be obtained a dispersion relation. A dispersion equation is achieved for neutral dust-particle systems where a generalized Jeans wave-number is obtained. This analysis gives rise to unstable modes not present in the standard Jeans analysis (derived assuming Newtonian gravity as weak filed limit of $f(R)=R$). In this perspective, we discuss several self-gravitating astrophysical systems whose dynamics could be fully addressed in the framework of $f(R)$-gravity.

Gravitational waves detected from well-localized inspiraling binaries would allow us to determine, directly and independently, binary luminosity and redshift. In this case, such systems could behave as "standard candles" providing an excellent probe of cosmic distances up to z <0.1 and complementing other indicators of cosmological distance ladder.

Gravitational waves detected from well-localized inspiraling binaries would allow to determine, directly and independently, both binary luminosity and redshift. In this case, such systems could behave as "standard candles" providing an excellent probe of cosmic distances up to z < 0.1 and thus complementing other indicators of cosmological distance ladder.

We discuss the possibility to obtain an electromagnetic emission accompanying the gravitational waves emitted in the coalescence of a compact binary system. Motivated by the existence of black hole configurations with open magnetic field lines along the rotation axis, we consider a magnetic dipole in the system, the evolution of which leads to (i) electromagnetic radiation, and (ii) a contribution to the gravitational radiation, the luminosity of both being evaluated. Starting from the observations on magnetars, we impose upper limits for both the electromagnetic emission and the contribution of the magnetic dipole to the gravitational wave emission. Adopting this model for the evolution of neutron star binaries leading to short gamma ray bursts, we compare the correction originated by the electromagnetic field to the gravitational waves emission, finding that they are comparable for particular values of the magnetic field and of the orbital radius of the binary system. Finally we calculate the electromagnetic and gravitational wave energy outputs which result comparable for some values of magnetic field and radius.

Coalescing binary systems, consisting of two collapsed objects, are among the most promising sources of high frequency gravitational waves signals detectable, in principle, by ground-based interferometers. Binary systems of Neutron Star or Black Hole/Neutron Star mergers should also give rise to short Gamma Ray Bursts, a subclass of Gamma Ray Bursts. Short-hard-Gamma Ray Bursts might thus provide a powerful way to infer the merger rate of two-collapsed object binaries. Under the hypothesis that most short Gamma Ray Bursts originate from binaries of Neutron Star or Black Hole/Neutron Star mergers, we outline here the possibility to associate short Gamma Ray Bursts as electromagnetic counterpart of coalescing binary systems.