We review some of our recent results about the Radial Acceleration Relation (RAR) and its interpretation as either a fundamental or an emergent law. The former interpretation is in agreement with a class of modified gravity theories that dismiss the need for dark matter in galaxies (MOND in particular). Our most recent analysis, which includes refinements on the priors and the Bayesian test for compatibility between the posteriors, confirms that the hypothesis of a fundamental RAR is rejected at more than 5$\sigma$ from the very same data that was used to infer the RAR.

General relativity (GR) extensions based on renormalization group (RG) flows may lead to scale-dependent couplings with nontrivial effects at large distance scales. Here we develop further the approach in which RG effects at large distance scales are fully encoded in an effective action and we apply it to cosmology. In order to evaluate the cosmological consequences, our main assumption is the use of a RG scale such that the (infrared) RG effects only appear at perturbative order (not at the background level). The emphasis here is on analytical results and qualitative understanding of the implied cosmology. We employ commonly used parametrizations for describing modified gravity in cosmology (as the slip parameter). From them, we describe the dynamics of the first order perturbations and estimate bounds on the single dimensionless parameter ($\nu$) introduced by this framework. Possible impacts on dark matter and dark energy are discussed. It is also shown here that the $\nu$ parameter effects to $f\sigma_8$ are stronger at low redshifts ($z<1.5$), while different values for $\nu$ do not appreciably change $f\sigma_8$ at higher redshifts, thus opening a window to alleviate an issue that is currently faced by $\Lambda$CDM.

We present the first direct-detection search for eV-to-GeV dark matter using a new ~2-gram high-resistivity Skipper-CCD from a dedicated fabrication batch that was optimized for dark-matter searches. Using 24 days of data acquired in the MINOS cavern at the Fermi National Accelerator Laboratory, we measure the lowest rates in silicon detectors of events containing one, two, three, or four electrons, and achieve world-leading sensitivity for a large range of sub-GeV dark matter masses. Data taken with different thicknesses of the detector shield suggest a correlation between the rate of high-energy tracks and the rate of single-electron events previously classified as "dark current." We detail key characteristics of the new Skipper-CCDs, which augur well for the planned construction of the ~100-gram SENSEI experiment at SNOLAB.

The Radial Acceleration Relation (RAR) shows a strong correlation between two accelerations associated to galaxy rotation curves. The relation between these accelerations is given by a nonlinear function which depends on an acceleration scale $a_\dagger$. Some have interpreted this as an evidence for a gravity model, such as Modified Newtonian Dynamics (MOND), which posits a fundamental acceleration scale $a_0$ common to all the galaxies. However, it was later shown, using Bayesian inference, that this seems not to be the case: the $a_0$ credible intervals for individual galaxies were not found to be compatible among themselves. This type of test is a fundamental test for MOND as a theory for gravity, since it directly evaluates its basic assumption and this using the data that most favor MOND: galaxy rotation curves. Here we improve upon the previous analyses by introducing a more robust method to assess the compatibility between the credible intervals, in particular without Gaussian approximations. We directly estimate, using a Monte Carlo simulation, that the existence of a fundamental acceleration is incompatible with the data at more than $5\sigma$. We also consider quality cuts in order to show that our results are robust against outliers. In conclusion, the new analysis further supports the claim that the acceleration scale found in the RAR is an emergent quantity.

Cameron {\it et al.}~2019 (hereafter C19) recommends a more cautious and rigorous approach to statistical analysis in astronomy. We welcome this particular side of their communication as it helps stimulating the effort towards the adoption of better statistical methods in galaxy rotation curves, an effort to which we contributed with Rodrigues {\it et al.}~2018 (hereafter R18). Indeed, R18 was the first work that, in order to conclude on the universality of the acceleration scale $a_0$, studied the posterior distributions on $a_0$ of a large set of galaxies. As C19 agrees, the credible intervals were found within the Bayesian framework, that is, the marginalized posteriors on $a_0$ were found using Bayes' theorem to update the priors in light of the observational data; this process was done without introducing any approximation. Considering R18, C19 also remarks that: i) better methods to select the nuisance parameters and the corresponding priors could be used; ii) a quality cut based on $\chi^2$ values should not be used, and iii) the compatibility of the posteriors should be assessed in a more robust way. In the following, after first clarifying the context of our work, we address these criticisms.

