We consider some observational consequences of replacing all black holes (BHs) with a class of non-singular solutions that mimic BHs but with Dark Energy (DE) interiors; GEneric Objects of Dark Energy (GEODEs). We focus on the BH mass function and chirp-mass redshift distribution of mergers visible to gravitational wave observatories. We incorporate the GEODE blueshift into an initially Salpeter stellar remnant distribution, and model the binary population by evolving synthesized binary remnant distributions, published before LIGO's first measurements. We find that a GEODE produced between $20 < z < 40$, and observed at $z \sim 7$, will have its initial mass amplified by $18-134\times$. This can relieve tension between accretion-only growth models and the inferred masses of BHs in quasars at $z \gtrsim 6$. Moreover, we find that merger rates of GEODE binaries increase by a factor of $\sim 2\times$ relative to classical BHs. The resulting GEODE mass function is consistent with the most recent LIGO constraints at $< 0.5\sigma$. In contrast, a Salpeter stellar distribution that evolves into classical remnants is in tension at $\gtrsim 2\sigma$. This agreement occurs without low-metallicity regions, abnormally large progenitor stars, novel binary formation channels, or primordial object formation at extreme rates. In particular, we find that solar metallicity progenitors, which produce $1.1-1.8\mathrm{M}_\odot$ remnants, overlap with many LIGO observations when evolved as GEODEs.

The hybrid Ricci scalar $R_{\mu\nu}\eta^{\mu\nu}$ formed from a Robertson-Walker derived Ricci tensor and a distinct flat space metric becomes a dark, pseudo-isothermal density under asymptotic constraint. While extension of this constraint to more general scenarios is unclear, it is model independent and suggestive of a radical premise: all galactic dark matter halos are sourced by very many distinct Robertson-Walker spacetimes. We explore this premise and present a cored galactic halo, which encloses finite mass over all space. We compute fit parameters to published data on the Milky Way dwarf spheroidal galaxies, and find consistency with a predicted reciprocal square scaling relation. Spherical symmetry permits treatment of each individual halo as a point particle interacting under a modified gravitational force. This force is Newtonian at large separation with a leading dipole correction, and naturally softens to zero at small separation. Taken as cold dark matter, ensembles of such halos present a possible resolution to tension between large-scale $N$-body simulations and small-scale astrophysical data.

Gravitational vacuum stars (gravastars) have become a viable theoretical alternative to the black hole (BH) end-stage of stellar evolution. These objects gravitate in vacuum like BHs, yet have no event horizons. In this paper, we present tests of the gravastar hypothesis within flat Friedmann cosmology. Such tests are complementary to optical and gravitational wave merger signatures, which have uncertainties dominated by the poorly constrained gravastar crust. We motivate our analysis directly from the action principle, and show that a population of gravastars must induce a time-dependent dark energy (DE). The possibility of such a contribution has been overlooked due to a subtlety in the de facto definition of the isotropic and homogeneous stress. Conservation of stress-energy directly relates the time evolution of this DE contribution to the measured BH comoving mass function and a gravastar population crust parameter $\eta$. We show that a population of gravastars formed between redshift $8 \lesssim z \lesssim 20$ readily produces the present-day DE density over a large range of $\eta$, providing a natural resolution to the coincidence problem. We compute an effective DE equation of state $w_\mathrm{eff}(a)$ as a function of present-epoch BH population observables and $\eta$. Using a BH population model developed from the cosmic star formation history, we obtain $w_0$ and $w_a$ consistent with Planck best-fit values. In summary, the gravastar hypothesis leads to an unexpected correlation between the BH population and the magnitude and time-evolution of DE.

We discuss an extension of the massively parallel cosmological simulation code GADGET-2, which now enables investigation of multiple and distinct gravitational force laws, provided they are dominated by a constant scaling of the Newtonian force. In addition to simplifying investigations of a universally modified force law, the ngravs extension provides a foundation for state-of-the-art collisionless cosmological simulations of exotic gravitational scenarios. We briefly review the algorithms used by GADGET-2, and present our extension to multiple gravities, highlighting additional features that facilitate consideration of exotic force laws. We discuss the accuracy and performance of the ngravs extension, both internally and with an unaltered GADGET-2, in the relevant operational modes. The ngravs extension is publicly released to the research community.