We present time-of-arrival measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of time-of-arrival (TOA) measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for fourteen pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al.; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision multi-band measurements in several pulsars indicate frequency-dependent DMs. The TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.

We describe the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) and its efforts to directly detect and study gravitational waves and other synergistic physics and astrophysics using radio timing observations of millisecond pulsars.

A major goal for the astronomy and astrophysics communities is the pursuit of diversity, equity, and inclusion (DEI) in all ranks, from students through professional scientific researchers. Large scientific collaborations - increasingly a primary place for both professional interactions and research opportunities - can play an important role in the DEI effort. Multimessenger astronomy, a new and growing field, is based on the principle that working collaboratively produces synergies, enabling advances that would not be possible without cooperation. The nascent Multimessenger Diversity Network (MDN) is extending this collaborative approach to include DEI initiatives. After we review of the current state of DEI in astronomy and astrophysics, we describe the strategies the MDN is developing and disseminating to support and increase DEI in the fields over the coming decade: provide opportunities (real and virtual) to share DEI knowledge and resources, include DEI in collaboration-level activities, including external reviews, and develop and implement ways to recognize the DEI work of collaboration members.

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.

Pulsar timing arrays are sensitive to gravitational wave perturbations produced by individual supermassive black hole binaries during their early inspiral phase. Modified gravity theories allow for the emission of gravitational dipole radiation, which is enhanced relative to the quadrupole contribution for low orbital velocities, making the early inspiral an ideal regime to test for the presence of modified gravity effects. Using a theory-agnostic description of modified gravity theories based on the parametrized post-Einsteinian framework, we explore the possibility of detecting deviations from General Relativity using simulated pulsar timing array data, and provide forecasts for the constraints that can be achieved. We generalize the {\tt enterprise} pulsar timing software to account for possible additional polarization states and modifications to the phase evolution, and study how accurately the parameters of simulated signals can be recovered. We find that while a pure dipole model can partially recover a pure quadrupole signal, there is little possibility for confusion when the full model with all polarization states is used. With no signal present, and using noise levels comparable to those seen in contemporary arrays, we produce forecasts for the upper limits that can be placed on the amplitudes of alternative polarization modes as a function of the sky location of the source.

Precision timing of large arrays (>50) of millisecond pulsars will detect the nanohertz gravitational-wave emission from supermassive binary black holes within the next ~3-7 years. We review the scientific opportunities of these detections, the requirements for success, and the synergies with electromagnetic instruments operating in the 2020s.

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $\sigma_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $\sigma_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $\sigma_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years.

We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; this allows a precise estimate of the transverse velocity, 123(5) km/s. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a full 3-D determination of the system's velocity; b) precise measurements of the rate of advance of periastron, which after subtraction of the contribution of the proper motion yields a total system mass of $1.6518^{+0.0033}_{-0.0035}$ solar masses; c) a Shapiro delay measurement, h_3 = $82 \pm 14$ ns despite the orbital inclination not being near 90 deg; combined with the measurement of the total mass, this yields a pulsar mass of $1.353^{+0.014}_{-0.017}$ solar masses and a companion mass of $0.298^{+0.015}_{-0.012}$ solar masses; d) we measure precisely the secular variation of the projected semi-major axis and detect significant annual orbital parallax; together these allow a determination of the full 3-D orbital geometry, including an unambiguous orbital inclination (i = $138.7^{+2.5}_{-2.2}$ deg) and a position angle for the line of nodes (Omega = $44^{+5}_{-4}$ deg). We discuss the component masses to investigate hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the full 3-D position, motion and orbital orientation of the system, plus the precisely measured pulsar and WD mass and the latter's optical detection make this system an unique test of our understanding of white dwarfs and their atmospheres.

Low-frequency gravitational-wave experiments require the highest timing precision from an array of the most stable millisecond pulsars. Several known sources of noise on short timescales in single radio-pulsar observations are well described by a simple model of three components: template-fitting from a finite signal-to-noise ratio, pulse phase/amplitude jitter from single-pulse stochasticity, and scintillation errors from short-timescale interstellar scattering variations. Currently template-fitting errors dominate, but as radio telescopes push towards higher signal-to-noise ratios, jitter becomes the next dominant term for most millisecond pulsars. Understanding the statistics of jitter becomes crucial for properly characterizing arrival-time uncertainties. We characterize the radio-frequency dependence of jitter using data on 48 pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing program. We detect significant jitter in 43 of the pulsars and test several functional forms for its frequency dependence; we find significant frequency dependence for 30 pulsars. We find moderate correlations of rms jitter with pulse width (R = 0.62) and number of profile components (R = 0.40); the single-pulse rms jitter is typically ~1% of pulse phase. The average frequency dependence for all pulsars using a power-law model has index -0.42. We investigate the jitter variations for the interpulse of PSR~B1937+21 and find no significant deviations from the main pulse rms jitter. We also test the time-variation of jitter in two pulsars and find that systematics likely bias the results for high-precision pulsars. Pulsar timing array analyses must properly model jitter as a significant component of the noise within the detector.