We investigate local environmental effects from dark matter (DM) on thermonuclear supernovae (SNe Ia) using publicly available archival data of 224 low-redshift events, in an attempt to shed light on the SN Ia progenitor systems. SNe Ia are explosions of carbon-oxygen (CO) white dwarfs (WDs) that have recently been shown to explode at sub-Chandrasekhar masses; the ignition mechanism remains, however, unknown. Recently, it has been shown that both weakly interacting massive particles (WIMPs) and macroscopic DM candidates such as primordial black holes (PBHs) are capable of triggering the ignition. Here, we present a method to estimate the DM density and velocity dispersion in the vicinity of SN Ia events and nearby WDs; we argue that (i) WIMP ignition is highly unlikely, and that (ii) DM in the form of PBHs distributed according to a (quasi-) log-normal mass distribution with peak $\log_{10}(m_0/1$g$)=24.9\pm 0.9$ and width $\sigma= 3.3\pm 1.0$ is consistent with SN Ia data, the nearby population of WDs and roughly consistent with other constraints from the literature.

This white paper summarizes the activities of the Brazilian community concerning dark matter physics and highlights the importance of financial support to Brazilian groups that are deeply involved in experimental endeavours. The flagships of the Brazilian dark matter program are the Cherenkov Telescope Array, DARKSIDE, SBN and LHC experiments, but we emphasize that smaller experiments such as DAMIC and CONNIE constitute important probes to dark sectors as well and should receive special attention. Small experimental projects showing the potential to probe new regions of parameter space of dark matter models are encouraged. On the theoretical and phenomenological side, some groups are devoted to astrophysical aspects such as the dark matter density profile while others explore the signature of dark matter models at colliders, direct and indirect detection experiments. In summary, the Brazilian dark matter community that was born not long ago has grown tremendously in the past years and now plays an important role in the hunt for a dark matter particle.

This correspondence answers the criticisms raised in McGaugh et al 2018 and Kroupa et al 2018.

The Radial Acceleration Relation confirms that a nontrivial acceleration scale can be found in the average internal dynamics of galaxies. The existence of such a scale is not obvious as far as the standard cosmological model is concerned, and it has been interpreted as a possible sign of modified gravity. The implications could be profound: it could in principle explain galactic dynamics without large amounts of yet-undetected dark matter and address issues that the standard cosmological model faces at galactic scales. Here, we consider 193 disk galaxies from the SPARC and THINGS databases and, using Bayesian inference, we show that the probability of existence of a fundamental acceleration that is common to all the galaxies is essentially zero: the $p$-value is smaller than $10^{-20}$ or, equivalently, the null hypothesis is rejected at more than 10$\sigma$. We conclude that the acceleration scale unveiled by the Radial Acceleration Relation is of emergent nature, possibly caused by a complex interplay between baryons and dark matter. In particular, the MOND theory, or any other theory that behaves like it at galactic scales, is ruled out as a fundamental theory for galaxies at more than 10$\sigma$.

We develop and apply new techniques in order to uncover galaxy rotation curves (RC) systematics. Considering that an ideal dark matter (DM) profile should yield RCs that have no bias towards any particular radius, we find that the Burkert DM profile satisfies the test, while the Navarro-Frenk-While (NFW) profile has a tendency of better fitting the region between one and two disc scale lengths than the inner disc scale length region. Our sample indicates that this behaviour happens to more than 75% of the galaxies fitted with an NFW halo. Also, this tendency does not weaken by considering "large" galaxies, for instance those with $M_*\gtrsim 10^{10} M_\odot$. Besides the tests on the homogeneity of the fits, we also use a sample of 62 galaxies of diverse types to perform tests on the quality of the overall fit of each galaxy, and to search for correlations with stellar mass, gas mass and the disc scale length. In particular, we find that only 13 galaxies are better fitted by the NFW halo; and that even for the galaxies with $M_* \gtrsim 10^{10} M_\odot$ the Burkert profile either fits as good as, or better than, the NFW profile. This result is relevant since different baryonic effects important for the smaller galaxies, like supernova feedback and dynamical friction from baryonic clumps, indicate that at such large stellar masses the NFW profile should be preferred over the Burkert profile. Hence, our results either suggest a new baryonic effect or a change of the dark matter physics.

