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.

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.

The Green Bank North Celestial Cap (GBNCC) pulsar survey will cover the entire northern sky ($\delta > -40^\circ$) at 350 MHz, and is one of the most uniform and sensitive all-sky pulsar surveys to date. We have created a pipeline to re-analyze GBNCC survey data to take a 350MHz census of all pulsars detected by the survey, regardless of their discovery survey. Of the 1413 pulsars in the survey region, we were able to recover 661. For these we present measured signal-to-noise ratios (S/N), flux densities, pulse widths, profiles, and where appropriate, refined dispersion measurements (647 out of 661) and new or improved spectral indices (276 out of 661 total, 15 new, 261 improved). Detection scans for several hundred sources were reanalyzed in order to inspect pulsars' single pulse behavior and 223 were confirmed to exhibit evidence of nulling. With a detailed analysis of measured and expected S/N values and the evolving radio frequency interference environment at 350MHz, we assess the GBNCC survey's sensitivity as a function of spin period, dispersion measure, and sky position. We find the sky-averaged limiting flux density of the survey to be 0.74mJy. Combining this analysis with PsrPopPy pulsar population simulations, we predict 60/5 non-recycled/millisecond pulsar discoveries in the survey's remaining 21,000 pointings, and we begin to place constraints on population model parameters.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

Despite its importance to our understanding of physics at supranuclear densities, the equation of state (EoS) of matter deep within neutron stars remains poorly understood. Millisecond pulsars (MSPs) are among the most useful astrophysical objects in the Universe for testing fundamental physics, and place some of the most stringent constraints on this high-density EoS. Pulsar timing - the process of accounting for every rotation of a pulsar over long time periods - can precisely measure a wide variety of physical phenomena, including those that allow the measurement of the masses of the components of a pulsar binary system (Lorimer & Kramer 2005). One of these, called relativistic Shapiro delay (Shapiro 1964), can yield precise masses for both an MSP and its companion; however, it is only easily observed in a small subset of high-precision, highly inclined (nearly edge-on) binary pulsar systems. By combining data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-year data set with recent orbital-phase-specific observations using the Green Bank Telescope, we have measured the mass of the MSP J0740+6620 to be $\mathbf{2.14^{+0.10}_{-0.09}}$ solar masses (68.3% credibility interval; 95.4% credibility interval is $\mathbf{2.14^{+0.20}_{-0.18}}$ solar masses). It is highly likely to be the most massive neutron star yet observed, and serves as a strong constraint on the neutron star interior EoS.

The Pulsar Search Collaboratory (PSC) is a collaboration between the Green Bank Observatory and West Virginia University, funded by the National Science Foundation. The PSC program is currently expanding nationwide and engages high school students, teachers, and undergraduate mentors in real-world research by searching for pulsars in data collected with the 100-m Green Bank Telescope. In the process, students learn about observational radio astronomy, radio frequency interference, pulsar timing, and data analysis procedures. The primary goals of the PSC are to stimulate student interest in Science, Technology, Engineering, and Mathematics (STEM) careers, to prepare teachers in implementing authentic research with students by training them within a professional scientific community, and to promote student use of information technologies through online activities and workshops. In this paper, we provide an overview of pulsar science and the data analysis students undertake, as well as a general overview of the program. We then discuss evaluation data collected from participants through a series of survey questions to determine if the program's initial goals were met. The program had a positive impact on the students according to multiple measures, in particular, on their understanding of the nature of scientific inquiry and motivation to pursue STEM career paths.

The new field of gravitational wave astrophysics requires a growing pool of students and researchers with unique, interdisciplinary skill sets. It also offers an opportunity to build a diverse, inclusive astronomy community from the ground up. We describe the efforts used by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) NSF Physics Frontiers Center to foster such growth by involving students at all levels in low-frequency gravitational wave astrophysics with pulsar timing arrays (PTAs) and establishing collaboration policies that ensure broad participation by diverse groups. We describe and illustrate the impact of these techniques on our collaboration as a case study for other distributed collaborations.