Observations have revealed that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, these SMBHs form SMBH binaries (SMBHBs) that emit low-frequency gravitational waves (GWs). The incoherent superposition of these sources produce a stochastic GW background (GWB) that can be observed by pulsar timing arrays (PTAs). The optimal statistic is a frequentist estimator of the amplitude of the GWB that specifically looks for the spatial correlations between pulsars induced by the GWB. In this paper, we introduce an improved method for computing the optimal statistic that marginalizes over the red noise in individual pulsars. We use simulations to demonstrate that this method more accurately determines the strength of the GWB, and we use the noise-marginalized optimal statistic to compare the significance of monopole, dipole, and Hellings-Downs (HD) spatial correlations and perform sky scrambles.

The frequency dependence of radio pulse arrival times provides a probe of structures in the intervening media. Demorest et al. 2013 was the first to show a short-term (~100-200 days) reduction in the electron content along the line of sight to PSR J1713+0747 in data from 2008 (approximately MJD 54750) based on an apparent dip in the dispersion measure of the pulsar. We report on a similar event in 2016 (approximately MJD 57510), with average residual pulse-arrival times of approximately 3.0,-1.3, and -0.7 microseconds at 820, 1400, and 2300 MHz, respectively. Timing analyses indicate possible departures from the standard nu^-2 dispersive-delay dependence. We discuss and rule out a wide variety of potential interpretations. We find the likeliest scenario to be lensing of the radio emission by some structure in the interstellar medium, which causes multiple frequency-dependent pulse arrival-time delays.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.

Previous optical and radio observations of the binary millisecond pulsar PSR J1640+2224 have come to inconsistent conclusions about the identity of its companion, with some observations suggesting the companion is a low-mass helium-core (He-core) white dwarf (WD), while others indicate it is most likely a high-mass carbon-oxygen (CO) WD. Binary evolution models predict PSR J1640+2224 most likely formed in a low-mass X-ray binary (LMXB) based on the pulsar's short spin period and long-period, low-eccentricity orbit, in which case its companion should be a He-core WD with mass about $0.35 - 0.39 \, M_\odot$, depending on metallicity. If it is instead a CO WD, that would suggest the system has an unusual formation history. In this paper we present the first astrometric parallax measurement for this system from observations made with the Very Long Baseline Array (VLBA), from which we determine the distance to be $1520^{+170}_{-150}\,\mathrm{pc}$. We use this distance and a reanalysis of archival optical observations originally taken in 1995 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST) in order to measure the WD's mass. We also incorporate improvements in calibration, extinction model, and WD cooling models. We find that the existing observations are not sufficient to tightly constrain the companion mass, but we conclude the WD mass is $>0.4\,M_\odot$ with $>90\%$ confidence. The limiting factor in our analysis is the low signal-to-noise ratio of the original HST observations.

One of the primary gravitational wave (GW) sources for pulsar timing arrays (PTAs) is the stochastic background formed by supermassive black holes binaries (SMBHBs). In this paper, we investigate how the environments of SMBHBs will effect the sensitivity of PTAs by deriving scaling laws for the signal-to-noise ratio (SNR) of the optimal cross-correlation statistic. The presence of gas and stars around SMBHBs will accelerate the merger at large distances, depleting the GW stochastic background at low frequencies. We show that environmental interactions may delay detection by a few years or more, depending on the PTA configuration and the frequency at which the dynamical evolution transitions from being dominated by environmental effects to GW-dominated.

A small body orbiting a black hole follows a trajectory that, at leading order, is a geodesic of the black hole spacetime. Much effort has gone into computing "self force" corrections to this motion, arising from the small body's own contributions to the system's spacetime. Another correction to the motion arises from coupling of the small body's spin to the black hole's spacetime curvature. Spin-curvature coupling drives a precession of the small body, and introduces a "force" (relative to the geodesic) which shifts the small body's worldline. These effects scale with the small body's spin at leading order. In this paper, we show that the equations which govern spin-curvature coupling can be analyzed with a frequency-domain decomposition, at least to leading order in the small body's spin. We show how to compute the frequency of precession along generic orbits, and how to describe the small body's precession and motion in the frequency domain. We illustrate this approach with a number of examples. This approach is likely to be useful for understanding spin coupling effects in the extreme mass ratio limit, and may provide insight into modeling spin effects in the strong field for non-extreme mass ratios.