The k-essence theory with a power-law function of $(\partial\phi)^2$ and Rastall's non-conservative theory of gravity with a scalar field are shown to have the same solutions for the metric under the assumption that both the metric and the scalar fields depend on a single coordinate. This equivalence (called k-R duality) holds for static configurations with various symmetries (spherical, plane, cylindrical, etc.) and all homogeneous cosmologies. In the presence of matter, Rastall's theory requires additional assumptions on how the stress-energy tensor non-conservation is distributed between different contributions. Two versions of such non-conservation are considered in the case of isotropic spatially flat cosmological models with a perfect fluid: one (R1) in which there is no coupling between the scalar field and the fluid, and another (R2) in which the fluid separately obeys the usual conservation law. In version R1 it is shown that k-R duality holds not only for the cosmological models themselves but also for their adiabatic perturbations. In version R2, among other results, a particular model is singled out that reproduces the same cosmological expansion history as the standard $\Lambda$CDM model but predicts different behaviors of small fluctuations in the k-essence and Rastall frameworks.

A Friedman cosmology is investigated based on scalar-tensor gravitation with general metric coupling and scalar potential functions. We show that for a broad class of such functions, the scalar field can be dynamically trapped using a recently suggested mechanism. The trapped scalar can drive inflation and accelerated cosmic expansion, compatible with standard requirements. The inflationary phase admits a natural exit with a value of the Hubble parameter dictated by the duration of inflation in a parameter independent manner. For inflationary duration consistent with the GUT description, the resulting Hubble parameter is found to be consistent with its observed value.

There are some approaches, either based on General Relativity (GR) or modified gravity, that use galaxy rotation curves to derive the matter density of the corresponding galaxy, and this procedure would either indicate a partial or a complete elimination of dark matter in galaxies. Here we review these approaches, clarify the difficulties on this inverted procedure, present a method for evaluating them, and use it to test two specific approaches that are based on GR: the Cooperstock-Tieu (CT) and the Balasin-Grumiller (BG) approaches. Using this new method, we find that neither of the tested approaches can satisfactorily fit the observational data without dark matter. The CT approach results can be significantly improved if some dark matter is considered, while for the BG approach no usual dark matter halo can improve its results.

Scalar-tensor (ST) gravity is considered in the case where the scalar is an external field. We show that General Relativity (GR) and usual ST gravity are particular cases of the External Scalar-Tensor (EST) gravity. It is shown with a particular cosmological example that it is possible to join a part of a GR solution to a part of a ST one such that the complete solution neither belongs to GR nor to ST, but fully satisfies the EST field equations. We argue that external fields may effectively work as a type of screening mechanism for ST theories.

We review some properties of black hole structures appearing in gravity with a massless scalar field, with both minimal and nonminimal coupling. The main properties of the resulting cold black holes are described. The study of black holes in scalar-gravity systems is extended to $k$-essence theories, and some examples are explicitly worked out. In these cases, even while the existence of horizons is possible, the metric regularity requirement on the horizon implies either a cold black type structure or a singular behavior of the scalar field.

The Brans-Dicke action is one of the most natural extensions of the Einstein-Hilbert action. It is based on the introduction of a fundamental scalar field that effectively incorporates a dynamics to the gravitational coupling $G$. In spite of the diverse motivations and the rich phenomenology that comes from its solutions, Solar System tests impose strong constraints on the Brans-Dicke theory, rendering it indistinguishable from General Relativity. In the present text, new perspectives for the Brans-Dicke theory are presented, based on the possibility that the scalar field presented in the BD theory can be external, as well as on the applications to black hole physics and the primordial universe.

