We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 seconds on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications.

The neutron star equation of state (EOS) is now being constrained from a diverse set of multi-messenger data, including gravitational waves from binary neutron star mergers, X-ray observations of the neutron star radius, and many types of laboratory nuclear experiments. These measurements are typically mapped to a common domain -- either to a corresponding radius or to a parametrized EOS using a Bayesian inference scheme -- for comparison with one another. We explore here the statistical biases that can arise when such multi-messenger data are mapped to a common domain for comparison. We find that placing Bayesian priors individually in each domain of measurement can transform to biased constraints in the domain of comparison. Using the first two binary neutron star mergers as an example, we show that a uniform prior in the tidal deformability can produce artificial evidence for large radii, which the data do not support. We present a new prescription for defining Bayesian priors in any domain of measurement, that will allow for minimally-biased constraints in the domain of comparison. Finally, using this new prescription, we provide a status update on multi-messenger EOS constraints on the neutron star radius.

The nuclear symmetry energy plays a role in determining both the nuclear properties of terrestrial matter as well as the astrophysical properties of neutron stars. The first measurement of the neutron star tidal deformability, from gravitational wave event GW170817, provides a new way of probing the symmetry energy. In this work, we report on new constraints on the symmetry energy from GW170817. We focus in particular on the low-order coefficients: namely, the value of the symmetry energy at the nuclear saturation density, S_0, and the slope of the symmetry energy, L_0. We find that the gravitational wave data are relatively insensitive to S_0, but that they depend strongly on L_0 and point to lower values of L_0 than have previously been reported, with a peak likelihood near L_0 ~ 20 MeV. Finally, we use the inferred posteriors on L_0 to derive new analytic constraints on higher-order nuclear terms.

The first detection of gravitational waves from a neutron star-neutron star merger, GW170817, has opened up a new avenue for constraining the ultradense-matter equation of state (EOS). The deviation of the observed waveform from a point-particle waveform is a sensitive probe of the EOS controlling the merging neutron stars' structure. In this topical review, I discuss the various constraints that have been made on the EOS in the year following the discovery of GW170817. In particular, I review the surprising relationship that has emerged between the effective tidal deformability of the binary system and the neutron star radius. I also report new results that make use of this relationship, finding that the radius inferred from GW170817 lies between 9.8 and 13.2 km at 90% confidence, with distinct likelihood peaks at 10.8 and 12.3 km. I compare these radii, as well as those inferred in the literature, to X-ray measurements of the neutron star radius. I also summarize the various maximum mass constraints, which point towards a maximum mass < 2.3 M_sun, depending on the fate of the remnant, and which can be used to additionally constrain the high-density EOS. I review the constraints on the EOS that have been performed directly, through Bayesian inference schemes. Finally, I comment on the importance of disentangling thermal effects in future EOS constraints from neutron star mergers.

Observations of isolated neutron stars place constraints on the equation of state (EOS) of cold, neutron-rich matter, while nuclear physics experiments probe the EOS of hot, symmetric matter. Many dynamical phenomena, such as core-collapse supernovae, the formation and cooling of proto-neutron stars, and neutron star mergers, lie between these two regimes and depend on the EOS at finite temperatures for matter with varying proton fractions. In this paper, we introduce a new framework to accurately calculate the thermal pressure of neutron-proton-electron matter at arbitrary density, temperature, and proton fraction. This framework can be expressed using a set of five physically-motivated parameters that span a narrow range of values for realistic EOS and are able to capture the leading-order effects of degenerate matter on the thermal pressure. We base two of these parameters on a new approximation of the Dirac effective mass, with which we reproduce the thermal pressure to within <~30% for a variety of realistic EOS at densities of interest. Three additional parameters, based on the behavior of the symmetry energy near the nuclear saturation density, allow for the extrapolation of any cold EOS in beta-equilibrium to arbitrary proton fractions. Our model thus allows a user to extend any cold nucleonic EOS, including piecewise-polytropes, to arbitrary temperature and proton fraction, for use in calculations and numerical simulations of astrophysical phenomena. We find that our formalism is able to reproduce realistic finite-temperature EOS with errors of <~20% and offers a 1-3 orders-of-magnitude improvement over existing ideal-fluid models.

Gravitational waves from the coalescence of two neutron stars were recently detected for the first time by the LIGO-Virgo collaboration, in event GW170817. This detection placed an upper limit on the effective tidal deformability of the two neutron stars and tightly constrained the chirp mass of the system. We report here on a new simplification that arises in the effective tidal deformability of the binary, when the chirp mass is specified. We find that, in this case, the effective tidal deformability of the binary is surprisingly independent of the component masses of the individual neutron stars, and instead depends primarily on the ratio of the chirp mass to the neutron star radius. Thus, a measurement of the effective tidal deformability can be used to directly measure the neutron star radius. We find that the upper limit on the effective tidal deformability from GW170817 implies that the radius cannot be larger than ~13km, at the 90% level, independent of the assumed masses for the component stars. The result can be applied generally, to probe the stellar radii in any neutron star-neutron star merger with a measured chirp mass. The approximate mass-independence disappears for neutron star-black hole mergers. Finally, we discuss a Bayesian inference of the equation of state that uses the measured chirp mass and tidal deformability from GW170817 combined with nuclear and astrophysical priors and discuss possible statistical biases in this inference.