We present a Chandra and XMM-Newton study of the SNR CTB 37B, along with archival radio observations. In radio wavelengths, CTB 37B is an incomplete shell showing bright emission from the eastern side while the X-ray morphology shows diffuse emission from regions surrounding the magnetar CXOU J171405.7-381031. We used archival HI absorption measurements to constrain the distance to the remnant as 9.8+/-1.5 kpc. The X-ray spectrum of the remnant is described by a thermal model in the 1-5 keV energy range, with a temperature of 1.3+/-0.1 keV. The abundances from the spectral fits are consistent with being solar or sub-solar. A small region of diffuse emission is seen to the southern side of the remnant, best fitted by a nonthermal spectrum with an unusually hard photon index of 1.3+/-0.3. Assuming a distance of 9.8 kpc to the SNR, we infer a shock velocity of 915+/-70 km/s and explosion energy of (1.8+/-0.6)e50 ergs. The overall imaging and spectral properties of CTB 37B favor the interpretation of a young SNR (6200 yr) propagating in a low-density medium, under the assumption of a Sedov evolutionary phase.

Enabled by the Fermi Large Area Telescope, we now know young and recycled pulsars fill the gamma-ray sky, and we are beginning to understand their emission mechanism and their distribution throughout the Galaxy. However, key questions remain: Is there a large population of pulsars near the Galactic center? Why do the most energetic pulsars shine so brightly in MeV gamma rays but not always at GeV energies? What is the source and nature of the pair plasma in pulsar magnetospheres, and what role does the polar cap accelerator play? Addressing these questions calls for a sensitive, wide-field MeV telescope, which can detect the population of MeV-peaked pulsars hinted at by Fermi and hard X-ray telescopes and characterize their spectral shape and polarization.

A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Energetic particles in our own galaxy, galactic cosmic rays (GCRs), engage in a complex interplay with the interstellar medium and magnetic fields in the galaxy, giving rise to many of its key characteristics. This White Paper is the first of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is galactic cosmic rays, supernova remnants, protostellar jets and superbubbles, and colliding wind binaries.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. This White Paper is the second of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is active galactic nuclei and gamma-ray bursts.

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.

The Pulsar Search Collaboratory (PSC) engages high school students and teachers in analyzing real data from the Robert C. Byrd Green Bank Telescope for the purpose of discovering exotic stars called pulsars. These cosmic clocks can be used as a galactic-scale detector of gravitational waves, ripples in space-time that have recently been directly detected from the mergers of stellar-mass black holes. Through immersing students in an authentic, positive learning environment to build a sense of belonging and competency, the goal of the PSC is to promote students' long-term interests in science and science careers. PSC students have discovered 7 pulsars since the start of the PSC in 2008. Originally targeted at teachers and students in West Virginia, over time the program has grown to 18 states. In a new effort to scale the PSC nationally, the PSC has developed an integrated online training program with both self-guided lectures and homework and real-time interactions with pulsar astronomers. Now, any high school student can join in the exciting search for pulsars and the discovery of a new type of gravitational waves.

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 present a Chandra Director's Discretionary Time observation of PSR J1119-6127 and its compact X-ray pulsar wind nebula (PWN) obtained on 27 October 2016, three months after the Fermi and Swift detection of millisecond bursts in hard X-rays, accompanied by >160 times increase in flux. This magnetar-like activity, the first observed from a rotation-powered radio pulsar, provides an important probe of the physical processes that differentiate radio pulsars from magnetars. The post-burst X-ray spectrum of the pulsar can be described by a single powerlaw model with a photon index of 2.0\pm0.2 and an unabsorbed flux of 5.7e-12 ergs/cm^2/s in the 0.5-7.0 keV energy range. At the time of Chandra observations, the pulsar was still brighter by a factor of ~22 in comparison with its quiescence. The X-ray images reveal a nebula brighter than in the pre-burst Chandra observations (from 2002 and 2004), with an unabsorbed flux of 2.2e-13 ergs/cm^-2/s. This implies a current X-ray efficiency of ~0.001 at a distance of 8.4 kpc. In addition, a faint torus-like structure is visible along the southeast-northwest direction and a jet-like feature perpendicular to the torus towards the southwest. The PWN is best fitted by an absorbed powerlaw with a photon index of 2.2\pm0.5 (post-burst). While the pulsar can still be energetically powered by rotation, the observed changes in PSR~J1119-6127 and its PWN following the magnetar-like bursts point to an additional source of energy powered by its high-magnetic field.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.