We study the properties of possible static, spherically symmetric configurations in k-essence theories with the Lagrangian functions of the form $F(X)$, $X \equiv \phi_{,\alpha} \phi^{,\alpha}$. A no-go theorem has been proved, claiming that a possible black-hole-like Killing horizon of finite radius cannot exist if the function $F(X)$ is required to have a finite derivative $dF/dX$. Two exact solutions are obtained for special cases of k-essence: one for $F(X) =F_0 X^{1/3}$, another for $F(X) = F_0 |X|^{1/2} - 2 \Lambda$, where $F_0$ and $\Lambda$ are constants. Both solutions contain horizons, are not asymptotically flat, and provide illustrations for the obtained no-go theorem. The first solution may be interpreted as describing a black hole in an asymptotically singular space-time, while in the second solution two horizons of infinite area are connected by a wormhole.

In this paper, we explore two different ways of implementing quantum effects in a classical structure. The first one is through an external field. The other one is modifying the classical conservation laws. In both cases, the consequences for the description of the evolution of the universe are discussed.

The N-body dark matter (DM) simulations point that DM density profiles, e.g. the NFW halo, should be cuspy in its center, but observations disfavour this kind of DM profile. Here we consider whether the observed rotation curves "close" to the galactic centre can favour modified gravity models in comparison to the NFW halo, and how to quantify such difference. Two explicit modified gravity models are considered, MOND and a more recent approach called RGGR (in reference to Renormalization Group effects in General Relativity). It is also the purpose of this work to significantly extend the sample on which RGGR has been tested in comparison to other approaches. By analysing 62 galaxies from five samples, we find that: i) there is a radius, given by half the disk scale length, below which RGGR and MOND can match the data about as well or better than NFW, albeit the formers have fewer free parameters; ii) considering the complete rotation curve data, RGGR could achieve fits with better agreement than MOND, and almost as good as a NFW halo with two free parameters (NFW and RGGR have respectively two and one more free parameters than MOND).

We study, by means of the model proposed in Del Popolo (2009), the effect of baryon physics on the small scale problems of the CDM model. We show that, using this model, the cusp/core problem, the missing satellite problem (MSP), the Too Big to Fail (TBTF) problem, and the angular momentum catastrophe can be reconciled with observations. Concerning the cusp/core problem, the interaction among dark matter (DM) and baryonic clumps of 1% the mass of the halo, through dynamical friction (DF), is able to flatten the inner cusp of the density profiles. We moreover assume that haloes form primarily through quiescent accretion, in agreement with the spherical collapse model (SCM)-secondary infall model (SIM) prescriptions. The results of this paper follow from the two assumptions above. Concerning the MSP and TBTF problem, applying to the Via Lactea II (VL2) subhaloes a series of corrections similar to those of Brooks et al. (2013), namely applying a Zolotov et al. (2012)-like correction obtained with our model, and further correcting for the UV heating and tidal stripping, we obtain that the number of massive, luminous satellites is in agreement with the number observed in the MW. The model also produces an angular momentum distribution in agreement with observations, that is with the distribution of the angular spin parameter and angular momentum of the dwarfs studied by van den Bosch, Burkert, & Swaters (2001). In conclusion, the small scale problems of the CDM model can all be solved by introducing baryon physics.

The Rastall's theory is a modification of General Relativity touching one of the cornestone of gravity theory: the conservation laws. In Rastall's theory, the energy-momentum tensor is not conserved anymore, depending now on the gradient of the Ricci curvature. In this sense, this theory can be seen as a classical implementation of quantum effects in a curved background space-time. We exploit this structure in order to reproduce some results of an effective theory of quantum loop cosmology. Later, we propose a model for the dark sector of the universe. In this case, the corresponding $\Lambda$CDM model appears as the only model consistent with observational data.

We discuss solutions of Vlasov-Einstein equation for collisionless dark matter particles in the context of a flat Friedmann universe. We show that, after decoupling from the primordial plasma, the dark matter phase-space density indicator Q remains constant during the expansion of the universe, prior to structure formation. This well known result is valid for non-relativistic particles and is not "observer dependent" as in solutions derived from the Vlasov-Poisson system. In the linear regime, the inclusion of velocity dispersion effects permits to define a physical Jeans length for collisionless matter as function of the primordial phase-space density indicator: \lambda_J = (5\pi/G)^(1/2)Q^(-1/3)\rho_dm^(-1/6). The comoving Jeans wavenumber at matter-radiation equality is smaller by a factor of 2-3 than the comoving wavenumber due to free-streaming, contributing to the cut-off of the density fluctuation power spectrum at the lowest scales. We discuss the physical differences between these two scales. For dark matter particles of mass equal to 200 GeV, the derived Jeans mass is 4.3 x 10^(-6) solar masses.

