Compact Object Mergers: Population Astrophysics and Statistics (COMPAS; https://compas.science) is a public rapid binary population synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parameterized assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It includes a number of tools for population processing in addition to the core binary evolution components. COMPAS is publicly available via the GitHub repository https://github.com/TeamCOMPAS/COMPAS/, and is designed to allow for flexible modifications as evolutionary models improve. This paper describes the methodology and implementation of COMPAS. It is a living document that will be updated as new features are added to COMPAS; the current document describes COMPAS v02.21.00.

Einstein Telescope (ET) is the European project for a gravitational-wave (GW) observatory of third-generation. In this paper we present a comprehensive discussion of its science objectives, providing state-of-the-art predictions for the capabilities of ET in both geometries currently under consideration, a single-site triangular configuration or two L-shaped detectors. We discuss the impact that ET will have on domains as broad and diverse as fundamental physics, cosmology, early Universe, astrophysics of compact objects, physics of matter in extreme conditions, and dynamics of stellar collapse. We discuss how the study of extreme astrophysical events will be enhanced by multi-messenger observations. We highlight the ET synergies with ground-based and space-borne GW observatories, including multi-band investigations of the same sources, improved parameter estimation, and complementary information on astrophysical or cosmological mechanisms obtained combining observations from different frequency bands. We present advancements in waveform modeling dedicated to third-generation observatories, along with open tools developed within the ET Collaboration for assessing the scientific potentials of different detector configurations. We finally discuss the data analysis challenges posed by third-generation observatories, which will enable access to large populations of sources and provide unprecedented precision.

As the era of gravitational-wave astronomy has well and truly begun, gravitational radiation from rotating neutron stars remains elusive. Rapidly spinning neutron stars are the main targets for continuous-wave searches since, according to general relativity, provided they are asymmetrically deformed, they will emit gravitational waves. It is believed that detecting such radiation will unlock the answer to why no pulsars have been observed to spin close to the break-up frequency. We review existing studies on the maximum mountain that a neutron star crust can support, critique the key assumptions and identify issues relating to boundary conditions that need to be resolved. In light of this discussion, we present a new scheme for modelling neutron star mountains. The crucial ingredient for this scheme is a description of the fiducial force which takes the star away from sphericity. We consider three examples: a source potential which is a solution to Laplace’s equation, another solution which does not act in the core of the star and a thermal pressure perturbation. For all the cases, we find that the largest quadrupoles are between a factor of a few to two orders of magnitude below previous estimates of the maximum-mountain size.

Rapidly spinning, deformed neutron stars have long been considered potential gravitational-wave emitters. However, so far only upper limits on the size of the involved quadrupole deformations have been obtained. For this reason, it is pertinent to ask how large a mountain can be before the neutron star crust fractures. This is the question we consider in this paper, which describes how mountains can be calculated in relativistic gravity. Formally, this is a perturbative calculation that requires a fiducial force to source the mountain. Therefore, we consider three simple examples and increase their deforming amplitudes until the crust yields. We demonstrate how the derived mountains depend on the equation of state by considering a range of models obtained from chiral effective field theory. We find that the largest mountains depend sensitively on both the mechanism that sources them and the nuclear-matter equation of state.

With the first detections of binary neutron star mergers by gravitational-wave detectors, it proves timely to consider how the internal structure of neutron stars affects the way in which they can be asymmetrically deformed. Such deformations may leave measurable imprints on gravitational-wave signals and can be sourced through tidal interactions or the formation of mountains. We detail the formalism that describes fully relativistic neutron star models with elastic crusts undergoing static perturbations. This formalism primes the problem for studies into a variety of mechanisms that can deform a neutron star. We present results for a barotropic equation of state and a realistic model for the elastic crust, which enables us to compute relevant quantities such as the tidal deformability parameter. We find that the inclusion of an elastic crust provides a very small correction to the tidal deformability. The results allow us to demonstrate when and where the crust starts to fail during a binary inspiral, and we find that the majority of the crust will remain intact up until merger.

The fastest-spinning neutron stars in low-mass X-ray binaries, despite having undergone millions of years of accretion, have been observed to spin well below the Keplerian break-up frequency. We simulate the spin evolution of synthetic populations of accreting neutron stars in order to assess whether gravitational waves can explain this behaviour and provide the distribution of spins that is observed. We model both persistent and transient accretion and consider two gravitational-wave-production mechanisms that could be present in these systems: thermal mountains and unstable rmodes. We consider the case of no gravitational-wave emission and observe that this does not match well with observation. We find evidence for gravitational waves being able to provide the observed spin distribution; the most promising mechanisms being a permanent quadrupole, thermal mountains, and unstable r modes. However, based on the resultant distributions alone, it is difficult to distinguish between the competing mechanisms.