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.

We have implemented prescriptions for modelling pulsars in the rapid binary population synthesis code COMPAS. We perform a detailed analysis of the double neutron star (DNS) population, accounting for radio survey selection effects. The surface magnetic field decay timescale (${\sim}1000$\,Myr) and mass scale (${\sim}0.02$\,M$_\odot$) are the dominant uncertainties in our model. Mass accretion during common envelope evolution plays a non-trivial role in recycling pulsars. We find a best-fit model that is in broad agreement with the observed Galactic DNS population. Though the pulsar parameters (period and period derivative) are strongly biased by radio selection effects, the observed orbital parameters (orbital period and eccentricity) closely represent the intrinsic distributions. The number of radio observable DNSs in the Milky Way at present is about 2500 in our model, corresponding to approximately 10\% of the predicted total number of DNSs in the galaxy. Using our model calibrated to the Galactic DNS population, we make predictions for DNS mergers observed in gravitational waves. The DNS chirp mass distribution varies from 1.1M$_\odot$ to 2.1M$_\odot$ and the median is found to be 1.14\,M$_\mathrm{\odot}$. The expected effective spin $\chi_\mathrm{eff}$ for isolated DNSs is $\lesssim$0.03 from our model. We predict that 34\% of the current Galactic isolated DNSs will merge within a Hubble time, and have a median total mass of 2.7\,M$_\mathrm{\odot}$. Finally, we discuss implications for fast radio bursts and post-merger remnant gravitational-waves.

Supernova remnants (SNRs) have a variety of overall morphology as well as rich structures over a wide range of scales. Quantitative study of these structures can potentially reveal fluctuations of density and magnetic field originating from the interaction with ambient medium and turbulence in the expanding ejecta. We have used $1.5$GHz (L band) and $5$GHz (C band) VLA data to estimate the angular power spectrum $C_{\ell}$ of the synchrotron emission fluctuations of the Kepler SNR. This is done using the novel, visibility based, Tapered Gridded Estimator of $C_{\ell}$. We have found that, for $\ell = (1.9 - 6.9) \times 10^{4}$, the power spectrum is a broken power law with a break at $\ell = 3.3 \times 10^{4}$, and power law index of $-2.84\pm 0.07$ and $-4.39\pm 0.04$ before and after the break respectively. The slope $-2.84$ is consistent with 2D Kolmogorov turbulence and earlier measurements for the Tycho SNR. We interpret the break to be related to the shell thickness of the SNR ($0.35 $ pc) which approximately matches $\ell = 3.3 \times 10^{4}$ (i.e., $0.48$ pc). However, for $\ell > 6.9 \times 10^{4}$, the estimated $C_{\ell}$ of L band is likely to have dominant contribution from the foregrounds while for C band the power law slope $-3.07\pm 0.02$ is roughly consistent with $3$D Kolmogorov turbulence like that observed at large $\ell$ for Cas A and Crab SNRs.