We report on the impact of a probabilistic prescription for compact remnant masses and spins on massive binary population synthesis. We find that this prescription populates the putative mass gap between neutron stars and black holes with low-mass black holes. However, evolutionary effects reduce the number of X-ray binary candidates with low-mass black holes, consistent with the dearth of such systems in the observed sample. We further find that this prescription is consistent with the formation of heavier binary neutron stars such as GW190425, but over-predicts the masses of Galactic double neutron stars. The revised natal kicks, particularly increased ultra-stripped supernova kicks, do not directly explain the observed Galactic double neutron star orbital period--eccentricity distribution. Finally, this prescription allows for the formation of systems similar to the recently discovered extreme mass ratio binary GW190814, but only if we allow for the survival of binaries in which the common envelope is initiated by a donor crossing the Hertzsprung gap, contrary to our standard model.

Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

The interaction of a supernova with a circumstellar medium (CSM) can dramatically increase the emitted luminosity by converting kinetic energy to thermal energy. In 'superluminous' supernovae (SLSNe) of Type IIn -- named for narrow hydrogen lines in their spectra -- the integrated emission can reach $\sim 10^{51}$ erg, attainable by thermalising most of the kinetic energy of a conventional SN. A few transients in the centres of active galaxies have shown similar spectra and even larger energies, but are difficult to distinguish from accretion onto the supermassive black hole. Here we present a new event, SN2016aps, offset from the centre of a low-mass galaxy, that radiated $\gtrsim 5 \times 10^{51}$ erg, necessitating a hyper-energetic supernova explosion. We find a total (SN ejecta $+$ CSM) mass likely exceeding 50-100 M$_\odot$, with energy $\gtrsim 10^{52}$ erg, consistent with some models of pair-instability supernovae (PISNe) or pulsational PISNe -- theoretically-predicted thermonuclear explosions from helium cores $>50$ M$_\odot$. Independent of the explosion mechanism, this event demonstrates the existence of extremely energetic stellar explosions, detectable at very high redshifts, and provides insight into dense CSM formation in the most massive stars.

As we enter the era of gravitational wave astronomy, we are beginning to collect observations which will enable us to explore aspects of astrophysics of massive stellar binaries which were previously beyond reach. In this paper we describe COMPAS (Compact Object Mergers: Population Astrophysics and Statistics), a new platform to allow us to deepen our understanding of isolated binary evolution and the formation of gravitational-wave sources. We describe the computational challenges associated with their exploration, and present preliminary results on overcoming them using Gaussian process regression as a simulation emulation technique.