We explore the phenomenology of nontrivial quantum effects on low-energy gravity. These effects come from the running of the gravitational coupling parameter G and the cosmological constant L in the Einstein-Hilbert action, as induced by the Renormalization Group (RG). The Renormalization Group corrected General Relativity (RGGR model) is used to parametrize these quantum effects, and it is assumed that the dominant dark matter-like effects inside galaxies is due to these nontrivial RG effects. Here we present additional details on the RGGR model application, in particular on the Poisson equation extension that defines the effective potential, also we re-analyse the ordinary elliptical galaxy NGC 4494 using a slightly different model for its baryonic contribution, and explicit solutions are presented for the running of G and L. The values of the NGC 4494 parameters as shown here have a better agreement with the general RGGR picture for galaxies, and suggest a larger radial anisotropy than the previously published result.

The paper is about possible effects of infrared quantum contributions to General Relativity on disk and elliptical galaxies. The Renormalization Group corrected General Relativity (RGGR model) is used to parametrize these quantum effects. The new RGGR results presented here concern the elliptical galaxy NGC 4374 and the dwarf disk galaxy DDO 47. Using the effective approach to Quantum Field Theory in curved background, one can argue that the proper RG energy scale, in the weak field limit, should be related to the Newtonian potential. In the context of galaxies, this led to a remarkably small variation of the gravitational coupling G, while also capable of generating galaxy rotation and dispersion curves of similar quality to the the best dark matter profiles (i.e., the profiles that have a core).

The renormalization group framework can be applied to Quantum Field Theory on curved space-time, but there is no proof whether the beta-function of the gravitational coupling indeed goes to zero in the far infrared or not. In a recent paper we have shown that the amount of dark matter inside spiral galaxies may be negligible if a small running of the General Relativity coupling G is present. Here we extend the proposed model to elliptical galaxies and present a detailed analysis on the modeling of NGC 4494 (an ordinary elliptical) and NGC 4374 (a giant elliptical). In order to compare our results to a well known alternative model to the standard dark matter picture, we also evaluate NGC 4374 with MOND. In this galaxy MOND leads to a significative discrepancy with the observed velocity dispersion curve and has a significative tendency towards tangential anisotropy. On the other hand, the approach based on the renormalization group and general relativity (RGGR) could be applied with good results to these elliptical galaxies and is compatible with lower mass-to-light ratios (of about the Kroupa IMF type).

We review the difficulties of the generalized Chaplygin gas model to fit observational data, due to the tension between background and perturbative tests. We argue that such issues may be circumvented by means of a self-interacting scalar field representation of the model. However, this proposal seems to be successful only if the self-interacting scalar field has a non-canonical form. The latter can be implemented in Rastall's theory of gravity, which is based on a modification of the usual matter conservation law. We show that, besides its application to the generalized Chaplygin gas model, other cosmological models based on Rastall's theory have many interesting and unexpected new features.

Rastall's theory is based on the non-conservation of the energy-momentum tensor. We show that, in this theory, if we introduce a two-fluid model, one component representing vacuum energy whereas the other pressureless matter (e.g. baryons plus cold dark matter), the cosmological scenario is the same as for the \Lambda CDM model, both at background and linear perturbative levels, except for one aspect: now dark energy may cluster. We speculate that this can lead to a possibility of distinguishing the models at the non-linear perturbative level.

The quantum contributions to the gravitational action are relatively easy to calculate in the higher derivative sector of the theory. However, the applications to the post-inflationary cosmology and astrophysics require the corrections to the Einstein-Hilbert action and to the cosmological constant, and those we can not derive yet in a consistent and safe way. At the same time, if we assume that these quantum terms are covariant and that they have relevant magnitude, their functional form can be defined up to a single free parameter, which can be defined on the phenomenological basis. It turns out that the quantum correction may lead, in principle, to surprisingly strong and interesting effects in astrophysics and cosmology.

