We present a family of modified Hermite integrators of arbitrary order possessing superior behaviour for the integration of Keplerian and near-Keplerian orbits. After recounting the derivation of Hermite N-body integrators of arbitrary order, we derive a corrector expression that minimises integrated errors in the argument of periapsis for any such integrator. In addition to providing an alternate derivation of the modified corrector for the 4th-order Hermite integrator, we focus on improved correctors for the 6th- and 8th-order integrators previously featured in the literature. We present a set of numerical examples and find that the higher-order schemes improve performance, even when considering their slightly higher cost in floating point operations. The algorithms presented herein hold promise for systems dominated by central potentials, such as planetary systems and the centres of galaxies. Existing Hermite integrators of any order can be modified to use the expressions presented here with minimal effort. Accordingly the schemes presented herein can be easily implemented on massively parallel architectures.

We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes.

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density.

Accretion disks around active galactic nuclei are potentially unstable to star formation at large radii. We note that when the compact objects formed from some of these stars spiral into the central supermassive black hole, there is no radiative feedback and therefore the accretion rate is not limited by radiation forces. Using a set of accretion disk models, we calculate the accretion rate onto the central supermassive black hole in both gas and compact objects. We find that the timescale for a supermassive black hole to double in mass can decrease by factors ranging from $\sim0.7$ to as low as $\sim0.1$ in extreme cases, compared to gas accretion alone. Our results suggest that the formation of extremely massive black holes at high redshift may occur without prolonged super-Eddington gas accretion or very massive seed black holes. We comment on potential observational signatures as well as implications for other observations of active galactic nuclei.

A small fraction of quasars show an unusually high nitrogen-to-carbon ratio (N/C) in their spectra. These "nitrogen-rich" (N-rich) quasars are a long-standing puzzle because their interstellar medium implies stellar populations with abnormally high metallicities. It has recently been proposed that N-rich quasars may result from tidal disruption events (TDEs) of stars by supermassive black holes. The rapid enhancement of nitrogen and the depletion of carbon due to the carbon--nitrogen--oxygen cycle in supersolar mass stars could naturally produce high N/C. However, the TDE hypothesis predicts that the N/C should change with time, which has never hitherto been observed. Here we report the discovery of the first N-rich quasar with rapid N/C variability that could be caused by a TDE. Two spectra separated by 1.7 years (rest-frame) show that the N III]\lambda 1750/C III]\lambda 1909 intensity ratio decayed by ~86%+/-14% (1 \sigma). Optical (rest-frame UV) light-curve and X-ray observations are qualitatively consistent with the TDE hypothesis; though, the time baseline falls short of a definitive proof. Putting the single-object discovery into context, statistical analyses of the ~80 known N-rich quasars with high-quality archival spectra show evidence (at a 5\sigma\ significance level) of a decrease in N/C on timescales of >1 year (rest-frame) and a constant level of ionization (indicated by the C III]\lambda 1909/C IV\lambda 1549 intensity ratio). If confirmed, our results demonstrate the method of identifying TDE candidates in quasars via abundance ratio variability, opening a new window of TDE observations at high redshift (z>2) with upcoming large-scale time-domain spectroscopic surveys.

In the summer of 2012, during a Pulsar Search Collaboratory workshop, two high-school students discovered J1930$-$1852, a pulsar in a double neutron star (DNS) system. Most DNS systems are characterized by short orbital periods, rapid spin periods and eccentric orbits. However, J1930$-$1852 has the longest spin period ($P_{\rm spin}\sim$185 ms) and orbital period ($P_{\rm b}\sim$45 days) yet measured among known, recycled pulsars in DNS systems, implying a shorter than average and/or inefficient recycling period before its companion went supernova. We measure the relativistic advance of periastron for J1930$-$1852, $\dot{\omega}=0.00078$(4) deg/yr, which implies a total mass (M$_{\rm{tot}}=2.59$(4) M$_{\odot}$) consistent with other DNS systems. The $2\sigma$ constraints on M$_{\rm{tot}}$ place limits on the pulsar and companion masses ($m_{\rm p}<1.32$ M$_{\odot}$ and $m_{\rm c}>1.30$ M$_{\odot}$ respectively). J1930$-$1852's spin and orbital parameters challenge current DNS population models and make J1930$-$1852 an important system for further investigation.