The mass distribution of compact objects provides a fossil record that can be studied to uncover information on the late stages of massive star evolution, the supernova explosion mechanism, and the dense matter equation of state. Observations of neutron star masses indicate a bimodal Gaussian distribution, while the observed black hole mass distribution decays exponentially for stellar-mass black holes. We use these observed distributions to directly confront the predictions of stellar evolution models and the neutrino-driven supernova simulations of Sukhbold et al. (2016). We find excellent agreement between the black hole and low-mass neutron star distributions created by these simulations and the observations. We show that a large fraction of the stellar envelope must be ejected, either during the formation of stellar-mass black holes or prior to the implosion through tidal stripping due to a binary companion, in order to reproduce the observed black hole mass distribution. We also determine the origins of the bimodal peaks of the neutron star mass distribution, finding that the low-mass peak (centered at ~1.4 M_sun) originates from progenitors with M_zams ~ 9-18 M_sun. The simulations fail to reproduce the observed peak of high-mass neutron stars (centered at ~1.8 M_sun) and we explore several possible explanations. We argue that the close agreement between the observed and predicted black hole and low-mass neutron star mass distributions provides new promising evidence that these stellar evolution and explosion models are accurately capturing the relevant stellar, nuclear, and explosion physics involved in the formation of compact objects.

One of the key goals of observing neutron stars is to infer the equation of state (EoS) of the cold, ultradense matter in their interiors. We present here a Bayesian statistical method of inferring the pressures at five fixed densities, from a sample of mock neutron star masses and radii. We show that while five polytropic segments are needed for maximum flexibility in the absence of any prior knowledge of the EoS, regularizers are also necessary to ensure that simple underlying EoS are not over-parametrized. For ideal data with small measurement uncertainties, we show that the pressure at roughly twice the nuclear saturation density, rho_sat, can be inferred to within 0.3 dex for many realizations of potential sources of uncertainties. The pressures of more complicated EoS with significant phase transitions can also be inferred to within ~30%. We also find that marginalizing the multi-dimensional parameter space of pressure to infer a mass-radius relation can lead to biases of nearly 1 km in radius, towards larger radii. Using the full, five-dimensional posterior likelihoods avoids this bias.

The increasing number and precision of measurements of neutron star masses, radii, and, in the near future, moments of inertia offer the possibility of precisely determining the neutron star equation of state. One way to facilitate the mapping of observables to the equation of state is through a parametrization of the latter. We present here a generic method for optimizing the parametrization of any physically allowed EoS. We use mock equations of state that incorporate physically diverse and extreme behavior to test how well our parametrization reproduces the global properties of the stars, by minimizing the errors in the observables mass, radius, and the moment of inertia. We find that using piecewise polytropes and sampling the EoS with five fiducial densities between ~1-8 times the nuclear saturation density results in optimal errors for the smallest number of parameters. Specifically, it recreates the radii of the assumed EoS to within less than 0.5 km for the extreme mock equations of state and to within less than 0.12 km for 95% of a sample of 42 proposed, physically-motivated equations of state. Such a parametrization is also able to reproduce the maximum mass to within 0.04 M_sun and the moment of inertia of a 1.338 M_sun neutron star to within less than 10% for 95% of the proposed sample of equations of state.

A precise moment of inertia measurement for PSR J0737-3039A in the double pulsar system is expected within the next five years. We present here a new method of mapping the anticipated measurement of the moment of inertia directly into the neutron star structure. We determine the maximum and minimum values possible for the moment of inertia of a neutron star of a given radius based on physical stability arguments, assuming knowledge of the equation of state only at densities below the nuclear saturation density. If the equation of state is trusted up to the nuclear saturation density, we find that a measurement of the moment of inertia will place absolute bounds on the radius of PSR J0737-3039A to within $\pm$1 km. The resulting combination of moment of inertia, mass, and radius measurements for a single source will allow for new, stringent constraints on the dense-matter equation of state.

We have detected sporadic, bright, short-duration radio pulses from PSR J0901$-$4624. These pulses are emitted simultaneously with persistent, periodic emission that dominates the flux density when averaging over many periods of the pulsar. The bright pulses have energies that are consistent with a power-law distribution. The integrated profile of PSR J0901$-$4624 is highly polarized and shows four distinct components. The bright pulses appear to originate near the magnetic pole of the pulsar and have polarization properties unlike that of the underlying emission at the same pulse phase. We conclude that the bright pulses represent a secondary giant-micropulse emission process, possibly from a different region in the pulsar magnetosphere.

High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broadband, pulse-shape variations intrinsic to the pulsar emission process. These variations induce jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short-term profile and timing stability of 22 millisecond pulsars that are part of the Parkes Pulsar Timing Array sample by examining intra-observation arrival time variability and single-pulse phenomenology. In 7 of the 22 pulsars, in the band centred at approximately 1400MHz, we find that the brightest observations are limited by intrinsic jitter. We find consistent results, either detections or upper limits, for jitter noise in other frequency bands. PSR J1909-3744 shows the lowest levels of jitter noise, which we estimate to contribute $\sim$10 ns root mean square error to the arrival times for hour-duration observations. Larger levels of jitter noise are found in pulsars with wider pulses and distributions of pulse intensities. The jitter noise in PSR J0437-4715 decorrelates over a bandwidth of $\sim$2 GHz. We show that the uncertainties associated with timing pulsar models can be improved by including physically motivated jitter uncertainties. Pulse-shape variations will limit the timing precision at future, more sensitive, telescopes; it is imperative to account for this noise when designing instrumentation and timing campaigns for these facilities.