We review our contribution to infrared Renormalization Group (RG) effects to General Relativity in the context of galaxies. Considering the effective action approach to Quantum Field Theory in curved background, we argued that the proper RG energy scale, in the weak field limit, should be related to the Newtonian potential. In the galaxy context, even without dark matter, this led to a remarkably small gravitational coupling G variation (about or less than 10^{-12} of its value per light-year), while also capable of generating galaxy rotation curves about as good as the best phenomenological dark matter profiles (considering both the rotation curve shape and the expected mass-to-light ratios). Here we also comment on related developments, open issues and perspectives.

We consider the application of quantum corrections computed using renormalization group arguments in the astrophysical domain and show that, for the most natural interpretation of the renormalization group scale parameter, a gravitational coupling parameter $G$ varying $10^{-7}$ of its value across a galaxy (which is roughly a variation of $10^{-12}$ per light-year) is sufficient to generate galaxy rotation curves in agreement with the observations. The quality of the resulting fit is similar to the Isothermal profile quality once both the shape of the rotation curve and the mass-to-light ratios are considered for evaluation. In order to perform the analysis, we use recent high quality data from nine regular disk galaxies. For the sake of comparison, the same set of data is modeled also for the Modified Newtonian Dynamics (MOND) and for the recently proposed Scalar Tensor Vector Gravity (STVG). At face value, the model based on quantum corrections clearly leads to better fits than these two alternative theories.

We show that a class of bi-gravity theories contain solutions describing dark matter. A particular member of this class is also shown to be equivalent to the Eddington-Born-Infeld gravity, recently proposed as a candidate for dark matter. Bigravity theories also have cosmological de Sitter backgrounds and we find solutions interpolating between matter and acceleration eras.

Recently a Born-Infeld action for dark energy and dark matter that uses additional affine connections was proposed. At background level, it was shown that the new proposal can mimic the standard cosmological evolution. In Bianchi cosmologies, contrary to the scalar field approach (e.g., Chaplygin gas), the new approach leads to anisotropic pressure, raising the issues of stability of the isotropic solution under anisotropic perturbations and, being it stable, how the anisotropies evolve. In this work, the Eddington-Born-Infeld proposal is extended to a Bianchi type I scenario and residual post-inflationary anisotropies are shown to decay in time. Moreover, it is shown that the shears decay following a damped oscillatory pattern, instead of the standard exponential-like decay. Allowing for some fine tuning on the initial conditions, standard theoretical bounds on the shears can be avoided.

There are evidences that the cosmic microwave background (CMB) large-angle anomalies imply a departure from statistical isotropy and hence from the standard cosmological model. We propose a LCDM model extension whose dark energy component preserves its nondynamical character but wield anisotropic vacuum pressure. Exact solutions for the cosmological scale factors are presented, upper bounds for the deformation parameter are evaluated and its value is estimated considering the elliptical universe proposal to solve the quadrupole anomaly. This model can be constructed from a Bianchi I cosmology with cosmological constant from two different ways: i) a straightforward anisotropic modification of the vacuum pressure consistently with energy-momentum conservation; ii) a Poisson structure deformation between canonical momenta such that the dynamics remain invariant under scale factors rescalings.

We present a new method for calculating the merger history of matter halos in hierarchical clustering cosmologies. The linear density field is smoothed on a range of scales, these are then ordered in decreasing density and a merger tree constructed. The method is similar in many respects to the block model of Cole \& Kaiser but has a number of advantages: (i) it retains information about the spatial correlations between halos, (ii) it uses a series of overlapping grids and is thereby much better at finding rare, high-mass halos, (iii) it is not limited to halos whose mass ratios are powers of two, and (iv) it is based on an actual realization of the density field and so can be tested against N-body simulations. The major disadvantages are (i) the minimum halo mass is eight times the unit cell with a corresponding loss of dynamic range, and (ii) occasionally the relative location of halos in the tree does not reflect the correct ordering of their collapse times, as computed from the mean halo density. We show that our model exhibits the required scaling behaviour when tested on power-law spectra of density perturbations, but that it predicts far more massive halos than does the Press-Schechter formalism for flat spectra. We suggest reasons why this should be so.