On June 20, 2020, the LIGO-Virgo collaboration announced the discovery of GW190814, a gravitational wave event originating from a binary system merger between a black hole of mass $M_1 = 23.2^{+1.1} _ {-1.0}M_\odot$ and an unidentified object with a mass of $M_2 = 2.59^{+0.08} _ {-0.09}M_\odot$. This second object would be either the heaviest neutron star or lightest black hole observed to date. Here we investigate the possibility of the $\sim 2.6M_\odot$ object being a primordial black hole (PBH). We find that a primordial black hole explanation to GW190814 is unlikely as it is limited by the formation rate of the primary stellar progenitor and the observed merger rates of $\mathcal{O}(20)M_\odot$ massive black hole pairs.

It has been suggested that late-universe dark matter decays can alleviate the tension between measurements of $H_0$ in the local universe and its value inferred from cosmic microwave background fluctuations. Decaying dark matter can potentially account for this discrepancy as it reshuffles the energy density between matter and radiation and as a result allows dark energy to become dominant at earlier times. We show that the low multipoles amplitude of the cosmic microwave background anisotropy power spectrum severely constrains the feasibility of late-time decays as a solution to the $H_0$ tension.

We study the cosmological effects of two-body dark matter decays where the products of the decay include a massless and a massive particle. We show that if the massive daughter particle is slightly warm it is possible to relieve the tension between distance ladder measurements of the present day Hubble parameter with measurements from the cosmic microwave background.

A thick dark matter disk is predicted in cold dark matter simulations as the outcome of the interaction between accreted satellites and the stellar disk in Milky Way sized halos. We study the effects of a self-interacting thick dark disk on the energetic neutrino flux from the Sun. We find that for particle masses between 100 GeV and 1 TeV and dark matter annihilation to heavy leptons either the self-interaction may not be strong enough to solve the small scale structure motivation or a dark disk cannot be present in the Milky Way.

We demonstrate that the Vlasov equation describing collisionless self-gravitating matter may be solved with the so-called Schr\"odinger method (ScM). With the ScM, one solves the Schr\"odinger-Poisson system of equations for a complex wave function in d dimensions, rather than the Vlasov equation for a 2d-dimensional phase space density. The ScM also allows calculating the d-dimensional cumulants directly through quasi-local manipulations of the wave function, avoiding the complexity of 2d-dimensional phase space. We perform for the first time a quantitive comparison of the ScM and a conventional Vlasov solver in d=2 dimensions. Our numerical tests were carried out using two types of cold cosmological initial conditions: the classic collapse of a sine wave and those of a gaussian random field as commonly used in cosmological cold dark matter N-body simulations. We compare the first three cumulants, that is, the density, velocity and velocity dispersion, to those obtained by solving the Vlasov equation using the publicly available code ColDICE. We find excellent qualitative and quantitative agreement between these codes, demonstrating the feasibility and advantages of the ScM as an alternative to N-body simulations. We discuss, the emergence of effective vorticity in the ScM through the winding number around the points where the wave function vanishes. As an application we evaluate the background pressure induced by the non-linearity of large scale structure formation, thereby estimating the magnitude of cosmological backreaction. We find that it is negligibly small and has time dependence and magnitude compatible with expectations from the effective field theory of large scale structure.