Model-independent distance constraints to binary millisecond pulsars (MSPs) are of great value to both the timing observations of the radio pulsars, and multiwavelength observations of their companion stars. Very Long Baseline Interferometry (VLBI) astrometry can be employed to provide these model-independent distances with very high precision via the detection of annual geometric parallax. Using the Very Long Baseline Array, we have observed two binary millisecond pulsars, PSR J1022+1001 and J2145-0750, over a two-year period and measured their distances to be 700 +14 -10 pc and 613 +16 -14 pc respectively. We use the well-calibrated distance in conjunction with revised analysis of optical photometry to tightly constrain the nature of their massive (M ~ 0.85 Msun) white dwarf companions. Finally, we show that several measurements of their parallax and proper motion of PSR J1022+1001 and PSR J2145-0750 obtained by pulsar timing array projects are incorrect, differing from the more precise VLBI values by up to 5 sigma. We investigate possible causes for the discrepancy, and find that imperfect modeling of the solar wind is a likely candidate for the timing model errors given the low ecliptic latitude of these two pulsars.

The extremely regular, periodic radio emission from millisecond pulsars makes them useful tools for studying neutron star astrophysics, general relativity, and low-frequency gravitational waves. These studies require that the observed pulse times of arrival be fit to complex timing models that describe numerous effects such as the astrometry of the source, the evolution of the pulsar's spin, the presence of a binary companion, and the propagation of the pulses through the interstellar medium. In this paper, we discuss the benefits of using Bayesian inference to obtain pulsar timing solutions. These benefits include the validation of linearized least-squares model fits when they are correct, and the proper characterization of parameter uncertainties when they are not; the incorporation of prior parameter information and of models of correlated noise; and the Bayesian comparison of alternative timing models. We describe our computational setup, which combines the timing models of Tempo2 with the nested-sampling integrator MultiNest. We compare the timing solutions generated using Bayesian inference and linearized least-squares for three pulsars: B1953+29, J2317+1439, and J1640+2224, which demonstrate a variety of the benefits that we posit.

General relativity predicts the existence of black holes, compact objects whose spacetimes depend on only their mass, spin, and charge in vacuum (the "no hair" theorem). As various observations probe deeper into the strong fields of black hole candidates, it is becoming possible to test this prediction. Previous work suggested that such tests can be performed by measuring whether the multipolar structure of black hole candidates has the form that general relativity demands, and introduced a family of "bumpy black hole" spacetimes to be used for making these measurements. These spacetimes have generalized multipoles, where the deviation from the Kerr metric depends on the spacetime's "bumpiness." In this paper, we show how to compute the Geroch-Hansen moments of a bumpy black hole, demonstrating that there is a clean mapping between the deviations used in the bumpy black hole formalism and the Geroch-Hansen moments. We also extend our previous results to define bumpy black holes whose {\it current} moments, analogous to magnetic moments of electrodynamics, deviate from the canonical Kerr value.

Our universe contains a great number of extremely compact and massive objects which are generally accepted to be black holes. Precise observations of orbital motion near candidate black holes have the potential to determine if they have the spacetime structure that general relativity demands. As a means of formulating measurements to test the black hole nature of these objects, Collins and Hughes introduced "bumpy black holes": objects that are almost, but not quite, general relativity's black holes. The spacetimes of these objects have multipoles that deviate slightly from the black hole solution, reducing to black holes when the deviation is zero. In this paper, we extend this work in two ways. First, we show how to introduce bumps which are smoother and lead to better behaved orbits than those in the original presentation. Second, we show how to make bumpy Kerr black holes -- objects which reduce to the Kerr solution when the deviation goes to zero. This greatly extends the astrophysical applicability of bumpy black holes. Using Hamilton-Jacobi techniques, we show how a spacetime's bumps are imprinted on orbital frequencies, and thus can be determined by measurements which coherently track a small orbiting body's orbital phase. We find that weak-field orbits of bumpy black holes are modified exactly as expected from a Newtonian analysis of a body with a prescribed multipolar structure, reproducing well-known results from the celestial mechanics literature. The impact of bumps on strong-field orbits is especially strong, suggesting that this framework will allow observations to set robust limits on the extent to which a spacetime's multipoles deviate from the black hole expectation.

We present strong observational evidence for a relationship between the direction of a pulsar's motion and its rotation axis. We show carefully calibrated polarization data for 25 pulsars, 20 of which display linearly polarized emission from the pulse longitude at closest approach to the magnetic pole. Such data allow determination of the position angle of the linear polarisation which in turn reflects the position angle of the rotation axis. Of these 20 pulsars, 10 show an offset between the velocity vector and the polarisation position angle which is either less than 10\degr or more than 80\degr, a fraction which is very unlikely by random chance. We believe that the bimodal nature of the distribution arises from the presence of orthogonal polarisation modes in the pulsar radio emission. In some cases this orthogonal ambiguity is resolved by observations at other wavelengths so that we conclude that the velocity vector and the rotation axis are aligned at birth. Strengthening the case is the fact that 4 of the 5 pulsars with ages less than 3 Myr show this relationship, including the Vela pulsar. We discuss the implications of these findings in the context of the Spruit & Phinney (1998)\nocite{sp98} model of pulsar birth-kicks. We point out that, contrary to claims in the literature, observations of double neutron star systems do not rule out aligned kick models and describe a possible observational test involving the double pulsar system.