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.

Close Double Neutron Stars (DNSs) have been observed as Galactic radio pulsars, while their mergers have been detected as gamma-ray bursts and gravitational-wave sources. They are believed to have experienced at least one common-envelope episode (CEE) during their evolution prior to DNS formation. In the last decades there have been numerous efforts to understand the details of the common-envelope phase, but its computational modelling remains challenging. We present and discuss the properties of the donor and the binary at the onset of the Roche-lobe overflow (RLOF) leading to these CEEs as predicted by rapid binary population synthesis models. These properties can be used as initial conditions for detailed simulations of the common-envelope phase. There are three distinctive populations, classified by the evolutionary stage of the donor at the moment of the onset of the RLOF: giant donors with fully-convective envelopes, cool donors with partially-convective envelopes, and hot donors with radiative envelopes. We also estimate that, for standard assumptions, tides would not circularise a large fraction of these systems by the onset of RLOF. This makes the study and understanding of eccentric mass transferring systems relevant for DNS populations.

A class of optical transients known as Luminous Red Novae (LRNe) have recently been associated with mass ejections from binary stars undergoing common-envelope evolution. We use the population synthesis code COMPAS to explore the impact of a range of assumptions about the physics of common-envelope evolution on the properties of LRNe. In particular, we investigate the influence of various models for the energetics of LRNe on the expected event rate and light curve characteristics, and compare with the existing sample. We find that the Galactic rate of LRNe is $\sim 0.2$ yr$^{-1}$, in agreement with the observed rate. In our models, the luminosity function of Galactic LRNe covers multiple decades in luminosity and is dominated by signals from stellar mergers, consistent with observational constraints from iPTF and the Galactic sample of LRNe. We discuss how observations of the brightest LRNe may provide indirect evidence for the existence of massive ($> 40$ M$_\odot$) red supergiants. Such LRNe could be markers along the evolutionary pathway leading to the formation of double compact objects. We make predictions for the population of LRNe observable in future transient surveys with the Large Synoptic Survey Telescope and the Zwicky Transient Facility. In all plausible circumstances, we predict a selection-limited observable population dominated by bright, long-duration events caused by common envelope ejections. We show that the Large Synoptic Survey Telescope will observe $20$--$750$ LRNe per year, quickly constraining the luminosity function of LRNe and probing the physics of common-envelope events.

Glitch size and waiting time probability density functions (PDFs) are estimated for the five pulsars that have glitched most using the nonparametric kernel density estimator. Two objects exhibit decreasing size and waiting time PDFs. Their activity is Poisson-like, and their size statistics are approximately scale-invariant. Three objects exhibit a statistically significant local maximum in the PDFs, including one (PSR J1341$-$6220) which was classified as Poisson-like in previous analyses. Their activity is quasiperiodic, although the dispersion in waiting times is relatively broad. The classification is robust: it is preserved across a wide range of bandwidth choices. There is no compelling evidence for multimodality, but this issue should be revisited when more data become available. The implications for superfluid vortex avalanche models of pulsar glitches are explored briefly.

Few statistically compelling correlations are found in pulsar timing data between the size of a rotational glitch and the time to the preceding glitch (backward waiting time) or the succeeding glitch (forward waiting time), except for a strong correlation between sizes and forward waiting times in PSR J0537-6910. This situation is counterintuitive, if glitches are threshold-triggered events, as in standard theories (e.g. starquakes, superfluid vortex avalanches). Here it is shown that the lack of correlation emerges naturally, when a threshold trigger is combined with secular stellar braking slower than a critical, calculable rate. The Pearson and Spearman correlation coefficients are computed and interpreted within the framework of a state-dependent Poisson process. Specific, falsifiable predictions are made regarding what objects currently targeted by long-term timing campaigns should develop strong size-waiting-time correlations, as more data are collected in the future.

Glitches are sudden jumps in the spin frequency of pulsars believed to originate in the superfluid interior of neutron stars. Superfluid flow in a model neutron star is simulated by solving the equations of motion of a two-component superfluid consisting of a viscous proton-electron plasma and an inviscid neutron condensate in a spherical Couette geometry. We examine the response of our model neutron star to glitches induced in three different ways: by instantaneous changes of the spin frequency of the inner and outer boundaries, and by instantaneous recoupling of the fluid components in the bulk. All simulations are performed with strong and weak mutual friction. It is found that the maximum size of a glitch that originates in the bulk decreases as the mutual friction strengthens. It is also found that mutual friction determines the fraction of the frequency jump which is later recovered, a quantity known as the 'healing parameter'. These behaviours may explain some of the diversity in observed glitch recoveries.